Next Article in Journal
A Chaotic Electromagnetic Field Optimization Algorithm Based on Fuzzy Entropy for Multilevel Thresholding Color Image Segmentation
Previous Article in Journal
Explicit Lump Solitary Wave of Certain Interesting (3+1)-Dimensional Waves in Physics via Some Recent Traveling Wave Methods
Previous Article in Special Issue
Evaluating the Transient Energy Dissipation in a Centrifugal Impeller under Rotor-Stator Interaction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model Development and Exergy Analysis of a Microreactor for the Steam Methane Reforming Process in a CFD Environment

1
Department of Chemical Engineering, National University of Sciences and Technology (NUST), H-12 Islamabad 46000, Pakistan
2
Department of Systems Science, Graduate School of Informatics, Kyoto University, Yoshida-Honmachi, Sakyo-Ku, Kyoto 606-8501, Japan
3
Department of Chemical and Petroleum Engineering, College of Engineering, United Arab Emirates University, Al-Ain 73000, UAE
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(4), 399; https://doi.org/10.3390/e21040399
Submission received: 18 February 2019 / Revised: 4 April 2019 / Accepted: 8 April 2019 / Published: 15 April 2019
(This article belongs to the Special Issue Entropy in Computational Fluid Dynamics II )

Abstract

:
Steam methane reforming (SMR) is a dominant technology for hydrogen production. For the highly energy-efficient operation, robust energy analysis is crucial. In particular, exergy analysis has received the attention of researchers due to its advantage over the conventional energy analysis. In this work, an exergy analysis based on the computational fluid dynamics (CFD)-based method was applied to a monolith microreactor of SMR. Initially, a CFD model of SMR was developed using literature data. Then, the design and operating conditions of the microreactor were optimized based on the developed CFD model to achieve higher conversion efficiency and shorter length. Exergy analysis of the optimized microreactor was performed using the custom field function (CFF) integrated with the CFD environment. The optimized catalytic monolith microreactor of SMR achieved higher conversion efficiency at a smaller consumption of energy, catalyst, and material of construction than the reactor reported in the literature. The exergy analysis algorithm helped in evaluating length-wise profiles of all three types of exergy, namely, physical exergy, chemical exergy, and mixing exergy, in the microreactor.

1. Introduction

The demand for hydrogen production is rapidly increasing due to its immense importance in the petroleum and chemical industry. Hydrogen is mainly used for the up-gradation of fossil fuels and ammonia production. The use of hydrogen as a fuel has also recently increased due to its high heating value and less-polluting nature [1]. There are a variety of other applications of hydrogen such as fuel cells, metal production, and power generation [2].
The steam reforming of methane gas is the most abundantly used process for the production of hydrogen. Saeidi et al. estimated that 48% of all hydrogen production came from the steam reforming of natural gas [3]. However, steam methane reforming (SMR) has several limitations, such as heat transfer limitation, mass transfer limitation, and catalyst deactivation [1]. To cope with these limitations, conventional SMR is being transformed into micro-reforming technology, which overcomes the heat transfer limitation and enhances the mass transfer; as a result, micro-reforming technology improves the conversion efficiency [4]. The intensified form of SMR (i.e., micro size plant) increases profit margins above 70% for the same capacity of H2 production [1].
Granlund et al. revealed that a microreactor having multiple air inlets produces a high yield of hydrogen as compared to the conventional monolith reactor by using computational fluid dynamics (CFD) [5]. Bhat et al. presented a review on various research directions in the area of process intensification of the SMR process [6]. Chen et al. developed the CFD-based model of a catalytic microcombustor and evaluated the combustion characteristics and stability of a methane–air mixture [7]. Yu et al. developed a microreactor for the steam reforming of methanol and analyzed the performance of the microreactor with optimum catalyst coating under different operating conditions [8]. Theampetch, A. et al. produced a hybrid system of SMR which transformed some part of methane to hydrogen by using the heat of exhaust gases of an internal combustion engine, and they showed through CFD simulation that the homogeneity of temperature in the longitudinal and radial directions was indispensable for high methane conversion efficiency [9]. Kansha, Y. et al. compared the energy requirement for self-heat recuperative systems and conventional thermal processes. They concluded that the recuperative thermal processes had a higher capability to save energy than the conventional thermal processes in the industries [10]. Hosseini and Wahid investigated the conventional and flameless combustion in a lab-scale furnace based on exergy analysis and concluded that the major irreversibilities were caused by the high temperature gradient present in the reactor chamber [11].
An et al. investigated different microreactor structures, such as parallel, oblique pin, pinhole, wavy, and coil, with CFD simulation to achieve a high-performance configuration with respect to heat transfer, reaction rates, and their flow characteristics. They concluded that the pinhole configuration could achieve better performance than the other structures [12]. Kashid et al. developed the CFD model of a slug flow microreactor to examine the impact of viscosity on the fluid flow streamlines within the slugs. They summarized that the variation in viscosity had no effect on the flow patterns within the slugs [13].
Since SMR is endothermic and energy intensive, a robust energy analysis of the micro-scale SMR is also needed [14]. The concept of exergy has attracted the attention of researchers and process designers because of its capability to evaluate the true thermodynamic potential of a process. Exergy analysis helps in the identification of the locations, causes, and true magnitude of irreversibility and losses in a system [15]. Exergy analysis of the process models in the CFD environment has been reported in the literature.
Boulenouar and Ouadha performed a CFD-based exergy analysis of flow in a supersonic steam ejector and found that the main source of irreversibility was the nozzle due to high levels of pressure and velocity gradients [16]. Mustafa et al. developed a CFD-based model for exergy analysis of naphtha reforming reactors and concluded that the total exergy of the stream increased along the reactor [17]. Alabi et al. compared experimental and CFD-based exergy analysis procedures for the design optimization of the airframe subsystem of aircraft. They observed that exergy-based analysis has a significant advantage over an energy-based approach [18]. Gunjo et al. developed a CFD model for a flat plate solar collector to perform exergy and energy analysis. It was observed that the performance of the collector was affected by inlet water temperature, ambient temperature, solar insolation, mass flow rate, working fluids, insulation thickness, and collector heat loss factor [19].
Debnath et al. performed CFD-based exergy efficiency analysis of air hydrogen detonation in a pulse detonation combustor. They concluded that the deflagration combustion process had higher exergy losses than the detonation combustion process [20]. Erguvan et al. analyzed energy and exergy of unsteady cross-flow overheated circular cylinders in the CFD environment. It was found that exergy efficiency could be increased by selecting specific pitch ratios for different Reynolds numbers [21]. Alabi et al. reported the use of a CFD-based exergy analysis procedure for an aircraft system. They presented different exergy calculations and compared empirical models of the airframe subsystem to CFD models [22]. Yong-an et al. performed an exergy analysis of exhaust gas in a kitchen burning liquefied gas [23]. Farmahini-Farahani investigated different geometrical parameters for the thermal stratification of cold water tanks by CFD-based exergy analysis and concluded that decreasing the diameter of inlet/outlet reduced initial mixing within the tank and thus thermal stratification happened sooner, leading to higher exergy in the tank [24].
In this work, exergy analysis of a monolith microreactor of SMR was performed by using a computational fluid dynamics (CFD) model. Initially, the CFD model was developed using literature data. The monolith catalyst bed was used by considering its large surface area, low-pressure drop, high durability, and high mechanical strength. The surface-based approach was used for modeling reaction rates. To incorporate reaction kinetics, the CHEMKIN mechanism was used. Size of the reactor and operation conditions were optimized to realize higher conversion efficiency and smaller capital and operational costs. To analyze true thermodynamic efficiency, its exergy analysis was performed by developing a custom field function (CFF)-based algorithm. The exergy analysis algorithm helps in evaluating length-wise profiles of all three types of exergy, namely, physical exergy, chemical exergy, and mixing exergy, in the microreactor.
In Section 2, the process description and concepts of exergy are discussed, followed by a description of the model development in Section 3. Section 4 describes the results and discussion, while Section 5 concludes the work.

