Dynamic Model Study and Analysis of DME Auto-Thermal Steam Reforming Reaction

Through the study of the kinetics of dimethyl ether steam reforming reaction, a two-dimensional model of the reactor is established. At the same time, rate equations of simplified elementary reactions of DME steam reforming reaction are deduced by the mechanism of Langmuir-Hinshelwood, and the rate constants are obtained by correcting the pre-exponential factor of Arrhenius equation. Finally, the analog data of DME steam reforming reaction under a variety of conditions are obtained, and making a simulation diagram. The two-dimensional model is substantially correct because of the analog data more in line with a variety of knowledge.


Introduction
Dimethyl ether has the advantage of economy, security, and the sources of diversity, bound from the research stage to the industrial stage [1].However, due to the engine of DME has a shortcoming of fire earlier, making it smaller work load range [2], while the hydrogen fuel cell and electric vehicles are being tested and toward industrialization in the world, but blocked in hydrogen storage technology [3], so automotive fuel reforming is an important way to solve the problem of the hydrogen energy [4] [5].
DME is an ideal vehicle fuel, but also has capability of chemical hydrogen storage.Therefore, many scholars have been studied to the DME reforming processes [6]- [8], confirming its value of commerce and scientific and finding some technical issues, such as the hydrogen yield.Generally speaking, the main method of DME reforming is the traditional methods and the plasma reforming method.The traditional methods include Auto-thermal reforming, steam reforming, partial oxidation reforming and so on.However, the study of kinetics model about DME reforming reaction is still relatively lacking in the current.Weiqi Qian et al. [9] used the quasi-steady method to establish a simplified chemical kinetic model for combustion of hydrocarbon fuels, effectively reducing the reaction components.Solution variables and differential equations are reduced when it is calculated in the flow field, improving the efficiency of calculation.Yumin Chen et al. [10] conducted kinetic analysis to the process of methane auto-thermal reforming reaction.Wang Feng et al. [11] conducted kinetic analysis to the process of methanol steam reforming in the micro-reactor.Dongmei Feng et al. [12] made the Kinetic Study of the process about dimethyl ether steam reforming reaction using of a catalyst CuO − ZnO − Al 2 O 3 + ZSM − 5.This article will use the catalyst

Experimental Methods
During the whole experiment, the space velocity of DME steam reforming reactor is 0.3 m/s.The molar ratio of water and DME is 1 to 3. The feed mole fraction of dimethyl ether is between 5 and 15.The flow rate of DME is

Reforming Reactor Model
The kinetic model of DME steam reforming reactor applies to a two-dimensional model.Although the kinetic parameters of the reactor are difficult to be determined in the process of reaction, it can improve the accuracy of the change of amount of substance and energy conversion.The two-dimensional model is non-continuous, axial symmetry, similar to the two-dimensional reactor of Groppi, G. et al. [13].The reaction of DME and H 2 O on the catalyst coating layer of the reactor produced H 2 and oxides of carbon.In the course of the reaction, the mass of reactants and product are equal, the mass conservation equation as follows: ( ) Because the reactants are gas and in the high-speed motion in the reactor, and so axial diffusion of the gas mass is negligible.Because the heat exchanger conduit wall can isolate gas, so the radial diffusion of the gas mass in the reactor is also negligible.Then the above equation can be simplified to as follows: ( ) Further, on the catalyst coating layer of the reactor, the mass conservation equation of each element in the course of the reaction and diffusion is as follows: ( ) In this reactor, if the catalyst coating layer is a uniformly thickness, the boundary conditions are shown in Table 1.The average thickness of the catalyst coating layer is about 35 microns, while the effective diffusion coefficient e D ( ) m s − ⋅ is estimated by random pore model [14].On the catalyst coating layer, according to the mass conservation of the elements, the average reaction rate js r can be obtained: Similar to the gas phase mass balance of (2), axial and radial dispersions were neglected in the gas-phase heat balance: Table 1.Boundary conditions for the 2D monolithic reactor wash-coat mass balance and solid heat balance.

Position
Boundary conditions External energy transfers to DME and water using of the heat conduction manner of the axial and radial through the reactor wall and the catalyst coating layer, using the equation: Wherein, t m -The total mass flow rate ( ) , w-Reactor wall, s-Catalyst coating surface, w-Mass fraction.Table 1 shows the boundary conditions of conservation of energy and the quality of the reactor model catalyst coating.

Kinetics Model of DME Steam Reforming Reaction
There are three major reactions in the process of DME steam reforming [15], which can be shown in Table 2, and DME oxidation reaction is the heat source.In the ideal case, the inlet of reactor is water, DME and O 2 , and the outlet is CO 2 and H 2 .When each elementary reaction is in a relatively stable state of the reaction rate, adsorption equilibrium constant of H 2 and CO 2 on the catalyst is relatively small, and it can be ignored as the case in the calculation process.The reaction rate equation takes Langmuir-Hinshelwood model of water molecules and other molecules competitive adsorption, and it is ideal for surface reaction.

Parameters of the Kinetic Model
Arrhenius equation is used to solve the rate constants, and the expression is: e Due to the range of experimental temperature is 200˚C to 400˚C, so temperature variation range is small, and the pre-exponential factor A is corrected and the activation energy not, then Arrhenius equation: e where, j k ( ) is the rate constant at temperature T.
Then, the activation energy A E ( ) 1 J mol − ⋅ can be considered as it of the initial temperature j s T (200˚C) and unchanged in the whole calculation process.The pre-exponential factor j s A is calculated as follows: The fitting parameter of Adsorption equilibrium constant ( ) j K comprises of mj K of average temperature (300˚C) and adsorption enthalpy ( ) adsj H ∆ , and its expression is as follows: The main elementary reaction rate constant can be obtained by Equation ( 8) and ( 9) using ofrate constant j s k and activation energy j A E in the Table 3 [16].While adsorption equilibrium constant ( ) j K of any tempera- ture can be known by formula (10)

Results and Discussion
The legend of all graphics is the same with the legend of Figure 1(a).

Conclusion
The two-dimensional model is substantially correct because of the analog data more in line with a variety of knowledge.For example, the ratio of water-DME and the ratio of oxygen-DME are important step of controlling the hydrogen yield.Finally, this article can provide important theoretical references to the kinetics of dimethyl ether steam reforming reaction.
studies of Dimethyl ether steam reforming reaction, establishing a two-dimensional model of the reactor, deducing the rate equation of simplified elementary reaction, learning the influence of the reaction rate and the mole fraction with temperature and time.

−
in the feed mouth.The pressure inside the reactor is 0.1 Mpa and the range of temperature is 200˚C to 400˚C.
⋅and H 2 to temperature is the maximum, indicating that their speed of response here has a jump, which can be proven from Figure1(b), Figure 1(c) and Figure 1(d).When the curve of Figure 2(a) and Figure 2(b) trend to steady, the curve of Figure 2(c) and Figure 2(d) are becoming more stable, indicating that the reactants at the inlet of the reactor achieves to conditions optimally.That is, if this state is changed, the reactants conditions of the reactor at the inlet need to change.

Table 3 .
Finally, the rate formula of major elementary reactions in Table2can be computed for getting a variety of data.

Table 2 .
Kinetic model overall reactions and rate expression.

Table 3 .
Fitted kinetic parameter values and 95% confidence intervals.