Mass Transfer Coefficient in Multi-Stage Reformer/Membrane Modules for Hydrogen Production

Hydrogen is a promising energy carrier, and is exploitable to extract energy from fossil fuels, biomasses, and intermittent renewable energy sources and its generation from fossil fuels, with CO2 separation at the source being one of the most promising pathways for fossil fuels’ utilization. This work focuses on a particular configuration called the Reformer and Membrane Module (RMM), which alternates between stages of Steam Reforming (SR) reactions with H2 separation stages to overcome the thermodynamic limit of the conventional SR. The configuration has numerous advantages with respect to the more widely studied and tested membrane reactors, and has been tested during a pilot-scale research project. Although numerous modelling works appeared in the literature, the design features of the material exchanger (in the so-called RMM architecture) of different geometrical configurations have not been developed, and the mass transfer correlations, capable of providing design tools useful for such membrane modules, are not available. The purpose of this work is therefore to apply a physical-mathematical model of the mass transfer, in three different geometries, considering both concentration polarization and membrane permeation, in order to: (i) simulate the cited experimental results; (ii) estimate the scaling-up correlations for the “material exchange modules”; and (iii) identify the mass transfer limiting regime in relation to the gas mass flow rate.


Introduction
Hydrogen is a promising energy carrier, and is exploitable to extract energy from fossil fuels, biomasses, and intermittent renewable energy sources [1,2]. Generating hydrogen from fossil fuels is one of the most promising alternatives in the framework of carbonaceous fuels' utilization, with simultaneous CO 2 sequestration [1,2]. The adoption of H 2 as an energy vector would represent a radical change in the energy sector, impacting its production, distribution, and consumption since it can be converted into both electrical power and heat using fuel cells or combustion engines without generating CO 2 emissions locally [3][4][5].
The current world hydrogen consumption is more than 50 Mton/year [1] and is consistently being devoted to the chemical and petrochemical sectors, such as for ammonia and methanol synthesis, The first experimental data collected at the Chieti Pilot plant during the testing phase, as well as the average H2 permeability at the operating conditions, can be found in our previous works [16][17][18][19]. The authors have also assessed and compared the benefits and drawbacks of the MMR configuration in relation to the MR [17][18][19][20]22]. Globally, the Research Project proved that the RMM presents a production cost which is 10% lower than those of the conventional H2 scheme, allows for implementation of a direct CO2 sequestration unit [22], and shows the following main advantages with respect to the MR arrangement: • The MR is mechanically complex and presents a large and unpractical heat transfer surface-in the MR, the concentric tube geometry yields an imbalance between the surfaces required for heat transfer (outer tube) and the available surface for mass transfer (of the inner membrane tube); • RMM enables the de-coupling of separation and reaction operating temperatures, increasing the stability and the durability of the membranes and enabling independent optimization of the reforming temperature; • RMM simplifies the mechanical design of membrane tubes compared with the one embedded in a catalyst tube, and a simple shell and tube geometry can be selected for the tubular separation module; • RMM simplifies maintenance of the Pd/Ag membrane modules and catalyst replacement.
With these premises, we analyzed multiple data related to the pilot plant experiments, and through the mathematical modeling dedicated to the membrane modules, we also estimated the mass transfer coefficients. By selecting several data at different flow-rates, gas compositions, temperatures, and different geometrical configurations, we calculated semi-empirical correlations useful for scaling up the RMM. The aforementioned mathematical model has been carefully validated through experimental data in order to extrapolate the values of the transport coefficients and the correlations necessary to scale-up the "material exchangers", which represents the core of the RMM architecture.
From a permeation modelling point of view, a first benchmark in the theoretical modelling of H2 permeation trough a Pd self-supported membrane has been developed by Ward et al. [23]. The model was utilized to individuate the rate-limiting processes among the fundamental kinetic steps: (i) external mass transfer (binary mixtures); (ii) surface adsorption and desorption; (iii) transitions to and from the bulk metal; and (iv) diffusion within the metal. On this basis, Caravella et al. modelled The first experimental data collected at the Chieti Pilot plant during the testing phase, as well as the average H 2 permeability at the operating conditions, can be found in our previous works [16][17][18][19]. The authors have also assessed and compared the benefits and drawbacks of the MMR configuration in relation to the MR [17][18][19][20]22]. Globally, the Research Project proved that the RMM presents a production cost which is 10% lower than those of the conventional H 2 scheme, allows for implementation of a direct CO 2 sequestration unit [22], and shows the following main advantages with respect to the MR arrangement:

