Next Article in Journal
Separation and Concentration of Nitrogen and Phosphorus in a Bipolar Membrane Electrodialysis System
Next Article in Special Issue
Effect of Packing Nonuniformity at the Fiber Bundle–Case Interface on Performance of Hollow Fiber Membrane Gas Separation Modules
Previous Article in Journal
Dense Electrode Layers-Supported Microtubular Oxygen Pump
Previous Article in Special Issue
Effect of Polystyrene Synthesis Method on Water Sorption and Glass Transition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Toward Minimal Complexity Models of Membrane Reactors for Hydrogen Production

by
Maria Anna Murmura
*,
Stefano Cerbelli
,
Ludovica Manozzi
and
Maria Cristina Annesini
Department of Chemical Engineering Materials and Environment, University of Rome “La Sapienza”, Via Eudossiana 18, 00184 Rome, Italy
*
Author to whom correspondence should be addressed.
Membranes 2022, 12(11), 1115; https://doi.org/10.3390/membranes12111115
Submission received: 19 October 2022 / Revised: 1 November 2022 / Accepted: 3 November 2022 / Published: 8 November 2022

Abstract

:
Membrane reactors are inherently two-dimensional systems that require complex models for an accurate description of the different transport phenomena involved. However, when their performance is limited by mass transport within the reactor rather than by the selective product permeation across the membrane, the 2D model may be significantly simplified. Here we extend results previously found for methane steam reforming membrane reactors to show that such simplified two-dimensional model admits either a straightforward analytical solution for the cross-section averaged concentration profile, or can be reduced to a 1D model with an enhanced Sherwood number, depending on the stoichiometry of the reaction considered. Interestingly, the stoichiometry does not affect the expression of the enhanced Sherwood number, indicating that a versatile tool has been developed for the determination of membrane reactor performance at an extremely low computational cost and good degree of accuracy.

1. Introduction

Membrane reactors represent an interesting solution for the decentralized production of high-purity hydrogen, as they allow the integrated production and separation of hydrogen in a single device. In addition, the removal of hydrogen as it is being produced allows to shift the equilibrium of the reaction toward the products, thereby increasing the overall conversion and allowing to carry out endothermic reactions at temperatures lower than those traditionally employed. This solution represents an advantage in terms of energy efficiency, but could also help avoid undesired reactions that take place at higher temperatures, such as coke formation in the presence of hydrocarbons. The use of these devices has been proposed and investigated for many applications, including the reforming of light hydrocarbons [1,2,3], the water-gas shift reaction [4,5], ammonia decomposition [6,7,8], and the dehydrogenation of alkanes [9,10,11]. The simplest reactor configuration envisaged is a tube-in-tube reactor, in which the reaction takes place in the annular volume filled with catalyst particles between the two concentric tubes. Heat may be provided through the outer wall of the tube. A hydrogen-permeable membrane is placed on the outer wall of the inner-most tube. Hydrogen is removed through the selective membrane and flows into the inner tube. A sweep gas is often employed in the permeate side to reduce the partial pressure of hydrogen and increase the driving force for permeation.
In the past years a significant evolution in modeling membrane reactors has been observed, from simple one-dimensional models [12,13] to more complex two- or three-dimensional models [14,15]. While the former models are relatively simple to develop and can be solved at a low computational cost, they could be inaccurate even in predicting integral quantities [16,17]. In addition one-dimensional models seem to be conceptually inadequate for the description of membrane reactors, in which the transport of hydrogen takes place in the direction transversal to that of the main gas flow. Multi-dimensional models, on the other hand, allow to simultaneously account for several transport phenomena and could potentially be very accurate; however, the numerous parameters included in the more complex models are generally difficult to determine with accuracy and the effect of changes in the operating conditions or geometry are more difficult to grasp intuitively from a small set of results, because a single parameter could affect more than one phenomenon.
In this context, the development of one-dimensional models in which transversal concentration gradients are accounted for through a mass transfer coefficient is of interest, as it could combine the benefits of both one- and multi-dimensional models. The difficulty lies in the derivation of an expression for a mass transfer coefficient capable of accounting for the effects not only of convection and dispersion, but also of the reaction and of hydrogen permeation across the selective membrane. Previous experience on modelling membrane reactors for methane steam reforming had shown that, when the resistance to hydrogen permeation across the membrane is negligible, the behavior of the system may be described by a simplified 2D model [18,19]. Here we show that, for reactions characterized by simple stoichiometries, a straightforward analytical solution to the problem exists; whereas for more complex stoichiometries, it is necessary to introduce a mass transfer coefficient. Interestingly, it is found that the stoichiometry of the reaction considered does not affect the expression of the mass transfer coefficient. The present article is divided as follows. In Section 2 the problem is stated in greater detail, the modelling approaches are described in Section 3, and the results are presented in Section 4.