2. Process and Exergy

In this section, the SMR process and exergy concepts are described.

2.1. Process Description

Hydrogen is produced mainly in four steps—natural gas pretreatment, reforming process, water–gas shift reaction, and gas purification, as shown in Figure 1 [25].
In the pretreatment step, the impurities present in the natural gas feed are separated to prevent the catalyst from poisoning. The zinc oxide bed captures the sulfur-containing compounds such as hydrogen sulfide and leaves the gas with sulfur content less than 1 ppm at 335 °C [26].
In the reforming process, the natural gas–steam mixture is catalytically converted into hydrogen, carbon monoxide, and carbon dioxide in the reformer according to the following reactions:
CH4 + H2O ↔ 3H2 + CO       ∆H°298 = +206 kJ/mol
CH4 + H2O ↔ 4H2 + CO2       ∆H°298 = +165 kJ/mol
The reforming reaction is highly endothermic; therefore, heat is provided through external burners to maintain the temperature and pressure at 1100–1500 °C and 1–5 atm, respectively. The side burners are operated to heat the furnaces which supply heat to the reformer tubes through forced convection and radiation. The reformer furnace normally consists of a number of tubes (i.e., in the range of 40–400) determined by the plant design capacity. Different shapes and arrangements of the monolithic support channels are shown in Figure 2 [4].
The third step is the water–gas shift reaction. The syngas exiting the reformer takes part in a water–gas shift (WGS) reaction. It converts CO and H2O in the syngas to CO2 and H2 using the following reaction:
CO + H2O ↔ CO2 + H2       ∆H°298 = −41 kJ/mol
In the fourth step, which is the purification of hydrogen stream, CO2 is removed through chemical absorption, and unreacted methane is separated by pressure swing adsorption (PSA).
In this work, a two-dimensional CFD model of a single reformer tube of monolith catalytic reactor was designed [27]. The single reformer tube is shown in Figure 3, and the model parameters are mentioned in Table 1.
The reactor is fed with a mixture of methane and water vapor. Inlet velocity is maintained at 0.45 m/s and a temperature of 800 °C. The mole fraction of methane and steam are 0.23 and 0.77, respectively. A heating jacket surrounds the monolithic reactor and provides the heat of reaction to maintain a constant temperature of 1477 °C at reactor walls. The inlet H2O/CH4 molar ratio is kept constant at 3, and the operating pressure is atmospheric.

2.2. Exergy Concept

The term exergy was first introduced by Zoran Rant, and it means the amount of work (-erg) that is released (ex-) [28]. Exergy is defined as the maximum amount of work obtainable when an energy carrier is brought from its initial state to a state of thermodynamic and chemical equilibrium with the environment, referred to as the dead state. Exergy can be viewed as a measure of energy usefulness or quality [17]. Exergy is consumed during real processes due to irreversibilities [29]. Conventionally, exergy of a stream is classified into three parts: physical exergy, chemical exergy, and mixing exergy. A general expression of exergy is given by the following:
E = Eph + Ech + Emix
where E represents total molar exergy of a stream, Eph is the molar physical exergy, Ech is the molar chemical exergy, and Emix is the molar mixture exergy [28].
Physical exergy represents the thermomechanical portion of the total exergy stream. It is the maximum obtainable amount of work when this stream is brought from actual conditions (T, P) to thermomechanical equilibrium at ambient temperature (To, Po) by reversible processes.
On a molar basis, physical exergy is given by the following:
E ph = R T o i   =   1 n l n p i p o +   i   =   1 n C p i mean ( T i T o T o ln ( T i T o ) )
C p i mean =   T 1 T 2 C p i d T
C p i ( j m o l . K ) =   a i +   b i T +   c i T 2 + d i T 3
where ai, bi, ci, and di are heat capacity coefficients and R is the ideal gas constant [28]. Pi and Ti represent the partial pressure and temperature of individual components, respectively, at each point in the reactor [17].
Chemical exergy is the maximum obtainable work from a material stream by taking it from a state of thermomechanical equilibrium to a state of thermomechanical and chemical equilibrium with the environment [17]. Chemical exergy of a material stream on a molar basis is given by Equation (8):
E ch = i   =   1 n ν i G ¯ i ( Reactants ) i   =   1 n ν i G ¯ i ( Products )
G ¯ i =   G f 0 + [ G ¯ i ( T , P )   G ¯ i ( T o , P o ) ]
where v i is the respective stoichiometric coefficients, G ¯ i is the molar Gibbs function of components i, and G f 0 is the molar Gibbs function of formation at a reference temperature and pressure [30].
Mixing exergy accounts for the mixing effect arising from the isothermal and isobaric mixing of pure components at process conditions [31]. It can be calculated by the following:
E m i x = i   =   1 n x i T o R ln x i ,
where x i is a mass fraction of component i. Mixing exergy is always a negative value because the mixing of different components decreases the exergy continuously along the length of the reactor [28]. It can also be written in the form of i   =   1 n x i R T o ln p i P o   , where Pi is the partial pressure of each component (Pi = χi Ptotal) according to Dalton’s law [32].

