In order to study the steam reforming process of dimethyl ether, a kind of reformer reactor with thermal jacket, heat pipe and catalytic reaction bed was designed. The effects of reaction gas temperature, molar r atio of water to ether and the structure parameters of the reactor on the conversion of dimethyl ether, the yield of hydrogen, the hydrogen concentration at the outlet of the reactor and the concentration of CO 2 were investigated experimentally. The mathematics and power of the reactor were established Learn model The COMSOL software was used to simulate it. The simulation results showed the temperature distribution, mass distribution, DME conversion and hydrogen production in DME steam reforming process. These simulation results will provide useful data for the design and operation of small scale catalytic reforming bed reactors.
With the rapid development of the world economy and the continuous improvement of the appointed living standard, the consumption of primary energy has rapidly risen. Oil is still the world’s leading fuel. The consumption of oil will inevitably bring many new challenges to the growth of the world economy, energy security and climate change. Therefore, the search for a new clean alternative energy is very necessary. As discussed by Chaubey et al. [
As discussed by Alves et al. [
Dimethyl ether as a reforming hydrogen source has many advantages, mainly in the following four points: 1) a wide range of sources of dimethyl ether, a new synthesis of dimethyl ether to ensure the economy of the source, is conducive to large-scale application; 2) dimethyl ether is non-toxic, non-corrosive, non-trigeminy (teratogenic, carcinogenic, mutagenic) and has no damage to the atmospheric ozone layer and is easily degraded in the troposphere, thus ensuring the utilization of dimethyl ether; 3) the physical properties of DME are similar to those of liquefied petroleum gas and natural gas, and the infrastructure can be compatible with it. As a result, the commercialization cost of dimethyl ether reforming hydrogen production is greatly reduced, meanwhile, it is also reduced This new fuel into the energy market threshold as discussed by elsewhere [
In this paper, the kinetics of the reaction rate of dimethyl ether reforming to hydrogen was studied. A monolithic channel reactor model was designed and its reaction process and operating parameters were simulated by COMSOL software. Water ether ratio and different reaction temperature on its hydrogen yield.
150 mm and the minimum element size is 7.5 mm. The gas phase in the wall is described by the non-slip boundary conditions. The reactants flow in from the reactor inlet and carry out the reforming and hydrogen production of dimethyl ether in a packed bed of a catalyst porous catalytic reaction bed.
The numerical simulation of DME reforming in COMSOL software is solved by the energy conservation control equation. The energy equation describes the temperature of the reaction gas in the reactor, as well as the thermal conductivity of the entire structure. Since temperature affects not only the reaction kinetics, but also the density and viscosity of the reactant gases, the energy equation is the reason that the heat pipes in the reactor structure are actually connected into a three-dimensional model.
In this reactor, the composition of the mixed gas is mainly composed of DME, H2O, CO2, H2 and CO. In addition, the volume flow is in the axial direction, whereas the mass transfer occurs predominantly in the lateral direction of the reactor wall. Therefore, the three-dimensional model is sufficient to deal with the reforming reaction.
The basic assumptions for this model are as follows:
1) the reactant gas forms a perfectly hypothetical laminar flow in the channel such that the average flow field proportionally crosses the pressure differential across the reactor;
2) the reaction rate is linear;
3) the gas flow in the direction of the channel to transmit the quality and energy only;
4) the reaction gas is incompressible, the ideal gas;
5) the reaction process is divided into three steps: dimethyl ether hydrolysis, methanol steam reforming and water gas shift reaction.
Assuming steady state, the mass balance equation for a plug flow reactor is:
d F i d V = R i (1)
where Fi is the species molar flow rate (mol/s), V represents the reactor volume (m3), and is Ri the species net reaction rate (mol/(m3・s)). The molar flow rate is related to the species concentrations, ci (mol/m3), through the volumetric flow rate, V (m3/s):
F i = v c i (2)
where the volumetric flow rate is given by the average flow velocity, u (SI unit: m/s), multiplied by the reactor cross-section A (m2):
v = u A (3)
The energy balance for the ideal reacting gas is:
∑ i F i C p , i d T d V = Q e x t + Q (4)
where C p , i is the species molar heat capacity (J/(mol・K)), and Q e x t is the heat added to the system per unit volume (J/(m3・s)). Q denotes the heat due to chemical reaction (J/(m3・s)).
Q = − ∑ i H j r j (5)
where Hj the heat of reaction (J/mol), and rj the reaction rate (mol/ (m3・s)).
Species concentrations are defined at the reactor inlet boundaries:
c = c i n (6)
At the outlet, use the Outflow condition:
n ⋅ ( − D ∇ c ) = 0 (7)
For the heat pipe in the reactor, heat can only be conducted by conduction:
− ∇ ⋅ ( k s ∇ T ) = 0 (8)
where ks (W/(m・K)) denotes the thermal conductivity of the heat pipe.
The temperature is specified at the reactor inlet boundaries:
T = T 0 (9)
At the outlet, use the Outflow condition:
n ⋅ ( k ∇ T ) = 0 (10)
DME steam reforming process is more complicated, this paper mainly studied three dimethyl ether hydrolysis, methanol steam reforming and water gas shift three reactions as discussed by Feng et al. [
DME hydrolysis reaction:
CH 3 OCH 3 + H 2 O ( g ) ⇔ 2 CH 3 OH Δ H = 36.6 kJ / mol (11)
MeOH steam reforming reaction:
CH 3 OH + H 2 O ⇔ CO 2 + 3 H 2 Δ H = 49.1 kJ / mol (12)
Water gas shift reaction:
CO + H 2 O ⇔ CO 2 + H 2 Δ H = − 47.17 kJ / mol (13)
The total reaction is:
CH 3 OCH 3 + 3 H 2 O ( g ) ⇔ 2 CO 2 + 6 H 2 Δ H = + 135 kJ / mol (14)
For the above reaction of dimethyl ether reformer hydrogenation, the reaction kinetic rate equation can draw as discussed by Oar-Arteta et al. [
r H = K H K M D , a K W , a ( P D P W − P M 2 / K H ) ( 1 + K M D , a ( P M + P D ) + K W , a P W ) 2 (15)
r S R = K S R K M , m K W , m ( P D P W − ( P C O 2 P H 2 3 / K S R ) ) 1 + K M , m P M + K W , m P W (16)
r r W G S = K r W G S ( P C O 2 P H 2 − ( P C O P W / K r W G S ) ) 1 + K M , m P M + K W , m P W (17)
In the above three formulas, K M D , a , K M , m respectively, the adsorption constants of the oxidant (MeOH + DME) at the acidic sites and the methanol sites at the metal sites, K W , a , K W , m represent the adsorption constants of water at the acidic site and the metal site, P represent the partial pressure of the chemical substances, K H , K S R , K r W G S represent are DME hydrolysis, methanol steam reforming and water gas shift reaction equilibrium constant.
Selected kinetic models have been used to model the reactor. Simulation allows the determination of a suitable range of operating conditions (temperature, reactor length and S/DME ratio) resulting in high H2 and low CO yields. The calculated reaction equation is as follows:
Dimethyl ether conversion:
X D M E = F D M E , i n − F D M E , o u t F D M E , i n × 100 % (18)
Yields of Hydrogen:
Y = F i F 0 ⋅ v i (19)
where F D M E , i n , F D M E , o u t denote the molar flow rate of DME at the inlet and outlet of the reactor and F i denote the molar flow rate of each product (H2, CO2) at the outlet of the reactor, which v i is the stoichiometric coefficient of formation of component i from methanol, CO2 is 2, H2 is 6.
parameter | K (623 K) | activation energy |
---|---|---|
K H | 3.26 (±1.31) (gcatalyst h atm)−1 | 145.1( ±33.7 ) KJ/mol |
K S R | 15.2 (±4.30) (gcatalyst h atm)−1 | 28.9( ±8.41 ) KJ/mol |
K r W G S | 20.1( ±6.7 ) (gcatalyst h atm)−1 | 63.3( ±18.4 ) KJ/mol |
K M D , a | 1.67( ±0.42 )atm−1 | 24.5( ±8.11 ) KJ/mol |
K M , m | 31.46( ±0.37 ) atm−1 | 40.2( ±10.7 ) KJ/mol |
K W , a | 6.15( ±1.89 ) atm−1 | 55.8( ±18.4 ) KJ/mol |
K W , m | 5.81( ±2.04 ) atm−1 | 33.7( ±8.11 ) KJ/mol |
The reaction conditions for the water ether molar ratio of 3,
temperature changes. Change the temperature values were 553 K, 573 K, 593 K, 613 K. It can be seen from
The molar ratio of steam/DME is one of the most important parameters that influence the steam reforming process.
water ether ether of 3, 3.5, 4, 4.5, 5, respectively, with the increasing molar ratio of water to ether and the temperature increases, dimethyl ether reforming hydrogen production reaction The reaction rate is also constantly accelerating. Similarly, we can also know that the hydrogen yield will be more and more as the reaction proceeds. Therefore, it can be inferred that excessive water and high temperature during steam reforming of DME may promote the reaction. This can also be called the concentration effect, that is, DME is a more expensive material than water vapor. In view of the reversible reaction rate and reaction conversion rate limited by the balance, in order to improve the DME utilization rate and to speed up the reaction rate, Very effective. However, it does not mean that the higher the water-to-water ratio, the better. Due to the relatively large heat capacity of water, excess water needs to consume more heat in the DME reforming hydrogen production reaction. Since the entire reforming reaction is endothermic, Higher water ether ratio is not conducive to DME conversion.
Steam and DME enter the gas mixer through the mass flow controller. The mixed gas flows into the reactor containing the CuO/ZnO/Al2O3 + ZSM5 catalyst. Product gas from the reactor flows into a dryer to separate the steam from the gas mixture. Finally, a portion of the gas product flows into the gas chromatograph (GC-9900) to detect the gas composition. Gas chromatograph connected to the computer, the composition of the product can be analyzed by the recorder. Thermocouples were used to measure the catalyst bed temperature in the reactor.
The experimental conditions in the reactor are exactly the same as those in the simulation. Inlet temperature is 553 K, exhaust temperature is 753 K. It is clear from
In this paper, a dimethyl ether steam reforming hydrogen production reactor was established by using COMSOL software, and the reactor performance was simulated with a suitable catalyst. The effects of different temperatures, molar ratios of water to ether and reactor geometry parameters impact on the reforming reactor. Simulation model can describe the whole reaction process well, the conclusion is as follows:
1) Higher water-ether ratio can promote the forward reaction and improve the conversion of dimethyl ether.
2) As the steam reforming reaction is generally endothermic reaction, higher
inlet temperature can provide more energy required for the reaction, thereby speeding up the reaction and increasing the hydrogen yield.
3) The mathematical model and physical model of the reactor can predict the performance of the reactor to a certain extent, and can improve the design of the catalytic reforming microreactor with certain reference value.
This work was supported by National Natural Science Foundation of China under Grant NO.51505275.
Guo, L. and Li, C. (2018) Modeling and Numerical Simulation of Hydrogen Production from Dimethyl Ether Steam Reforming. Open Access Library Journal, 5: e4531. https://doi.org/10.4236/oalib.1104531