2. Statement of the Problem

Previous work on the description of membrane reactors for the production of hydrogen through the steam reforming of methane had highlighted that, depending on the operating conditions, the behavior of such systems could be limited by either radial mass transport within the packed bed or hydrogen permeation across the selective membrane. At a given inlet velocity, the transition between the two regimes was found to depend on the operating pressure. This is due to the fact that while the rate of the convective transport increases almost proportionally with pressure, the rate of permeation depends on the square root of pressure. In fact, hydrogen-selective membranes are generally Pd-based and the rate of permeation depends on the membrane permeability and the difference between the square roots of the hydrogen partial pressures in the retentate and permeate sides according to Sieverts’ law [20]
J h = P m p H 2 r p H 2 p
where J h is the molar hydrogen flux across the membrane, P m is the membrane permeability, p H 2 is the hydrogen partial pressure, and the superscripts r and p indicate the retentate and permeate sides, respectively. Under conditions in which the resistance offered by the membrane is negligible, the behavior of the reactor may be described through a simplified two-dimensional model. A more detailed explanation of these observations may be found in [19].
In [21], an expression for an enhanced Sherwood number for the evaluation of the hydrogen concentration profiles in steam reforming membrane reactors was developed, to further simplify the description of the problem to a one-dimensional model in which concentration gradients are accounted for through a mass transfer coefficient, or Sherwood number when the problem is analyzed in its dimensionless form. It is well known that the Sherwood number is a dimensionless group commonly employed to evaluate the mass transport coefficient, k y , between the bulk of a fluid, b, and a given surface, s
S h = y i r s u r f a c e y i , s y i , b = k y l c D i
and that it expresses the ratio between the rate of mass transport by convection and by diffusion (or dispersion). The Sherwood number is generally evaluated as a function of the Reynolds and Schmidt numbers; however, in reactive systems, the chemical reactions alter the local composition of the mixtures, thereby affecting the species concentration profile and hence the overall transport coefficient. In membrane reactors, there is the additional effect of the selective membrane. Consequently, correlations traditionally employed for the evaluation of the mass transfer coefficient lose significance.
The aim of the present work has been to verify the possibility of extending the findings established in the conetext of steam reforming of hydrocarbons to different reactive systems, specifically to the dehydrogenation of propane. The range of operating conditions in which it is possible to describe the behavior of membrane reactors through this simplified approach with a good degree of accuracy was also investigated. In addition, the possibility of obtaining a simple analytical expression for the simplified 2D model was investigated for reactions characterized by a simple stoichiometry.
The problem was tackled through the following procedure
1.
A simplified 2D model, derived in [19], was applied to the propane dehydrogenation reaction;
2.
an analytical solution to the average hydrogen concentration profile along the length of the reactor was obtained;
3.
an analytical expression for the enhanced Sherwood number was derived and employed in a 1D model of the membrane reactor;
The model was developed for a tube-in-tube catalytic reactor configuration, in which the propane dehydrogenation reaction takes place
C 3 H 8 H 2 + C 3 H 6 Δ H 0 = 124.3 kJ / mol
Given its endothermicity, this reaction requires high temperatures and is prone to coke formation. To reduce both the energy cost and the risk of catalyst deactivation, it would be useful to lower the operating temperature; on the other hand, this would cause a reduction in the equilibrium conversion of propane. The problem could be solved by employing membrane reactors, in which hydrogen is selectively removed from the reacting volume as it is being produced, thereby shifting the equilibrium of the reaction.