3. Model Development

In Section 3.1, details regarding geometry and mesh preparation are provided, while in Section 3.2, boundary and cell zone conditions are described. A set of CFD-based conservative equations are provided in Section 3.3. The kinetic models describing the surface catalytic reactions are presented in Section 3.4, and numerical schemes to solve those equations are given in Section 3.5. The following assumptions are taken for model development:
  • The equilibrium state is reached and the maximum yield is achieved.
  • There is negligible heat loss from the reactor wall to the surroundings.
  • The catalytic wall is isothermal in condition.
  • The flow regime is laminar in the reaction channel, and steady-state operation is reached.
  • The gas mixture is treated as an incompressible ideal gas. The density of the mixture is constant as calculated from ideal gas law.
  • The gas-phase non-catalytic reaction can be neglected so that only the surface reactions are modeled.

3.1. Geometry and Meshing

ANSYS Design Modeler® 16.0 was used to create geometry and mesh. The reactor considered in this study is cylindrical, with a height of 1 mm and length of 14 mm. The reactor consists of an isothermal catalytic (rhodium catalyst) wall as shown in Figure 3 [33]. The simulation was carried out in an axisymmetric mesh; for effective visualization, the axisymmetric mesh was mirrored around its axis. The zoomed version of the mesh is shown in Figure 4. The computational mesh consists of 7200 cells and 7525 nodes. Three important parameters to evaluate the mesh quality are minimum orthogonal quality, maximum ortho-skew, and maximum aspect ratio; in this study, they are 1.0, 0.0, and 1.4, respectively. Table 2 presents the mesh properties.

3.2. Boundary and Cell Zone Conditions

A single tube of monolith reactor is shown in Figure 3. Boundary conditions are defined at the inlet and outlet of the reactor. The inlet boundary condition, that is, the velocity, temperature, pressure, and the composition of inlet gas mixture, in each channel is set at uniform values. The catalytic wall of the reactor tube is kept isothermal. Each reforming tube is tightly packed with rhodium (Rh) catalyst particles. It facilitates the formation of hydrogen fuel from steam and methane through the highly endothermic SMR reactions. It also plays a role as an intermediate medium to enhance the rate of heat transfer to the tube-side gas mixture in the simulations [25]. Considering the Knudsen number of the flow, the flow regime is a continuum, and hence the no-slip condition is used on the microreactor walls. The outflow mixture is discharged to atmospheric pressure. Table 1 shows the model parameters.

3.3. CFD Model Describing the Flow Field

A two-dimensional steady-state model is employed by incorporating it with the detailed homogeneous and heterogeneous reaction schemes in CHEMKIN and Surface-CHEMKIN format. Steady-state continuity, momentum, energy, and species equations in the fluid phase are as follows [34].
Continuity equation
ρ t + ( ρ u x ) x + ( ρ u y ) y =   S m
where ρ , u x , and u y are the density, and the velocity in x and y-directions, respectively. Sm represents a mass addition to the continuous phase, which is zero in this case. The first term of the left-hand side of the equation shows the local derivative which is physically the time rate of change of density at the fixed point. The second and third terms show the convective derivative of the density of the gas mixture [35].
Momentum equations in x and y-directions
( ρ u x ) t + ( ρ u x u x ) x + ( ρ u x u y ) y = P x + x [ 4 3 µ u x x 2 3 µ u y y ] +   y [ µ ( u y x + u x y ) ]
( ρ u y ) t + ( ρ u y u x ) x + ( ρ u y u y ) y = P y + y [ 4 3 µ y y 2 3 µ u x x ] + x [ µ ( u y x + u x y ) ]
where P represents the stream pressure and µ is the viscosity of the gas stream. The first term of the right-hand side of each equation shows pressure forces. The second and third terms of the right-hand side of each equation show the viscous forces [36].
Energy equation
ρ d e d t = P   d i v   u + d i v   ( k   g r a d   T ) + ϕ + S
d i v   u = u x x + u y y
g r a d   T = T x + T y
ι ˙ = C v T
P = ρ R T
ϕ = µ { 2 [ ( u x x ) 2 + ( u y y ) 2 ] + ( u y + v x ) 2 }
where e is the specific internal energy, C v is specific heat constant, and k is the thermal conductivity. ϕ shows the rate of dissipation energy per unit volume and s denotes the work done per unit volume by body forces. The first term of the right-hand side of Equation (14) is the rate of work done per unit volume, and the second term is the rate of heat transfer per unit volume through conduction [37].
Species transport equation
( ρ Y i ) t + ( ρ u x Y i ) x + ( ρ u y Y i ) y = [ J i , x x + J i , y y ] + S ˙ i
J i , x = ρ D i Y i x
J i , y = ρ D i Y i y    
where Y i is a mass fraction, D i is a diffusion coefficient, J i is a mass flux of component i, and S ˙ i is the net production rate of species through chemical reactions [36].

3.4. Kinetic Models Describing the Surface Catalytic Reactions

