Next Article in Journal
Characterization and Assessment of a Novel Plate and Frame MD Module for Single Pass Wastewater Concentration–FEED Gap Air Gap Membrane Distillation
Next Article in Special Issue
Hydrogen Production via Steam Reforming: A Critical Analysis of MR and RMM Technologies
Previous Article in Journal
Highly Selective Oxygen/Nitrogen Separation Membrane Engineered Using a Porphyrin-Based Oxygen Carrier
Previous Article in Special Issue
Theoretical Evaluation of Graphene Membrane Performance for Hydrogen Separation Using Molecular Dynamic Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Techno-Economic Assessment in a Fluidized Bed Membrane Reactor for Small-Scale H2 Production: Effect of Membrane Support Thickness

by
Gioele Di Marcoberardino
1,
Jasper Knijff
2,
Marco Binotti
3,
Fausto Gallucci
2 and
Giampaolo Manzolini
3,*
1
Department of Mechanical and Industrial Engineering, University of Brescia, 25123 Brescia, Italy
2
Inorganic Membranes and Membrane Reactors, Department of Chemical Engineering and Chemistry, Eindhoven University of Technology, 5612 AZ Eindhoven, The Netherlands
3
Politecnico di Milano, Department of Energy, 20156 Milano, Italy
*
Author to whom correspondence should be addressed.
Membranes 2019, 9(9), 116; https://doi.org/10.3390/membranes9090116
Submission received: 6 August 2019 / Revised: 3 September 2019 / Accepted: 3 September 2019 / Published: 6 September 2019

Abstract

:
This paper investigates the influence of the support material and its thickness on the hydrogen flux in Palladium membranes in the presence of sweep gas in fluidized bed membrane reactors. The analysis is performed considering both ceramic and metallic supports with different properties. In general, ceramic supports are cheaper but suffer sealing problems, while metallic ones are more expensive but with much less sealing problems. Firstly, a preliminary analysis is performed to assess the impact of the support in the permeation flux, which shows that the membrane permeance can be halved when the H2 diffusion through the support is considered. The most relevant parameter which affects the permeation is the porosity over tortuosity ratio of the porous support. Afterward, the different supports are compared from an economic point of view when applied to a membrane reactor designed for 100 kg/day of hydrogen, using biogas as feedstock. The stainless steel supports have lower impact on the hydrogen permeation so the required membrane surface area is 2.6 m2 compared to 3.6 m2 of the best ceramic support. This ends up as 5.6 €/kg H2@20bar and 6.6 €/kg H2@700bar for the best stainless steel support, which is 3% lower than the price calculated for the best ceramic support.

1. Introduction

The development of innovative technologies based on renewable energy is becoming essential for a more sustainable future. Among alternative energy carriers, hydrogen is certainly the most popular because of the high energy density on mass base and the absence of CO2 emissions at the point of use. The global hydrogen production is around 600–720 Billion Nm3/year [1], increasing at a rate of 5–6% per year. About 96% of the hydrogen is produced from fossil fuels [2], mostly natural gas (NG) as feedstock accounting for 50% of the world hydrogen production, and only 3–4% is produced without using fossil fuels [3]. The hydrogen produced by fossil fuels cannot be considered CO2 neutral, as conventional hydrogen production plants emit large CO2 emissions into the atmosphere (about 380–420 gCO2/Nm3H2) [4]. The production of H2 from renewables is usually performed via water electrolysis using the electricity generated by photovoltaic or wind power stations. The resulting price of H2 production is around 5 $/kg [2], above the target of 2.4 $/kg (original reference reports the cost in 2 €/kg [5]), and about three times the H2 cost from a large-scale steam reforming plant equal to 1.28 $/kg [6]. An alternative option consists of producing the H2 from biogas: The conversion process can be similar to the steam reforming one (biogas composition is mainly CH4 diluted in CO2), but this has to be applied at a much smaller scale. The typical size of existing BG plants corresponds to around 0.20–25 Nm3/s of biogas (with a methane content of 40–60%mol), which is 100 times smaller than the large-scale reforming plant, but there is a significantly higher number of plants (by the end of 2017, in Europe there were over 17,700 biogas plants for an overall installed capacity of 10.5 MWel) [7]. Considering the small-scale of biogas plants, it is mandatory to adopt an efficient system for biogas conversion to hydrogen which is slightly affected by the size of the components. At these scales, the hydrogen price is much higher with conventional systems, as heat recovery is more complicated, and export of steam/electricity is not possible.
Membrane reactors have been always foreseen as a promising technology because of their capability to simultaneously produce and separate the hydrogen [8,9,10,11]. Hydrogen production is performed via steam-reforming reaction and hydrogen separation through H2 selective membranes [12]. The problem of integrating membranes in a steam reformer reactor relies in the heat management: The reaction is very endothermic, so heat must be supplied to the reaction side while avoiding damaging the membrane [13]. Fluidized membrane reactors with auto-thermal reforming were demonstrated to solve this problem, as well as reducing the mass transfer limitation of the hydrogen from the bulk of the reaction to the membrane wall [14,15,16].
Among the different types of hydrogen perm-selective membranes, the Pd-based type has always been identified as the most suitable because of the operating temperature range (up to 500–600 °C) similar to the steam-reforming reaction [11,17]. Moreover, they offer the best compromise between permeance, selectivity, and chemical stability [16]. As a drawback, one of the most relevant Pd membrane limitations, which is beyond technological development, is the high cost of palladium [18]. For this reason, Pd membranes can be deposited on a support to reduce the thickness with cost benefits and to increase the mechanical stability with respect to the self-supported one. More resistant membranes can work with larger total pressure differences between the feed and permeate sides, with benefits for the hydrogen permeating across the membrane according to permeation law [19].
Hence, the membrane layer support must have specific properties to allow the deposition of the Pd layer on top of this [20]. One of the most important properties is the pore diameter: On one hand, it shall be significantly lower than the Pd layer, which is between 2 to 5 µm [20]; on the other hand, the smaller the support pores, the higher the resistance to hydrogen permeation and; therefore, the impact of the support to the overall flux [16,21]. An option that partly overcomes this issue is the adoption of asymmetric supports: The substrate where the Pd layer is deposited consists of small pores, while the majority of the support consists of larger pore sizes [18].
Concerning the reactor design, as anticipated, the feed gas pressure has a direct effect on the permeation driving force through the membrane, but also the pressure on the permeate side is relevant. The typical configuration adopted to increase the H2 flux consists of decreasing the hydrogen partial pressure in the permeate side by using a vacuum or a sweep gas [22]. The former option leads to simpler plant configurations but lower performance. Using sweep gas would reduce the overall cost as no vacuum pump had to be used on the permeate side [23].
However, some studies suggest that using sweep gas increases the diffusion resistance of the support, thus reducing the driving force across the membrane [24,25,26,27]. Furthermore Zeng et al. and Li et al. discovered that the H2 transport out of the ceramic support was slowed down by the counter diffusion of N2 [28,29]. According to Pinacci and Drago [22], an increase in sweep gas flow rate decreased the H2 flux instead of increasing it, due to the resistance in mass transfer through the support [22]. The main limitation of this work was the absence of the polarization effects in the retentate and the permeate side, as pure hydrogen was used as feedstock to implement the dusty-gas model (DGM). Nordio et al. [18] performed a study on the effect of sweep gas on the supported Pd membranes with concentration polarization and diluted hydrogen. The authors compared three different membranes with different support thickness and with different sweep gas flow. The experimental campaign revealed that the support thickness has a direct effect on the pressure drop, thus on the mass transfer limitation. The mass transfer limitation is inversely proportional to the hydrogen binary diffusivity in the mixture, so higher partial pressure leads to high mass transfer resistance. In addition, a high pressure drop over the porous support gives a negative impact on the interface pressure of the palladium surface increasing the pressure at the interface with penalties on the flux [18]. The experimental results confirmed the importance of accounting for the diffusion barrier of the porous support, so mathematical models have been developed [26,27] or innovative types of membranes as duplex (the Pd layer is on both sides of the tube) have been manufactured [24].
To sum up, the sweep gas benefits can be counterbalanced by the support geometry (thickness) and material properties. Therefore, the identification of the optimal sweep flow requires detailed information on the type of support and its properties and specific modelling studies. In general, ceramic (alumina, α-Al2O3) and metallic (porous stainless steel (PSS)) materials are used as membrane support. Whilst the first can be easily made to obtain small pores and/or different pore size distribution, but has low mechanical strength, the latter do not suffer the leakage problem thanks to the easier connection to the reactor shell, but shows low hydrogen selectivity due to large pores. Further comparison between the two are given below in Table 1.
Typically, thicker ceramic supports (up to 3.5 mm) are the most common support used in Pd membranes [21,30,31], being cheaper than the metallic ones and ensuring high selectivity.
This works aims at giving more insight on the effect of the support type integrated in a fluidized membrane reactor using sweep gas on the amount of hydrogen permeated at different operating pressures and sweep gas flow rates. A detailed comparison in terms of hydrogen flux and costs between ceramic and porous stainless steel (PSS) supports is performed. The analysis starts with the performance assessment of the reactor over a wide range of operating conditions (feed and permeate pressures, sweep flows) assuming methane as feedstock in the membrane reactor. This analysis points out the support properties that mostly affect the hydrogen permeation. Then, the impact of the support on the hydrogen permeation are applied to biogas as feedstock, so to evaluate the more convenient type of membrane from both thermodynamic and economic point of views. The analysis is carried out using Aspen and integrates the details of the hydrogen permeation phenomena with the overall plant performance. To the authors’ knowledge, this is the first time that this kind of comparison is performed, in particular identifying the most convenient type of support from an economic point of view. This research is performed within the BIONICO European project, where a fluidized membrane reactor capable of producing up to 100 kg/d of H2 is manufactured and tested in a real biogas plant site.
The paper is organized as follows: Section 2 describes the mathematical model of the membrane reactor, Section 3 reports the membrane reactor performance as a function of the type of support and its properties. The case study is described in Section 4, presenting the calculated cost of hydrogen for the different considered supports. Finally, the conclusions of the work and future works are outlined in the last section.