3. Modelling

The 2D simplified model was developed under the following simplifying assumptions:
constant temperature;
negligible axial dispersion;
negligible radial convection;
vanishing hydrogen pressure in the permeate side;
negligible resistance to hydrogen transport across the membrane;
constant density;
local equilibrium conditions
With these assumptions, the mass balance equations read
U y i z + D r r 1 r y i r + 2 y i r 2 + ν i r p = 0
where y i is the molar fraction of the i-th component; U is the axial velocity, which is constant following the assumption of constant density; D r r is the effective dispersion coefficient in the radial direction, considered to be the same for all components; ν i is the stoichiometric coefficient of the i-th component; r p is the volumetric rate of propane consumption; and the radial and axial coordinates are represented by r and z, respectively.
The rate of propane dehydrogenation was expressed according to the correlation proposed by Sheintuch et al. [22] while neglecting the inhibiting effect on kinetics by propylene adsorption on the catalyst
r p = k p p C 3 H 8 1 η
where k p is the kinetic constant, which, for a given catalyst, depends on temperature only; p C 3 H 8 is the partial pressure of propane, whereas the term 1 η accounts for the distance from chemical equilibrium conditions, with values ranging between 0 and 1
η = p H 2 p C 3 H 6 p C 3 H 8 1 K e q
The boundary conditions are of fixed inlet composition, namely the equilibrium composition under the temperature and pressure conditions considered
y i = y i 0 ; in z = 0
impermeability to the flux of all components on the outer wall
y i r = 0 ; in r = R 2
impermeability of the inner wall to propane and propylene
y i r = 0 ; in r = R 1 for i H 2
and hydrogen permeation according to Sieverts’ law through the membrane placed on the outer wall of the inner tube
P R T D r r y H 2 r = P m p H 2 ; in r = R 1
In its dimensionless formulation, Equation (4) becomes
y i z ˜ + D ˜ r r P e 1 r ˜ y i r ˜ + 2 y i r ˜ 2 + ν i D a y C 3 H 8 1 η = 0
where D ˜ r r is the dimensionless ratio between the effective dispersion coefficient in the radial direction and the molecular diffusion coefficient, D m , the Peclet number
P e = U l c D m
represents the ratio between the characteristic times of diffusion and convection, and the Damkholer number
D a = R T k l c U
The dimensionless axial and radial coordinates were both defined on the basis of the characteristic length, l c , defined as the difference between the outer and inner reactor radii
l c = R 2 R 1
The boundary conditions remain virtually unchanged, with the exception of the one for hydrogen on the membrane wall, which becomes
D ˜ r r P e y i r = γ P ˜ y H 2
in which a dimensionless parameter representing the ratio between the characteristic times of convection and permeation, γ , is introduced
γ = P m R T P a t m 1 / 2 U
having made the operating pressure dimensionless by defining it in terms of the atmospheric pressure
P ˜ = P P a t m
For thin annular channels it is possible to describe the problem in terms of Cartesian coordinates
y i z ˜ + D ˜ r r P e 2 y i x ˜ 2 + ν i D a y C 3 H 8 1 η = 0
From here on, the dimensionless axial and transversal coordinates will be written as z and x, respectively, for easier readability.
As described in greater detail elsewhere [23] the problem described above may be simplified by considering an auxiliary variable, defined as a linear combination of the molar fractions of hydrogen and of the species from which it is produced, propane in the present case
Y = ν H 2 y C 3 H 8 ν C 3 H 8 y H 2
for which the following balance equation holds
Y z ε 2 Y x 2 = 0
with
ε = D ˜ r r P e
and
x = r R 1 R 2 R 1
To be solved with the boundary conditions
Y z = 0 = Y 0 = ν H 2 y C 3 H 8 0 ν C 3 H 8 y H 2 0
Y x x = 1 = 0
Y x = 0 = 0
The latter boundary condition is a consequence of having considered negligible resistance to the transport of hydrogen across the membrane, from which one can set a y H 2 = 0 on the membrane wall, and local chemical equilibrium, a condition that can be met only if the partial pressures of hydrogen and propane both go to zero when approaching the membrane wall. The above homogeneous problem admits the following analytical solution
Y = 2 Y 0 l = 0 1 λ l exp λ l 2 ε z sin λ l x
where λ l = π 2 ( 2 l + 1 ) .
Keeping in mind that, given the impermeability of the membrane to all components other than hydrogen,
Y x x = 0 = y H 2 x x = 0
the Sherwood number, defined in Equation (2), is given by
S h ( z ) = 2 Y 0 l = 0 exp λ l 2 ε z y ¯ H 2 ( z )
where
Y 0 = ν H 2 y C 3 H 8 0 ν C 3 H 8 y H 2 0
and where y ¯ H 2 ( z ) is the average molar fraction of hydrogen and may be determined by solving the following set of equations
Y ¯ ( z ) = 2 Y 0 l = 0 1 λ l 2 e x p λ l 2 ε z
Y ¯ ( z ) = ν H 2 y ¯ C 3 H 8 ( z ) ν C 3 H 8 y ¯ H 2 ( z )
K e q = P 1 y ¯ C 3 H 8 y ¯ H 2 y ¯ H 2 y ¯ C 3 H 8
Given the simple stoichiometry of the reaction, the above set of equations admits the analytical solution
y ¯ H 2 = K e q Y ¯ ( z ) P ν H 2 Y ¯ ( z ) K e q ν C 3 H 8
Having obtained the above solution, it is interesting to determine how such an average concentration profile would differ from the one obtained with a 1D model with an enhanced Sherwood number and to compare the expression of S h ( z ) obtained from the propane dehydrogenation reaction with the one previously obtained for methane steam reforming.