For the steam reforming of methane over rhodium, the detailed heterogeneous reaction scheme proposed by Karakaya et al. is employed [38]. It involves 44 elementary reactions with 6 gaseous and 13 surface species. The rhodium catalyst having surface site density ΓRh = 2.72 × 10−9 mol/cm2 is used. The total molar production rate of species i-th by surface reactions is given by the following [39]:
S i ˙ = n = 1 n s v i , n k f , r j = 1 N g + N s C j v j , n ,
where S i ˙ is the overall rate expression of i-th species for the gas phase or surface phase, n s is the number of surface elementary reactions, v i , n and v j , n are the stoichiometric coefficients, and N g and N s are the number of gas-phase and surface-phase species, respectively. The concentrations C j of adsorbed species are given in mol/m2.
As binding states of the adsorption of all species vary with the surface coverages, the temperature dependence of the reaction rate coefficients is determined using the modified Arrhenius expression [40]:
K f , r = A k r T β k r e x p ( E a r R T ) θ i μ i r e x p [ ε i r θ i R T ]
where k f , r is a forward rate coefficient, A k r is the pre-exponential factor, β k r is the temperature exponent, E a r is the activation energy of the reaction r, θi is the surface coverage with adsorbed species, and coefficients μ i r and ε i r describe the dependence of the rate coefficients on the surface coverage of i-th species. Sticking coefficients are commonly used for adsorption reactions and are converted to conventional rate coefficients according to [41]:
K f , r a d s = S i o Γ τ R T 2 π M i
where K f ,   r a d s are the reaction rate constant for adsorption reactions, S i o is the initial sticking coefficient, Γ is the surface site density, τ is the number of sites occupied by the adsorbing species, and M i is the molar mass of i-th species.

3.5. Computational Schemes

The governing equations (i.e., momentum, energy, continuity, and species conservation) were discretized by using the finite volume method and were solved numerically by Fluent 16.0®. The second-order upwind scheme was used to discretize the mathematical model. The semi-implicit method for pressure-linked equations (SIMPLE) algorithm was used. Exergy analysis was performed by developing a custom field function (CFF)-based algorithm. An under-relaxation factor was used to slow down the rate of change. The default reference frame was used for velocity initialization. The Mach number, the ratio of flow velocity to the speed of sound, was used to determine whether the flow was compressible or incompressible. The Mach number was less than 0.3 in our model due to the low velocity of the gas stream; consequently, compressibility effects were ignored, and a pressure-based solver was used. Furthermore, the laminar finite rate model was used for calculating the rate of reactions. Computations are very intensive, and the convergence of CFD simulations was evaluated based on the residuals of all governing equations. The governing equations in our CFD model were converged at 2196, the number of iterations, as shown in Figure 5.
A schematic flowchart of the methodology adopted in this study is shown in Figure 6. Model development starts with geometry preceded by mesh preparation and identification of the boundary and zone cell conditions. Then, the ANSYS Fluent simulator reads the mesh and its properties. As temperature variations occur due to reaction kinetics, the energy equation was enabled in ANSYS Fluent to analyze the effect of temperature. Reaction kinetics and reaction thermodynamics were imported through the CHEMKIN file. The formulas for physical, chemical, mixing, and total exergy were generated and imported to the ANSYS Fluent software through the CFF. Changes in a mole fraction of species and exergy profiles were evaluated.

4. Results and Discussion

The results and discussion encompass contours and profiles of temperature, pressure, mole fractions of reactants and products, and exergy profiles of the reactor model. Contours of temperature, pressure, mole fractions of reactants, and products are shown in Figure 7, whereas the contours of three types of exergy and total exergy are shown in Figure 8. The profiles of temperature and pressure are shown in Figure 9 and Figure 10, respectively, whereas the mole fractions of methane, steam, hydrogen, and carbon monoxide are shown in Figure 11. The mixing exergy profile is demonstrated in Figure 12, followed by the physical, chemical, and total exergy profiles in Figure 13.
The temperature contours and profile shown in Figure 7 and Figure 9, respectively, portray the change in temperature in the longitudinal direction. A constant amount of heat is supplied through the wall of the reactor where the temperature drops in the first part of the reactor (i.e., up to 8 mm) due to the consumption of heat in an endothermic reaction. On the completion of the reactions, in the later part of the reactor (i.e., up to 14 mm), less amount of heat is consumed, and the temperature remains high and constant.
Figure 7 and Figure 10 show the pressure contours and their graph, respectively, revealing the decreasing trend linearly along the length of the reactor. The continuous pressure drop is due to the increase of fluid velocity along the length of the reactor. The increase in velocity causes the pressure drops inversely in fluids. The same effect is stated by the Bernoulli equation.
The contours and the graph of methane conversion are shown in Figure 7 and Figure 11, respectively, portraying that the mole fraction is decreasing rapidly from 0.23 to 0.075 (i.e., up to 12 mm), then remains constant. Methane, a limiting reactant, is converted up to 67.4%. In a previous study reported in the literature, 60% conversion of methane was achieved in a monolith reactor. The conversion of excess reactant (i.e., steam) is shown in Figure 7 and Figure 11, respectively, where the mole fraction of steam drops from 0.77 to 0.56.
The formation of hydrogen is shown in Figure 7 and Figure 11, respectively, demonstrating that the concentration of hydrogen increases along the length of the reactor up to 9 mm. A slight decrease was then noticed up to 12 mm. The decrease in concentration may be caused by the formation of intermediates and reversible reactions, as shown in Table 3. Similarly, the concentration profile of carbon monoxide is shown in Figure 7 and Figure 11, respectively, where the concentration of carbon monoxide is continuously increasing along the length of the reactor.
The performance of the proposed model is compared with a model reported by Cao et al. in the literature [36]; see Table 4. The proposed model achieved a 7.4% higher conversion with 76.6 % smaller surface area compared to the reported work [36]. The higher conversion at the shorter reactor is due to a higher wall temperature, namely, 1477 °C), in the proposed work compared to 900 °C in the reported work. However, total heat consumption in the proposed model is 65.0% lower than the reported work due to the smaller surface area. Furthermore, the reduction in the size of the reactor resulted in a 76.7% lesser requirement of the catalyst in comparison with the reported work [36].
The mixing exergy always has a negative value as exergy of pure components is higher than the components in the mixed form [28]. Decrease in total exergy due to mixing is demonstrated in Figure 8 and Figure 12, respectively. The high conversion rate in the reactor produces new product species at a faster rate. The high-speed molecules intensify the mixing effects. These effects create major irreversibility and contribute significantly to the overall exergy destruction in the reactor. This trend of mixing exergy is validated by another study reported in the literature [17].
The contour and profile of physical exergy are shown in Figure 8 and Figure 13, respectively. At the start of the reactor, the quantity of physical exergy is low which then increases rapidly up to 7 mm. The low quantity of physical exergy at the start of the reactor is due to the low temperature as shown in Figure 7 and Figure 9, respectively. The low temperature is due to the consumption of heat by endothermic reactions of steam methane reforming. The temperature drop causes irreversibility, which results in a decrease in physical exergy at the start of the reactor. A slight decrease in the physical exergy from 7 mm to 12 mm is observed, which is caused by the continued decrease of pressure from start to end of the reactor; physical exergy at this portion only depends on the pressure because its temperature remains constant.
The contours and profile of chemical exergy increase from left to right of the reactor as shown in Figure 8 and Figure 13, respectively. Chemical reactions take place in the catalytic bed and produce new species. These species have high chemical potential which increases the total chemical exergy.
Total exergy is the summation of physical, mixing, and chemical exergy, which is shown in Figure 8 and Figure 13, respectively. The total exergy increases from the start of the reactor till the end of the reactor. The increase in total exergy is due to the combined increasing effect of physical and chemical exergies. The increase in total exergy results in higher work potential of the product, syngas.