2. Materials and Methods

The membrane reactor model is developed in Aspen Custom Modeler (ACM) which allows the development of proprietary components to be integrated in Aspen Plus and Aspen Plus Dynamics for steady-state and dynamic simulations of complex chemical/energy systems. The ACM adopted in this work was developed and validated in previous studies performed by Nordio et al. and Foresti et al. [16,18]. The model is based on the dusty-gas model [26,36]. Figure 1 shows the scheme of the fluidized membrane reactor outlining the H2 partial pressure trend in presence (left side) or not (right side) of the sweep gas. This model works in a fluidized bed regime, and the bed is divided into three regions. The first region is used for the reaction of steam methane reforming (SMR), water gas shift (WGS), and oxidation (OX). The oxidation zone is considered because an auto-thermal reactor is adopted to simplify the heat management and reactor design [37]. The second region, also known as the membrane region, is used to separate the produced hydrogen. The third region is used as overhead for the expansion of the bed. The three reactions that are happening in the bottom region (SMR, WGS, and OX) are modelled by the kinetic laws from Trimm and Lam, and Numaguchi and Kikuchi [38,39]. The reactor is designed as isothermal, which corresponds to preliminary tests running on the prototype of this reactor [31]. The performance is represented by the hydrogen recovery factor (HRF). As shown in Equation (1), this is defined as the moles of hydrogen permeated over the theoretical recoverable number of hydrogen moles of the complete reforming of the feed (i.e., 4 per each mole of methane), subtracting the methane combusted to support the reaction [16].
HRF = n ˙ H 2 , perm 4 · n CH 4 , ret + n CO , ret + n H 2 , ret + n H 2 , perm = n ˙ H 2 , perm 4 · n ˙ CH 4 , in 2 n ˙ O 2 , in   .
The model can include also the concentration polarization (CP) losses which, as reported in a previous work [16,18], are described by a parameter to be tuned according to the operating conditions (pressure, temperature, compositions). Since for the considered case, this information was not available, the CP losses were neglected; therefore, the partial pressure of H2 in the emulsion phase is the same as on the Pd membrane surface, so to only focus on the effect of the membrane support on the flux.
As mentioned in the introduction, the sweep gas plays an important role. The diluent (e.g., N2 or steam) is often used as the H2 carrier. In this model, steam is used as it is easier to be produced on site and separated. Concerning the support type, the main parameters, which determines the hydrogen permeation, are the porosity, tortuosity, and pore diameter. In Table 2 and Table 3 the characteristics of different types of supports available in literature for both PSS and ceramic are reported. It can be noted that values cover a wide range; therefore, the analysis was limited to selected specific cases.

Modelling of Pd Membranes:

The permeation is modelled (see Figure 2) assuming a two-step process: Firstly the hydrogen permeates across the Pd membrane, and then the support resistance is considered assuming the sweep gas on the permeate side. The flux through the Pd membrane is described as:
J s   ( mol / m 2 s ) = Pe t ( P Feed , H 2 n P Perm , H 2 n )   .
Based on experimental data, the H2 flux is interpolated and according to Equation (2) the n-value is determined. The membrane permeability is temperature dependent. For this reason, it is described as an Arrhenius temperature dependency Equation (3).
Pe = Pe o exp ( E A RT )   .
For the activation energy, the tests reported in [18] are used for the PSS and ceramic supports. The H2 diffusion across the Pd membrane is described by Equations (2) and (3), while the DGM [45] (see Equation (4)) is used to account for the diffusion of the H2 across the porous support.
j = 1 j i n x i N j x j N i pD i , j eff N i pD i , K eff = 1 RT x i r + x i pRT ( B 0 p μ D i , K eff + 1 ) p r ,
where B0 is the viscous permeability of the support, Di,K is the Knudsen diffusivity for the i species, and N is the flux of either hydrogen or steam. In this model, the Knudsen diffusion flow (first term), viscous flow (second term), and the binary diffusion (third term) are accounted for [25].
The effective Maxwell–Stefan gas–gas diffusivity [18,46] is given by:
D i , j eff = ε τ D i , j   .
The effective Knudsen diffusivity is given by:
D i , K eff = ε τ d pore 3   8 RT π M i   .
The viscous permeability is a function of the pore diameter, ε and τ. This formula [18,46] is used if convection through the support is modelled as laminar and for cylindrical pores:
B 0 = ε τ d pore 2 32   .

