1. Introduction
The development of mathematical models of chemical and technological processes has been an actual problem for many decades. Mathematical modeling allows to avoid full-scale experiments and has proven itself as an effective method for solving many technological problems [
1,
2].
An important stage of industrial oil refining is the process of cleaning the catalyst from the coke sediments that accumulate on the surface of the layer and inside the catalyst grains [
3]. One of the least expensive cleaning methods is oxidative regeneration—burning of coke sediments with oxygen-containing gas. An actual problem is to ensure such conditions for the flow of oxidative regeneration, in which coke burns out quickly enough, but at the same time there is no overheating of the catalyst grain.
Currently, the oil and gas industry is focused on the maximum selectivity of raw materials. This is often reached by using complex technological schemes of reactors with significant gas-dynamic resistance. The problem of gas-dynamic resistance can be solved by developing catalysts of various geometric shapes and sizes. The variety of catalyst types leads to the need to develop classes of mathematical models that allow calculating the main characteristics of the oxidative regeneration process for various grain forms, sediment compositions and catalyst materials.
The article [
4] presents a study of processes inside a spherical catalyst grain in the presence of coke sediments. It is noted that coke in the pores of the catalyst significantly reduces diffusion. It is also noted that diffusion models are quite complex; in most cases, they require a numerical solution.
The article [
5] is devoted to the study of the influence of coke on heterogeneous catalysis. As in [
4], a strong influence of coke in the catalyst pores on internal diffusion was noted. The Monte Carlo method was used to simulate the deposition of coke in the pores. Similarly, the removal of coke from the pores is modeled. This shows that coke can be removed unevenly over time.
The paper [
6] presents a two-dimensional (bed length and pellet radius coordinates), two-temperature (gas and catalyst phases) catalyst bed model used for simulation of the oxidative regeneration of coked hydrofluorination CrF3/MgF2 catalyst in the adiabatic reactor. The study discusses the rather complicated influence of process parameters on the main characteristics of the regeneration procedure. In particular, it was shown that feeding of the regenerating gas with temperatures lower than the coke ignition point may lead to the formation of a superadiabatic combustion wave with extremely high temperatures, even under low inlet oxygen concentration.
The article [
7] presents a model of adsorption–catalytic removal of volatile organic compounds. This model also gives the successful example of the adiabatic catalyst bed model, where the interplay of external and internal mass transfer limitations with exothermic reaction may also give the rather complicated and nontrivial behavior of the oxidative regeneration process.
These models are rather detailed, but they all consider the catalyst pellets as spheres or propose the reduction of other shapes to pseudospheres. Such an approach may seriously worsen the simulation adequacy and accuracy in cases of nonspherical particles (cylinders, rings, multilobes, etc.).
The article [
8] contains a coked naphtha reforming catalyst (Pt/Re-Al2O3) regeneration model. Modeling is carried out at the level of individual particles and a fixed catalyst bed, the rate of coke burnout is investigated with varying some characteristics of the catalyst. This work is related to cylinder-shaped pellets, though using the simplified model assuming the absence of axial coke and oxygen gradients in the pellet, which is correct only for cylinders with some fixed diameter and infinite length; this model is not applicable to real shapes, where the diameter is commensurable with length.
Therefore, the correct simulation of catalyst pellets with nonsphere shapes is an important problem, definitely lacking the corresponding research, especially in the case of transient processes such as oxidative regeneration.
A common problem that arises in models that combine heat and mass transfer, chemical reactions and diffusion is the complexity of systems of equations describing these models. Well-known computational methods and algorithms are often ineffective for solving such problems and do not allow numerous computational experiments to be carried out quickly. A possible solution seems to be the reduction of complex nonlinear models to simpler ones that reflect some of the main characteristics [
9].
The examples of such an approach in modeling the oxidation of organic substances is the model of an unresponsive core and the model of progressive transformation.
The unreacted core model is used for modeling of nonporous environments. If the catalyst grain is nonporous, then reactions can only take place on its surface, which leads to a significant simplification of the model [
10]. This model is used for layer-by-layer burnout of the reagent with a decrease in the reaction zone. It takes into account kinetics and diffusion. However, many researchers note that adequate modeling of some processes requires consideration of flows [
11]. Moreover, this model is characterized by the presence of a movable boundary, which adequately describes the process of particle reduction over time [
12] but does not correspond to the process of oxidative regeneration.
The progressive conversion model can be used in the case of a high level of porosity of the catalyst grain [
13]. This model adequately describes the processes accompanied by rapid diffusion of substances into the grain, however, it does not take into account the change in the reaction surface area that inevitably occurs during coke burning. Moreover, a necessary condition for the application is the absence of the influence of solid concentrations on the chemical reactions rate, which is impossible in the presence of sorption processes, as well as the assumption of quasi-stationarity of the process [
14].
Previously, the authors obtained a mathematical model for burning coke sediments from a catalyst spherical grain [
15,
16], developed on the basis of the well-known model [
17,
18]. It is described by a system of parabolic equations corresponding to the thermal, diffusion and chemical processes accompanying oxidative regeneration.
The present model uses a similar multistage kinetic scheme that takes into account the gas and solid phases of the reaction. Regeneration of the catalyst based on
is considered. Kinetic characteristics are taken from the work [
17].
Effective coefficients of diffusion, thermal conductivity and heat capacity, depending on temperature, are used to calculate the main characteristics of the regeneration process. The dependencies of heat capacity and thermal conductivity are approximated by polynomials based on experimental data, and the diffusion coefficient is obtained based on the Knudsen diffusion coefficient.
The grain geometry makes its own adjustments to the dynamics of the main processes’ characteristics, despite the uniform nature of the processes occurring in spherical and cylindrical grains. The material balance equations are supplemented with a convective term corresponding to mass transfer along the axis of the cylinder to account for the processes in a cylindrical grain. Similar changes have been made to thermal conductivity—it is necessary to take into account its propagation along the cylindrical axis in addition to the propagation of heat along the radius.
The calculation of the Stefan flow velocity has been corrected due to the complexity of the diffusion processes description, which in turn leads to the complexity of its discrete analogue. The flow in the spatial cell with the number is integrated using the values of grain concentrations and temperature in the cell , despite that only the first-order derivative is used in the model. In addition, the transition from a one-dimensional problem in the case of central symmetry to a two-dimensional axisymmetric problem leads to a significant increase in computational complexity for refinement of the spatial grid. All this leads to the use of parallel computing as an effective means of reducing time costs.
The aim of this work is to construct a mathematical model of the oxidative regeneration of solid cylindrical catalyst grains, as well as to develop and implement an effective numerical algorithm for calculating concentrations of substances involved in chemical reactions and the temperature of the catalyst grain.
2. Development of a Mathematical Model of Processes in a Cylindrical Catalyst Grain
The most important part of the mathematical description of any process involving chemical reactions is their kinetic model. This paper presents a description of the burning of coke sediments, presented in the works [
17,
18,
19,
20], as a kinetic scheme:
Here,
,
—rates of chemical interaction stages, dimension
,
—mole/(m
s),
and
—g/(m
c);
—catalyst grain temperature, K;
,
—rate constants of chemical interaction stages, the dimension of
coincides with the dimension of
;
,
—the degree of covering of the coke surface with various carbon complexes (
—hydrogen-carbon complex,
—oxygen–carbon complex and
—free carbon surface);
—oxygen concentration in the gas phase in molar fractions;
and
—concentrations of hydrogen and oxygen in the coke layer in mass fractions;
and
—the amount of hydrogen and oxygen adsorbed by coke relative to the current state of the surface of coke sediments;
and
are the density (g/m
) and the average radius of the pellets (m) of coke. Meanwhile,
The first reaction in the scheme (
1) is a multistage dissociative adsorption of oxygen, leading to the formation of an oxygen–carbon complex on the surface of the catalyst grain. The second stage of the scheme (
1) is the oxidation of the oxygen–carbon complex with the release of carbon dioxide. The third stage describes the destruction of the carbon–oxygen complex with the release of a free carbon surface and the release of carbon monoxide [
19]. The fourth reaction describes the multistage oxidation of a hydrogen–carbon complex with the formation of an oxygen–carbon complex and water vapor. The fifth stage corresponds to the recombination of the oxygen–carbon complex with the release of a free carbon surface. The concomitant release of carbon dioxide is the result of oxygen desorption, which, without creating special conditions, is possible only with reaction products [
21].
The determination of the rates of the first five stages of oxidative regeneration is subject to the law of conservation of active masses. This assumption allows us to obtain expressions for reaction rates consistent with experimental data.
The sixth and seventh reactions in the (
1) scheme are heterogeneous processes that lead to the release of the free carbon surface of the coke pellet. Taking into account the density and radius is necessary due to the fact that the coke structure, and hence its behavior in the oxidation reaction, depend on the density and size of the coke granules. The expressions for the reaction rates
and
are derived from empirical data. The reaction scheme and its detailed description are given in [
17].
The constants
are given in the works [
17] as the constants of reaction rates at temperature
K,
. In this case, it is necessary to take into account the change in the rate constants of reactions with temperature changes, expressed from the Arrhenius equation:
Here, —dimensionless catalyst grain temperature obtained by dividing the current temperature by ; —activation energies, J/mole; J/mole/K.
It should be noted that, as coke sediments are burned during oxidative regeneration, the radius and area of coke pellets decrease. The following dependencies are given in [
17]:
where
and
are the current and initial mass fraction of coke sediments on the catalyst grain, respectively;
is the initial average radius of coke granules on the catalyst grain, m; and
—the initial specific surface area of coke sediments, which plays the role of the reaction surface in the burning process, m
/g.
The flow of catalytic processes in the catalyst grain is classically described by the equations of thermal and material balance. The developed model is called diffusion, and it is based on the laws of conservation of mass and energy.
The reduction of the three-dimensional formulation of the problem to its two-dimensional analogue, the axisymmetric problem, has proven itself well. This approach assumes a transition from a three-dimensional Cartesian rectangular coordinate system to a cylindrical coordinate system with two spatial variables: r corresponds to a cylindrical radius, and z corresponds to the semiaxis of the cylinder.
The thermal balance implies taking into account the distribution of heat through the catalyst grain and the thermal effect of the reaction. It can be described using the thermal conductivity equation, while the role of an inhomogeneous term in its composition is played by the heating of the grain due to the release of heat as a result of exothermic reactions. The equation of thermal conductivity will have the form:
where
is the effective heat capacity coefficient of the catalyst, J/m
/K;
—effective thermal conductivity coefficient of the catalyst, W/m/K;
—bulk density of the catalyst, g/m
; and
,
—thermal effect of the
j-th reaction stage in the (
1) scheme.
The effective heat capacity coefficient calculated by the formula [
22] is taken as
:
where
is the density of the catalyst material;
,
,
and
,
,
are the coefficients of the polynomial temperature dependence of the specific heat capacity of the catalyst material and coke, respectively; and
is the porosity of the catalyst grain.
The effective coefficient of thermal conductivity
is calculated by the formula [
23]:
where
,
,
are the coefficients of the polynomial temperature dependence of the thermal conductivity of the catalyst material.
The source term
for the thermal conductivity equation is derived from accounting for chemical transformations on grain and takes into account its heating due to the exothermic nature of reactions in the (
1) scheme.
The material balance of the catalyst grain is described on the basis of Fick’s law:
where
,
—concentrations of substances in grain pores in molar fractions;
—effective diffusion coefficient, m
/s;
—the velocity of the Stefan flow, m/s;
—molar gas density, mole/m
;
,
,
—stoichiometric coefficients of substances.
The detailed kinetic scheme assumes heterogeneous reactions occurring with a change in the reaction volume, which leads to the loss of some volume of the reaction gas phase. It follows from the law of mass conservation that the Stefan flow plays an important role in the regeneration process in addition to the diffusion flow. The calculation of its velocity is implemented by summing the equations of the material balance and equating to zero the changes in the sum of the substance concentrations of the gas phase in time and space. It is necessary to take into account two components of the Stefan flow—along the cylindrical radius and along the axis of the cylinder—since the assumption of a uniform grain structure is accepted in the work. The velocity of the Stefan flow does not depend on time, but it depends on space and is a scalar quantity:
It is necessary to take into account the consumption and formation of substances due to chemical reactions in addition to the diffusion and Stefan flow—the source term in the diffusion equation.
The material balance also needs to be supplemented with equations describing changes in the substances’ concentrations in the solid phase of the reaction:
where
is the molar mass of coke, g/mole.
Various types of cylindrical catalysts are used in industrial processes, which can be conditionally divided into two large groups: whole and hollow cylinders. The processes occurring in whole and hollow cylinders are the same, however, in terms of mathematical physics, boundary conditions for the mathematical model play an important role, implying for this task the heat and mass transfer of grain with the environment.
In the case of a whole cylindrical catalyst grain, heat and mass transfer occurs through the cylindrical wall and the base of the cylinder. Initial boundary conditions have the form:
Thus, the mathematical model of oxidative regeneration of the cylindrical catalyst grain has the form:
with initial boundary conditions (
11).
7. Conclusions
The article develops the mathematical model of the oxidative regeneration of a cylindrical catalyst grain using a diffusion approach—diffusion, Stefan flow, thermal conductivity of the grain and chemical reactions are taken into account.
The model is presented in an axisymmetric formulation. The calculation of the main characteristics was conducted for whole catalyst grains. This corresponds by the boundary conditions: for a whole grain, heat and mass transfer occurs through the cylindrical wall and the ends of the cylinder.
The calculation of the thermal balance of the catalyst grain was carried out using the effective coefficients of heat capacity and thermal conductivity of the catalyst grain, taking into account its porosity. The dependence of these coefficients on temperature is restored on the basis of experimental data.
The efficient computational algorithm is constructed, taking into account the very different times of diffusion processes and chemical reactions. Splitting by physical processes with sequential calculation of chemical reactions, diffusion and thermal conductivity were used to reduce the complexity of the algorithm. The resulting difference scheme formed the basis of the code that allows calculations of the main characteristics of the process.
Technological parameters corresponding to the actual process of burning coke sediments are used for the computational experiment.
The result of the program is the distribution patterns of substances and temperature by the catalyst grain.