•
The MR is mechanically complex and presents a large and unpractical heat transfer surface-in the MR, the concentric tube geometry yields an imbalance between the surfaces required for heat transfer (outer tube) and the available surface for mass transfer (of the inner membrane tube); • RMM enables the de-coupling of separation and reaction operating temperatures, increasing the stability and the durability of the membranes and enabling independent optimization of the reforming temperature; • RMM simplifies the mechanical design of membrane tubes compared with the one embedded in a catalyst tube, and a simple shell and tube geometry can be selected for the tubular separation module; • RMM simplifies maintenance of the Pd/Ag membrane modules and catalyst replacement.
With these premises, we analyzed multiple data related to the pilot plant experiments, and through the mathematical modeling dedicated to the membrane modules, we also estimated the mass transfer coefficients. By selecting several data at different flow-rates, gas compositions, temperatures, and different geometrical configurations, we calculated semi-empirical correlations useful for scaling up the RMM. The aforementioned mathematical model has been carefully validated through experimental data in order to extrapolate the values of the transport coefficients and the correlations necessary to scale-up the "material exchangers", which represents the core of the RMM architecture.
From a permeation modelling point of view, a first benchmark in the theoretical modelling of H 2 permeation trough a Pd self-supported membrane has been developed by Ward et al. [23]. The model was utilized to individuate the rate-limiting processes among the fundamental kinetic steps: (i) external mass transfer (binary mixtures); (ii) surface adsorption and desorption; (iii) transitions to and from the bulk metal; and (iv) diffusion within the metal. On this basis, Caravella et al. modelled the transport in Pd-alloy membranes, considering both the non-ideal hydrogen transportation in membrane (influences of support, inhibition by CO by the Sievert-Langmuir equation, and the effect of membrane polarization) and external mass transfer resistance on hydrogen permeation [24][25][26]. Their results are reported as an "adjusted" pressure exponent of a Sieverts-type empirical law, and it is function of several factors, such as temperature, total pressure, membrane thickness, and non-ideal behaviours. The theoretical meaning of the adjusted-Sieverts exponent can be explained by dividing the overall hydrogen permeation into several elementary steps (adsorption, desorption, diffusion in the Pd-based layer, and the two transitions, surface-to-bulk and bulk-to surface) [26,27]. Sarti et al. tested a Pd 80 -Ag 20 (NGK) in a shell and tube configuration with different mixtures and operating conditions. The experimental tests were performed with a mixture of H 2 /N 2 /CH 4 . At the same time, a theoretical model based on the previous model of Ward et al. was developed to analyze the competitive adsorption of hydrogen and carbon monoxide molecules. The experimental tests show the existence of a concentration polarization phenomena due to non-permeable species [28][29][30][31]. In the first experimental setup [28], the Sherwood number followed a boundary-layer type of correlation, whereas in the second [30], a linear correlation between the Sherwood and Péclet numbers was found. Globally, the concentration polarization has been extensively discussed in the literature, and has been included in several models and at different conditions, mainly by modifying the exponent for the pressure dependence in the Sievert law equation [26][27][28][29][30][31][32].
Other authors discussed the transport-reaction-permeation regimes, also addressing the competitive adsorption limiting the overall H 2 permeation [14,15]. Barbieri et al. [33] interpreted the observed decrease in hydrogen flux through a palladium-silver membrane over time with a CO inhibition (by up to 2 bars) in terms of a Sieverts-Langmuir model, assuming a linear correlation between the decrease in hydrogen permeance and surface coverage by carbon monoxide. Consequently, they accounted for the membrane surface fraction not available for hydrogen permeation using a Langmuir affinity constant for carbon monoxide and a temperature-dependent "permeance reduction factor" [33].
To the best of our knowledge, although numerous modelling works have appeared in the literature, the design features of the material exchanger (in the so-called RMM architecture) of different geometrical configurations have not been developed, and the mass transfer correlations, capable of providing design tools useful for such modules, are also not available.
The purpose of this work is therefore to apply a physical-mathematical mass-transfer model at three different geometries, and by considering both concentration polarization and membrane permeation, to simulate the cited experimental results, thereby estimating the scaling-up correlations for the "material exchange modules".

Materials and Methods
The Pilot Plant, realized by Tecnimont-KT Kinetics Technologies in Chieti (Italy), includes two reaction zones (R1, R2) and two separation zones, as depicted in the block diagram of Figure 2. The reaction zone is 15 m high, and consists of a radiant area (with burners and catalytic tubes) and a convective zone, where superheated steam is produced. The separation zone consists of two membrane separators-the first presenting two membranes (M-01 and M-02) working in parallel, and the second including a single module (M-03). The natural gas (NG) was provided by the NG network at 12 barg and desulfurized in the HDS unit. Steam was added to the feed and the mixture was preheated in the convective zone. The stream enters the first reformer where the reactions take place at 550-680 • C [16][17][18][19]. The effluent (syngas) is cooled in an air cooler to 450 • C, and routed to the first stage of separation. The depleted syngas (30-35% of H 2 is removed) passes to the next reaction and separation zone. Permeates of both membranes are mixed and sent to the flare or to the cooling system in case of sweeping (water vapor condensation with cooling water in the closed cycle). Finally, the retentate stream is routed to the flare. The tests have been realized for 4 months to investigate the long-term stability of the structured catalyst and membrane modules, to carry out a comparison between structured catalyst and traditional Ni-based pellets, and to test the effectiveness of the RMM architecture integrated with the CPO reactor and to prove the improvement obtained by the RMM enhancement. Further details on the pilot plant can be found in the cited literature [16][17][18][19]. of the RMM architecture integrated with the CPO reactor and to prove the improvement obtained by the RMM enhancement. Further details on the pilot plant can be found in the cited literature [16][17][18][19]. More than 70 tests, each running for around 10 h, were performed by varying the main operating conditions, such as temperature and pressure, as well as the steam-to-carbon ratio in the reforming section and the flowrate in the membrane modules. At the beginning of the test runs, the heating was realized by steam and nitrogen to reach a temperature above 300 °C (heating rate under 2.5-3 °C/min). Due to the endothermicity of the reaction, a slight temperature drop is observed after the feed introduction. This paper focused on the results related to the permeation zone, highlighted in Figure 2 and presented in Section 4. The main independent variable considered for the scaling-up correlations was the syngas flowrate (out of the reforming) that generates different mass transfer coefficients.
The main geometrical characteristics of the separation units are summarized in Table 1. The first membrane module, M-01, contained 13 tubular membranes of palladium supported on alumina ( Figure 2a). The M-02 module ( Figure 3b) consisted of five plate membranes, with each side consisting of a dense Pd-Ag layer deposited on the external surface of an α-alumina support (each side formed by two Pd-Ag panels welded in series). Both membranes can be used with sweeping gas. The last membrane module, M-03, used in the present work, was similar to M-01 but included three tubular membranes ( Figure 3c).
The permeability value that accounts for the diffusion and solubility of hydrogen in palladium and palladium-silver alloys, was calculated for the different materials based on the theoretical and experimental works of Holleck [34] and Sarti et al. [28] according to the following general expression that diffusion is an energetically activated process and that Sievert's constant represents the equilibrium reaction constant for hydrogen dissociation: where P is the permeability pre-exponential factor, and Ea is the activation energy for hydrogen permeability which contains contributions from the activation energy for the diffusion of hydrogen atoms, the standard enthalpy of the surface dissociation reaction, as well as the entropy change of the dissociation reaction. These estimated values are reported in Table 1. More than 70 tests, each running for around 10 h, were performed by varying the main operating conditions, such as temperature and pressure, as well as the steam-to-carbon ratio in the reforming section and the flowrate in the membrane modules. At the beginning of the test runs, the heating was realized by steam and nitrogen to reach a temperature above 300 • C (heating rate under 2.5-3 • C/min). Due to the endothermicity of the reaction, a slight temperature drop is observed after the feed introduction. This paper focused on the results related to the permeation zone, highlighted in Figure 2 and presented in Section 4. The main independent variable considered for the scaling-up correlations was the syngas flowrate (out of the reforming) that generates different mass transfer coefficients.
The main geometrical characteristics of the separation units are summarized in Table 1. The first membrane module, M-01, contained 13 tubular membranes of palladium supported on alumina (Figure 2a). The M-02 module (Figure 3b) consisted of five plate membranes, with each side consisting of a dense Pd-Ag layer deposited on the external surface of an α-alumina support (each side formed by two Pd-Ag panels welded in series). Both membranes can be used with sweeping gas. The last membrane module, M-03, used in the present work, was similar to M-01 but included three tubular membranes (Figure 3c).
The permeability value that accounts for the diffusion and solubility of hydrogen in palladium and palladium-silver alloys, was calculated for the different materials based on the theoretical and experimental works of Holleck [34] and Sarti et al. [28] according to the following general expression that diffusion is an energetically activated process and that Sievert's constant represents the equilibrium reaction constant for hydrogen dissociation: where P 0 H 2 is the permeability pre-exponential factor, and E a is the activation energy for hydrogen permeability which contains contributions from the activation energy for the diffusion of hydrogen atoms, the standard enthalpy of the surface dissociation reaction, as well as the entropy change of the dissociation reaction. These estimated values are reported in Table 1.

Mathematical Modelling
The present section describes the differential equations used to model the permeation phenomena and the procedure to estimate the mass transfer coefficient from the experimental data.
The main variables of the mathematical model are represented in Figure 4, which shows the longitudinal flux of the syngas (in the external tube) flowing in counter-current to the permeated flux of hydrogen and sweep gas (in the inner tube). The schematic representation is based on the experiments related to the shell and tube configuration (M-01 and M-03) with the H 2 flux occurring from the outside to the inside of the tubes, where OD t is the external diameter of the inner tube and ID s is the diameter of the shell. The same configuration can be used for plate membranes (M-02) where the permeation channel is realized between two interspaces.

Mathematical Modelling
The present section describes the differential equations used to model the permeation phenomena and the procedure to estimate the mass transfer coefficient from the experimental data.
The main variables of the mathematical model are represented in Figure 4, which shows the longitudinal flux of the syngas (in the external tube) flowing in counter-current to the permeated flux of hydrogen and sweep gas (in the inner tube). The schematic representation is based on the experiments related to the shell and tube configuration (M-01 and M-03) with the H2 flux occurring from the outside to the inside of the tubes, where ODt is the external diameter of the inner tube and IDs is the diameter of the shell. The same configuration can be used for plate membranes (M-02) where the permeation channel is realized between two interspaces. The membrane is simulated as an isotherm and isobar material exchanger (considering both permeance and molar fraction in the retentate side as constants) enabling the selective hydrogen permeation. Referring to the membrane's infinitesimal volume of a length dz, the variation of the mass flow Fj can be written as in Equations (2) and (3), respectively, for the permeate and the retentate side (with the minus sign if the flow is counter-current).
The Hydrogen flux JH2 through the membrane was calculated by adopting the schematization of the film theory represented in Figure 5. The mass transfer resistances considered in the diffusive process were concentrated from the bulk to the membrane wall on both the side of the membrane (retentate/permeate side), as well as through the palladium membrane (internal diffusion according to the Sievert-Fick law). In the absence of sweeping gas, the resistance in the permeate side is null and = . The fluxes , , , respectively related to the driving force in the three abovementioned zones, are expressed by the Equations (4)-(6). For the experimental tests performed in absence of vapour sweeping, the resistance in the permeate side is neglected and the analysis focuses on the mass transport coefficient on the retentate side, explicated by Equation (7), that takes into account the concentration polarization phenomena with the mass transfer coefficient * , depending on the concentration at the interface [35]. The membrane is simulated as an isotherm and isobar material exchanger (considering both permeance and molar fraction in the retentate side as constants) enabling the selective hydrogen permeation. Referring to the membrane's infinitesimal volume of a length dz, the variation of the mass flow F j can be written as in Equations (2) and (3), respectively, for the permeate and the retentate side (with the minus sign if the flow is counter-current).
The Hydrogen flux J H2 through the membrane was calculated by adopting the schematization of the film theory represented in Figure 5. The mass transfer resistances considered in the diffusive process were concentrated from the bulk to the membrane wall on both the side of the membrane (retentate/permeate side), as well as through the palladium membrane (internal diffusion according to the Sievert-Fick law). In the absence of sweeping gas, the resistance in the permeate side is null and p RI H2 = p P H2 . The fluxes J R H2 , J M H2 , J P H2 , respectively related to the driving force in the three abovementioned zones, are expressed by the Equations (4)-(6). For the experimental tests performed in absence of vapour sweeping, the resistance in the permeate side is neglected and the analysis focuses on the mass transport coefficient on the retentate side, explicated by Equation (7), that takes into account the concentration polarization phenomena with the mass transfer coefficient F R * G , depending on the concentration at the interface [35].
A mathematical algorithm was implemented to estimate the mass transfer coefficients in the retentate side by regarding the observed experimental compositions, with partial pressure, as the unknown variable alongside the membrane. On the other hand, the membrane permeance was fixed on the basis of Equation (1) (calculated in agreement with the literature results [28,29,34]) by dividing the permeability by the membrane thickness δ, and thus does not represent a dependent variable of the mathematical model. The calculation procedure starts by assuming an initial value of FG. At each step, the algorithm calculates the molar fraction of hydrogen on the retentate side by solving the nonlinear Equation (7) (by means of the Levenberg-Marquardt algorithm).
The Equations (8) and (9) show the linearization of the driving force of Sievert's Law, here implemented to calculate the overall mass transfer coefficient FOG.
A mathematical algorithm was implemented to estimate the mass transfer coefficients in the retentate side F R G by regarding the observed experimental compositions, with partial pressure, as the unknown variable alongside the membrane. On the other hand, the membrane permeance was fixed on the basis of Equation (1) (calculated in agreement with the literature results [28,29,34]) by dividing the permeability P H 2 by the membrane thickness δ, and thus does not represent a dependent variable of the mathematical model. The calculation procedure starts by assuming an initial value of F G . At each step, the algorithm calculates the molar fraction of hydrogen on the retentate side y LI H 2 by solving the nonlinear Equation (7) (by means of the Levenberg-Marquardt algorithm).
The Equations (9) and (10) show the linearization of the driving force of Sievert's Law, here implemented to calculate the overall mass transfer coefficient F OG . The numerical integration was performed in Matlab by dividing the domain into 100 elements. The algorithm stopped when the difference between the flow of hydrogen on the permeate side (at the inlet of the membrane) was simulated by the mathematical model and the experimental value was minimized. In particular, an error of 1% between the calculated and experimental values was accepted.
Eventually, the estimated mass transfer coefficient F R G could be correlated with the geometrical and operating conditions, keeping in mind that the experiments were realized with different membranes, at different temperature and flow rates by means of the non-dimensional numbers Re, Sc, and Sh, also depending on the physical properties of the mixture on the retentate side. The analysis of the competitive adsorption, well-characterized in the abovementioned papers [13][14][15], is kept outside of the model due to the low concentration of CO. The operation and geometrical parameters assumed to calculate these numbers are reported in Table 2, for the tubular (M-01 and M-02) and the plate (M-03) modules, respectively.

Results and Discussion
As described in Section 2, 70 test runs were selected for the purpose of this work. The experimental data obtained from a singular test (as a clarifying example) are reported in Figure 6. The inlet and outlet flowrates and temperatures were available as experimental data for the coefficient estimation.
The figure shows the gas concentrations both at the outlet of the first reactor (entrance to the membrane) and at the outlet of the membrane. The example is related to the first installed membrane module (M-01). The operating conditions are also reported, referring to the process diagram of Figure 2. The table also reports the characterization of the natural gas supplied to the overall RMM process.
At this point, to perform the mass transfer coefficient estimation, the experimental data were grouped in different levels of flowrate to find a possible correlation between F G and the Reynolds number relative to the three tested membranes: five values for M-01, three values for M-02, and three values for M-03. These average experimental data (e.g., composition, hydrogen recovery, flow-rate, and permeability) are summarized in Tables 3-5.    Tables 3-5. At this point, to perform the mass transfer coefficient estimation, the experimental data were grouped in different levels of flowrate to find a possible correlation between FG and the Reynolds number relative to the three tested membranes: five values for M-01, three values for M-02, and three values for M-03. These average experimental data (e.g., composition, hydrogen recovery, flow-rate, and permeability) are summarized in Tables 3-5.  Tables 3-5.  By implementing the described procedure with the independent variables fixed at the values in  Tables 3-5, the mass transfer coefficient F G was estimated. Figure 7 reports the F G estimation results versus the Reynolds number characterizing each of the experiments. Each point represents the average value at each Reynolds number with error bars confined in the range of ±15%. The values globally follow a good linear correlation with the Reynolds number in the semi-log diagram, qualitatively confirming the observations made by Catalano et al. [28].
By adopting the approach of non-dimensional analysis, the mass transfer dimensionless group j D (Equation (11)) [35] can be found to correlate with the Reynolds number in Figure 8.
This representation produced a correlation equation (Equation (12)) which was valid in the range of tested conditions for the different geometries.
j D = 2.172 × 10 −10 ∆Re 2.79 (12) Furthermore, to test again the reliability of the procedure, the obtained correlation (Equation (12)) was tested by reproducing the experimental results by implementing our model with a feedforward approach where, in this case, both the permeability and the F G were known parameters and the latter was calculated though Equation (12). The modelling results, presented again in comparison with the observed ones, are in the form of calculated compositions at the exit of the membrane modules. These results are compared against the experimental values in the Parity diagram of Figure 9, finding an acceptable agreement. This last validation result confirms the goodness of the correlation obtained for the gas-side coefficient.     On this basis, referring to the film theory schematized in Figure 5, it is possible to calculate the relative weight of the mass transfer resistance R of the elementary steps (in-series) reported in Equations (12) and (13), respectively for the gas-phase in the retentate zone R cp (including the concentration polarization) and the one related to the permeation through the membrane R M . On this basis, referring to the film theory schematized in Figure 5, it is possible to calculate the relative weight of the mass transfer resistance R of the elementary steps (in-series) reported in Equations (13) and (14), respectively for the gas-phase in the retentate zone R cp (including the concentration polarization) and the one related to the permeation through the membrane R M . Figure 10 shows the calculated resistance for the experimental conditions versus the Reynolds number. The resistance relative to the permeation step through the membrane is confined in a band (red) between the calculated values for the tubular (M-01,02) and plane (M-03) modules, having different thicknesses and different Pd-Ag ratios. The concentration polarization resistance R cp is an order of magnitude higher than the resistance R M due to the permeation in the membrane layer (according to the mechanisms postulated by Sieverts). It is worth noting that these data were obtained for membranes of different materials and different configurations. Therefore, in this work we enabled the calculation of the mass transfer coefficients in these peculiar "material exchangers" highlighting the two regimes of Figure 10, practically extending the work of Catalano et al. [27,28], being realized at very low numbers of Reynolds.

Conclusions
The extensive use of hydrogen as a carrier would be a solution to the current conflict between economic expansion and pollution. It is also the only way to decarbonize the conversion processes of fossil fuels, waste, and biomass. In this paper, we gave new strength to the research on multistage In the left side of the graph, the gas-side transport appears to be the limiting step; in the second, for Reynolds numbers greater than 25,000, the two resistances are comparable. As the number of Reynolds increases, the controlling step becomes the transport across the membrane, characterized by Sievert's law.

Conclusions
The extensive use of hydrogen as a carrier would be a solution to the current conflict between economic expansion and pollution. It is also the only way to decarbonize the conversion processes of fossil fuels, waste, and biomass. In this paper, we gave new strength to the research on multistage reactors with intermediate hydrogen separation (so-called "RMM architecture"), focusing in particular on the "material exchanger" design, basically on the scale-up of these arrangements of metal membrane devices. The correlations now available in the literature focused on the estimation of the permeation direct from the Sievert Law (with or without changes of the exponent relating to partial pressures) resulting in a partial or inaccurate approximation of the phenomenology. Our approach, although still affected by simplifications, allows for a more accurate estimation of the transport coefficients in the membrane material exchangers for H 2 separation, tested for different membrane layers and composition, as well as different operating conditions straddling two zones characterized by two different rate-limiting steps.
In this work, we reported the experimental tests obtained at the pilot plant, designed during the R&D Project, entitled "Pure hydrogen from natural gas reforming up to total conversion obtained by integrating chemical reaction and membrane separation", and constructed by Tecnimont-KT Kinetics Technologies in Chieti (Italy). Many test runs have been collected and organized in relation to the flow rates and the membrane type.
A good correlation between our model and the experimental results validates the estimation of the mass transfer coefficient. Furthermore, it was possible to find reliable scaling-up correlations, including the whole set of data from different membrane configuration and operation conditions. The proposed correlation also allowed us to show that in the operating conditions, the mass transfer resistance due to concentration polarization can limit the hydrogen flux.
Author Contributions: Conceptualization and methodology, D.B. and M.C.; resources, data curation and validation, G.F., all authors contributed to the remaining activities regarding the paper production.
Funding: Part of this work was carried out within the framework of the project "Pure hydrogen from natural gas reforming up to total conversion obtained by integrating chemical reaction and membrane separation", financially supported by MIUR (FISR DM 17/12/2002)-Italy.

Conflicts of Interest:
The authors declare no conflict of interest.