4. Results

Figure 1a,b show the Sherwood number along the length of the reactor, evaluated according to Equation (28) at atmospheric pressure, temperatures between 400 and 700 ° C and at Pe values of 10 and 100, respectively. With regards to the results shown in Figure 1b, it is interesting to note that at 400 ° C, the Sherwood number remains almost constant along the entire reactor length. This is mainly attributable to the fact that at low temperatures, the equilibrium conversion of propane is lower than 10%, meaning that the hydrogen concentration is low and the flux permeating across the membrane is not sufficiently high to significantly increase the degree of propane conversion.
The results show that, as in the case of the methane steam reforming reaction, the variation of the Sherwood number along the reactor axis may be described through an expression of the kind
S h ( z ) = S h 0 z n
at high values of the Peclet number, whereas the trend for Sh changes along the length of the reactor when the Peclet number is low. The different behavior observed depending on the value of P e was expected because, when all other conditions are kept constant, increasing the value of Pe corresponds to a decrease in the dispersion coefficient; consequently, the performance of the reactor is increasingly limited by mass transport within the packed bed, i.e., the main assumption under which the simplified 2D model was developed.
We now move on to consider the same approximated expression for S h ( z ) as the one employed for the methane steam reforming reaction, namely [21]
S h ( z ) = S h 0 z z * 0.45
with
S h 0 = 2 Ω 0 1 ν p f i n W h l = 0 exp λ l 2 ε z y h ( z * )
and z * 1 , and determine the average hydrogen concentration profiles by solving the following problem
d F ˜ H 2 d z ˜ = D a P ˜ F ˜ C 3 H 8 i = 1 3 F ˜ i ( 1 η ) S h P e 2 R 1 R 2 + R 1 P ˜ F ˜ H 2 i = 1 3 F ˜ i
d F ˜ C 3 H 8 d z ˜ = D a P ˜ F ˜ C 3 H 8 i = 1 3 F ˜ i ( 1 η )
d F ˜ C 3 H 6 d z ˜ = D a P ˜ F ˜ C 3 H 8 i = 1 3 F ˜ i ( 1 η )
Figure 2 shows a comparison between the values of the Sherwood number along the length of the reactor, evaluated from the simplified 2D model and from Equation (33). The agreement is very good, particularly if taking into account that the approximate expression was simply predicted from results previously obtained for the methane steam reforming reaction.
Figure 3 shows the average hydrogen concentration profiles obtained at P e = 100 under different temperature and pressure conditions by solving the 1D model (black dashed curves) of Equations (35a–c) and the simplified 2D model (Equation (31), solid red curves).
It is clear that the agreement between the trends predicted with the two models is excellent, even though the expression for the Sherwood number was not developed specifically for this reactive system. The results also show that there is an error in the evaluation of the hydrogen molar fraction close to the reactor inlet, which causes the discrepancy between the two sets of results. On the other hand, the hydrogen recovery, defined as the ratio between the flowrate of hydrogen permeated across the membrane and the inlet flow rate of hydrogen, evaluated from the two models coincides, as shown from the results reported in Figure 4. This indicates that the proposed solution can be successfully employed for the evaluation of integral quantities at a very low computational cost.

5. Conclusions

A simplified 2D model, initially developed to describe the behavior of a membrane reactor for the steam reforming of methane was applied to the propane dehydrogenation reaction in the same type of device. It was shown that, due to the simple stoichiometry of the reaction considered, it is possible to obtain a straightforward analytical expression for the average hydrogen concentration profile along the length of the reactor. Although, in comparison to two-dimensional models, this solution does not allow to have detailed information on the transversal concentration profiles formed within the reactor, significant information is obtained at almost no computational cost. In addition, an attempt was made to describe the performance of the system through a 1D model making use of the same expression for an enhanced Sherwood number previously obtained for the methane steam reforming reactor, with the aim of validating the previously proposed approach. An excellent agreement was obtained between the trends of the average hydrogen concentration profiles along the reactor under conditions in which the behavior of device is limited by mass transport across the reactor bed. The results show that it is possible to evaluate integral quantities with a good degree of accuracy at very low computational cost. The possibility of employing the same expression for S h ( z ) regardless of the reactive system indicates that the tool developed is extremely versatile.

Author Contributions

Conceptualization: M.A.M., S.C. and M.C.A.; Methodology: M.A.M. and S.C.; Software: L.M.; Validation M.A.M. and L.M.; Data curation: L.M.; Formal Analysis: M.A.M., S.C. and M.C.A.; Original draft preparation: M.A.M.; Writing—Review & Editing, M.A.M. and M.C.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Habib, M.A.; Harale, A.; Paglieri, S.; Alrashed, F.S.; Al-Sayoud, A.; Rao, M.V.; Nemitallah, M.A.; Hossain, S.; Hussien, M.; Ali, A. Palladium-alloy membrane reactors for fuel reforming and hydrogen production: A review. Energy Fuels 2021, 35, 5558–5593. [Google Scholar] [CrossRef]
  2. Zhao, Q.; Su, B.; Wang, H.; He, A.; He, R.; Kong, H.; Hu, X. Mid/low-temperature solar hydrogen generation via dry reforming of methane enhanced in a membrane reactor. Energy Convers. Manag. 2021, 240, 114254. [Google Scholar] [CrossRef]
  3. Saidi, M.; Moradi, P. Conversion of biodiesel synthesis waste to hydrogen in membrane reactor: Theoretical study of glycerol steam reforming. Int. J. Hydrog. Energy 2020, 45, 8715–8726. [Google Scholar] [CrossRef]
  4. Pati, S.; Jangam, A.; Wang, Z.; Dewangan, N.; Wai, M.H.; Kawi, S. Catalytic Pd0.77Ag0.23 alloy membrane reactor for high temperature water-gas shift reaction: Methane suppression. Chem. Eng. J. 2019, 362, 116–125. [Google Scholar] [CrossRef]
  5. Basile, A.; Curcio, S.; Bagnato, G.; Liguori, S.; Jokar, S.M.; Iulianelli, A. Water gas shift reaction in membrane reactors: Theoretical investigation by artificial neural networks model and experimental validation. Int. J. Hydrogen Energy 2015, 40, 5897–5906. [Google Scholar] [CrossRef]
  6. Wang, B.; Kong, H.; Wang, H.; Wang, Y.; Hu, X. Kinetic and thermodynamic analyses of mid/low-temperature ammonia decomposition in solar-driven hydrogen permeation membrane reactor. Int. J. Hydrogen Energy 2019, 44, 26874–26887. [Google Scholar] [CrossRef]
  7. Itoh, N.; Kikuchi, Y.; Furusawa, T.; Sato, T. Tube-wall catalytic membrane reactor for hydrogen production by low-temperature ammonia decomposition. Int. J. Hydrogen Energy 2021, 46, 20257–20265. [Google Scholar] [CrossRef]
  8. Cechetto, V.; Di Felice, L.; Medrano, J.A.; Makhloufi, C.; Zuniga, J.; Gallucci, F. H2 production via ammonia decomposition in a catalytic membrane reactor. Fuel Process. Technol. 2021, 216, 106772. [Google Scholar] [CrossRef]
  9. Shelepova, E.V.; Vedyagin, A.A. Intensification of the dehydrogenation process of different hydrocarbons in a catalytic membrane reactor. Chem. Eng. Process.-Process Intensif. 2020, 155, 108072. [Google Scholar] [CrossRef]
  10. He, R.; Wang, Y.; Wang, H.; Lundin, S.B.; Wang, B.; Kong, H.; Lu, X.; Wang, J.; Li, W. A mid/low-temperature solar-driven integrated membrane reactor for the dehydrogenation of propane—A thermodynamic assessment. Appl. Therm. Eng. 2021, 193, 116952. [Google Scholar] [CrossRef]
  11. Dittrich, C.J. The role of heat transfer on the feasibility of a packed-bed membrane reactor for propane dehydrogenation. Chem. Eng. J. 2020, 381, 122492. [Google Scholar] [CrossRef]
  12. Patrascu, M.; Sheintuch, M. On-site pure hydrogen production by methane steam reforming in high flux membrane reactor: Experimental validation, model predictions and membrane inhibition. Chem. Eng. J. 2015, 262, 862–874. [Google Scholar] [CrossRef]
  13. Di Marcoberardino, G.; Sosio, F.; Manzolini, G.; Campanari, S. Fixed bed membrane reactor for hydrogen production from steam methane reforming: Experimental and modeling approach. Int. J. Hydrogen Energy 2015, 40, 7559–7567. [Google Scholar] [CrossRef]
  14. Ma, R.; Castro-Dominguez, B.; Dixon, A.G.; Ma, Y.H. CFD study of heat and mass transfer in ethanol steam reforming in a catalytic membrane reactor. Int. J. Hydrogen Energy 2018, 43, 7662–7674. [Google Scholar] [CrossRef]
  15. Bian, Z.; Xia, H.; Wang, Z.; Jiang, B.; Yu, Y.; Yu, K.; Zhong, W.; Kawi, S. CFD simulation of a hydrogen-permeable membrane reactor for CO2 reforming of CH4: The interplay of the reaction and hydrogen permeation. Energy Fuels 2020, 34, 12366–12378. [Google Scholar] [CrossRef]
  16. Mori, N.; Nakamura, T.; Noda, K.; Sakai, O.; Takahashi, A.; Ogawa, N.; Sakai, H.; Iwamoto, Y.; Hattori, T. Reactor configuration and concentration polarization in methane steam reforming by a membrane reactor with a highly hydrogen-permeable membrane. Ind. Eng. Chem. Res. 2007, 46, 1952–1958. [Google Scholar] [CrossRef]
  17. Bi, Y.; Xu, H.; Li, W.; Goldbach, A. Water–gas shift reaction in a Pd membrane reactor over Pt/Ce0.6Zr0.4O2 catalyst. Int. J. Hydrogen Energy 2009, 34, 2965–2971. [Google Scholar] [CrossRef]
  18. Murmura, M.A.; Cerbelli, S.; Annesini, M.C. Transport-reaction-permeation regimes in catalytic membrane reactors for hydrogen production. The steam reforming of methane as a case study. Chem. Eng. Sci. 2017, 162, 88–103. [Google Scholar]
  19. Murmura, M.A.; Cerbelli, S.; Annesini, M.C. An equilibrium theory for catalytic steam reforming in membrane reactors. Chem. Eng. Sci. 2017, 160, 291–303. [Google Scholar] [CrossRef]
  20. Ward, T.; Dao, T. Model of hydrogen permeation behavior in palladium membranes. J. Membr. Sci. 1999, 153, 211–231. [Google Scholar] [CrossRef]
  21. Murmura, M.A.; Cerbelli, S.; Annesini, M.C.; Sheintuch, M. Derivation of an enhanced Sherwood number accounting for reaction rate in membrane reactors. Steam reforming of methane as case study. Catal. Today 2021, 364, 285–293. [Google Scholar] [CrossRef]
  22. Sheintuch, M.; Liron, O.; Ricca, A.; Palma, V. Propane dehydrogenation kinetics on supported Pt catalyst. Appl. Catal. A Gen. 2016, 516, 17–29. [Google Scholar] [CrossRef]
  23. Murmura, M.A.; Roccehtti, C.; Annesini, M.C. An Enhanced Sherwood Number to Model the Hydrogen Transport in Membrane Steam Reformers. Membranes 2021, 11, 805. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Sherwood number along the length of the reactor at temperatures of 400, 500, 600, and 700 ° C and 1 atm evaluated from Equation (28) at Pe = 10 (a) and 100 (b).