5. Conclusions

In this work, the computational fluid dynamics (CFD)-based method was adopted to perform an exergy analysis of the monolith microreactor of the steam methane reforming (SMR) process. Initially, the CFD model of SMR was developed using literature data. In order to incorporate reaction kinetics, CHEMKIN was used. By optimizing the size and the operation conditions, the optimal SMR microreactor achieved a 7.4% higher conversion with 76.6 % smaller surface area compared to the reported work. The higher conversion achieved by the shorter reactor is due to a higher wall temperature, namely, 1477 °C, in the proposed work compared to 900 °C in the reported work. Although the temperature used in the proposed work is higher than the reported work, total heat consumption in the proposed work is 65.0% lower due to the smaller surface area. Furthermore, the reduction in the size of the reactor resulted in a 76.7% reduction in the catalyst requirement.
The exergy analysis was performed by developing a custom field function (CFF)-based algorithm. The exergy analysis helped in evaluating length-wise profiles of all three types of exergy, namely, physical exergy, chemical exergy, and mixing exergy, in the microreactor. The results show that the physical and chemical exergies increase due to the increase in temperature and high chemical potential of product species, respectively. On the other hand, the mixing exergy decreases due to the high rate of mixing effects that causes irreversibility.
In future work, sensitivity analysis and uncertainty analysis will be performed to achieve a further optimization of the process conditions. The sensitivity analysis helps in evaluating the individual impact of process conditions on its outcome, whereas the uncertainty analysis is used to quantify the collective impact of variation in process conditions on its outcome.

Author Contributions

Conceptualization, I.A. and Z.u.R.; Formal analysis, Z.u.R.; Methodology, Z.u.R. and J.M.; Supervision, I.A.; Writing—review and editing, Z.u.R., I.A. and M.K.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CFDcomputational fluid dynamics
SMRsteam methane reforming
PSApressure swing adsorption
Ephmolar physical exergy
Emixmolar mixing exergy
Echmolar chemical exergy
G ¯ i molar Gibbs function
RIdeal gas constant
C v specific heat constant
ϕ dissipation energy per unit volume
Y i mass fraction
D i diffusion coefficient
J i mass flux of component i
S ˙ i production rate
k f , r forward rate coefficient
A k r pre-exponential factor of the reaction r
β k temperature exponent
E a r activation energy of the reaction r
Γ surface site density
M i molar mass of i-th species
RhRhodium catalyst