3. Results

The impact of membrane support on hydrogen permeation is firstly evaluated using natural gas as feedstock. The selection of natural gas is due to the greater availability of experimental measurements and works available in literature to be compared against [16,18,23]. The membrane reactor characteristics and operating conditions are summarized in Table 4. The feed pressure is varied from 8 to 20 bars. Table 4 summarizes the feed streams and reactor conditions for the sweep gas case. Four different cases (see with several support thickness (from 0.5 to 3 mm)) are proposed to fully describe the effect of the pore size and ε/τ and cover a wide range of membrane support types.
The characteristics of the membrane supports investigated are reported Table 5. Figure 2 shows the Pd composite membrane where the three points are, respectively, the emulsion phase, interface, and permeate partial pressure of H2. Since, the length of the membrane is kept equal to 0.4 m, the H2 flux will decrease along the length.
Figure 3 shows the different pressure gradients for the cases considered with support thicknesses of 0.5 and 3 mm. Even though the initial conditions are the same (limited differences can be noted due to the solid circulation, which is modified by the model depending on the permeated hydrogen), significant differences in the pressure profiles can be noted: Assuming a fixed length, lower hydrogen permeation ends up in higher partial pressure at the feed side. Other considerations are:
  • Case 2 is the one that has the highest penalties in terms of partial pressures with thicker membranes, while Case 4 has very limited differences; Case 2 have roughly the same pore radius of Case 4, while 10 times lower ε/τ;
  • Case 1 and Case 3 have pretty similar with limited impact of the support thickness on the partial pressure trend. Case 1 and Case 3 share similar ε/τ with larger pore diameters for the alumina support (Case 1);
The considerations on the pressures can be translated into the HRF parameter, as shown in Figure 4.
  • Case 4 shows the highest HRF thanks to the lower resistance introduced by the support, hence a lower driving force requested. Case 1 and Case 3 have almost the same HRF and trend;
  • The support thickness can affect the HRF by more than 20%, ranging from 40% of Case 2 to 60% of Case 4;
  • The permeate pressure is an important factor to justify the higher HRF, with 3 mm supports compared to the one for 2.5 mm supports. As the calculations are performed at constant outer diameter and sweep inlet pressure, the thicker the membrane, the lower the final permeate pressure as a consequence of the higher pressure drops in the permeate side. In the 20 bar case, the permeate pressure can go even below the atmospheric one (this is not a design condition for a real plant, but it serves for the purpose of this analysis).
For sake of completeness, the partial pressure difference for the different cases as a function of the support thickness is reported in the Appendix A.
A further analysis is carried out to elucidate the different contributions of the pore diameter and ε/τ on the hydrogen permeation (see Figure 5) assuming a feed pressure of 15 bar and pore diameter of 0.36 µm. The HRF is affected by any ε/τ value for thicker supports (3 mm), while only for the very low values (below 0.4) for thicknesses of 0.5 mm. Moving to the pore size impact, a similar trend to the ε/τ parameter can be pointed out. The main difference is in the magnitude of the impact: The ε/τ can change the HRF from 0.35 to 0.55, while the pore size has half the impact as the HRF ranges from 0.32 to 0.46. To sum up, a higher ratio of ε/τ will diminish the effect of the support, since less collisions with the wall occurs during H2 permeation and a higher Knudsen diffusion occurs. For small pores, the thickness of the support becomes more important. Increasing the pore size, the limitations are minimized. Hence, a turning point occurs in the HRF because less collisions to the wall happens.
As anticipated, the effect of varied sweep flow also plays an important role in maximum-achieved HRF according to the two studies [18,22,27].The calculations performed in this work at varying sweep flow rates are reported in Figure 6, using the two thickness’ based on mechanical stability and sealing issues equal to 0.5 and 2 mm [31] for 316 L and α-Al2O3, respectively. On the x-axis, a sweep flow to H2 permeated ratio is used because it allows implementation of the results in different models. A clear asymptotic trend can be seen confirming the limited advantages of using large sweep flow rates. This graph also shows the HRF when diffusion effects are neglected, whereby the stainless steel Case 4 (i.e., high ε/τ) comes close to the ideal case.
This preliminary analysis confirmed the importance of accounting for the diffusion limitation in the support, as it reduces the hydrogen permeated. In general, the metallic and the ceramic supports show similar performance when they have similar thicknesses and ε/τ. It seems that metallic membranes can have higher ε/τ, with benefits for the hydrogen permeation. Certainly, the thickness of the ceramic support shall be higher than the metallic ones to limit the leakages in correspondence to the sealings.
Finally, a new parameter called equivalent permeance is introduced to quantify the impact of the support in the permeation process. The equivalent permeance represents the value of the Pd membrane permeance without any support, which will determine the same amount of hydrogen permeated of the actual Pd membrane with the support. It is determined changing the permeance value in the model so as to achieve the same results above reported.
The values of the equivalent permeance determined for the different cases as a function of the support thickness are reported in Figure 7. The effect of the support can even half the initial permeance value of the Pd membrane in the case of thick support and high ε/τ.
In the next section, a case study will be introduced so as to evaluate both the thermodynamic performance and the economic costs of the different supports considered.

4. Case Study

The different supports are integrated in a membrane reactor for hydrogen production in the BIONICO reactor. The BIONICO reactor is designed to produce 100 kg/day of pure hydrogen, starting from a landfill type of biogas, with the following molar composition: CH4 equal to 44.2%, CO2 equal to 34.0%, N2 equal to 16.0%, and O2 equal to 2.7%. The plant layout (Figure 8), which was introduced in a previous work [21], is here briefly discussed. Table 6 summarizes the thermodynamic properties of the main streams involved for one of the simulation case with landfill biogas [21]. The biogas is compressed, then preheated up to 300 °C in HX-3. It is mixed with steam and water before final heating in HX-4 prior to the inlet of the ATR-MR. At the ATR-MR outlet, the retentate (point 3) and hydrogen (point 6) leave from the top section of the reactor. The retentate stream, mainly consisting of H2O, CO2, and N2 with the remaining fraction of CO, CH4, and H2 is firstly cooled in a series of heat exchangers down to 200 °C (point 4) and then sent to a catalytic burner for a complete combustion of the remaining fuels with an additional air stream. The thermal power released during retentate, exhaust gases, and permeate cooling covers the heat duty for the steam production (steam for the sweep gas and for the reforming reaction). In particular, the sweep gas steam is evaporated, cooling the exhaust gas and superheated using the permeate stream. The geometry adopted for the heat transfer corresponds to counter-current configuration with hot gas flow externally to staggered tube bundles and cold streams in tube. Such configuration provides compact design while keeping a high heat transfer coefficient. The resulting exchange area is then used for the economic analysis. The separators are simulated as simple adiabatic flash units. The main parameters values adopted for the system auxiliary components (compressors, fans, pumps, heat exchangers, etc.) are taken from previous works [20].
The plant layout is implemented in Aspen Plus® [47], where mass and energy balances are solved. The methodology adopted is consistent with previous works [16,21,48]. The Peng–Robinson cubic equation of state [49] is used for all thermodynamic properties, except for liquid molar volume evaluation where the Rackett model [50] is used, and for steam properties where NBS/NRC steam tables [51] are adopted. To save computational time, the achievements of the previous analysis (in particular the equivalent permeance) are included in the model implemented. The number of membranes and, consequently, the reactor design are determined so as to achieve the HRF that closes the thermal balance of the system avoiding any additional biogas input.
The membrane reactor characteristics are reported in Table 7, pointing out the significant area variation (i.e., number of tubes) between the different cases requested for producing the same amount of hydrogen. As a term of reference, the area calculated in a previous work for the BIONICO project with simplified calculation of the support diffusion resistance is also added.