Figure 1. Sherwood number along the length of the reactor at temperatures of 400, 500, 600, and 700 ° C and 1 atm evaluated from Equation (28) at Pe = 10 (a) and 100 (b).
Membranes 12 01115 g001
Figure 2. Sherwood number along the length of the reactor evaluated from the simplified 2D model (black curves) and the approximated correlation of Equation (33) (red curves) at Pe = 100, 700 ° C, and pressures of 1, 3, and 10 bar.
Figure 2. Sherwood number along the length of the reactor evaluated from the simplified 2D model (black curves) and the approximated correlation of Equation (33) (red curves) at Pe = 100, 700 ° C, and pressures of 1, 3, and 10 bar.
Membranes 12 01115 g002
Figure 3. Average hydrogen concentration profiles evaluated from the simplified 2D model (Equation (31), solid red curves) and from the 1D model (Equation (35a–c), dashed black curves) at Pe = 100, and (a) T = 500 ° C at 1 and 10 bar; (b) P = 10 bar at 500, 600, and 700 ° C.
Figure 3. Average hydrogen concentration profiles evaluated from the simplified 2D model (Equation (31), solid red curves) and from the 1D model (Equation (35a–c), dashed black curves) at Pe = 100, and (a) T = 500 ° C at 1 and 10 bar; (b) P = 10 bar at 500, 600, and 700 ° C.
Membranes 12 01115 g003
Figure 4. Hydrogen recoveries evaluated from the simplified 2D (solid curves) and 1D (dashed curves) models at Pe = 100, P = 10 atm and temperatures of 600 and 700 ° C.
Figure 4. Hydrogen recoveries evaluated from the simplified 2D (solid curves) and 1D (dashed curves) models at Pe = 100, P = 10 atm and temperatures of 600 and 700 ° C.
Membranes 12 01115 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Murmura, M.A.; Cerbelli, S.; Manozzi, L.; Annesini, M.C. Toward Minimal Complexity Models of Membrane Reactors for Hydrogen Production. Membranes 2022, 12, 1115. https://doi.org/10.3390/membranes12111115

AMA Style

Murmura MA, Cerbelli S, Manozzi L, Annesini MC. Toward Minimal Complexity Models of Membrane Reactors for Hydrogen Production. Membranes. 2022; 12(11):1115. https://doi.org/10.3390/membranes12111115

Chicago/Turabian Style

Murmura, Maria Anna, Stefano Cerbelli, Ludovica Manozzi, and Maria Cristina Annesini. 2022. "Toward Minimal Complexity Models of Membrane Reactors for Hydrogen Production" Membranes 12, no. 11: 1115. https://doi.org/10.3390/membranes12111115

APA Style

Murmura, M. A., Cerbelli, S., Manozzi, L., & Annesini, M. C. (2022). Toward Minimal Complexity Models of Membrane Reactors for Hydrogen Production. Membranes, 12(11), 1115. https://doi.org/10.3390/membranes12111115

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