References

  1. Zhai, X.; Ding, S.; Cheng, Y.; Jin, Y.; Cheng, Y. CFD simulation with detailed chemistry of steam reforming of methane for hydrogen production in an integrated micro-reactor. Int. J. Hydrogen Energy 2010, 35, 5383–5392. [Google Scholar] [CrossRef]
  2. Irani, M.; Alizadehdakhel, A.; Pour, A.N.; Hoseini, N.; Adinehnia, M. CFD modeling of hydrogen production using steam reforming of methane in monolith reactors: Surface or volume-base reaction model? Int. J. Hydrogen Energy 2011, 36, 15602–15610. [Google Scholar]
  3. Saeidi, S.; Fazlollahi, F.; Najari, S.; Iranshahi, D.; Klemeš, J.J.; Baxter, L.L. Hydrogen production: Perspectives, separation with special emphasis on kinetics of WGS reaction: A state-of-the-art review. J. Ind. Eng. Chem. 2017, 49, 1–25. [Google Scholar] [CrossRef]
  4. Baharudin, L.; Watson, M.J. Monolithic substrate support catalyst design considerations for steam methane reforming operation. Rev. Chem. Eng. 2018, 34, 481–501. [Google Scholar] [CrossRef]
  5. Granlund, M.Z.; Görke, O.; Pfeifer, P.; Pettersson, L.J. Comparison between a micro reactor with multiple air inlets and a monolith reactor for oxidative steam reforming of diesel. Int. J. Hydrogen Energy 2014, 39, 18037–18045. [Google Scholar] [CrossRef]
  6. Bhat, S.A.; Sadhukhan, J. Process intensification aspects for steam methane reforming: An overview. AIChE J. 2009, 55, 408–422. [Google Scholar] [CrossRef]
  7. Chen, J.; Yan, L.; Song, W.; Xu, D. Effect of heat and mass transfer on the combustion stability in catalytic micro-combustors. Appl. Therm. Eng. 2018, 131, 750–765. [Google Scholar] [CrossRef]
  8. Yu, X.; Tu, S.T.; Wang, Z.; Qi, Y. Development of a microchannel reactor concerning steam reforming of methanol. Chem. Eng. J. 2006, 116, 123–132. [Google Scholar] [CrossRef]
  9. Theampetch, A.; Chaiwang, W.; Jermkwan, N.; Narataruksa, P.; Sornchamni, T.; Prapainainar, C. Design of Microreactor Flow Channel for Fischer Tropsch Synthesis Using Computational Fluid Dynamic. Energy Procedia 2016, 100, 439–447. [Google Scholar] [CrossRef]
  10. Kansha, Y.; Kotani, Y.; Aziz, M.; Kishimoto, A.; Tsutsumi, A. Evaluation of a self-heat recuperative thermal process based on thermodynamic irreversibility and exergy. J. Chem. Eng. Jpn. 2012, 12we084. [Google Scholar] [CrossRef]
  11. Hosseini, S.E.; Wahid, M.A. Enhancement of exergy efficiency in combustion systems using flameless mode. Energy Convers. Manag. 2014, 86, 1154–1163. [Google Scholar] [CrossRef]
  12. An, H.; Li, A.; Sasmito, A.P.; Kurnia, J.C.; Jangam, S.V.; Mujumdar, A.S. Computational fluid dynamics (CFD) analysis of micro-reactor performance: Effect of various configurations. Chem. Eng. Sci. 2012, 75, 85–95. [Google Scholar] [CrossRef]
  13. Kashid, M.N.; Agar, D.W.; Turek, S. CFD modeling of mass transfer with and without chemical reaction in the liquid–liquid slug flow microreactor. Chem. Eng. Sci. 2007, 62, 5102–5109. [Google Scholar] [CrossRef]
  14. Tzanetis, K.; Martavaltzi, C.; Lemonidou, A. Comparative exergy analysis of sorption enhanced and conventional methane steam reforming. Int. J. Hydrogen Energy 2012, 37, 16308–16320. [Google Scholar] [CrossRef]
  15. Simpson, A.P.; Lutz, A.E. Exergy analysis of hydrogen production via steam methane reforming. Int. J. Hydrogen Energy 2007, 32, 4811–4820. [Google Scholar] [CrossRef]
  16. Boulenouar, M.; Ouadha, A. CFD-Exergy analysis of the flow in a supersonic steam ejector. J. Phys. Conf. Ser. 2015. [Google Scholar] [CrossRef]
  17. Mustafa, J.; Ahmad, I.; Ahsan, M.; Kano, M. Computational fluid dynamics based model development and exergy analysis of naphtha reforming reactors. Int. J. Exergy 2017, 24, 344–363. [Google Scholar] [CrossRef]
  18. Alabi, K.; von Spakovskyρ, M.; Ladeinde, F.; Moorhouse, D.; Camberos, J. A comparison of empirical and CFD-based exergy modeling for the airframe subsystem of aircraft design. In Proceedings of the 25th International Congress of the Aeronautical Sciences, Hamburg, Germany, 3–8 September 2006. [Google Scholar]
  19. Gunjo, D.G.; Mahanta, P.; Robi, P.S. Exergy and energy analysis of a novel type solar collector under steady state condition: Experimental and CFD analysis. Renew. Energy 2017, 114, 655–669. [Google Scholar] [CrossRef]
  20. Debnath, P.; Pandey, K. Exergetic efficiency analysis of hydrogen–air detonation in pulse detonation combustor using computational fluid dynamics. Int. J. Spray Combust. Dyn. 2017, 9, 44–54. [Google Scholar] [CrossRef]
  21. Erguvan, M.; MacPhee, D. Energy and Exergy Analyses of Tube Banks in Waste Heat Recovery Applications. Energies 2018, 11, 2094. [Google Scholar] [CrossRef]
  22. Alabi, K.; Ladeinde, F. Utilizing CFD-based exergy calculations in the design/optimization of a complete aircraft system. In Proceedings of the 45th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 8–11 January 2007; p. 1130. [Google Scholar]
  23. Ao, Y.-A.; Gao, X.-Q.; Shen, L.; Wang, Y.-R.; Feng, G.-H. Exergy analysis of exhaust-gas of burning liquefied-gas in a Chinese kitchen. In Proceedings of the 2009 International Conference on Energy and Environment Technology, Guilin, China, 16–18 October 2009; pp. 40–43. [Google Scholar]
  24. Farmahini-Farahani, M. Investigation of four geometrical parameters on thermal stratification of cold water tanks by exergy analysis. Int. J. Exergy 2012, 10, 332–345. [Google Scholar] [CrossRef]
  25. Lao, L.; Aguirre, A.; Tran, A.; Wu, Z.; Durand, H.; Christofides, P.D. CFD modeling and control of a steam methane reforming reactor. Chem. Eng. Sci. 2016, 148, 78–92. [Google Scholar] [CrossRef]
  26. Hajjaji, N.; Pons, M.N.; Houas, A.; Renaudin, V. Exergy analysis: An efficient tool for understanding and improving hydrogen production via the steam methane reforming process. Energy Policy 2012, 42, 392–399. [Google Scholar] [CrossRef]
  27. Palma, V.; Miccio, M.; Ricca, A.; Meloni, E.; Ciambelli, P. Monolithic catalysts for methane steam reforming intensification: Experimental and numerical investigations. Fuel 2014, 138, 80–90. [Google Scholar] [CrossRef]
  28. Sato, N. Chemical Energy and Exergy: An Introduction to Chemical Thermodynamics for Engineers; Elsevier: Amsterdam, The Netherlands, 2004. [Google Scholar]
  29. Rosen, M.A. Exergy as a tool for sustainability. In Proceedings of the 3rd IASME/WSEAS International Conference on Energy & Environment, Cambridge, UK, 23–25 February 2008; pp. 90–98. [Google Scholar]
  30. Moran, M.J.; Shapiro, H.N.; Boettner, D.D.; Bailey, M.B. Fundamentals of Engineering Thermodynamics; John Wiley & Sons: Hoboken, NJ, USA, 2010. [Google Scholar]
  31. Hinderink, A.; Kerkhof, F.; Lie, A.; Arons, J.D.S.; Van Der Kooi, H. Exergy analysis with a flowsheeting simulator—I. Theory; calculating exergies of material streams. Chem. Eng. Sci. 1996, 51, 4693–4700. [Google Scholar] [CrossRef]
  32. Querol, E.; Gonzalez-Regueral, B.; Perez-Benedito, J.L. Exergy concept and determination. In Practical Approach to Exergy and Thermoeconomic Analyses of Industrial Processes; Springer: Berlin, Germany, 2013; pp. 9–28. [Google Scholar]
  33. Deutschmann, O.; Schmidt, L.D. Two-dimensional modeling of partial oxidation of methane on rhodium in a short contact time reactor. Symp. (Int.) Combust. 1998, 27, 2283–2291. [Google Scholar] [CrossRef]
  34. Chen, J.; Yan, L.; Song, W.; Xu, D. Methane steam reforming thermally coupled with catalytic combustion in catalytic microreactors for hydrogen production. Int. J. Hydrogen Energy 2017, 42, 664–680. [Google Scholar] [CrossRef]
  35. Menon, K.G.; Patnaikuni, V.S. CFD simulation of fuel reactor for chemical looping combustion of Indian coal. Fuel 2017, 203, 90–101. [Google Scholar] [CrossRef]
  36. Cao, C.; Zhang, N.; Chen, X.; Cheng, Y. A comparative study of Rh and Ni coated microchannel reactor for steam methane reforming using CFD with detailed chemistry. Chem. Eng. Sci. 2015, 137, 276–286. [Google Scholar] [CrossRef]
  37. Versteeg, H.K.; Malalasekera, W. An Introduction to Computational Fluid Dynamics: The Finite Volume Method; Pearson Education: London, UK, 2007. [Google Scholar]
  38. Karakaya, C.; Maier, L.; Deutschmann, O. Surface reaction kinetics of the oxidation and reforming of CH4 over Rh/Al2O3 catalysts. Int. J. Chem. Kinet. 2016, 48, 144–160. [Google Scholar] [CrossRef]
  39. Thormann, J.; Maier, L.; Pfeifer, P.; Kunz, U.; Deutschmann, O.; Schubert, K. Steam reforming of hexadecane over a Rh/CeO2 catalyst in microchannels: Experimental and numerical investigation. Int. J. Hydrogen Energy 2009, 34, 5108–5120. [Google Scholar] [CrossRef]
  40. Mundhwa, M.; Thurgood, C.P. Numerical study of methane steam reforming and methane combustion over the segmented and continuously coated layers of catalysts in a plate reactor. Fuel Process. Technol. 2017, 158, 57–72. [Google Scholar] [CrossRef]
  41. Deutschmann, O.; Schwiedemoch, R.; Maier, L.I.; Chatterjee, D. Natural gas conversion in monolithic catalysts: Interaction of chemical reactions and transport phenomena. In Studies in Surface Science and Catalysis; Elsevier: Amsterdam, The Netherlands, 2001; Volume 136, pp. 251–258. [Google Scholar]