5. Economic Analysis and Discussion

Starting from the results determined in the previous section, an economic analysis is carried out to identify the most cost effective configuration between the different reactor designs as a function of the support type. This analysis is important as stainless steel supports are much more expensive than ceramic ones [52]. This analysis has an accuracy of ±50% [53] because of the development stage of the components. However, it will provide a good indication on the most suitable configuration to be adopted.
The parameter of merit adopted for comparing the different cases is the cost of hydrogen production (LCOH) [2], defined as:
LCOH = ( TPC CCF ) + C O & M fix + ( C O & M var h eq ) kg H 2   ,
where TOC is the total overnight cost, CCF is the capital charge factor, C O & M fix is the fix operating and maintenance (O&M) costs, while C O & m var is the O&M variable costs. The fixed cost includes insurance, maintenance, and labor wages, while the variable cost accounts for water, fuel, and the cost of electricity consumed. The equivalent operating hours is equal to 7500 h and the CCF is equal to 0.13 (corresponding to a plant lifetime of 15 years was considered with discount rate of 10%), consistent with [48,54].
Table 8 summarizes the methodology to calculate the TOC. The equipment cost is determined as in Equation (9).
C i , 2018 = ( C i , 0 ( S i S i , 0 ) f ) y × CEPCI 2018 CEPCI y   .
The component cost actualized to 2018 is scaled and actualized by this formula to account for the different sizes (Si) with respect to the reference size (Si,0) and cost (Ci,0) by a scale factor f. Moreover, to actualize the component cost the chemical engineering plant cost index (CEPCI) of the year 2018 is used, where this value is equal to 603.1 [55]. The corresponding CEPCI of the reference equipment cost should be applied.
The bare equipment costs for the different plant components is summarized in, they are taken from [48] (Table 9).
Table 10 reports the O&M fixed and variable costs (i.e., catalyst, biogas, water, auxiliaries, maintenance, insurance, and operators costs). For the cost of ceramic and stainless steel, data analysis carried out in Ferret project is used [52]. This value includes the cost of the support, membrane manufacturing, and the cost of the selective membrane layer. The used catalyst is the BIONICO catalyst, whose quantity is defined by the ACM model (50% catalyst and 50% filler particles).
Table 11 summarizes the CAPEX and OPEX for the different cases considered. As expected, Case 2 has the highest cost of the membrane reactor as a consequence of the larger membrane area/tubes (+60%), though the costs increase is quite limited (+15%). Case 4 has the lowest reactor cost because of the smaller membrane surface area; however, Case 1 has the cheapest membrane because of the lower price of the ceramic support over stainless steel. From an overall point of view and considering the cost assessment accuracy, Case 1, Case 2, and Case 4 share the same TOC.
The fixed OPEX is mainly determined as a fixed percentage of the TOC; therefore, there is a proportionality between the two costs and no significant difference between Cases 1, 3, and 4 can be pointed out (as it was for CAPEX).
The variable cost (the cost of process water, biogas, and electricity) remains the same for all cases because they all have the same HRF and efficiency. However, the membrane replacement cost and the catalyst is very significant with impacts on O&Mvar. The cost reduction is mainly in the amount of catalyst needed, proportional to the size of the reactor, while limited differences can be pointed out for the membrane replacement costs (as it was in the CAPEX, because the cheaper support balances the higher membrane area). Variable O&M indicates a limited preference for the stainless steel support with respect to the ceramic one
Finally, the LCOH is reported in Figure 9 considering a total amount of H2 produced per year equal to 31,250 kg/y. Case 4 has the lowest LCOH, equal to 5.3 €/kg H2@20bar and 6.3 €/kg H2@700bar. Cases 1 and 3 have the same LCOH in the range of 5.4 €/kg H2@20bar and 6.4 €/kg H2@700bar, so with differences below 3% with respect to the optimal case. Case 4 has a LCOH higher by 10% with respect to Case 1. The accuracy of the economic assessment does not allow the drawing of any final conclusions, but certainly the higher price of stainless steel tubes can be balanced by the higher permeability leading to the competitive cost of LCOH.
Considering the limited accuracy of the economic analysis a sensitivity analysis on the main cost assumption was performed and the main outcomes reported in Figure 10. Results show that the overall capital cost and electricity price have the highest impact on LCOH leading to a variation around 10%. The other parameters are less relevant and, in particular, reducing the cost of the membrane by 50% results in only a 5% LCOH cost reduction.
The calculate LCOH is significantly higher than the target H2 price, which is 2 €/kg [5]. The higher price is due to the smaller-scale production (which does not allow export of steam/electricity), high cost of biogas, and the need for further development in the membrane technology focusing on higher permeability, lower manufacturing costs, and increased lifetime. However, the sensitivity analysis pointed out that the considered layout is penalized because of the small-scale and the utilization of biogas as feedstock, with respect to the adoption of natural gas.

6. Conclusions

This paper performed a detailed techno-economic assessment on support thickness for hydrogen production in a membrane reactor, using Palladium membranes in fluidized bed conditions in the presence of sweep gas. The advantage of stainless steel is the easy integration, but the disadvantage if the high cost. Ceramic supports are cheaper but need to be thicker since integration is more difficult at higher temperatures. The effect of the support used has been investigated in detail to identify the most relevant characteristics. It was found that the porosity over tortuosity ratio is more relevant than the pore diameter, as it significantly affects gas diffusion. To assess the impact of the support on the permeance, an equivalent permeance is introduced. For the best case scenario, stainless steel has a reduction of 5% with respect to a pure Pd membrane. For the ceramic cases, this parameter shows a reduction ranging between 30% and 60% compared to the reference condition.
These results are applied to a fluidized membrane reactor operating with biogas for hydrogen production equal to 100 kg/day. This reactor is currently under testing in the BIONICO project (summer 2019). The analysis showed that, for the same biogas and hydrogen production, the best stainless steel support requires 2.6 m2 compared to 3.6 m2 of the best ceramic support, meaning that the reactor is also larger. The resulting cost of hydrogen production resulted equal to 5.6 €/kg H2@20bar and 6.6 €/kg H2@700bar for the best stainless steel support, which is 3% lower than the price calculated for the best ceramic support. The main difference is related to the variable O&M costs as a consequence of the larger membrane reactor. The accuracy of the economic assessment does not allow for the drawing of any final conclusions, but certainly the higher price of stainless steel tubes can be balanced by the higher permeability leading to the competitive cost of LCOH.

Author Contributions

Conceptualization, G.M. and G.D.M.; methodology, G.M., M.B. and G.D.M.; software, J.K. and G.D.M.; validation, J.K. and G.D.M.; writing—original draft preparation, J.K. and G.M.; writing—review and editing, G.D.M., M.B. and F.G.; supervision, G.M. and F.G.

Funding

This research was funded by BIONICO project grant agreement No. 671459.

Acknowledgments

The BIONICO project leading to this application has received funding from the Fuel Cells and Hydrogen 2 Joint Undertaking under grant agreement No 671459. This joint undertaking receives support from the European Union’s Horizon 2020 research and innovation programme, Hydrogen Europe, and Hydrogen Europe Research.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

AmemMembrane area, m2
BViscous permeability
DKnudsen diffusion
dDiameter
EaEnergy activation, kJ/mol
k0Pre exponential factor, mol/smPan
MMolar mass
nExponential factor, -
NMolar flux
PePermeability, mol/sm2Pan
pPressure, bar
p i Partial pressure, bar
rradius
RUniversal gas constant
tSupport thickness [m]
TTemperature, °C
Acronyms
ACMAspen Custom Modeler
ATRAuto-thermal reformer
ATR-MRAuto-thermal membrane reformer
BGBiogas
BSRBiogas steam reforming
CAPEXCapital expenditure
C&OCOwner’s and contingencies costs
CCFCapital charge factor
CPConcentration polarization
DGMDusty-gas model
EUEurope
HRFHydrogen recovery factor, -
HXHeat exchanger
ICIndirect costs, €
ID/ODInner/outer diameter, mm
LCOHLevelized cost of hydrogen, €/kg
LFLandfill
LHVLow heating value
LTLow temperature
NBS/NRCNational Bureau of Standards/Nuclear Regulatory Commission
NGNatural gas
O&MOperation and maintenance costs, €
OPEXOperating expenditure
OXOxidation
PSSPorous stainless steel
S/CSteam to carbon molar ratio, -
SMRSteam methane reforming
TECTotal equipment cost, €
TICInstallation costs, €
TPCTotal plant cost, €
WGSWater gas shift
Subscripts
Ffeed
permATR-MR permeate side
retATR-MR retentate side
memMembrane
Greek letters
εPorosity
τtortuosity
ηel,refAverage electric efficiency of the power generating park, %

Appendix A

Additional results for the different types of supports considered.
Membranes 09 00116 g0a1

References

  1. Ball, M.; Wietschel, M. The Hydrogen Economy: Opportunities and Challenges; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  2. IRENA. Hydrogen From Renewable Power: Technology Outlook for the Energy Transition; IRENA: Abu Dhabi, United Arab Emirates, 2018. [Google Scholar]
  3. Sumbramani, V.; Basile, A.; Verizoglu, N.T. Compendium of Hydrogen Energy: Hydrgoen Production and Purification; Woodhead Publishing: Sawston, UK, 2015. [Google Scholar]
  4. Dincer, I.; Acar, C. Review and evaluation of hydrogen production methods for better sustainability. Int. J. Hydrogen Energy 2015, 40, 11094–11111. [Google Scholar] [CrossRef]
  5. Hydrogen Price Target. Available online: https://www.iea.org/hydrogen2019/ (accessed on 20 July 2019).
  6. Parkinson, B.; Balcombe, P.; Speirs, J.F.; Hawkes, A.D.; Hellgardt, K. Levelized cost of CO2 mitigation from hydrogen production routes. Energy Environ. Sci. 2019, 12, 19–40. [Google Scholar] [CrossRef]
  7. European Biogas Association. EBA Annual Report 2018; European Biogas Association: Brussels, Belgium, 2019. [Google Scholar]
  8. Uemiya, S.; Sato, N.; Ando, H.; Matsuda, T. Steam reforming of methane in a hydrogen-permeable membrane reactor. Appl. Catal. A Gen. 1991, 67, 223–230. [Google Scholar] [CrossRef]
  9. Adris, A.M.; Grace, J.R. Characteristics of Fluidized-Bed Membrane Reactors: Scale-up and Practical Issues. Ind. Eng. Chem. Res. 1997, 36, 4549–4556. [Google Scholar] [CrossRef]
  10. Lu, G.Q.; Diniz da Costa, J.C.; Duke, M.; Giessler, S.; Socolow, R.; Williams, R.H.; Kreutz, T. Inorganic membranes for hydrogen production and purification: A critical review and perspective. J. Colloid Interface Sci. 2007, 314, 589–603. [Google Scholar] [CrossRef] [PubMed]
  11. Gallucci, F.; Fernandez, E.; Corengia, P.; van Sint Annaland, M. Recent advances on membranes and membrane reactors for hydrogen production. Chem. Eng. Sci. 2013, 92, 40–66. [Google Scholar] [CrossRef]
  12. Tosti, S.; Basile, A.; Bettinali, L.; Borgognoni, F.; Gallucci, F.; Rizzello, C. Design and process study of Pd membrane reactors. Int. J. Hydrogen Energy 2008, 33, 5098–5105. [Google Scholar] [CrossRef]
  13. Roses, L.; Gallucci, F.; Manzolini, G.; Campanari, S.; van Sint Annaland, M. Comparison between fixed bed and fluidized bed membrane reactor configurations for PEM based micro-cogeneration systems. Chem. Eng. J. 2011, 171, 1415–1427. [Google Scholar] [CrossRef]
  14. Gallucci, F.; van Sint Annaland, M.; Kuipers, J.A.M. Autothermal reforming of methane with integrated CO2 capture in novel fluidized bed membrane reactors. Asia-Pac. J. Chem. Eng. 2009, 4, 334–344. [Google Scholar] [CrossRef]
  15. Deshmukh, S.A.R.K.; Volkers, S.; van Sint Annaland, M.; Kuipers, J.A.M. Heat Transfer in a Membrane Assisted Bubbling Fluidized Bed with Immersed Horizontal tubes. Int. J. Chem. React. Eng. 2005. [Google Scholar] [CrossRef]
  16. Foresti, S.; Di Marcoberardino, G.; Manzolini, G.; De Nooijer, N.; Gallucci, F.; van Sint Annaland, M. A comprehensive model of a fluidized bed membrane reactor for small-scale hydrogen production. Chem. Eng. Process. Process Intensif. 2018, 127, 136–144. [Google Scholar] [CrossRef]
  17. Fernandez, E.; Coenen, K.; Helmi, A.; Melendez, J.; Zuñiga, J.; Pacheco Tanaka, D.A.; Van Sint Annaland, M.; Gallucci, F. Preparation and characterization of thin-film Pd-Ag supported membranes for high-temperature applications. Int. J. Hydrogen Energy 2015, 40, 13463–13478. [Google Scholar] [CrossRef]
  18. Nordio, M.; Soresi, S.; Manzolini, G.; Melendez, J.; Van Sint Annaland, M.; Pacheco Tanaka, D.A.; Gallucci, F. Effect of sweep gas on hydrogen permeation of supported Pd membranes: Experimental and modeling. Int. J. Hydrogen Energy 2019, 44, 4228–4239. [Google Scholar] [CrossRef]
  19. Basile, A.; Campanari, S.; Manzolini, G.; Iulianelli, A.; Longo, T.; Liguori, S.; De Falco, M.; Piemonte, V. Methane steam reforming in a Pd-Ag membrane reformer: An experimental study on reaction pressure influence at middle temperature. Int. J. Hydrogen Energy 2011, 36, 1531–1539. [Google Scholar] [CrossRef]
  20. Fernandez, E.; Helmi, A.; Coenen, K.; Melendez, J.; Viviente, J.L.; Pacheco Tanaka, D.A.; van Sint Annaland, M.; Gallucci, F. Development of thin Pd–Ag supported membranes for fluidized bed membrane reactors including WGS related gases. Int. J. Hydrogen Energy 2015, 40, 3506–3519. [Google Scholar] [CrossRef]
  21. Di Marcoberardino, G.; Foresti, S.; Binotti, M.; Manzolini, G. Potentiality of a biogas membrane reformer for decentralized hydrogen production. Chem. Eng. Process. Process Intensif. 2018, 129, 131–141. [Google Scholar] [CrossRef]
  22. Pinacci, P.; Drago, F. Influence of the support on permeation of palladium composite membranes in presence of sweep gas. Catal. Today 2012, 193, 186–193. [Google Scholar] [CrossRef]
  23. Marcoberardino, G.D.; Gallucci, F.; Manzolini, G.; van Sint Annaland, M. Definition of validated membrane reactor model for 5 kW power output CHP system for different natural gas compositions. Int. J. Hydrogen Energy 2016, 41, 19141–19153. [Google Scholar] [CrossRef]
  24. Zhao, C.; Xu, H.; Goldbach, A. Duplex Pd/ceramic/Pd composite membrane for sweep gas-enhanced CO2 capture. J. Membr. Sci. 2018, 563, 388–397. [Google Scholar] [CrossRef]
  25. Boon, J.; Pieterse, J.A.Z.; Dijkstra, J.W.; van Sint Annaland, M. Modelling and systematic experimental investigation of mass transfer in supported palladium-based membrane separators. Int. J. Greenh. Gas Control 2012, 11, 122–129. [Google Scholar] [CrossRef]
  26. Li, X.; Li, A.; Lim, C.J.; Grace, J.R. Hydrogen permeation through Pd-based composite membranes: Effects of porous substrate, diffusion barrier and sweep gas. J. Membr. Sci. 2016, 499, 143–155. [Google Scholar] [CrossRef]
  27. Zhao, C.; Caravella, A.; Xu, H.; Brunetti, A.; Barbieri, G.; Goldbach, A. Support mass transfer resistance of Pd/ceramic composite membranes in the presence of sweep gas. J. Membr. Sci. 2018, 550, 365–376. [Google Scholar] [CrossRef]
  28. Zeng, G.; Goldbach, A.; Xu, H. Impact of support mass flow resistance on low-temperature H2 permeation characteristics of a Pd95Ag5/Al2O3 composite membrane. J. Membr. Sci. 2009, 326, 681–687. [Google Scholar] [CrossRef]
  29. Li, H.; Dijkstra, J.W.; Pieterse, J.A.Z.; Boon, J.; Van Den Brink, R.W.; Jansen, D. WGS-mixture separation and WGS reaction test in a bench-scale multi-tubular membrane reactor. Energy Procedia 2011, 4, 666–673. [Google Scholar] [CrossRef] [Green Version]
  30. Fernandez, E.; Helmi, A.; Medrano, J.A.; Coenen, K.; Arratibel, A.; Melendez, J.; de Nooijer, N.C.A.; Spallina, V.; Viviente, J.L.; Zuñiga, J.; et al. Palladium based membranes and membrane reactors for hydrogen production and purification: An overview of research activities at Tecnalia and TU/e. Int. J. Hydrogen Energy 2017, 42, 13763–13776. [Google Scholar] [CrossRef] [Green Version]
  31. Nooijer, N.; Arratibel Plazaola, A.; Meléndez Rey, J.; Fernandez, E.; Pacheco Tanaka, D.; van Sint Annaland, M.; Gallucci, F. Long-Term Stability of Thin-Film Pd-Based Supported Membranes. Processes 2019, 7, 106. [Google Scholar] [CrossRef]
  32. Nayebossadri, S.; Fletcher, S.; Speight, J.D.; Book, D. Hydrogen permeation through porous stainless steel for palladium-based composite porous membranes. J. Membr. Sci. 2016, 515, 22–28. [Google Scholar] [CrossRef]
  33. Zhu, B.; Duke, M.; Dumée, L.F.; Merenda, A.; des Ligneris, E.; Kong, L.; Hodgson, P.D.; Gray, S. Short review on porous metal membranes—Fabrication, commercial products, and applications. Membranes 2018, 8, 83. [Google Scholar] [CrossRef] [PubMed]
  34. Peters, T.A.; Tucho, W.M.; Ramachandran, A.; Stange, M.; Walmsley, J.C.; Holmestad, R.; Borg, A.; Bredesen, R. Thin Pd–23%Ag/stainless steel composite membranes: Long-term stability, life-time estimation and post-process characterisation. J. Membr. Sci. 2009, 326, 572–581. [Google Scholar] [CrossRef]
  35. Hwang, K.R.; Oh, D.K.; Lee, S.W.; Park, J.S.; Song, M.H.; Rhee, W.H. Porous stainless steel support for hydrogen separation Pd membrane; fabrication by metal injection molding and simple surface modification. Int. J. Hydrogen Energy 2017, 42, 14583–14592. [Google Scholar] [CrossRef]
  36. Caravella, A. Dusty-Gas Model (DGM). In Encyclopedia of Membranes; Springer: Berlin, Germany, 2016. [Google Scholar]
  37. Gallucci, F.; Van Sint Annaland, M.; Kuipers, J.A.M. Pure hydrogen production via autothermal reforming of ethanol in a fluidized bed membrane reactor: A simulation study. Int. J. Hydrogen Energy 2010, 35, 1659–1668. [Google Scholar] [CrossRef]
  38. Trimm, D.L.; Lam, C.-W. The Combustion of Methane on Platinum-Alumina Fibre Catalysts-1. Chem. Eng. Sci. 1979, 35, 1405–1413. [Google Scholar] [CrossRef]
  39. Numaguchi, T.; Kikuchi, K. Intrinsic Kinetics and Design Simulation in a Complex Reaction Network; Steam-Methane Reforming. Chem. Eng. Sci. 1988, 43, 2295–2301. [Google Scholar] [CrossRef]
  40. Nam, S.E.; Lee, K.H. Hydrogen separation by Pd alloy composite membranes: Introduction of diffusion barrier. J. Membr. Sci. 2001, 192, 177–185. [Google Scholar] [CrossRef]
  41. Peters, T.A.; Stange, M.; Klette, H.; Bredesen, R. High pressure performance of thin Pd-23%Ag/stainless steel composite membranes in water gas shift gas mixtures; influence of dilution, mass transfer and surface effects on the hydrogen flux. J. Membr. Sci. 2008, 316, 119–127. [Google Scholar] [CrossRef]
  42. Rothenberger, K.S.; Cugini, A.V.; Howard, B.H.; Killmeyer, R.P.; Ciocco, M.; Morreale, B.; Enick, R.; Bustamante, F.; Mardilovich, I.; Ma, Y. High pressure hydrogen permeance of porous stainless steel coated with a thin palladium film via electroless plating. J. Membr. Sci. 2004, 244, 55–68. [Google Scholar] [CrossRef]
  43. Mardilovich, P.P.; She, Y.; Ma, Y.H.; Rei, M.-H. Defect-free palladium membranes on porous stainless-steel support. AIChE J. 1998, 44, 310–322. [Google Scholar] [CrossRef]
  44. Huang, Y.; Dittmeyer, R. Preparation of thin palladium membranes on a porous support with rough surface. J. Membr. Sci. 2007, 302, 160–170. [Google Scholar] [CrossRef]
  45. Veldsink, J.W.; van Damme, R.M.J.; Versteeg, G.F.; van Swaaij, W.P.M. The use of the dusty-gas model for description of mass transport with chemical reaction in porous media. Chem. Eng. J. 1995, 115–125. [Google Scholar] [CrossRef]
  46. Nooijer, N.D.; Gallucci, F.; Pellizzari, E.; Melendez, J.; Tanaka, D.A.P.; Manzolini, G.; van Sint Annaland, M. On concentration polarisation in a fluidized bed membrane reactor for biogas steam reforming: Modelling and experimental validation. Chem. Eng. J. 2018, 348, 232–243. [Google Scholar] [CrossRef]
  47. AspenTech Aspen Plus. Available online: https://www.aspentech.com/en/products/engineering/aspen-plus (accessed on 20 July 2019).
  48. Di Marcoberardino, G.; Vitali, D.; Spinelli, F.; Binotti, M.; Manzolini, G. Green Hydrogen Production from Raw Biogas: A Techno-Economic Investigation of Conventional Processes Using Pressure Swing Adsorption Unit. Processes 2018, 6, 19. [Google Scholar] [CrossRef]
  49. Peng, D.Y.; Robinson, D.B. New Two-Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59–64. [Google Scholar] [CrossRef]
  50. Rackett, H.G. Equation of state for saturated liquids. J. Chem. Eng. Data 1970, 15, 514–517. [Google Scholar] [CrossRef]
  51. Haar, L.; Gallagher, J.; Kell, G.S. NBS/NRC Steam Tables: Thermodynamic and Transport Properties and Computer Programs for Vapor and Liquid States of Water in S.I. Units; Hemisphere Publishing Corporation: Washington, DC, USA, 1984. [Google Scholar]
  52. Tecnalia FERRET WP3—Public Report on Production Costs and Scale Up of the Membrane Production, 2017.
  53. Manzolini, G.; Sanchez Fernandez, E.; Rezvani, S.; Macchi, E.; Goetheer, E.L.V.; Vlugt, T.J.H. Economic assessment of novel amine based CO2 capture technologies integrated in power plants based on European Benchmarking Task Force methodology. Appl. Energy 2015, 138, 546–558. [Google Scholar] [CrossRef]
  54. Di Marcoberardino, G.; Liao, X.; Dauriat, A.; Binotti, M.; Manzolini, G. Life Cycle Assessment and Economic Analysis of an Innovative Biogas Membrane Reformer for Hydrogen Production. Processes 2019, 7, 86. [Google Scholar] [CrossRef]
  55. Jenkins, S. Chemical Engineering Plant Cost INdex: 2018 Annual Value. Available online: https://www.chemengonline.com/2019-cepci-updates-january-prelim-and-december-2018-final/ (accessed on 20 July 2019).
  56. Spallina, V.; Pandolfo, D.; Battistella, A.; Romano, M.C.; Van Sint Annaland, M.; Gallucci, F. Techno-economic assessment of membrane assisted fluidized bed reactors for pure H2 production with CO2 capture. Energy Convers. Manag. 2016, 120, 257–273. [Google Scholar] [CrossRef]
Figure 1. Scheme of the fluidized bed membrane reactor developed in ACM with and without gas for PSS and ceramic supports [16].
Figure 1. Scheme of the fluidized bed membrane reactor developed in ACM with and without gas for PSS and ceramic supports [16].
Membranes 09 00116 g001
Figure 2. Model of the Pd composite membrane [22].
Figure 2. Model of the Pd composite membrane [22].
Membranes 09 00116 g002
Figure 3. Pressure gradients comparison for the considered cases.
Figure 3. Pressure gradients comparison for the considered cases.
Membranes 09 00116 g003
Figure 4. HRF comparison for the considered cases.
Figure 4. HRF comparison for the considered cases.
Membranes 09 00116 g004
Figure 5. Porosity over tortuosity (left side) and pore size (right side) effects on HRF.
Figure 5. Porosity over tortuosity (left side) and pore size (right side) effects on HRF.
Membranes 09 00116 g005
Figure 6. Sweep flow rate impact on the HRF for the considered cases.
Figure 6. Sweep flow rate impact on the HRF for the considered cases.
Membranes 09 00116 g006
Figure 7. Equivalent permeance calculated at 12 and 20 bar feed pressure for different support thicknesses.
Figure 7. Equivalent permeance calculated at 12 and 20 bar feed pressure for different support thicknesses.
Membranes 09 00116 g007
Figure 8. Layout of BIONICO system using sweep gas.
Figure 8. Layout of BIONICO system using sweep gas.
Membranes 09 00116 g008
Figure 9. LCOH comparison for the considered supports.
Figure 9. LCOH comparison for the considered supports.
Membranes 09 00116 g009
Figure 10. Sensitivity analysis on main economic assumptions for Case A (left) and Case B (right). LCOH cost calculated assuming H2 at 700 bar.
Figure 10. Sensitivity analysis on main economic assumptions for Case A (left) and Case B (right). LCOH cost calculated assuming H2 at 700 bar.
Membranes 09 00116 g010
Table 1. Difference between ceramic and porous metallic supports [32,33,34,35].
Table 1. Difference between ceramic and porous metallic supports [32,33,34,35].
Support MaterialAdvantagesDisadvantages
Ceramic supports
  • Variety of materials
  • Pore structures
  • Qualities
  • Brittle
  • Integration is difficult at high temperature
Metallic supports
  • Easy integration
  • Mechanical stability
  • Same thermal coefficient
  • Large open pores
  • Wide size distribution
  • High surface roughness
Table 2. Porous metallic support data.
Table 2. Porous metallic support data.
Pore SizePorosity/TortuositySupport GeometrySurface AreaManufacturerReference
0.5 μm-Thickness 1 mm7.07 cm2Mott metallurgical corporation[40]
2 μm-Thickness 0.48 mm6.8 cm2AccuSep PALL corporation USA[41]
0.2 μmε = 20–23%OD 6.4 mm
ID 3.2 mm
9 cm2Mott Metallurgical Corporation[42]
3 μmε = 17%OD 15.9 and 12.7 mm
Thickness 1.6 mm
-Mott Metallurgical Corporation[43]
0.5 μm 1–2 μm-OD 10 mmID 6 mm-GKN Sinter Metal Filters GmbH
Radevormwald
[44]
0.45 μmε/τ = 1.12OD 21 mm
Thickness 1 mm
-Mott Metallurgical Corporation[32]
0.30 μmε/τ = 2.89
0.87 μmε/τ = 0.42
0.24 μmε/τ = 2.60
Table 3. Ceramic support data.
Table 3. Ceramic support data.
Pore SizePorosityTortuosityThicknessSurface AreaManufacturer Reference
0.16 μmε = 35%τ = 1.2514 mm OD
10 mm ID
21 cm2Provided by ECN[27]
0.16 μmε = 35%τ = 1.2514 mm OD
10 mm ID
21 cm2Provided by ECN[27]
3.4 μmε = 43%τ = 1.2514 mm OD
10 mm ID
21 cm2Provided by ECN[27]
0.36 μmε/τ = 0.151.5 mm29/36.4 cm2-[16]
Table 4. Feed stream and reactor conditions, membrane properties [16,27,32].
Table 4. Feed stream and reactor conditions, membrane properties [16,27,32].
Feed Stream ConditionsValueMembrane CharacteristicsValue
Natural gas feed (kmol/h)0.055Membrane thickness (μm)5
Steam feed (kmol/h)0.165Permeance (kmol s−1 m−2 Pa−n)1.3 × 10−5
Air feed (kmol/h)0.062Ea (J/mol)10,171
Sweep gas (kmol/h)0.058n0.5
Temperature (°C)400Membrane length (m)0.4
Pressure (bar)8–20Outer diameter support (mm)10
S/C ratio3Inner diameter support range (mm)4–9
Membrane area (m2)0.14
Reactor conditions
Temperature (°C)500
Pressure (vessel/permeate side) (bar)8–20/1.3
Table 5. Characteristics of the different cases investigated in this work.
Table 5. Characteristics of the different cases investigated in this work.
CasesSupport TypeOD [mm]ID Range [mm]ε/τrp [μm]
Case 1 [27]α-Al2O3104–90.3443.4
Case 2 [16]α-Al2O3104–90.150.34
Case 3 [32]316 L104–90.420.87
Case 4 [32]316 L104–92.60.24
Table 6. Stream properties at in/exit of reactors (@pressure and T).
Table 6. Stream properties at in/exit of reactors (@pressure and T).
StreamFlowT (°C)p (bar)Composition (% Molar Basis)
Molar
(mol/s)
Mass
(g/s)
CH4H2COCO2H2OO2N2
11.1627.355352024.9--19.245.41.59.0
20.329.3252020-----2179
31.2435.51550206.45.42.032.924.3-29.1
40.928.2930.1208.67.42.741.10.2-39.9
52.3270.763351.1---20.29.64.765.5
61.1511.565501.1-50.0--50.0--
Table 7. New reactor dimensioning compared to BIONICO.
Table 7. New reactor dimensioning compared to BIONICO.
ParametersCase 1:
(rpore = 3.4 μm)
ε/τ = 0.344
Case 2:
(rpore = 0.36 μm)
ε/τ = 0.15
Case 3:
(rpore = 0.87 μm)
ε/τ = 0.42
Case 4:
(rpore = 0.24 μm)
ε/τ = 2.60
Thickness (mm)220.50.5
Nmem (–)267438223203
Amem (m2)3.365.502.802.55
Vvessel (m3)0.250.410.220.20
mcat (kg)165175135120
HRF (%)62.4062.8062.4062.50
Table 8. Methodology for the TOC calculation [56].
Table 8. Methodology for the TOC calculation [56].
Plant ComponentCost (M €)
CompressorA
Heat exchangerB
ReactorC
Bare erected cost (BEC)A + B + C
Direct costs as percentage of BEC
Includes piping/valves, civil works, instrumentation, steel structure, erections, etc.
Total installation cost (TIC)80% BEC
Total direct plant cost (TDPC)BEC + TIC
Indirect costs (IC)14% TDPC
Engineering procurement and construction (EPC)TDPC + IC
Contingencies and owner’s costs (C&OC)
Contingency5% EPC
Owner’s cost10% EPC
Total C&OC15% EPC
Total plant cost (TPC)EPC + C&OC
Table 9. Cost assumption for plant components.
Table 9. Cost assumption for plant components.
ComponentsAmountScaling ParameterS0C0 (k€)fYear CostCEPCI
Reactor1Weight (lb)130,00070.320.32007525.4
Heat Exchanger8Exchange area (m2)215.50.592007525.4
Biogas compressor1Power (kW)53.30.822006499.6
Air compressor1Power (MW)0.683.420.672009521.9
Water demineralizer1Water flow rate (lH2O/h)902.10.682011585.7
Water pump2Water flow rate (lH2O/h)901.20.72011585.7
Table 10. Assumptions to calculate the O&M costs [48,52,56].
Table 10. Assumptions to calculate the O&M costs [48,52,56].
O&M—Fixed
Labor costs60,000 €
Maintenance costs2.5% TOC
Insurance2.0% TOC
O&M—Variable
Catalyst cost258 €/kg/y
Filler particles12 €/kg/y
Membrane replacement ceramic360 €/m2/y
Membrane replacement SS2040 €/m2/y
Deionization Resin90 €/y
Lifetime5 Years
Process water0.35 €/m3
Biogas cost1.50 €/GJLHV
Electricity cost0.12 €/kWh
Table 11. CAPEX and OPEX for the different cases considered.
Table 11. CAPEX and OPEX for the different cases considered.
ComponentsCeramicStainless Steel
Case 1:
(rpore = 3.4 μm)
ε/τ = 0.34
Case 2:
(rpore = 0.36 μm)
ε/τ = 0.15
Case 3:
(rpore = 0.87 μm)
ε/τ = 0.42
Case 4:
(rpore = 0.24 μm)
ε/τ = 2.60
CAPEX
Reactor cost (k€)29.534.228.427.6
Membranes (k€)22.837.428.626.0
Heat exchangers (k€)73.973.973.973.9
Biogas compressors (k€)3.93.93.93.9
Balance of plant (k€)0.30.30.30.3
H2 compr @20 bar (k€)4.44.44.44.4
TPC @20 bar (k€)318.2363.7329.0321.1
H2 compr @700 bar (k€)22.122.122.122.1
TPC @700 bar (k€)370.3415.8381.2373.3
OPEX
Catalyst + filler22.223.618.216.2
Biogas9.19.19.19.1
Water cost0.20.20.20.2
Electricity @20 bar18.518.518.518.5
Membranes4.67.55.75.2
Deionization resin0.450.450.450.45
O&Mvar. total @ 20 bar55.059.352.149.6
Electric energy @ 700 bar21.621.621.621.6
O&Mvar. total @ 700 bar76.781.073.871.3
O&Mfix @ 20 bar74.3276.3674.8174.45
O&Mfix @ 700 bar76.6678.7177.1576.80

Share and Cite

MDPI and ACS Style

Di Marcoberardino, G.; Knijff, J.; Binotti, M.; Gallucci, F.; Manzolini, G. Techno-Economic Assessment in a Fluidized Bed Membrane Reactor for Small-Scale H2 Production: Effect of Membrane Support Thickness. Membranes 2019, 9, 116. https://doi.org/10.3390/membranes9090116

AMA Style

Di Marcoberardino G, Knijff J, Binotti M, Gallucci F, Manzolini G. Techno-Economic Assessment in a Fluidized Bed Membrane Reactor for Small-Scale H2 Production: Effect of Membrane Support Thickness. Membranes. 2019; 9(9):116. https://doi.org/10.3390/membranes9090116

Chicago/Turabian Style

Di Marcoberardino, Gioele, Jasper Knijff, Marco Binotti, Fausto Gallucci, and Giampaolo Manzolini. 2019. "Techno-Economic Assessment in a Fluidized Bed Membrane Reactor for Small-Scale H2 Production: Effect of Membrane Support Thickness" Membranes 9, no. 9: 116. https://doi.org/10.3390/membranes9090116

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