Next Article in Journal
Engineering of MSC-Derived Exosomes: A Promising Cell-Free Therapy for Osteoarthritis
Next Article in Special Issue
Classification of Nanomaterials and the Effect of Graphene Oxide (GO) and Recently Developed Nanoparticles on the Ultrafiltration Membrane and Their Applications: A Review
Previous Article in Journal
Anti-Pectin Fouling Performance of Dopamine and (3-Aminopropy) Triethoxysilane-Coated PVDF Ultrafiltration Membrane
Previous Article in Special Issue
Tribological Performance of Diamond Films with Different Roughnesses of Silicon Nitride Substrates and Carbon Source Concentrations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model-Based Performance Analysis of Membrane Reactor with Ethanol Steam Reforming over a Monolith

by
Ludmilla Bobrova
,
Nadezhda Vernikovskaya
,
Nikita Eremeev
and
Vladislav Sadykov
*
Federal Research Center Boreskov Institute of Catalysis SB RAS, 630090 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
Membranes 2022, 12(8), 741; https://doi.org/10.3390/membranes12080741
Submission received: 30 June 2022 / Revised: 21 July 2022 / Accepted: 26 July 2022 / Published: 28 July 2022
(This article belongs to the Collection Feature Papers in Inorganic Membranes and Hybrid Materials)

Abstract

:
Membrane reactors (MR) with an appropriate catalyst are considered to be an innovative and intensified technology for converting a fuel into the hydrogen-rich gas with the simultaneous recovery of high-quality hydrogen. Characteristics of an asymmetric membrane disk module consisting of a gas-tight nanocomposite functional coating (Ni + Cu/Nd5.5WO11.25-δ mixed proton-electron conducting nanocomposite) deposited on a gas-permeable functionally graded substrate has previously been extensively studied at lab-scale using MRs, containing the catalyst in a packed bed and in the form of a monolith. The catalytic monolith consisted of a FeCrAl substrate with a washcoat and an Ni + Ru/Pr0.35Ce0.35Zr0.35O2 active component. It has been shown that the driving potential for hydrogen permeation across the same membrane in a monolithic catalyst –assisted MR is greater compared to the packed bed catalyst. This paper presents results of the study where a one-dimensional isothermal model was used to interrelate catalytic and permeation phenomena in a MR with ethanol steam reforming over the monolith, operating at atmospheric pressure and in the temperature range of 700–900 °C. The developed mathematical reaction–transport model for the constituent layers of the catalyst-asymmetric membrane assembly together with a Sieverts’ equation for the functional dense layer, taking also into account the effect of boundary layers, was implemented in a COMSOL Multiphysics environment. Good agreement with the experimental data of the lab-scale MR with reasonable parameters values is provided. In numerical experiments, concentration profiles along the reactor axis were obtained, showing the effect of the emerging concentration gradient in the boundary layer adjacent to the membrane. Studies have shown that a MR with a catalytic monolith along with appropriate organization of a stagnant feed flow between the monolith and the membrane surface may enhance production and flux of hydrogen, as well as the efficiency characteristics of the reactor compared to a reactor with packed beds.

1. Introduction

The use of an expensive hydrogen-fueling network motivated by concerns about environmental issues would lead to high costs in the fuel delivery system. On site hydrogen generation from a hydrocarbon feedstock is considered to be preferable [1,2].
Membrane-based reactors for hydrogen production are among the leading process intensification technologies for small or medium scale applications that open new pathways for both material chemistry and process engineering. The term membrane reactor first began to appear in the literature on chemical technology around 1980 [3]. Although there is no generally accepted definition of a membrane reactor, this term usually refers to membrane devices in which chemical conversion is simultaneously carried out under conditions in which the unique contacting and separating properties of membranes are used. Basically, two configurations of the membrane reactor system are distinguished: in the first case, the reactor and the membrane separation equipment are simply connected in series, while in the second case, the real membrane reactor concept combines the membrane separation process with chemical or biochemical reactions into a single unit. The combination of a membrane separation and a catalytic reaction facilitates process miniaturization, continuous operation and energy saving [4,5].
Membrane reactors with a proper catalyst capable of providing a higher fuel conversion and a membrane possessing permselectivity toward a product gas allow very pure hydrogen to be produced, reducing the down-stream purification load [6,7,8,9,10]. Extensive research has focused on the application of membrane reactors for hydrogen production by the reforming of bio-ethanol, which is an important candidate as a chemical carrier of hydrogen [8,11,12,13,14,15]. The main reaction products are H2 and carbon-containing species (CO2, CO, CH4 and/or C) depending on the catalyst and ethanol-to-water ratio [16,17]. The largest amount of hydrogen is obtained by steam reforming of ethanol, so theoretically 6 moles of hydrogen are formed per mole of ethanol in the feedstock. According to the stoichiometry of the reaction, the molar ratio of water vapor to ethanol (S/E) for the complete conversion of the main reagents into carbon dioxide and hydrogen is 3. According to the thermodynamic study performed by Sun et al. [18], for the reforming reaction with the H2, CO, CO2 and CH4 productions, at about 1000 K more than 5.1 mol of H2 (close to 90% of equilibrium yield) can be obtained at S/E > 8. The CO2 yield maximizes at about 900 K with S/E > 8. Yet, increasing the S/E molar ratio above 6 has been found to decrease the hydrogen recovery, and more energy is required for feed preparation [15,19].
It has been well studied that the conversion of ethanol and selectivity with respect to various products largely depend on the physico-chemical properties of catalysts, active metal and catalysts supports. Noble metals in catalysts providing a high selectivity to hydrogen are Rh, Ru, Pd and Ir, while among transition metals the best performance ensure Ni, Co and Cu [20]. The order of activity for these metals is Ru > Rh > Ni ~ Ir > Pt > Pd [21]. Despite the high activities and the low tendency of the noble metals to carbon deposition, costs and limited availability of these metals prevent widespread industrial use. Transition metals have shown themselves to be more promising due to their relatively low cost. However, transition metals are more susceptible to coking and deactivate faster than noble metals. On the other hand, the introduction of rare earth oxides into the catalyst support can improve both activity and stability. It is also known that CeO2 is effectively used in the development of catalysts for steam reforming of ethanol to produce hydrogen [22,23,24,25,26,27]. Moreover, lanthanum oxide has been reported as an efficient additive to supports providing to catalysts a higher selectivity and stability [28]. Other catalyst types based on doped zirconia [29], alumina [25], zeolites [30], etc. supports [31,32,33] are used in ethanol steam reforming as well.
In the practical design and operating decisions, such multifunctional reactor configuration consists of two reactor volumes, which are the reaction (feed-side) compartment followed by a permeate sweep-side zone. A functionally selective diffusion barrier separating reactor compartments enables hydrogen produced by a catalytic reaction to migrate into the side of a lower hydrogen concentration and be carried out by a sweeping gas, while the retentate stream leaves the feed-side compartment. Mass transport across a membrane module may occur under induced driving forces. The driving force for hydrogen transmembrane permeation is a hydrogen partial pressure gradient at the sides of the membrane. A variety of membrane separation processes are often classified according to their driving forces [34,35,36].
In developing MRs with an integrated catalytic process, the effective characteristics and architecture of the catalyst and the membrane module under identical operating conditions are of utmost importance. Characteristics of a catalyst, such as its design and levels of intrinsic activity and selectivity sufficient to compete with membrane efficiency, are essential for creating a high hydrogen gradient all over the membrane module. The competition between the different phenomena determining the performance of the system, namely mass transport, reaction in the catalytic bed and hydrogen permeation across the membrane may occur. Morever, the catalyst and membrane need to be adapted to provide a low overall resistance to mass transport of hydrogen. Among factors influencing the net driving force of hydrogen transport is the fluid dynamics, which is highly dependent on the chosen experimental reactor configuration and operating conditions [37,38,39].
Clearly, a thorough understanding of all factors that could affect interrelated phenomena of mass transport, catalytic activity and permeation is required when designing a MR. Moreover, these types of reactors may inherently have serious drawbacks such as deactivation of catalyst and membrane due to coke formation, limitations for mass transfer and difficulty in replacing inside materials, etc. [40,41]. A fixed bed of catalyst is the most common approach to the design of reforming reactors for hydrogen production, whereas it suffers from poor heat and mass transfer behavior, high-temperature gradient, catalyst sintering, coke deposition and dust jamming. Catalytic packed beds consisting of particles (typical sizes 1–10 mm) make possible some limitations to hydrogen transport between the bulk of the catalytic bed (where hydrogen -rich gas is produced) and the membrane surface, thus lowering the driving potential for the transmembrane permeation [10,42,43].
Previously, a packed bed membrane reactor for production of hydrogen through ethanol steam reforming has been successfully tested at ambient pressure and the temperature range of 700–900 °C at a laboratory scale [44]. The asymmetric membrane disk module consisted of a gas-tight nanocomposite functional coating (Ni + Cu/Nd5.5WO11.25-δ mixed proton-electron conducting nanocomposite) deposited on a gas-permeable functionally graded foam substrates. A packed bed consisting of spherical catalyst particles of 1 mm diameter was placed on the top surface of the membrane module. The best results have been obtained at 900 °C and feed of ethanol/H2O mixture in Ar at steam-to-ethanol ratio of four. The overall hydrogen flux was achieved to be about 1.31 Nml cm−2 min−1.
By using the one-dimensional isothermal reaction–transport model for the constituent layers in the experimental reactor, numerical experiments had been performed to elucidate an impact of the structural parameters of composite membranes on the interrelated catalytic and permeation phenomena in the MR with a packed bed. It was shown that the asymmetric support contributes up to 70% to the overall resistances across the membrane module at the feed gas flow rates of 3 and 10 Nl h−1. Transmembrane transport was mainly controlled by chemistry in the catalytic layer located on the surface of the membrane [45].
It is known that structural catalysts (honeycomb monoliths as well as microchannel plates) with a thin catalytically active layer of about 0.01 mm thickness washcoated over metallic substrates provide both a low pressure drop and a low pore diffusion resistance, hence, a high activity. The mini-channels with characteristic diameters between 400 μm and about 1 mm have large surface-to-volume ratios (catalytically active surface area per the catalyst unit volume), which result in superior transfer properties and, consequently, offer a compact and modular solution for the devices. The excellent thermal conductivity of metallic monoliths provides more uniform temperature profiles along the catalysts length/diameters. Advantages such as the ability to intensify catalytic processes by increasing heat and mass transfer or by more precise control of contact times can also improve process efficiency, arising from a specific flow regime that occurs in small channels—nearly plug flow behavior [46,47].
The same membrane module was applied in a lab-scale MR, in which a structural catalytic monolith was installed. In order to evaluate the efficient production of hydrogen and the permeability of the membrane, intensive studies were conducted under various operating conditions concerning temperature, fuel concentration (ethanol and water in Ar) and various molar ratios between ethanol and water [48]. The catalytic monolith consisted of a FeCrAl substrate with corundum protective coating and an Ni + Ru/Pr0.35Ce0.35Zr0.35O2 active component. The catalytic monolith was placed above the membrane module and separated from the dense functional layer by a gap junction. Functionally, the gap allowed for the feed stream to flow in a direction perpendicular to the membrane surface, while retentate flow was discharged tangentially along the surface. The experiments were carried out with constant flow rates of 5 Nl h−1 for the feed gas (ethanol and water mixture in argon) and of 10 Nl h−1 for the Ar sweeping gas. The concentrations of both ethanol and water were systematically varied from 3 to 30 vol. % and from 13 to 80 vol. %, respectively, at the steam-to-ethanol molar ratios S/E = 2, 4, 6. It has been shown that the monolithic catalyst-assisted MR is capable of increasing the driving potential for hydrogen permeation through the same membrane as compared with that of the packed bed catalyst. The maximum flux of 3.03–3.5 Nml H2 cm−2min−1 was obtained at 900 °C, which was nearly independent of the S/E ratio while in the membrane reactor with a catalytic packed bed; this characteristic was calculated to be about 1.31 Nml cm−2 min−1 [44]. Additionally, an increase in the fuel concentration (ethanol and steam) did not result in improving reactor performance characteristics. Thus, in this mode of operation, hydrogen permeation was limited by passage through the membrane itself. Namely, the diffusion flux through the membrane was no longer governed by the driving partial pressure force, but rather by the ability of the membrane itself to hydrogen mass transfer. The observable phenomenon is a final result of the interplay between catalysis, kinetics, transfer effects and design and operation of the reactor. It makes us curious to understand in detail unobservable critical elements that have important roles to play in the effectiveness of the reaction–separation process in the MR.
The research concept of this paper, based on the intensive experimental study, is to investigate, by conducting a model-based performance analysis, all the main phenomena that determine the performance of the monolithic catalyst-assisted MR, namely design parameters, chemistry, transfer effects and hydrogen permeation across the membrane layers, as well as boundary layer effects. As far as we know, we have not found such a study in the available literature. Simulations and numerical methods are the powerful tools for scientific explanations and predictions, being also an alternative to the techniques of experimental science and observation, especially in cases when phenomena are not observable or when measurements are impractical or too expensive [49].
The computational one-dimensional reaction-transport isothermal model has been implemented in COMSOL Multiphysics. The mathematical model with reasonable parameter values was then verified with experimental data and shown to provide good agreements. Our findings also expand our understanding of the role of gap junction in the reactor, showing it is essential for hydrogen transmitted to the sweep-side compartment. Most significant results are detailed hereafter.

2. Reactor Model

Accurate and reliable modeling of a MR calls for a thorough understanding of the processing phenomena that are specific to the particular reactor configuration. Subsequently, in order to reflect the principal features in the mathematical description as accurately as possible the following important information should be provided: the size and general configuration of the reactor and the more important functions and dimensions of its internal structures, operating conditions and the composition of the product emerging therefrom. In this study, we used the following approach: fluxes through all constituent layers were coupled by first defining the concentrations at the interface and the flows on the boundaries inside all assemblies, and then by combining these interface concentrations. Molar flow change due to reaction and membrane transport is accounted for. The developed model is validated by confronting its predictions with data from experimental studies [48].
Moreover, the experimental data have been used to obtain and estimate all necessary parameters to quantify the characteristics for the developed mathematical model.
To capture concentration distributions in axial direction of the experimental MR, a computationally efficient pseudo-homogeneous, one-dimensional isothermal reactor model was developed consisting of the total gas-phase continuity and differential molar balances on species in the constituent layers.

2.1. Configuration of Reactor

Detailed information about the design of the experiments, including the experimental setup, the purpose of experimentation, the morphology and the operation of the developed asymmetric supported membrane, consisting of a thin Ni–Cu alloy–Nd tungstate nanocomposite dense permselective layer deposited on a hierarchically structured asymmetric support, has been reported in previous papers [44,45]. The main assumptions performed in the simulations of the MR with the monolith catalyst are summarized below: design, kinetics, and transfer effects.
Detailed information about the experimental setup and disk-shaped reactor with axis-symmetric flows with ethanol steam reforming over a metal honeycomb catalyst has been reported in detail in the previous paper [48]. The catalytic monolith is also, as in the case of the packed bed membrane reactor [45], located above the membrane module. A sketch of the MR is shown in Figure 1, showing all essential constructive details of the experimental reactor that is under present numerical study.
The use of a well-insulated cylindrical furnace ensures constant temperature in the reactor system. In this reactor configuration, the monolith and membrane surface are linked via a gap junction. It is the main difference between the monolith bed and the packed bed membrane reactors. The constructive gap, which separates catalyst and membrane, allows for the hydrogen-rich gas from the catalyst to flow in a direction perpendicular to the membrane surface, while retentate flow is discharged tangentially along the surface out of the reactor system. Hence, the experimental reactor utilizes a stagnation flow configuration in the constructive axial distance between the honeycomb catalyst and frontal surface of the membrane.
The operation at the water/ethanol molar ratio of 4 in the fuel mixture was found to be beneficial for ethanol conversion and hydrogen separation within an operating temperature range of 700–900 °C. These operation parameters were taken as a reference case in the simulations.
To verify a developed model and simulate the reactor behavior for different operational cases, it is necessary to have a reference case. The general parameters for the reactor and experimental conditions (reference case) are given in Table 1.

2.1.1. Catalyst Bed Configuration

The performance of the catalytic monolith is to a large extent determined by its structure, morphology and porosity. A metal support is formed by stacking flat and corrugated foil (Kanthal FeCrAl alloy) bands (each approximately 120 µm in thickness) and winding them into an Arkhimedes spiral (Figure 2). The corrugated foil had a wave height of 0.7 mm at a pitch of 2.5 mm, so the entire cylindrical honeycomb substrate with outer diameter of 24 mm had 12 layers of flat and corrugated foils in the cross section. The support was coated by γ-Al2O3 (10 wt. %) followed by the catalytically active composition, 10 wt. % Pr0.35Ce0.35Zr0.35O2 + 5 wt. % Ni + 1 wt. % Ru as described elsewhere [48]. The geometrical properties of the catalytic monolith together with the feed gas composition for the reference case are listed in Table 1.

2.1.2. Design of Membrane Module

The morphology and geometry of membrane module were discussed in detail in the previous paper [45]. The asymmetric supported membrane consists of a gas-tight Ni–Cu alloy–Nd tungstate nanocomposite dense functional layer deposited on a gas-permeable Ni-Al hierarchically structured asymmetric substrate (Figure 3) [44]. The Nd5.5WO11.25-δ (NW) powder synthesized by the mechanical activation as described elsewhere [50] was characterized as a single-phase defect fluorite. The Ni–Cu alloy was synthesized by modified Pechini route using a fluidized bed reactor for subsequent thermal treatment of products of polymeric precursors decomposition in Ar + H2 streams. The obtained nanoparticles were put into isopropanol to avoid oxidation. The Ni–Cu (30 wt. %)/Nd5.5WO11.25-δ nanocomposite was prepared by ultrasonic dispersion of Ni–Cu alloy and NW powder suspension in isopropanol with the addition of polyvinyl butyral. The membrane parameters are given in the Table 2.
Powder Metallurgy Institute (Minsk, Belarus) has provided functionally graded substrates. The method for making gas-permeable structures involves applying two thin low-porosity layers onto a thick Ni/Al foam. The tough foam is made from alumina-silica ceramics. In this process, 30 PPI polymeric foam is used as a template. Macroporous ceramics is obtained by 3-times impregnation of foam polyurethane with alumosilicate suspension followed by centrifugation, drying and sintering at 1350 °C. In addition, Ni covering of alumina-silica ceramics is done by electrolysis followed by drying and sintering at 1000 °C in cracked ammonia. Ni:Al2O3 weight ratio is obtained to be 2-2.2:1. Formation of two thin layers onto the Ni/Al substrate have been done by the use of dual doctor blades in series. The compaction of layers has been performed via sintering of the specimen at 1000 °C in cracked ammonia followed by in-pack aluminizing as described in Ref. [51]. The resulting gas-permeable substrate has a composite structure, in which the upper porous layer has pores of a smaller diameter than the intermediate porous layer on the foam substrate. In terms of pore size distributions and porosity, the three-layer structure of the gas-pearmable support is clearly quantified and designated in Table 1. Quantification of the constituent particles and pore systems has been done via the image analysis with a Joyce–Loebl Mini Magiscan (Joyce-Loebl, Ltd., Gateshead, UK) computerized image analysis system [52,53], and with a vector program, by applying an appropriate spatial calibration in a series of parallel and perpendicular lines along the image axes in the reference area.

2.1.3. Gap Flow Characterization

The reactor configuration introduces a space between the catalytic monolith and the membrane dense surface, in which case the hydrogen-rich product stream from the catalyst layer must pass across the gap before it reaches on the membrane surface. That is, the dense surface of the membrane is perpendicular to the flow, thereby creating a stagnation region, in which the feed flow velocity vanishes, and the gas continues flowing out radially outward to the exhaust through an annular space around, causing the formation of a boundary layer adjacent to the stagnate surface. By convention, the value of the boundary layer thickness is defined as a height needed to obtain 99% of the mainstream velocity.
It is known that the stagnation region encounters the highest pressure, the highest transfer rate and the highest rates of mass exchange and deposition to a flat surface [54,55]. Under uniform distribution and effective and stable contact with the membrane surface, hydrogen from the product gas selectively permeates into the sweep-side compartment. The axisymmetric stagnation flow geometry is a desirable flow scenario for controlled measurements [56,57], as it ensures that an active surface sees at least the same gas-phase composition. The constructive and the operating parameters in the stagnant flow arrangement, affecting the experimental information, were investigated quantitatively by F. Zanier et al. [58]. It was shown that the flow structure may even become the controlling factor determining the global activity of disk-shaped, flat active surfaces.
Since the gap is full of stagnant hydrogen-rich gas, there should be a reduction in the mass transfer resistance and, thus, a higher hydrogen flux compared to membrane reactors with packed beds. Understanding whether our experimental evidence supports the conclusion is a goal of the present paper.
In the modeling approach, a steady similarity numerical solution has been applied to describe the axisymmetric stagnation fluid flow toward the membrane.

2.2. Model Description

The mathematical model of the membrane reactor utilizes a stagnation flow configuration in the gap between the catalytic monolith and dense surface of the membrane module. The molar balances comprise axial mass transfer by convection and axial dispersion as well as mass transfer through the phase interfacial surface area.
The following assumptions have been made in the formulation of the mathematical model:
  • Both the reaction and transport are conducted isothermally. This assumption for the endothermic reaction of ethanol reforming signifies that sufficient heat is supplied to maintain an approximately constant reaction temperature.
  • The internal mass transfer resistance of the catalytic washcoat over the metallic monolith is neglected.
  • Negligible pressure drop occurs across the membrane reactor.
  • The permeation of hydrogen through the dense layer of the membrane module follows Sievert’s law.
  • Molecular and Knudsen diffusions control in the hydrogen and argon transport through the gas-permeable layers of the membrane support.
  • A perfect-gas equation of state relates the density, pressure, temperature and composition.
Important constitutive equations for the reaction kinetics, membrane flux and dispersion coefficients are discussed in following sections. A set of equations for the feed-side is coupled through the hydrogen permeation flux depending on the partial pressures on both sides of the dense permselective layer of the membrane module. As there are not significant pressure gradients, the gas flow has been modeled considering a constant gas density in each layer of the membrane reactor (Table 1).

2.2.1. Feed-Side Compartment Model

The model equations and the requisite boundary conditions have been listed in Table 3 and Table 4.
It is believed [59,60] that in reactors with the separation option, reduced Reynolds numbers (Re) are desired, indicating highly ordered laminar flow and a low Péclet number (Pe of less than or around 1), when diffusive transports between phases are preponderant. Such is the case with the experimental MR operating conditions for the reference case (Table 4). Thus, at T = 900 °C, the Reynolds number is about 1800, while the Péclet number is 0.66, assuring laminar flow conditions. Chemical reactions result in production or consumption of species (Ri), which is modeled as a molar source or sink for the i-th specie. A continuity equation for species molar balance in the feed-side compartment contains partial derivatives of molar flow rates and species concentrations with respect to the axial position.

Reaction Description

The steam reforming of ethanol is assumed to include a stepwise reaction scheme, see Table 4 [61]. The observed reaction rates depend on the concentration of the reactants, the temperature and the catalyst external surface area. Effective parameters and variables have been evaluated in an earlier experimental study of ethanol steam reforming performed over the 5 wt.% Ni + 1 wt. % Ru/Sm0.15Pr0.15Ce0.35Zr0.3O2-δ catalyst [44].

Dispersion of Mass

It is well known [59] that the complexity of the molecular transport processes does not allow a purely theoretical fundamental approach in analysis of the diffusive transport. In the numerical simulations, transport coefficients (the diffusion coefficients, viscosity coefficients etc.) are calculated from the transport coefficients for the individual species (Table 5). For laminar flow in tubular pipes (Re < 2100), Taylor [62] and Aris [63] established a relation between the molecular diffusion coefficient, the hydraulic diameter and the average velocity in the channel. For the contribution of molecular diffusion, an effective diffusion coefficient was calculated with the Wilke equation [64], where binary diffusion coefficients Dij ([m2/s]) (Table 5) were calculated by using Fuller-the Schettler-Giddings empirical Eq. [65]. There is considerable uncertainty about hydrodynamics in the feed-side and effective diffusivities of species, which warrants the use of fitting parameters for dispersion coefficients, when agreement with experimental data is desired.

Stagnate Flow in the Gap

The monolithic catalytic structure imposes a well-defined plug-flow condition for the product gas at the outlet F f | z = h m t h + , which stagnates on the surface of the dense layer. The thickness δ h , m of the laminar boundary layer developed over a flat dense surface is the result of the interplay between convective and diffuse fluxes in the gap. That is, diffusive effects only become important in a thin region near the membrane surface. The change of hydrogen concentrations on each side of the stagnant boundary layer is influenced not only by driving force, but also transport properties in the membrane itself. Hence, the feed composition is changed in this intermediate region resulting from the diffusion of hydrogen to the membrane surface across the boundary layer, and the concentrations decrease from their feed stream values to that of the retentate composition-controlled values.
The retentate product F r depends on the net production rates of the species produced ( F f ) z = h m t h over the monolithic catalyst and hydrogen flux permeated ( F H 2 , p e r m ) through the membrane (Table 3). When a concentration gradient exists, the species tend to flow in a direction such as to reduce the concentration gradient. The mass transfer coefficient β, also known as the transfer velocity, is used to characterize transfer of a substance through another on a molecular scale due to concentration difference or gradient. Mass transfer rate is dependent on the thickness of the boundary layer.
The interfacial effect on diffusion resistance in the stagnation region due to changing geometrical and operational parameters should be introduced obviously. Hiemenz [66] first studied the steady flow in the neighborhood of a semi-infinite wall. For the axisymmetric stagnation fluid flow, a steady similarity numerical solution has been applied quite often [54,55,67]. A measure of a change of flow conditions near the stagnation surface is determined as a stagnation velocity gradient u [ 1 s ] [68,69]. Expression for the boundary layer thickness in a stagnation-flow given in Table 5 as the square root of the reciprocal value of the velocity gradient [70] takes into account a change velocity with the respect of the gap geometry and gas flow rate.
The mechanism by which hydrogen diffuses through the boundary layer in the gas mixture due to the concentration gradient is exactly analogous to that by which heat diffuses. By adopting the approach of non-dimensional analysis, the mass transfer coefficient β is defined in a manner analogous to the heat transfer coefficient in terms of the Reynolds and Schmidt numbers. A theoretical correlation for predicting the heat transfer near the forward stagnation point assuming laminar, incompressible, low-speed flow is used [68]. When used in connection with mass transfer, the Prandtl number is replaced by the Schmidt number, Sc, which expresses the ratio of the momentum diffusivity to the mass diffusivity. The stagnation Reynolds number, Re, was based on the diameter of the impingement membrane surface, dm, and the stagnation velocity gradient.

2.2.2. Sweep-Side Reactor Model

Dense Layer

Mass transport through dense membranes is a purely diffusive process [71]. Permselectivity is the intrinsic property of the active material, while the driving force for hydrogen transport across the dense membrane is caused by the difference of hydrogen partial pressure on each side. Sieverts’ Law, which identifies the difference of the square roots of the partial pressure of the hydrogen as the the driving force of the permeation, is a temperature-activated phenomena. It is valid with an underlying assumption that the surface coverage is low when interfacial equilibrium is achieved and the rate-limiting step is an atom diffusion through the membrane. These criteria are satisfied in most cases when the temperature is relatively high [72]. Obviously, the steady state of the hydrogen flux through the dense layer of a thickness hdm being described by Sieverts’ law is equal to the permeate flux through the adjacent boundary layer (parity–flux equation in Table 6). The permeation rate is limited by the concentration of hydrogen on the gas phase /dense layer interphase x H 2 , d m . An explicit solution for x H 2 , d m is attained from the parity–flux equation (Table 6).
It is known that values for the hydrogen flux and the hydrogen permeation coefficient can be estimated more precisely only for given conditions. Therefore, for predicting the hydrogen flux, the apparent gas permeability is useful to quantify for certain operating regimes [73,74]. Corresponding parameters in the expression for the hydrogen permeation flow rate ( F H 2 , p e r m ) , the apparent activation energy, and the permeability were obtained from the temperature-dependent conductivity study [75] and the fitting of this equation to the reference case experimental data [48].

Gas-Permeable Layers

The mathematical model for the porosity graded substrate was described in detail in a previous article published by the author [45]. The gaseous mixture inside the gas-permeable layers is a binary mixture of argon and hydrogen. The hydrogen gas transport inside is governed by Fick’s diffusion with effective binary diffusion coefficients. The linkage between the layers requires the equality of the fluxes at each boundary (Table 6). The effective pore diffusivity in the foam layer is taken as a function of the structural parameters of the heterogeneous porous medium. Hydrogen diffusion resistance accounted for with the specific mass transfer correlation for the foam configurations [76,77,78].

Sweep Compartment

Using the ideal gas law and assuming a perfectly mixed gas, the molar flow rates for the species H2 and Ar in the sweep gas compartment are related to the volumetric gas flow accounting for a change in the total molar flow rate of the sweep gas at a given temperature (Table 6). The unknown variable x H 2 , s w , the molar fraction of the hydrogen in the permeate flow leaving the sweep compartment, is yielded at steady state to be the real solution of a quadratic.

2.3. Numerical Solution Strategy

The scaled reactor model equations for the constituent layers of the experimental membrane reactor equipped with the catalytic monolith (Table 1) with the boundary conditions have been implemented in COMSOL Multiphysics 4.3a software. The numerical model has been executed and solved simultaneously in the domains along the full length. The finite-element method for numerical solutions of differential equations being used employs a uniform (fine) mesh and error control in the domains. A numerical problem arises from the initial values in the feed, which generates a division by zero in the reaction rate equations. This problem has been sorted out by using very small values for the mole fractions of the generated species at the inlet. The solver, an implicit time-stepping scheme, is well suited for solving stiff and nonstiff nonlinear boundary value problems. At the end of the solving process, species concentrations and fluxes are known at each axial location within every reactor’s domains. The physical properties are determined by the local composition and temperature, which are affected by the chemical reactions and the permeation flux. A fast and unconditionally stable solution is provided, when solving the species molar balance.

3. Results and Discussion

3.1. Base, Reference Case

3.1.1. Verification of the Model Formulation and Computational Simulations

The feasibility of performing ethanol steam reforming over a monolith in the lab-scale MR in isothermal mode has been investigated extensively at different operating conditions regarding the temperature, the fuel concentration and the different molar ratios between ethanol and water [48]. An isothermal condition in a reactor is the most ideal mode of operation, because the constant temperature along both a catalyst and a membrane is advantageous for the life duration and stability. The reactor operation in the temperature range of 700–900 °C with the water/ethanol molar ratio of 4 in the feed was found to be quite sufficient in terms of the ethanol transformation into syngas and CO2 as well as the amount of hydrogen produced from 1 mol of ethanol. These operation parameters were taken as a reference case in the simulations.
The most important question in any simulation is the reliability and validity of an assessment tool and whether we can get more out of our model by tuning the model parameters. The applicability of the one-dimensional model formulation and computational simulations are verified and validated by direct comparison of model results with existing measurements of interest. In Figure 4, the results of the simulation for the hydrogen flux and concentrations of species are compared to the known values in the experimental study of the reference case. The operating conditions and general parameters for the reactor internals for the base (reference case) have been listed in Table 1.
Under the examined operating condition, correction factors were applied to the theoretically calculated beta values (feed-side mass transfer coefficients in Figure 5) in order to match satisfactorily the experimental data in the studied temperature range with the same magnitude of permeability ( Θ = 1.4801 × 10 2 mol m−1s−1atm−0.5) used in modeling.
Consistent with the previous results of modeling the MR with a packed bed catalyst, our simulations are in reasonably good agreement with the experimentally deduced species concentration in both the feed- and sweep-side and the hydrogen permeation flux. In particular, by considering the effect of the temperature shown in Figure 4, the average error between the experimental and the modeling data for the hydrogen fluxes is about 5%, while a highest deviation of about 8% is observed for hydrogen concentration in the retentate gas at 700 °C.
The concentration distribution along the reactor internals affects the overall hydrogen permeation flux, which, in turn, depends on mass transport parameters: reaction kinetics, Péclet number and effective diffusivity in the constituent layers [79,80]. Understanding all these phenomena is important when studying factors that direct the hydrogen removal.

3.1.2. Performance Analysis

Generally, comparative analysis of MRs is performed by using characteristics of recovery, yield of hydrogen, fuel conversion and resistances to internal mass transfer of a permeating component. The details of introducing a structural catalyst in the MR reactor instead of a packed bed on the overall efficiency of the integrated reaction–separation process were discussed in detail in the previous paper [48]. It was shown that the reactor operating with the catalytic monolith showed better performance in terms of both hydrogen recovery and the yield with respect to the packed bed catalyst. Thus, a comparative study showed that in the case of a packed bed, an increase in temperature from 700 °C to 900 °C leads to an increase in the yield of hydrogen by almost 40% (from 24% to 33%), while the yield for the monolith practically does not change and averaged 37%. Hydrogen recovery for the monolith was about 10% higher in all operating temperatures.
The resistance concept proposed by J.M. Henis and M.K. Tripodi [81] gives insight and understanding about which of the reactor internals is controlling the total flux. Irrespective of the transport mechanisms, a total transport resistance consists of the layers resistances in series and includes also resistances by the gas film on both sides of the membrane assembly [45,82,83,84]. A layer with the highest concentration gradient would provide the major mass transfer resistance to the overall permeation process.
Typical concentration profiles, as those calculated for the reference case at T = 800 °C, Figure 6a, evolve throughout the reactor length at isothermal operation. First, the hydrogen concentration increases up to a maximum value (29.8 vol. %) due to the high ethanol reforming reaction rate and then, after its maximum value, the concentration gradually decreases down towards the membrane surface up to 23 vol. %. The concentration gradient of about 3.6 vol. % is established across the boundary layer adjacent to the membrane. The percentage values of hydrogen recovery and yield estimated by experiments at 800 °C give 47% and 44%, respectively.
The hydrogen removal results in an increase in the overall mole fraction of other species and a depletion of hydrogen in the retentate gas. This causes a concentration gradient build up in the boundary layer. The polarization reduces the overall efficiency of a membrane separation [79]. The polarization resistance of the membrane assembly at different temperatures in the reference case are shown in Figure 6b. The slope of the curves is mostly dependent on the temperature and on the permeability of the membrane layers. High temperature is beneficial to syngas formation and membrane permeability. When the operating temperature increases from 700 up to 900 °C, the hydrogen flux increases by ~35% while the overall resistance to the hydrogen mass transfer decreases by ~38% as the temperature increases. The values with high impact on the overall efficiency are at least in the same order of magnitude. It means that resistance values have to be estimated in a simulation procedure. The derived model allows describing resistances to the hydrogen mass transfer as a function of the operating conditions and the structural properties of each single layer.
The MR components that contribute to the hydrogen permeation resistance performance are boundary layers from both sides of the membrane module and its constituents —dense, powder, intermediate and foam layers. The contribution to the overall resistance of the components is illustrated in Figure 7.
It can be seen that the hierarchically structured asymmetric support, which is required for desirable industrial applications of membranes, demonstrates a similar extent of resistance to hydrogen transport in the constituent gas-permeable layers for both the configuration of the catalyst and for the different feed-flow rates, while the variable contribution of the permselective layer is observed. In the stagnant flow condition, molecular hydrogen is able to pass easily to the dense layer of the membrane. However, the transmembrane potential, which can be characterized by the concentration gradient across the membrane and by the intrinsic permeability, fails to permeate more flux. Indeed, depending on the conditions at the entry (feed side) and exit (sweep side) faces of the membrane, a concentration gradient establishing through the membrane module gives rise to the membrane permeation. Yet, the maximum flux was detected at 900 °C to be about 3.03–3.5 Nml H2 cm−2 min−1 [48]. Thereby, the stagnation flow mode is supposed to elicit a limit of the membrane permeation, when the diffusion flux through the membrane is no longer governed by the driving partial pressure force, but rather by the intrinsic features of its constituents. Thus, the permeation rates are higher for the monolith-assisted MR than for the cases of the packed bed and limited by the permeation ability of the dense membrane layer.

3.2. Parametric Study

As it is drawn in Figure 6a, the nonvanishing gradients of concentrations in the gap are restricted to the boundary layer. This is due to the convective nature of the mass transport in the gap. The gas streams toward the dense layer and the permeation process at the surface have no influence on the fields far away. Species move with the flow, which dominates over any non-convective transport. In this flow configuration, only when the axial velocity approaches zero near the surface, the diffusion becomes of similar size as the convective transport [58]. The remaining retentate products are transported away sufficiently quickly.
We modeled various distances between the monolithic catalyst and the membrane surface, as well as feed rates, to determine their effect on permeability characteristics. As expected, the thickness of the boundary layer affects the hydrogen diffusion rate through it. Figure 8a, left panel, shows that by enlarging the gap above the membrane surface from 4 to 15 mm, the thickness of the boundary layer increases from 1.28 mm to 2.48 mm at 700 °C, and from 1.46 mm to 2.83 mm at 900 °C. This leads to a decrease in the mass transfer velocity of hydrogen β (at 700 °C from 0.35 to 0.09 m/s, at 900 °C from 0.81 to 0.21 m/s) with a corresponding decrease in the hydrogen flux (about 3–4%). Indeed, a smaller distance between the monolith and the membrane can lead to compressing the boundary layer and, thereby, can provide better diffusion of hydrogen into the dense layer.
Higher feed flow rates have the same effect on the boundary layer: the thickness of the boundary layer on feed-side decreases. However, this leads to lower values of both the hydrogen flux and hydrogen concentration in the sweep compartment, as shown in Figure 8b, due to an increase in the hydrogen concentration in the retentate flow coming out of the gap space. The higher the feed rate, the higher the negative effect of a larger gap height above the membrane surface can be.
The higher hydrogen flux in the case of the MR with a monolith may also be due to the influence of the boundary layer phenomenon after the impingent of the feed stream (product gas from catalytic monolith), which affects the corresponding driving forces of interfacial hydrogen transfer rate in the membrane. The effect of the operating temperatures on the driving force, expressed here as the difference between the concentrations on each side of a boundary layer adjacent to the membrane surface and the magnitude of the corresponding permeation fluxes, is shown in Figure 9. It can be seen that in the case of the stagnate boundary layer, the driving force is about an order of magnitude higher, and accordingly the permeation flux is 30–40% greater. It is obvious that the hydrogen flux through the membrane is dictated by the number of molecules impinging and interacting with the interfacial surface under the boundary condition being formed.
At the same time, the stagnate boundary layer phenomenon must be properly taken into account for the specific reactor configuration and operating conditions; otherwise incorrect conclusions about the technological parameters will be drawn. For example, if the feed flow is moving at a very high speed, it may have some kinetic energy, and this may lead to undesirable effects on the penetration rate.

4. Conclusions

In this study, the conceptual feasibility of a membrane reactor equipped with a catalytic monolith for ethanol steam reforming has been investigated by a detailed numerical simulation using a 1-D reactor model. The formulated mathematical description takes into account the design parameters of the MR, the chemistry and transfer effects in the monolith, the hydrogen permeation across the membrane constituent layers together with a Sieverts’ equation for the functional dense layer as well as the boundary layer effects. The scaled reactor model equations were implemented in a COMSOL Multiphysics environment. Good agreement with the experimental data of the lab-scale MR with reasonable parameters values is provided. In the numerical experiments for the MR operating at atmospheric pressure and in the temperature range of 700–900 °C, the concentration profiles along the reactor axis were obtained, showing the effect of the emerging concentration gradient in the boundary layer adjacent to the membrane.
A parametric study was carried out to determine the influence of the distance between the catalyst monolith and the membrane surface, as well as the feed rate, in order to determine their influence on the permeability characteristics. Polarization resistance, which reduces the overall membrane separation efficiency, has also been elucidated for the reference case and compared with the packed bed MRs.
The main findings are as follows:
  • the use of a catalytic monolith with stagnant flow between the catalyst and membrane surface may increase both the production and flux of hydrogen, as well as reactor efficiency characteristics compared to a packed bed MR;
  • the stagnate feed flow configuration must be adequately taken into account for the specific reactor design and the operating conditions; otherwise incorrect conclusions about the technological parameters will be derived. Thus, increasing the distance between the monolith and the membrane surface from 4 to 15 mm increases the thickness of the boundary layer from 1.28 mm to 2.48 mm at 700 °C, and from 1.46 mm to 2.83 mm at 900 °C with a corresponding reduction in the hydrogen flux (about 3–4%);
  • at higher feed flow rates, the thickness of the boundary layer on the feed-side of the membrane decreases. However, in this case, lower values of both the hydrogen flux and hydrogen concentration in the sweep compartment are obtained, due to an increase in the hydrogen concentration in the retentate flow leaving the gap space. The higher the feed rate, the higher the negative effect of a greater gap height above the membrane surface can be;
  • the hierarchically structured asymmetric membrane support, which is often necessary for industrial applications, demonstrates a similar extent of resistance to hydrogen transport in the constituent gas-permeable layers of the membrane used in the study of MRs with monolith and packed bed catalysts and for different feed-flow rates. A variable contribution of the permselective layer is observed, and the permeation rate in the MR with the monolith is limited by the permeability of the dense membrane layer.
This study proved that catalytic monoliths could be effectively used inside MRs to produce hydrogen. Comparative studies of hydrogen mass transfer through the same membrane module using an in series model resistance showed a stronger effect for the monolith-assisted MR with stagnant flow over the membrane when compared to a packed bed. In this case, all mass transfer resistance will lie in the membrane module itself, especially in the permselective layer, and almost not be in the boundary layers.

Author Contributions

L.B. designed the model, carried out simulations and wrote the paper; N.E. contributed analysis of samples and experiments; N.V. scaled mathematical model and implemented in COMSOL Multiphysics 4.3a software; V.S. supervised this research study and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The studies were supported by the Ministry of Science and Higher Education of the Russian Federation within the governmental order for Boreskov Institute of Catalysis (projects AAAA-A21-121011390007-7, AAAA-A21-121011390009-1).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Ameffective area of the permeating surface (m2):
Amthcatalytic monolith cross section area (m2);
dcellcell diameter of the foam layer in the membrane (m);
dpore,pl(il)hydraulic pore diameter in the powder (intermediate) layer (m);
dp,foamaverage diameter of the windows which connect the cells (m);
ds,avgaverage cylindrical strut size (m);
Dijbinary diffusion coefficients (m2 s−1);
dh,mthequivalent channel diameter of the catalytic monolith (m);
F f 0 inlet feed molar gas flow rate (mol s−1);
F r retentate product flow (mol s1);
hdmthickness of the dense layer in the membrane module (m);
hfoamthickness of the foam layer in the membrane module (m);
hgapheight of the gap with stagnate flow (m);
hilthickness of the intermediate layer in the membrane module (m);
hmthheight of the catalytic monolith (m);
hplthickness of the powder layer in the membrane module (m);
J H 2 hydrogen permeate flux through the asymmetric membrane (mol m−2 s−1);
Mimolecular mass of i-th component (g mol−1);
Mfaverage molecular mass (g mol−1);
P H 2 , d m hydrogen partial pressure at the catalyst/membrane interface (Pa);
Qdmmembrane permeance (mol m−2 s−1);
ReReynolds number;
ScSchmidt number;
SV,foamfoam layer of the membrane volumetric surface area (m2 m−3);
SV,ilintermediate layer of the membrane volumetric surface area (m2 m−3);
SV,mthcatalytic monolith volumetric surface area (m2 m−3);
SV,plpowder layer of the membrane volumetric surface area (m2 m−3);
Toperation temperature (K or °C);
ustagnation velocity gradient (s−1);
x H 2 , d m molar fraction of hydrogen the catalyst/dense layer interphase;
x H 2 , f molar fraction of hydrogen in the feed-side (-);
xi,fmolar fraction of gas species in the feed-side (-);
Greek letters
βeffective mass transfer coefficient (m s−1);
εfoamporosity of the foam layer in the membrane module (-);
εilporosity of the intermediate layer in the membrane (-);
εmthporosity of the catalytic monolith (-);
εplporosity of the powder layer in the membrane module (-);
Θpermeability (mol m−1 s−1 Pa−0.5);
mfdynamic viscosity of the feed mixture (kg s−1 m−1) or (Pa s);
mfoamdynamic viscosity of the gas mixture in foam layer (kg s−1 m−1) or (Pa s);
mIdynamic viscosities of a specie (kg s−1 m−1) or (Pa s);
ρffluid density in the feed side (kg m−3);
ρtottotal molar density (mol m−3);
τfoamtortuosity of the foam layer in the membrane module (-);
τiltortuosity of the intermediate layer in the membrane module (-);
τpltortuosity of the powder layer in the membrane module (-).

References

  1. Osman, A.I.; Mehta, N.; Elgarahy, A.M.; Hefny, M.; Al-Hinai, A.; Al-Muhtaseb, A.H.; Rooney, D.W. Hydrogen production, storage, utilisation and environmental impacts: A review. Environ. Chem. Lett. 2022, 20, 153–188. [Google Scholar] [CrossRef]
  2. Agyekum, E.B.; Nutakor, C.; Agwa, A.M.; Kamel, S. A critical review of renewable hydrogen production methods: Factors affecting their scale-up and its role in future energy generation. Membranes 2022, 12, 173. [Google Scholar] [CrossRef]
  3. Matson, S.L.; Quinn, J.A. Membrane reactors. In Membrane Handbook; Ho, W.S., Sirkar, K.K., Eds.; Van Nostrand Reinhold: New York, NY, USA, 1992; pp. 809–832. [Google Scholar]
  4. Nagy, E. Basic Equations of Mass Transport through a Membrane Layer, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2019. [Google Scholar]
  5. Algieri, C.; Coppola, G.; Mukherjee, D.; Shammas, M.I.; Calabro, V.; Curcio, S.; Chakraborty, S. Catalytic membrane reactors: The industrial applications perspective. Catalysts 2021, 11, 691. [Google Scholar] [CrossRef]
  6. Drioli, E.; Brunetti, A.; Profio, G.D.; Barbieri, G. Process intensification strategies and membrane engineering. Green Chem. 2012, 14, 1561–1572. [Google Scholar] [CrossRef]
  7. Amiri, T.Y.; Ghasemzageh, K.; Iulianelli, I.A. Membrane reactors for sustainable hydrogen production through steam reforming of hydrocarbons: A review. Chem. Eng. Proces. Process Intensif. 2020, 157, 108148. [Google Scholar] [CrossRef]
  8. Alique, D.; Bruni, G.; Sanz, R.; Calles, J.A.; Tosti, S. Ultra-pure hydrogen via co-valorization of olive mill wastewater and bioethanol in Pd-membrane reactors. Processes 2020, 8, 219. [Google Scholar] [CrossRef] [Green Version]
  9. Di Marcoberardino, G.; Binotti, M.; Manzolini, G.; Viviente, J.L.; Arratibel, A.; Roses, L.; Gallucci, F. Achievements of European projects on membrane reactor for hydrogen production. J. Clean Prod. 2017, 161, 1442–1450. [Google Scholar] [CrossRef] [Green Version]
  10. Taghizadeh, M.; Aghili, F. Recent advances in membrane reactors for hydrogen production by steam reforming of ethanol as a renewable resource. Rev. Chem. Eng. 2019, 35, 377–392. [Google Scholar] [CrossRef]
  11. Iulianelli, A.; Palma, V.; Bagnato, G.; Ruocco, C.; Huang, Y.; Veziroğlu, N.T.; Basile, A. From bioethanol exploitation to high grade hydrogen generation: Steam reforming promoted by a Co-Pt catalyst in a Pd-based membrane reactor. Renew. Energy 2018, 119, 834–843. [Google Scholar] [CrossRef]
  12. Basile, A.; Pinacci, P.; Iulianelli, A.; Broglia, M.; Drago, F.; Liguori, S.; Longo, T.; Calabrò, V. Ethanol steam reforming reaction in a porous stainless steel supported palladium membrane reactor. Int. J. Hydrogen Energy 2011, 36, 2029–2037. [Google Scholar] [CrossRef]
  13. Iulianelli, A.; Liguori, S.; Vita, A.; Italiano, C.; Fabiano, C.; Huang, Y.; Basile, A. The oncoming energy vector: Hydrogen produced in Pd-composite membrane reactor via bioethanol reforming over Ni/CeO2 catalyst. Catal. Today 2016, 259, 368–375. [Google Scholar] [CrossRef]
  14. Iulianelli, A.; Basile, A. Hydrogen production from ethanol via inorganic membrane reactors technology: A review. Catal. Sci. Technol. 2011, 1, 366–379. [Google Scholar] [CrossRef]
  15. Chen, W.-H.; Li, S.-C.; Lim, S.; Chen, Z.-Y.; Juan, J.C. Reaction and hydrogen production phenomena of ethanol steam reforming in a catalytic membrane reactor. Energy 2021, 220, 119737. [Google Scholar] [CrossRef]
  16. Mattos, L.V.; Jacobs, G.; Davis, B.H.; Noronha, F.B. Production of hydrogen from ethanol: Review of reaction mechanism and catalyst deactivation. Chem Rev. 2012, 112, 4094–4123. [Google Scholar] [CrossRef]
  17. Andreev, D.V.; Makarshin, L.L.; Gribovskii, A.G.; Zaikovskii, V.I.; Parmon, V.N. Ethanol steam reforming over Co/ZnO and Rh/Al2O3 catalysts in a microchannel catalytic reactor. Kinet Catal. 2014, 55, 798–808. [Google Scholar] [CrossRef]
  18. Sun, S.; Yan, W.; Sun, P.; Chen, J. Thermodynamic analysis of ethanol reforming for hydrogen production. Energy 2012, 44, 911–924. [Google Scholar] [CrossRef]
  19. Saidi, M.; Jahangiri, A. Theoretical study of hydrogen production by ethanol steam reforming: Technical evaluation and performance analysis of catalytic membrane reactor. Int. J. Hydrogen Energy 2018, 43, 15306–15320. [Google Scholar] [CrossRef]
  20. Contreras, J.L.; Salmones, J.; Colín-Luna, J.A.; Nuco, L.; Quintana, B.; Cyrdova, I.; Zeifert, B.; Tapia, C.; Fuentes, G.A. Catalysts for H2 production using the ethanol steam reforming (a review). Int. J. Hydrogen Energy 2014, 39, 18835–18853. [Google Scholar] [CrossRef]
  21. Hou, T.; Zhang, S.; Chen, Y.; Wang, D.; Cai, W. Hydrogen production from ethanol reforming: Catalysts and reaction mechanism. Renew. Sustain. Energy Rev. 2015, 44, 132–148. [Google Scholar] [CrossRef]
  22. Soykal, I.I.; Sohn, H.; Ozkan, U.S. Effect of support particle size in steam reforming of ethanol over Co/CeO2 catalysts. ACS Catal. 2012, 2, 2335–2348. [Google Scholar] [CrossRef]
  23. Liu, Z.Y.; Duchon, T.; Wang, H.R.; Peterson, E.W.; Zhou, Y.H.; Luo, S.; Zhou, J.; Matolin, V.; Stacchiola, D.J.; Rodriguez, J.A.; et al. Mechanistic insights of ethanol steam reforming over Ni-CeOx(111): The importance of hydroxyl groups for suppressing coke formation. J. Phys. Chem. C 2015, 119, 18248–18256. [Google Scholar] [CrossRef]
  24. Mondal, T.; Pant, K.K.; Dalai, A.K. Mechanistic kinetic modeling of oxidative steam reforming of bioethanol for hydrogen production over Rh-Ni/CeO2-ZrO2 catalyst. Ind. Eng. Chem. Res. 2016, 55, 86–98. [Google Scholar] [CrossRef]
  25. Anil, S.; Indraja, S.; Singh, R.; Appari, S.; Roy, B. A review on ethanol steam reforming for hydrogen production over Ni/Al2O3 and Ni/CeO2 based catalyst powders. Int. J. Hydrogen Energy 2022, 47, 8177–8213. [Google Scholar] [CrossRef]
  26. Słowik, G.; Greluk, M. The Influence of active phase composition and reaction temperature on the catalytic properties of K-promoted Co–Ni/CeO2 catalysts in the steam reforming of ethanol. Catal. Lett. 2022, in press. [Google Scholar] [CrossRef]
  27. Ferreira, G.R.; Nogueira, F.G.; Lucrédio, A.F.; Assaf, E.M. Ethanol steam reforming by Ni catalysts for H2 production: Evaluation of Gd effect in CeO2 support. Catal. Lett. 2022, in press. [Google Scholar] [CrossRef]
  28. Mulewa, W.; Tahir, M.; Amin, N.A.S. Ethanol steam reforming for renewable hydrogen production over La-modified TiO2 catalyst. Chem. Eng. Trans. 2017, 56, 349–354. [Google Scholar]
  29. Rajabi, Z.; Martinelli, M.; Upton, G.F.; Watson, C.D.; Cronauer, D.C.; Kropf, A.J.; Jacobs, G. Low temperature ethanol steam reforming: Selectivity control with lithium doping of Pt/m-ZrO2. Catal. Today 2022, in press. [Google Scholar] [CrossRef]
  30. Dalena, F.; Giglio, E.; Marino, A.; Aloise, A.; Giorgianni, G.; Migliori, M.; Giordano, G. Steam Reforming of Bioethanol Using Metallic Catalysts on Zeolitic Supports: An Overview. Catalysts 2022, 12, 617. [Google Scholar] [CrossRef]
  31. Bepari, S.; Sarkar, J.J.; Pradhan, N.C. Kinetics of ethanol steam reforming over Ni/Olivine catalyst. Int. J. Hydrogen Energy 2022, in press. [Google Scholar] [CrossRef]
  32. Elharati, M.A.; Lee, K.M.; Hwang, S.; Hussain, A.M.; Miura, Y.; Dong, S.; Fukuyama, Y.; Dale, N.; Saunders, S.; Kim, T.; et al. The effect of silica oxide support on the catalytic activity of nickel-molybdenum bimetallic catalyst toward ethanol steam reforming for hydrogen production. Chem. Eng. J. 2022, 441, 135916. [Google Scholar] [CrossRef]
  33. Chen, B.; Xu, Y.; Xiong, C.; Rickard, S.; Boscoboinik, J.A.; Jiang, D.E.; Kidder, M.; Savara, A. Mechanism for acetone and crotonaldehyde production during steam reforming of ethanol over La0.7Sr0.3MnO3–x perovskite: Evidence for a shared C4 aldol addition intermediate. ACS Catal. 2022, 12, 4358–4374. [Google Scholar] [CrossRef]
  34. Li, X.; Li, J. Fluxes and driving forces in membrane separation processes. In Encyclopedia of Membranes; Drioli, E., Giorno, L., Eds.; Springer: Berlin/Heidelberg, Germany, 2015; pp. 1–3. [Google Scholar]
  35. Matsuura, T. Membrane separation technologies. In Wastewater Recycle, Reuse, and Reclamation; Vigneswaran, S., Ed.; EOLSS: Abu Dhabi, United Arab Emirates, 2009; Volume 1, pp. 98–135. [Google Scholar]
  36. Manohar. Development & characterization of ceramic membranes. Int. J. Mod. Eng. Res. 2021, 2, 1492–1506. [Google Scholar]
  37. Gozόlvez-Zafrilla, J.M.; Santafé-Moros, A.; Escolástico, S.; Serra, J.M. Fluid dynamic modeling of oxygen permeation through mixed ionic–electronic conducting membranes. J. Membr. Sci. 2011, 378, 290–300. [Google Scholar] [CrossRef]
  38. Hong, J.; Kirchen, P.; Ghoniem, A. Numerical simulation of ion transport membrane reactors: Oxygen permeation and transport and fuel conversion. J. Membr. Sci. 2012, 407–408, 71–85. [Google Scholar] [CrossRef]
  39. Khodaei, M.; Darabi, Z.; Omidifar, M. Modeling and simulation of the fluid dynamic and performance of the Pd-based membrane by CFD for hydrogen separation. Iran. J. Hydrog. Fuel Cell 2022, 9, 19–25. [Google Scholar]
  40. Sheintuch, M.; Nekhamkina, O. Architecture alternative for propane dehydrogenation in a membrane reactor. Chem. Eng. J. 2018, 347, 900–912. [Google Scholar] [CrossRef]
  41. Ghasemzadeh, K.; Aghaeinejad-Meybodi, A.; Iulianelli, A.; Basile, A. Theoretical performance evaluation of inorganic (non Pd-based) membranes for hydrogen separation. J. Membr. Sci. Res. 2018, 4, 198–203. [Google Scholar]
  42. Arratibel Plazaola, A.; Pacheco Tanaka, D.A.; van Sint Annaland, M.; Gallucci, F. Recent advances in Pd-based membranes for membrane reactors. Molecules 2017, 22, 51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Miachon, S.; Dalmon, J.-A. Catalysis in membrane reactors: What about the catalyst? Top. Catal. 2004, 29, 59–65. [Google Scholar] [CrossRef]
  44. Sadykov, V.A.; Eremeev, N.F.; Fedorova, Y.E.; Krasnov, A.V.; Bobrova, L.N.; Bespalko, Y.N.; Lukashevich, A.I.; Skriabin, P.I.; Smorygo, O.L.; Van Veen, A.C. Design and performance of asymmetric supported membranes for oxygen and hydrogen separation. Int. J. Hydrogen Energy 2021, 46, 20222–20239. [Google Scholar] [CrossRef]
  45. Bobrova, L.; Eremeev, N.; Vernikovskaya, N.; Sadykov, V.; Smorygo, O. Effect of asymmetric membrane structure on hydrogen transport resistance and performance of a catalytic membrane reactor for ethanol steam reforming. Membranes 2021, 11, 332. [Google Scholar] [CrossRef] [PubMed]
  46. Montebelli, A.; Visconti, C.G.; Groppi, G.; Tronconi, E.; Cristiani, C.; Ferreira, C.; Kohler, S. Methods for the catalytic activation of metallic structured substrates. Catal. Sci. Technol. 2014, 4, 2846–2870. [Google Scholar] [CrossRef]
  47. Haase, S.; Murzin, D.Y.; Salmi, T. Review on hydrodynamics and mass transfer in minichannel wall reactors with gas–liquid Taylor flow. Chem. Eng. Res. Des. 2016, 113, 304–329. [Google Scholar] [CrossRef]
  48. Eremeev, N.; Krasnov, A.; Bespalko, Y.; Bobrova, L.; Smorygo, O.; Sadykov, V. An experimental performance study of a catalytic membrane reactor for ethanol steam reforming over a metal honeycomb catalyst. Membranes 2021, 11, 790. [Google Scholar] [CrossRef]
  49. Boon, M. Scientific concepts in the engineering sciences epistemic tools for creating and intervening with phenomena. In Scientific Concepts and Investigative Practice; Feest, J., Steinle, F., Eds.; De Gruyter: Berlin, Germany, 2012; pp. 219–243. [Google Scholar]
  50. Bespalko, Y.; Eremeev, N.; Skryabin, P.; Krieger, T.; Chesalov, Y.; Lapina, O.; Khabibulin, D.; Ulihin, A.; Uvarov, N.; Sadykov, V. Structural and transport properties of neodymium tungstates prepared via mechanochemical activation. Ceram. Int. 2019, 45, 9529–9536. [Google Scholar] [CrossRef]
  51. Leonov, A.; Romashko, A. Laminar Ni-NixAl-Al2O3 foam material with high oxidation resistance. In Proceedings of the International Conference on Cellular metals and Metal Foaming Technology, Berlin, Germany, 23–25 June 2003; Banhart, J., Fleck, N.A., Mortensen, A., Eds.; MIT-Verlag: Berlin, Germany, 2003; pp. 271–274. [Google Scholar]
  52. Grelf.Net—Joyce Loebl Magiscan Image Analyser. Available online: https://www.grelf.net/magiscan.html (accessed on 20 July 2022).
  53. Reedy, C.L.; Anderson, J.; Reedy, T.J.; Liu, Y. Image analysis in quantitative particle studies of archaeological ceramic thin sections. Adv. Archaeol. Pract. 2014, 2, 252–268. [Google Scholar] [CrossRef] [Green Version]
  54. Wang, C.Y. Similarity stagnation point solutions of the Navier-Stokes equations—Review and extension. Eur. J. Mech. B Fluids 2008, 27, 678–683. [Google Scholar] [CrossRef]
  55. Homann, F. Der Einfluss grosser Zahigkeit bei der Stromung um den Zylinder und um die Kugel. Z. Angew. Math. Mech. 1936, 16, 153–164. [Google Scholar] [CrossRef]
  56. Matera, S.; Reuter, K. Transport limitations and bistability for in situ CO oxidation at RuO2 (110): First-principles based multiscale modeling. Phys. Rev. B Condens. Matter Mater. Phys. 2010, 82, 085446. [Google Scholar] [CrossRef] [Green Version]
  57. Matera, S.; Reuter, K. When atomic-scale resolution is not enough: Spatial effects on in situ model catalyst studies. J. Catal. 2012, 295, 261–268. [Google Scholar] [CrossRef] [Green Version]
  58. Zanier, F.; Michelon, N.; Canu, P. Design and characterization of a stagnation flow reactor for heterogeneous microkinetic studies. Chem. Eng. J. 2017, 315, 67–82. [Google Scholar] [CrossRef]
  59. Levenspiel, O. Chemical Reaction Engineering, 3rd ed.; John Wiley & Sons: New York, NY, USA, 1999. [Google Scholar]
  60. Incropera, F.P.; Bergman, T.L.; Lavine, A.S.; DeWitt, D.P. Fundamentals of Heat and Mass Transfer, 8th ed.; Wiley: Hoboken, NJ, USA, 2018; 992p. [Google Scholar]
  61. López, E.; Divins, N.J.; Anzola, A.; Schbib, S.; Borio, D.; Llorca, J. Ethanol steam reforming for hydrogen generation over structured catalysts. Int. J. Hydrogen Energy 2013, 38, 4418–4428. [Google Scholar] [CrossRef]
  62. Taylor, G. Dispersion of soluble matter in solvent flowing slowly through a tube. Proc. R. Soc. Lond. Math. Phys. Eng. Sci. 1953, 219, 186–203. [Google Scholar]
  63. Aris, R. On the dispersion of a solute in a fluid flowing through a tube. Proc. R.Soc. Lond. Math. Phys. Eng. Sci. 1956, 235, 67–77. [Google Scholar]
  64. Reid, R.C.; Sherwood, T.K.; Prausnitz, J. The Properties of Gases and Liquids, 3rd ed.; McGraw-Hill: New York, NY, USA, 1977. [Google Scholar]
  65. Fuller, E.N.; Schettler, P.D.; Giddings, J.C. A new method for prediction of binary gas-phase diffusion coefficients. Ind. Eng. Chem. 1966, 58, 19–27. [Google Scholar] [CrossRef]
  66. Hiemenz, K. Die Grenzschicht an einem in den gleichformingen Flussigkeits-strom einge-tauchten graden Kreiszylinder. Dinglers Polytech. J. 1911, 326, 321–324. [Google Scholar]
  67. Seader, J.D.; Henley, E.J.; Roper, K. Separation Process Principles: Chemical and Biochemical Operations, 3rd ed.; Wiley: Hoboken, NJ, USA, 2011. [Google Scholar]
  68. Siblukin, M. Heat transfer near the forward stagnation point of a body of revolution. J. Aeronaut. Sci. 1952, 19, 570–571. [Google Scholar]
  69. Truitt, R.W. Fundamentals of Aerodynamic Heating; Ronald Press: New York, NY, USA, 1960. [Google Scholar]
  70. Schlichting, H. Boundary Layer Theory, 4th ed.; McGraw-Hill Book Company: New York, NY, USA; London, UK, 1960. [Google Scholar]
  71. Nagy, E. Basic Equations of the Mass Transport Through a Membrane Layer, 1st ed.; Elsevier Insights, Elsevier: Amsterdam, The Netherlands; Boston, MA, USA, 2012. [Google Scholar]
  72. Deveau, N.D.; Ma, Y.H.; Datta, R. Beyond Sieverts’ law: A comprehensive microkinetic model of hydrogen permeation in dense metal membranes. J. Membr. Sci. 2013, 437, 298–311. [Google Scholar] [CrossRef]
  73. Suzuki, A.; Yukawa, H.; Nambu, T.; Matsumoto, Y.; Murata, Y. Consistent description of hydrogen permeability through metal membrane based on hydrogen chemical potential. Int. J. Hydrogen Energy 2014, 39, 7919–7924. [Google Scholar] [CrossRef]
  74. Suzuki, A.; Yukawa, H. A Review for consistent analysis of hydrogen permeability through dense metallic membranes. Membranes 2020, 10, 120. [Google Scholar] [CrossRef] [PubMed]
  75. Uvarov, N.F.; Ulichin, A.C.; Bespalko, Y.N.; Eremeeev, N.F.; Krasnov, A.V.; Skriabin, P.I.; Sadykov, V.A. Study of proton conductivity of composite metal-ceramic materials based on neodimium tugstates using a four-electrode technique with ionic probes. Int. J. Hydrogen Energy 2018, 43, 19521–19527. [Google Scholar] [CrossRef]
  76. Ambrosetti, M.; Bracconi, M.; Groppi, G.; Tronconi, E. Analytical geometrical model of open cell foams with detailed description of strut-node intersection. Chem. Ing. Tech. 2017, 89, 915–925. [Google Scholar] [CrossRef]
  77. Inayat, A.; Klumpp, M.; Lämmermann, M.; Freund, H.; Schwieger, W. Development of a new pressure drop correlation for open-cell foams. Chem. Eng. J. 2016, 287, 70–719. [Google Scholar] [CrossRef]
  78. Bracconi, M.; Ambrosetti, M.; Maestri, M.; Groppi, G.; Tronconi, E. A fundamental investigation of gas/solid mass transfer in open-cell foams using a combined experimental and CFD approach. Chem. Eng. J. 2018, 352, 558–571. [Google Scholar] [CrossRef]
  79. Bhattacharya, S.; Hwang, S.-T. Concentration polarization, separation factor, and Peclet number in membrane processes. J. Membr. Sci. 1997, 132, 73–90. [Google Scholar] [CrossRef]
  80. Nagy, E.; Vitai, M. Analysis of mass transport through anisotropic, catalytic/bio-catalytic membrane reactors. Catalysts 2019, 9, 358. [Google Scholar] [CrossRef] [Green Version]
  81. Henis, J.M.; Tripodi, M.K. Composite hollow fiber membranes for gas separation: The resistance model approach. J. Membr. Sci. 1981, 8, 233–246. [Google Scholar] [CrossRef]
  82. Huang, T.C.; Wie, M.C.; Chen, H.I. Permeation of hydrogen through palladium/alumina composite membrane. Sep. Sci. Technol. 2001, 36, 199–222. [Google Scholar] [CrossRef]
  83. Iwuchukwu, I.J.; Sheth, A. Mathematical modeling of high temperature and high-pressure dense membrane separation of hydrogen from gasification. Chem. Eng. Proc. 2008, 47, 1292–1304. [Google Scholar] [CrossRef]
  84. Cardoso, S.P.; Lin, Z.; Portugal, I.; Rodrigues, A.E.; Silva, C.M. Synthesis, dynamic characterization, and modeling studies of an AM-3 membrane for light gases separation. Microporous Mesoporous Mater. 2018, 261, 170–180. [Google Scholar] [CrossRef]
Figure 1. Sketch of a membrane reactor with catalytic monolith: 1—asymmetric supported hydrogen separation membrane, 2—catalytic monolith, 3—thermocouples pockets, 4—feeding tube, 5—Ar sweep gas tube, 6—permeate outlet tube, 7—retentate gas tube.
Figure 1. Sketch of a membrane reactor with catalytic monolith: 1—asymmetric supported hydrogen separation membrane, 2—catalytic monolith, 3—thermocouples pockets, 4—feeding tube, 5—Ar sweep gas tube, 6—permeate outlet tube, 7—retentate gas tube.
Membranes 12 00741 g001
Figure 2. Pictures of the experimental catalytic monolith: (a) a general view, (b) a fragment with dimensions.
Figure 2. Pictures of the experimental catalytic monolith: (a) a general view, (b) a fragment with dimensions.
Membranes 12 00741 g002
Figure 3. Cross-sectional images of the layer-by-layer assembled membrane module (a) and SEM microgragh of the inlet area (b).
Figure 3. Cross-sectional images of the layer-by-layer assembled membrane module (a) and SEM microgragh of the inlet area (b).
Membranes 12 00741 g003
Figure 4. Effect of the operating temperature (a) on species concentration (dry basis) and (b) on hydrogen permeation flux at the retentate and sweep streams. Data points correspond to the measurements of gas phase average composition.
Figure 4. Effect of the operating temperature (a) on species concentration (dry basis) and (b) on hydrogen permeation flux at the retentate and sweep streams. Data points correspond to the measurements of gas phase average composition.
Membranes 12 00741 g004
Figure 5. Temperature curves of feed-side mass transfer coefficients: a theoretically calculated one and those with application of correction factors.
Figure 5. Temperature curves of feed-side mass transfer coefficients: a theoretically calculated one and those with application of correction factors.
Membranes 12 00741 g005
Figure 6. Concentration (a) and hydrogen mass transfer resistance (b) profiles along the monolith-gap-membrane assembly in the MR at Pe = 0.76; T = 800 °C; Table 1.
Figure 6. Concentration (a) and hydrogen mass transfer resistance (b) profiles along the monolith-gap-membrane assembly in the MR at Pe = 0.76; T = 800 °C; Table 1.
Membranes 12 00741 g006
Figure 7. Comparison of the MRs with packed bed and monolith at different feed flow rates in contribution of each individual layers to the overall hydrogen mass transfer resistance.
Figure 7. Comparison of the MRs with packed bed and monolith at different feed flow rates in contribution of each individual layers to the overall hydrogen mass transfer resistance.
Membranes 12 00741 g007
Figure 8. Effect of the distance above the membrane surface (a) and the feed flow rate (b) on the permeation phenomena.
Figure 8. Effect of the distance above the membrane surface (a) and the feed flow rate (b) on the permeation phenomena.
Membranes 12 00741 g008
Figure 9. Effect of the temperature on the driving force (difference between concentrations on each side of a boundary layer adjacent to the membrane surface) and hydrogen permeation flux in the MRs with monolith (Feed: 5 Nl h−1) or the one with the packed bed (Feed: 3 Nl h−1).
Figure 9. Effect of the temperature on the driving force (difference between concentrations on each side of a boundary layer adjacent to the membrane surface) and hydrogen permeation flux in the MRs with monolith (Feed: 5 Nl h−1) or the one with the packed bed (Feed: 3 Nl h−1).
Membranes 12 00741 g009
Table 1. Structural parameters of the constituent layers of the experimental membrane reactor assisted with the catalytic monolith and experimental conditions (reference case) used in simulations.
Table 1. Structural parameters of the constituent layers of the experimental membrane reactor assisted with the catalytic monolith and experimental conditions (reference case) used in simulations.
VariableUnitsValue
Feed-side
Catalytic monolith:
Height (hmth)mm22
Diameter (dmth)mm24
Cross section area (Amth)mm2452.16
Equivalent channel diameter (dh,mth)mm0.6926
Porosity (εmth)(-)0.58
Volumetric surface area (SV,mth)m2 m−33355
Flow rate of ethanolNl h−10.6
Flow rate of steamNl h−12.4
Flow rate of argonNl h−12.0
Gap with stagnate flow:
height (hgap)mm4
Sweep-side
Membrane module:
Dense layer
Thickness (hdm)mm0.15
Powder layer
Thickness (hpl)mm0.4
Particle size (dpl)mm0.072
Hydraulic pore diameter (dpore,pl)mm0.012
Porosity (εpl)(-)0.2
Tortuosity(τpl)(-)4.2
Volumetric surface area (SV,pl)m2 m−366,667
Intermediate layer
Thickness (hil)mm0.6
Particle size (dil)mm0.061
Hydraulic pore diameter(dpore,il)mm0.027
Porosity (εil)(-)0.4
Tortuosity (τil)(-)3.4
Volumetric surface area (SV,il)m2 m−359,259
Foam layer
Thickness (hfoam)mm4.5
Cell diameter (dcell)mm2.2
Hydraulic pore diameter (dp,foam)mm1.006
Porosity (εfoam)(-)0.75
Tortuosity (τfoam)(-)1.42
Volumetric surface area (SV,foam)m2 m−31395.4
Flow rate of argonNl h−110
Table 2. Morphological and structural characteristics of the prepared asymmetric membrane module. Reprinted from Ref. [45] under the CC BY 4.0 License.
Table 2. Morphological and structural characteristics of the prepared asymmetric membrane module. Reprinted from Ref. [45] under the CC BY 4.0 License.
LayerCompositionThickness (µm)True Density (g cm−3)Particle Size b (µm)Pore Diameter b (µm)Porosity c
(%)
Dense layerNi–Cu/Nd5.5WO11.25-δ93.3–115 (center);
194–256 (edge)
6.60.045 for Ni–Cu,
0.1–1 for Nd5.5WO11.25-δ
15 (x)
42 (y)
~4
Powder layerNi-Al380–440~765 (x)
81 (y)
12 (x)
11 (y)
12–14
Intermediate layerNi-Al400–13005.3445 (x)
50 (y)
27 (x)
27 (y)
27–32
Foam layerAl2O3-SiO2 foam with Ni-Al coating4500–50004.632400(x) a
1800 (y) a
1000 (x)
1100 (y)
38–40
83 d
a Cell diameter; b x axis is parallel and y axis is perpendicular to the membrane surface; c Quantification with vector program and by image analysis; d Overall porosity including pores in Al2O3.
Table 3. Governing equations for the feed-side reactor.
Table 3. Governing equations for the feed-side reactor.
Monolith bed
Component molar balance
ε m t h ρ t o t x i , f t + 1 A m t h ( F f x i , f ) z D i , m t h ρ t o t 2 x i , f z 2 = S V , m t h R i , where R i = j = 1 4 ν i j r j
Conservation equation for the change in the total molar flow rate d F f d z = i F f , i d z = 2 A m t h S V , f ( r 1 + r 2 + r 4 ) .
Boundary conditions
Inlet (z = 0): x i , f | z = 0 = x i , f 0 F f = F f 0
Outlet (z = hmth):
i H 2
( F f x i , f A m t h D i , m t h ρ t o t x i , f z ) | z = h m t h = 0
i = H 2 ( F f x H 2 , f A m t h D H 2 , m t h ρ t o t x H 2 , f z ) | z = h m t h = β A m t h ρ t o t x H 2 , f + β A m ρ t o t x H 2 , d m
Gap with stagnate flow
Retentate product flow
F r
F r = F f | z = h m t h F H 2 , p e r m ,   where   F H 2 , p e r m = J d m , H 2 A m = β ρ t o t A m ( A m t h A m x H 2 , f | z = h m t h x H 2 , d m )
Mole fractions in the retentate gas flow
i H 2 x i , r = F f x i , f | z = h m t h F f | z = h m t h F H 2 , p e r m
i = H 2 x H 2 , r = F f x H 2 , f | z = h m t h F H 2 , p e r m F f | z = h m t h F H 2 , p e r m
Table 4. Summary of the kinetic expressions used in the simulation.
Table 4. Summary of the kinetic expressions used in the simulation.
ReactionsRate Equations
C 2 H 5 OH     CH 4   +   H 2   +   CO Δ r H 298   K o   =   49.0   kJ/mol r 1 = k f , 1 x Eth
CH 4   +   H 2 O     3 H 2   +   CO Δ r H 298   K o =   206.3   kJ/mol r 2 = k f , 2 x CH 4 α 2 x H 2 O β 2 [ 1 Q 2 K e q , 2 ]
CO   +   H 2 O     CO 2   +   H 2 Δ r H 298   K o =   41.2   kJ/mol r 3 = k f , I I I x CO α 3 x H 2 O β 3 [ 1 Q 3 K e q , 3 ]
CH 4   +   2 H 2 O     4 H 2   +   CO 2 Δ r H 298   K o =   165.1   kJ/mol r 4 = k f , I V x CH 4 α 4 x H 2 O β 4 [ 1 Q 4 K e q , 4 ]
Kinetic rate constant
k f , j = A j exp ( E j R T )
Equilibrium constant
K e q , j = exp ( Δ G j R T )
Reaction quotient for a reaction
c   C + d   D     a   A + b   B Q j = [ p A ] a [ p B ] b [ p C ] c [ p D ] d
Parameters of the rate equations for the reactions
A j ( mol m 2 s ) E j ( kJ mol ) α j β j K e q , j Unit
r 1 1.4 × 10451--
r 2 1.86 × 1057212 1.167 × 10 13 exp ( 26830 T ) ( atm 2 )
r 3 4.08 × 1045211 1.767 × 10 2 exp ( 4400 T ) ( atm 0 )
r 4 1.408 × 1048111.25 2.062 × 10 11 exp ( 22430 T ) ( atm 2 )
Table 5. Constitutive equations for the feed-side reactor model.
Table 5. Constitutive equations for the feed-side reactor model.
Effective axial dispersion of mass in the monolith
D i , m t h = D i m i x + u m t h 2 d h , m t h 2 192 D i m i x ,   where   D i m i x = 1 x i , f j = 1 , j i 7 x j , f D i j , D i j = 10 7 T 1.75 [ ( 1 M i + 1 M j ) ] 1 2 P [ ( υ ) i 1 3 + ( υ ) j 1 3 ] 2
Hydrogen mass transfer coefficient at feed-side of the membrane
β = S h D i m i x d h , m ,   where   S h = 0.763 R e 0.5 S c 0.4 ;   δ h , m = 2.4 ν [ m 2 s ] u [ 1 s ] = 2.4 μ f V g a p ρ t o t ρ f F f
R e = ρ f u d m 2 μ f and   S c = μ f ρ f D H 2 m i x ,
μ f = i = 1 n = 7 x i μ i i = 1 n = 7 x j Φ i j ,
μ i = μ i 0 T 0 + C i T + C i ( T T 0 ) 3 2 ,
Φ i j = [ 1 + ( μ i μ j ) 1 2 ( M j M i ) 1 4 ] 2 [ 8 ( 1 + M i M j ) ] 1 2 , Φ j i = Φ i j μ j μ i M i M j .
The reference viscosity at reference temperature and Sutherland’s temperature for gaseous substances
SubstanceC, Sutherland’s temperature (K)T0, reference temperature (K) μ i 0 reference viscosity (kg s−1 m−1)
H2O673873.163.09 × 10−5
CH4164873.162.46 × 10−5
CO2240873.163.61 × 10−5
CO102873.163.63 × 10−5
H272873.161.83 × 10−5
Ar142873.164.87 × 10−5
Table 6. Governing and constitutive equations for the sweep-side reactor.
Table 6. Governing and constitutive equations for the sweep-side reactor.
Dense layer of a thickness hdm in the membrane module
Hydrogen flow at steady state
F H 2 , p e r m = A m Q d m ( P H 2 , d m P H 2 , p l ) ,   where   P H 2 , d m = P x d m   and   P H 2 , p l = P x p l
Parity   flux   equation   at   the   feed - side   boundary   A m t h β ρ t o t x H 2 , f A m β ρ t o t x H 2 , d m = A m Q d m P ( x H 2 , d m x H 2 , p l ) .
Real solution of a quadratic from the parity equation for x d m
x H 2 , d m = Q d m P + P Q d m 2 + 4 β 2 ρ t o t 2 A m t h A m x H 2 , f | z = h m t h + 4 β ρ t o t Q d m P x H 2 , p l 2 β ρ t o t
Permeance   ( mol   m 2 s 1 atm 0.5 ) Q d m = Θ h d m exp ( E d m R T ) . Permeability   ( mol   m 1 s 1 atm 0.5 ) Θ = 1.4801 × 10 2 Activation energy (J mol−1)
60,000
Powder layer of the membrane module
ε p l ρ t o t x H 2 , p l t ρ t o t D H 2 Ar , p l e f f 2 x H 2 , p l z 2 = 0   and   x Ar , p l = 1 x H 2 , p l .
Boundary conditions:
at z = hmth + hgap + hdm:
ρ t o t D H 2 A r , p l e f f x H 2 , p l z | z = h m t h + h g a p + h d m = Q d m P ( x H 2 , d m x H 2 , p l ) ,
at z = hmth + hgap + hdm + hpl: ρ t o t D H 2 Ar , p l e f f x i , p l z = ρ t o t D H 2 Ar , i l e f f x i , i l z .
Intermediate layer of the membrane module
ε i l ρ t o t x H 2 , i l t = ρ t o t D H 2 Ar , i l e f f 2 x H 2 , i l z 2   and   x H 2 , i l = 1 x Ar , i l .
Boundary conditions:
at z = hmth + hgap + hdm + hpl:
ρ t o t D H 2 Ar , i l e f f x H 2 , i l z = ρ t o t D H 2 Ar , p l e f f x H 2 , p l z ,
at z = hmth + hgap + hdm + hpl + hil: ρ t o t D H 2 Ar , i l e f f x H 2 , i l z = ρ t o t D H 2 Ar , f o a m e f f x H 2 , f o a m z .
D H 2 Ar , p l ( i l ) e f f = D Ar H 2 , p l ( i l e f f = ε p l ( i l ) τ p l ( i l ) 1 2 ( 1 1 / D H 2 , p l ( i l ) k n + 1 / D H 2 Ar + 1 1 / D Ar , p l ( i l ) k n + 1 / D H 2 Ar ) .
D i , p l ( i l ) k n = d p o r e , p l ( i l ) 3 ( 8 R T π M i ) 1 / 2 = 48.5 d p o r e , p l ( i l ) T M i . d p o r e , p l ( i l ) = 4 ε p l ( i l ) S V , p l ( i l ) . d p o r e , p l ( i l ) = 4 ε p l ( i l ) S V , p l ( i l ) .
Foam layer in the membrane module
ε f o a m ρ t o t x H 2 , f o a m t = ρ t o t D H 2 Ar , f o a m e f f 2 x H 2 , f o a m z 2 ,   and   x Ar , f o a m = 1 x H 2 , f o a m ,
D H 2 Ar , f o a m e f f = ε f o a m ( D H 2 Ar τ f o a m + 0.5 d p , f o a m u f o a m ) ,   where pore diameter
d p , f o a m = d c e l l ( 3.7033 2.5516 ε f o a m + 0.7054 ε f o a m 2 ) ,
Tortuosity
τ f o a m = 1 + 4.867 [ 1 0.971 ( 1 ε f o a m ) 0.5 ] 4 ε f o a m ( 1 ε f o a m ) 0.5 ( 1 ε f o a m ) .
Boundary conditions:
at z = hmth + hgap + hdm + hpl + hil:
at z = hmth + hgap + hdm + hpl + hil + hfoam:
ρ t o t D H 2 Ar , f o a m e f f x H 2 , f o a m z = ρ t o t D H 2 Ar , f o a m e f f x H 2 , i l z , ρ t o t D H 2 Ar , f o a m e f f x H 2 , f o a m z = β f o a m ρ t o t ( x H 2 f o a m x H 2 , s w ) .
Hydrogen effective mass transfer coefficient at the sweep-side
S h d s , a v g = β f o a m d s , a v g D H 2 Ar , f o a m e f f ,   where   d s , a v g 2.85 ( 1 ε f o a m ) S V , f o a m ,   and   S h d s , a v g = ε f o a m 2 ( 0.566 R e d s , a v g 0.33 + 0.039 R e d s , a v g 0.8 ) S c f o a m 1 / 3
R e s , a v g = ρ f o a m u f o a m d s , a v g μ f o a m , S c f o a m = μ f o a m ρ f o a m D H 2 Ar , f o a m e f f , S V , f o a m = ( 2 3 π d c e l l ) ( 1 ε f o a m ) 0.5 .
Sweep compartment
Volumetric gas flow rate G s w = G s w 0 + β f o a m A m ( x H 2 , f o a m | z = h m t h + h g a p + h d m + h p l + h i l + h f o a m x H 2 , s w ) and x H 2 , s w + x Ar , sw = 1
x H 2 , s w = ( β f o a m + G s w 0 A m + β f o a m ) ( β f o a m + G s w 0 A m + β f o a m x H 2 , f o a m ) 2 4 β f o a m x H 2 , f o a m 2 β f o a m .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bobrova, L.; Vernikovskaya, N.; Eremeev, N.; Sadykov, V. Model-Based Performance Analysis of Membrane Reactor with Ethanol Steam Reforming over a Monolith. Membranes 2022, 12, 741. https://doi.org/10.3390/membranes12080741

AMA Style

Bobrova L, Vernikovskaya N, Eremeev N, Sadykov V. Model-Based Performance Analysis of Membrane Reactor with Ethanol Steam Reforming over a Monolith. Membranes. 2022; 12(8):741. https://doi.org/10.3390/membranes12080741

Chicago/Turabian Style

Bobrova, Ludmilla, Nadezhda Vernikovskaya, Nikita Eremeev, and Vladislav Sadykov. 2022. "Model-Based Performance Analysis of Membrane Reactor with Ethanol Steam Reforming over a Monolith" Membranes 12, no. 8: 741. https://doi.org/10.3390/membranes12080741

APA Style

Bobrova, L., Vernikovskaya, N., Eremeev, N., & Sadykov, V. (2022). Model-Based Performance Analysis of Membrane Reactor with Ethanol Steam Reforming over a Monolith. Membranes, 12(8), 741. https://doi.org/10.3390/membranes12080741

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