Figure 1. Flow diagram of the conventional steam methane reforming process [22].
Figure 1. Flow diagram of the conventional steam methane reforming process [22].
Entropy 21 00399 g001
Figure 2. Circular, square, and trigonal shape arrangement of monolith reactors [4].
Figure 2. Circular, square, and trigonal shape arrangement of monolith reactors [4].
Entropy 21 00399 g002
Figure 3. Tubular reactor model representing a single pore of a monolithic catalyst.
Figure 3. Tubular reactor model representing a single pore of a monolithic catalyst.
Entropy 21 00399 g003
Figure 4. Computational grid of monolith reactor and zoom portion of the grid.
Figure 4. Computational grid of monolith reactor and zoom portion of the grid.
Entropy 21 00399 g004
Figure 5. Convergence residual of our computational fluid dynamics (CFD) model.
Figure 5. Convergence residual of our computational fluid dynamics (CFD) model.
Entropy 21 00399 g005
Figure 6. Schematic of model development process.
Figure 6. Schematic of model development process.
Entropy 21 00399 g006
Figure 7. Total contours of simulation results along the length of the reactor (0–14 mm).
Figure 7. Total contours of simulation results along the length of the reactor (0–14 mm).
Entropy 21 00399 g007
Figure 8. Contours of three types of exergies and total exergy, along the length of the reactor (0–14 mm).
Figure 8. Contours of three types of exergies and total exergy, along the length of the reactor (0–14 mm).
Entropy 21 00399 g008
Figure 9. Temperature profile along the length of the reactor.
Figure 9. Temperature profile along the length of the reactor.
Entropy 21 00399 g009
Figure 10. Pressure profile along the length of the reactor.
Figure 10. Pressure profile along the length of the reactor.
Entropy 21 00399 g010
Figure 11. Conversion profiles of the reactants and formation of the products along the length of the reactor.
Figure 11. Conversion profiles of the reactants and formation of the products along the length of the reactor.
Entropy 21 00399 g011
Figure 12. Mixing exergy profile along the length of the reactor.
Figure 12. Mixing exergy profile along the length of the reactor.
Entropy 21 00399 g012
Figure 13. Profiles of two types of exergies and total exergy along the length of the reactor.
Figure 13. Profiles of two types of exergies and total exergy along the length of the reactor.
Entropy 21 00399 g013
Table 1. Model parameters.
Table 1. Model parameters.
ParametersSymbolValue
Channel lengthL14 mm
Channel diameterD1 mm
Inlet temperatureTin800 °C
Isothermal wall temperatureTw1477 °C
PressureP1 atm
Ratio of steam to methaneS/C3:1
Inlet velocityVin0.45m/s
Number of catalyst active sites (catalyst density)Γ2.7 × 10−9 mol/cm2
Table 2. Mesh properties.
Table 2. Mesh properties.
Properties of MeshValues
Orthogonal quality (minimum)1.0
Ortho skew (maximum)0.0
Aspect ratio (maximum)1.4
The number of nodes7525
The number of cells7200
Minimum volume (m3)2.4 × 10−13
Maximum volume (m3)2.3 × 10−11
Total volume (m3)4.2 × 10−8
Minimum face area (m2)4.1 × 10−5
Maximum face area (m2)4.5 × 10−5
Table 3. Elementary reactions and their kinetics for steam methane reforming on the Rh catalyst [33].
Table 3. Elementary reactions and their kinetics for steam methane reforming on the Rh catalyst [33].
Serial No.ReactionsAkr (mol, cm, s, K) S i o ( )   Ear [KJ/mol]
1.H2 + 2Rh(s)→2H(s) 0.010.0
2.2H(s)→H2 + 2Rh(s)3.0 × 1021 77.8
3.O2 + 2Rh(s)→2O(s)
4.2O(s)→O2 + 2Rh(s)1.33 × 1022 355.2
5.CH4 + Rh(s)→CH4(s) 8.0 × 10−30.0
6.CH4(s)→CH4 + Rh(s)2.0 × 1014 25.1
7.H2O + Rh(s)→H2O(s) 0.10.0
8.H2O(s)→H2O + Rh(s)6.0 × 1013 45
9.CO2 + Rh(s)→CO2(s) 1.0 × 10−50.0
10.CO2(s)→CO2 + Rh(s)3.0 × 108 21.7
11.CO + Rh(s)→CO(s) 5 × 10−10.0
12.CO(s)→CO + Rh(s)1.0 × 1013 133.4
13.H(s) + O(s)→OH(s) + Rh(s)5.0 × 1022 83.7
14.OH(s) + Rh(s)→ H(s) + O(s)3.0 × 1020 37.7
15.H(s) + OH(s)→H2O(s) + Rh(s)3.0 × 1020 33.5
16.H2O(s) + Rh(s)→ H(s) + OH(s)5.0 × 1022 106.4
17.2OH(s)→H2O(s) + O(s)3.0 × 1021 100.8
18.H2O(s) + O(s)→ 2OH(s)3.0 × 1021 171.8
19.C(s) + O(s)→CO(s) + Rh(s)5.0 × 1023 97.9
20.CO(s) + Rh(s)→C(s) + O(s)3.7 × 1021 169.0
21.CO(s) + O(s)→CO2(s) + Rh(s)1.0 × 1019 121.6
22.CO2 + Rh(s)→ CO(s) + O(s)5.0 × 1021 115.3
23.CO(s) + H(s)→HCO(s) + Rh(s)5.0 × 1019 108.9
24.HCO(s) + Rh(s)→ CO(s) + H(s)3.7 × 1021 0.0
25.HCO(s) + Rh(s)→CH(s) + O(s)8.0 × 1023 59.5
26.CH(s) + O(s)→ HCO(s) + Rh(s)3.7 × 1021 167.5
27.CH4(s) + Rh(s)→CH3 + H(s)5.5 × 1020 61.0
28.CH3 + H(s)→ CH4(s) + Rh(s)3.7 × 1021 51.0
29.CH3(s) + Rh(s)→CH2(s) + H(s)3.7 × 1021 103.0
30.CH2(s) + H(s)→ CH3(s) + Rh(s)3.7 × 1021 44.0
31.CH2(s) + Rh(s)→CH(s) + Rh(s)3.7 × 1034 100.0
32.CH(s) + Rh(s)→ CH2(s) + Rh(s)3.7 × 1034 68.0
33.CH(s) + Rh(s)→C(s) + H(s)3.7 × 1021 21.0
34.C(s) + H(s)→ CH(s) + Rh(s)3.7 × 1021 172.8
35.CH4(s) + O(s)→CH3(s) + OH(s)1.7 × 1024 80.3
36.CH3(s) + OH(s)→CH4(s) + O(s)3.7 × 1021 24.3
37.CH3(s) + O(s)→CH2(s) + OH(s)3.7 × 1024 120.3
38.CH2(s) + OH(s)→CH3(s) + O(s)3.7 × 1021 15.1
39.CH2(s) + O(s)→CH(s) + OH(s)3.7 × 1024 114.5
40.CH(s) + OH(s)→CH2(s) + O(s)3.7 × 1021 36.8
41.CH(s) + O(s)→C(s) + OH(s)3.7 × 1021 30.1
42.C(s) + OH(s)→CH(s) + O(s)3.7 × 1021 136.0
43.CO(s) + H(s)→C(s) + OH(s)3.7 × 1021 142.0
44.C(s) + OH(s)→CO(s) + H(s)3.7 × 1020 25.5
Table 4. Comparison of the proposed model with the model reported in the literature [36].
Table 4. Comparison of the proposed model with the model reported in the literature [36].
ParametersLiterature Data [36]Proposed Model Data
Feed temperature800 °C800 °C
Wall temperature900 °C1477 °C
Pressure1 atm1 atm
Steam to methane ratio3:13:1
Inlet velocity0.45 m/s0.45m/s
Length6.0 × 10−2 m1.4 × 10−2 m
Surface area1.88 × 10−4 m24.4 × 10−5 m2
Heat requirement22.05 kW7.7 kW
Conversion60 %67.4 %
Catalyst requirement per length5.076 × 10−9 moles1.18 × 10−9 moles

Share and Cite

MDPI and ACS Style

Rahman, Z.u.; Ahmad, I.; Kano, M.; Mustafa, J. Model Development and Exergy Analysis of a Microreactor for the Steam Methane Reforming Process in a CFD Environment. Entropy 2019, 21, 399. https://doi.org/10.3390/e21040399

AMA Style

Rahman Zu, Ahmad I, Kano M, Mustafa J. Model Development and Exergy Analysis of a Microreactor for the Steam Methane Reforming Process in a CFD Environment. Entropy. 2019; 21(4):399. https://doi.org/10.3390/e21040399

Chicago/Turabian Style

Rahman, Zia ur, Iftikhar Ahmad, Manabu Kano, and Jawad Mustafa. 2019. "Model Development and Exergy Analysis of a Microreactor for the Steam Methane Reforming Process in a CFD Environment" Entropy 21, no. 4: 399. https://doi.org/10.3390/e21040399

APA Style

Rahman, Z. u., Ahmad, I., Kano, M., & Mustafa, J. (2019). Model Development and Exergy Analysis of a Microreactor for the Steam Methane Reforming Process in a CFD Environment. Entropy, 21(4), 399. https://doi.org/10.3390/e21040399

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop