Next Article in Journal
Assessment of the Potential of Using Nanofiltration Polymeric and Ceramic Membranes to Treat Refinery Spent Caustic Effluents
Next Article in Special Issue
Polybenzimidazole Membrane Crosslinked with Epoxy-Containing Inorganic Networks for Organic Solvent Nanofiltration and Aqueous Nanofiltration under Extreme Basic Conditions
Previous Article in Journal
Low-Temperature Hydrophilic Pervaporation of Lactic Acid Esterification Reaction Media
Previous Article in Special Issue
Tailoring the Stabilization and Pyrolysis Processes of Carbon Molecular Sieve Membrane Derived from Polyacrylonitrile for Ethylene/Ethane Separation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation and Optimization of 4-Component LDG Separation in the Steelmaking Industry Using Polysulfone Hollow Fiber Membranes

C1 Gas & Carbon Convergent Research Center, Korea Research Institute of Chemical Technology, 141 Gajeongro, Yuseong-gu, Daejeon 34114, Korea
*
Author to whom correspondence should be addressed.
Membranes 2022, 12(1), 97; https://doi.org/10.3390/membranes12010097
Submission received: 4 December 2021 / Revised: 4 January 2022 / Accepted: 11 January 2022 / Published: 17 January 2022
(This article belongs to the Special Issue Recent Membrane Research and Development in Korea)

Abstract

:
A general finite element model and a new solution method were developed to simulate the permeances of Lintz Donawiz converter gas (LDG) components and the performance of a polysulfone membrane separation unit. The permeances at eight bars of CO, N2, and H2 in LDG simulated using the developed model equations employing the experimental mixed gas data were obtained by controlling the finite element numbers and comparing them with pure gas permeation data. At the optimal finite element numbers (s = 15, n = 1), the gas permeances under the mixed-gas condition were 6.3% (CO), 3.9% (N2), and 7.2% (H2) larger than those of the pure gases, On the other hand, the mixed-gas permeance of CO2 was 4.5% smaller than that of pure gas. These differences were attributed to the plasticization phenomenon of the polysulfone membrane used by CO2. The newly adopted solution method for the stiff nonlinear model functions enabled the simulation of the performance (in terms of gas recovery, concentration, and flow rate) of the first-stage membrane within two seconds under most gas flow conditions. The performance of a first-stage membrane unit separating LDG could be predicted by the developed model with a small error of <2.1%. These model and solution methods could be utilized effectively for simulating gas permeances of the membrane that is plasticized severely by the permeating gas and the separation performance of two- or multi-stage membrane processes.

1. Introduction

Emissions of greenhouse gases into the atmosphere are still growing despite the enormous efforts to mitigate them. In Korea, CO2 emissions from the combustion of fossil fuels increased from 594 to 704 million tons from 2010 to 2018 [1]. The Intergovernmental Panel on Climate Change has reported that anthropogenic CO2 emissions should be reduced to net-zero by 2050 relative to 2010 levels to limit global warming to 1.5 °C above pre-industrial levels [2]. We focused on the separation and utilization of CO2 and CO in the Lintz Donawiz converter gas (LDG, 64% CO, 18% CO2, 16% N2, and 2% H2) in the steelmaking industry while a variety of attempts have been devoted toward the reduction of CO2 emissions. In Korea, LDG is combusted to generate electricity with other byproduct gases in the steel industry, such as coke oven gas (COG) and blast furnace gas (BFG), causing a large amount of CO2 emissions. The main components of LDG are CO and CO2. Therefore, CO2 emissions to the atmosphere can be reduced by recovering CO and CO2 and utilizing them as raw materials to produce various chemical products.
Several separation technologies, such as cryogenic distillation, absorption, and adsorption, have been considered for industrial applications to separate CO2 and CO concentrations from the LDG mixed gases. CO2 separation using separation technologies has been extensively studied and commercialized. On the other hand, CO separation using adsorbents or absorbents using copper transition metals was developed and commercialized in 2000, including the chemical absorption process by R.C. Costello & Associate Inc., USA (Cosorb or Copure process) and pressure swing adsorption (PSA) process by Kobe Steel Co., Japan and Pioneer Technology Inc., China. However, these methods require high separation energy in the processes due to the strong binding of the adsorbents or adsorbents with CO. Moreover, their unsteady operation behavior due to low CO/CO2 selectivity of Cu(I) instability and impurity poisoning from hydrogen sulfide or water vapor has severely limited their scale-up and wide commercial applications.
To solve these problems, we have designed a new hybrid-type CO purification process employing membrane separation to remove CO2 and H2 selectively in LDG and recover high-purity CO from N2 residues using a typical PSA, which is fed to membrane modules to separate it into a residue stream containing mostly CO and N2 and a permeate stream consisting of CO2 and H2. Then, the residue CO/N2 mixture is sent to a PSA or chemical absorption process unit to separate it into high purity CO above 99% and the residue N2. Ag(I)-based facilitated transport membrane processes could also be an alternative to the typical PSA to separate the residue CO/N2 mixture. To accomplish this task effectively, a high-performance membrane that selectively permeates CO2 and H2 in LDG, the process simulation, and design techniques must be developed in which the membrane separation of LDG uses a hollow fiber membrane. This includes the mathematical modeling and simulation required for the scale-up of the membrane separation process.
Since the early 1970s, various mathematical models and solution methods for membrane gas separation processes have been reported in the literature [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. The models can be classified into two categories according to their equation forms: differential and finite element models. The differential models, where the governing equations of motion are expressed in the form of differential equations, can be conveniently applied to binary systems because all equations can be solved analytically. The differential models are still useful for multi-component concurrent systems where feed and permeate gases flow concurrently because their governing equations can be solved step by step by numerical integrations without any assumptions. However, in the case of multicomponent countercurrent systems, integrating differential equations numerically iteratively (trial-and-error shooting method) [28,29,30,31], the flow rates and concentrations of each component in the permeate leaving membrane module are very cumbersome and time-consuming because the nonlinear equations of membrane permeation are very stiff. To avoid integration, various methods that convert the differential equations into algebraic forms, such as series approximation [32,33,34], discretization [35], and orthogonal collocation [36,37] approaches, have been suggested.
The derived governing equations of the finite element models have algebraic forms. These models regarded the membrane as being composed of a number of finite elements having the same area in which both permeate, and residue side gases are in a perfectly mixed state. Therefore, the gas concentrations are all the same at any position within a finite element. The governing equations obtained from the mass balance around the finite elements are composed only of algebraic variables. These models can be applied effectively to hollow fiber membrane modules where gas mixing occurs compared to differential models that assume, in general, plug-flow on both permeate and residue sides. With the extent of gas mixing, the number of finite elements can be adjusted identically or differently on both sides, which is an advantage not found in other models. In a hollow fiber membrane module, if gas mixing takes place, it would be more likely on the shell side than on the bore side because the shell side has more free space for gas to move. In this case, it is desirable that the bore and shell sides have different finite element numbers from each other.
In this study, we investigated the general finite element model with an arbitrary ratio n of the bore to the shell side finite element number to determine the value of n at which the model agrees best with the experimental data. A set of derived algebraic equations were solved by the NEQNF subroutine in the IMSL library, which is a method of solving a system of nonlinear equations using the modified Powell hybrid algorithm and a finite-difference approximation to Jacobian. Using the subroutine, we can solve nonlinear equations much faster and easier compared to other methods suggested in previous studies. Based on the model, two Fortran programs have been developed for the simulation, one for calculating the permeances of LDG mixed gas using the experimental data obtained under one experimental condition and the other for predicting the performance of a used membrane (flow rates and concentrations of the LDG components in residue and permeate) using the permeance data calculated using the former program.
Generally, the permeances of mixed gas used in simulations of membrane separation processes are those of pure gases [29,30,36], those calculated using linear cross-flow models [38,39], or those obtained by experiments at a low stage-cut of less than 1% [40,41,42] with the assumption that the concentrations of feed and residue gases are the same or at a stage-cut of less than 20% [10,43,44,45]. The latter permeances were calculated using the logarithmic mean concentrations of feed and residue gases as residue-side concentrations. However, these are all approximate values that are not acquired under the real operating conditions of the membrane modules. Therefore, in this study, we calculated the permeances of the LDG mixed gas using the data obtained under real operating conditions of the membrane module.

2. Materials and Methods

2.1. Mathematical Model

2.1.1. Bore Side Feed

Katoh et al. [46] developed a finite element model for a hollow fiber membrane module called the “tank-in-series model.” Here, the value of n was fixed at 2. We deal with the most general case of this model, where n is not specified at a certain value for an asymmetric hollow fiber membrane. For each finite element of the shell side, the countercurrent hollow fiber membrane is divided into finite elements on the shell side and h finite elements on the bore side, as shown in Figure 1. For convenience, the k-th element is verified.
The following equations are obtained by mass balance and fluid dynamics for the k-th element of the bore side feed mode (Figure 1a) under isothermal and steady-state conditions if the membrane does not deform under pressure and the pressure drop of the bore side follows the Hagen–Poiseuille law.
(1)
The permeation rate of the i-component in the k,j-th element of the bore side is
g k , j , i = Q i A j ( P k , j x k , j , i p y k , i ) .
Summing g k , j , i for all the components (c) gives
G k , j = i = 1 c g k , j , i
The summation of Equation (1) for the i-component on the bore side is:
G C k , i = j = 1 h g k , j , i
The summation of the permeation rates of all the components in the bore side is
G T k = i = 1 c G C k , i
(2)
Component and total mass balances
  • The k,j-th element on the bore side is expressed as
    u k , j 1 x k , j 1 , i u k , j x k , j , i g k , j , i = 0
    u k , j 1 u k , j G k , j = 0
  • Meanwhile, the k,j-th element in the permeate side is given by
    v k + 1 y k + 1 , i v k y k , i + G C k , i = 0
    v k + 1 v k + G T k = 0
The pressure drop between the (k,j − 1)-th and k,j-th elements is
P k , j 1 P k , j = 128 u k , j μ m L j π D o 2 ,
where μm is the viscosity of the residue side gas calculated using the Wike equation [47].

2.1.2. Shell Side Feed

The lengthwise and backmixing of gases is very restrictive within the porous supporting layer of the asymmetric hollow fiber membrane. Therefore, gases can pass through the layer in the cross-flow type. When gases flow in this manner, a cross-flow model can be set up for the shell side feed mode. However, for the bore side feed, this model is difficult to realize mathematically. The concentrations of permeate are different at the wall of the selective layer, in which high permeability gases permeate selectively, and in the bulk permeate stream of the bore side. If the concentration of the i-component is yi at the selective layer wall and y ¯ i in the permeate bulk stream (Figure 1b) the governing equations for the cross-flow model are as follows:
(1)
The permeation rates in the k-th element:
g k , j , i = Q i A j ( p x k , i P k , j y k , j , i )
G k , j = i = 1 c g k , j , i
G C k , i = j = 1 h g k , j , i
G T k = i = 1 c G C k , i
(2)
The component and total mass balances in the k-th and k,j-th elements:
v k , j + 1 y ¯ k , j + 1 , i v k , j y ¯ k , j , i + g k , j , i = 0
v k , j + 1 v k , j + G k , j = 0
u k 1 x k 1 , i u k x k , i G C k , i = 0
u k 1 u k G T k = 0
(3)
The concentration yk,j,i at the k,j-th element:
y k , j , i = g k , j , i G k , j = Q i A j ( p x k , i P k , j y k , j , i ) i = 1 c Q i A j ( p x k , i P k , j y k , j , i )
(4)
The pressure drop between the (k,j − 1)-th and k,j-th elements:
P k , j + 1 P k , j = 128 v k , j μ m L j π D o 2 ,
where μm is the viscosity of the permeate side gas mixture.

2.2. Solution of the Model Equation

Katoh et al. [46] solved the nonlinear equations of their finite element model using the relaxation method, which has been applied to distillation processes. On the other hand, Coker et al. [48,49] used the Thomas algorithm after converting a set of flow rate equations into a tridiagonal matrix form. However, as mentioned previously, nonlinear equations can be solved more conveniently and rapidly by the NEQNF subroutine than by other methods. The subroutine requires only rough initial values for all variables to start the calculation. Given the initial values, it finds the final values of all the variables satisfying a given error tolerance, iteratively modifying the Jacobians of coded equations. Before coding, nonlinear equations must be rearranged such that one side is equal to zero, as shown in Equations (5)–(8). Subsequently, they must be set to variables F(x), which then becomes an error in each equation. The NEQNF subroutine provides the sum of F(x)2 instead of F(x), which can have negative real values. From F ( x ) 2 , the root mean square error of the NEQNF method can be calculated as
F ( x ) 2 x ,
where x is the total number of coded equations.
For the bore side feed operation, there are a total of 11 unknown variables, except for the feed compositions, feed rate, bore side pressure, and permeances of the gas components, when s = 1, n = 1, and c = 4. In the case of the shell side feed, there are four additional unknown variables (i.e., the permeate concentrations at the wall of the selective layer yi). We can find the values of the unknown variables (concentrations, flow rate, and pressure) by solving 11 and 15 equations from Equations (5)–(9) and (14)–(19) for the bore side feed and shell side feed, respectively. In contrast, if there are experimental data for the concentrations and flow rates of the residue and permeate leaving the membrane, we can calculate the mixed gas permeances of the four components. The number of unknown variables is five at s = 1 and n =1. Therefore, the permeances can be calculated using five equations from Equations (9) and (5) or Equation (7) for the bore-side feed. For the general case (i.e., s = s, h = h, and c = 4), there are a total of (6sh + 5s) and (6sh − 1) variables for the performance-simulating model that predicts concentrations and flow rates and permeance-simulating model for the bore side feed, respectively. The performance and gas permeance of the membrane can be obtained using as many equations from Equations (5)–(9) as the variables.

2.3. Experimental

Figure 2 shows a schematic of the membrane separation unit installed for LDG separation. To develop the model, polysulfone (PS) hollow fiber modules (Synopex Inc., Ltd., Hwaseong, Korea) were used as the membrane, whose specifications are listed in Table 1. In the LDG model gas, the gas mixture (64% CO, 18% CO2, 16% N2, and 2% H2 by volume) was made by flowing pure gases with mass flow controllers and fed to the bore side of the membrane module. The mass flow controllers were recalibrated with a standard wet gas meter (Shinagawa Co., Tokyo, Japan) before use. The membrane module was maintained at 20 °C by a heat exchanger that controlled the temperature of the feed gas and an air conditioner installed in the equipment room. The gas lines were connected to ensure that the membrane operated in countercurrent mode. The flow rates and concentrations of the residue and permeate leaving the membrane were measured with wet gas meters and an Agilent 6890 GC (Agilent Technologies, Inc., Santa Clara, CA, USA) using CH4 as an external standard gas. For precise measures of flow rates, the wet gas meters were installed in a constant temperature oven kept at 3 °C to suppress as much as possible the vaporization of their working fluid (water) into the process gas streams, which increased the gas flow rates. After the analysis, all waste gases were burned in a flare stack.

3. Results and Discussion

Glassy polymers such as PS and PI are plasticized by penetrant CO2 [38,50,51,52,53,54,55]. The permeabilities of other gases increase if a membrane is swollen by CO2. Therefore, their permeances differ when permeating with CO2 and as pure gases. The permeances of pure gases are not suitable for simulating CO2-containing mixed-gas separation processes. Here, we calculated the permeances of LDG mixed gas using the model with experimental data obtained at one feed flow rate condition. Then, the effectiveness of the approach was verified through experiments under other conditions. Two Fortran programs were prepared, one for calculating the gas permeances and the other for simulating the performance of the PS membrane.
Table 2 presents the experimental results at a feed flow rate of 10 L/min, feed pressure of eight bars, and permeate pressure of one bar.
Using these data, the permeances of the LDG components were calculated using the permeance-simulating program. Table 3 summarizes the simulated results while varying the finite element numbers of the bore and shell sides. In addition, the permeances and selectivities of the pure gases are also included in the table for comparison. The results show that the permeances of the more permeable gas CO2 are lower than those of pure CO2 in the mixed gas state, while those of the less permeable gases CO and N2 are slightly larger than those of the pure gases. In addition, the CO2 permeances decreased by 0.2-10.3% at the investigated finite element numbers, while those of CO and N2 increased by 5.9–9.67% and 3.5–4.7%, respectively, in comparison to those of pure gases. Considering that the permeance differences between the pure and mixed gases are not large, the plasticization effect is not significant in the used PS membrane. In our CO recovery process, the CO2/CO and CO/N2 selectivities are important properties. The CO selectivity for CO2 and N2 increased when permeated with CO2. Accordingly, the actual CO recovery is lower than predicted with the selectivity data of their pure gases. This means that gas mixing occurs over a wider range in the membrane as the number of finite elements decreases. In contrast, the gas flow pattern in the membrane module becomes closer to the plug-flow type as the number of finite elements increases. In general, the performance of the membrane is better in the plug-flow mode than in the mixed-flow mode. According to the simulation results, the permeances of the more permeable gases CO2 and H2 decrease with an increase in the finite element number, and those of the less permeable gases CO and N2 exhibit the opposite trend. This is due to the input performance data required for the simulation to be fixed at the values listed in Table 2. The performance of the membrane would increase and the permeances of more permeable gases should decrease if the number of finite elements increases to satisfy the input performance data.
Table 4 presents the experimental results at feed rates of 5, 20, and 30 L/min, a feed pressure of 8 bars, and a permeate pressure of one bar. We investigated the finite numbers at which the model best fits the experimental data presented in Table 3 and Table 4. In the LDG separation process, the most important performance measures of the membrane are CO recovery (the percentage of the CO flow rate in the residue to the flow rate in the feed gas) and CO and CO2 concentrations in the residue. Table 5 shows the simulation errors with respect to the experimental data for the performance measures of the PS membrane with the variation in feed flow rates. The simulations were performed with permeance data of s = 15 and n = 1, as shown in Table 3. Based on the calculated percent errors, the root mean square percent errors (RMSPE) for the experimental results of the four feed flow rates were calculated as follows:
( p e r c e n t   e r r o r ) 2 4 .
Table 5 lists the results. In this manner, the RMSPEs were calculated at other values of s and n. We can see that the model agrees very well with the experimental data. The RMSPEs for CO recovery and CO and CO2 concentrations in the residue were 1.42%, 0.16%, and 2.12%, respectively, in the feed flow rate range of 5–30 L/min, respectively. The RMSPE for the CO2 concentration in the residue was relatively large because the absolute values of the CO2 concentrations in the residue were small. When the gas concentration is small, a small experimental error can significantly affect the percent error.
Figure 3 shows the graphs of RMSPE against s at n = 1. The variation of RMSPEs for CO recovery and CO concentration in the residue was very small in the simulated range of s = 10–30. However, the RMSPE for the concentration of the more permeable gas CO2 varies considerably according to the s value. The lowest RMSPE for the concentration of CO2 in the residue was attained at the number of shell side finite elements s = 15. Overall, at s = 15, the RMSPEs for the major performance measures of the membrane were the lowest. As shown in Figure 4, at s = 15, the RMSPE for the CO2 concentration increases significantly with an increase in n, while the RMSPEs for CO recovery and CO concentration are nearly constant. From these results, the optimal values of s and n, where the RMSPEs are the lowest, are 15 and 1, respectively. In conclusion, the finite element model agrees best with the experimental data when there are 15 finite elements on both the shell and bore sides. It was expected that the model would fit better at an n value greater than one because the shell side has more free space for gases to move. However, it fits best at the same finite number of 15 on both sides. The RMSPEs of the three major performance measures of the membrane were below 2.1% for the four feed flow rates, as shown in Table 5. The agreement between the model and the experimental data is also excellent for the concentrations of the components in the permeate, as shown in Figure 5.
Table A1 (Appendix A) presents the RMSPEs for the major performance measures calculated similarly to what is outlined above using the permeance data of pure gases at various s. The RMSPEs for the CO concentrations in the residues were below 2.24% in the range of s = 2–60. Although the RMSPE for the CO2 concentrations in the residue was the lowest (2.26%) at s = 15, it significantly increased up to 171% at other values of s. The RMSPEs for CO recovery are nearly constant at approximately 3.06%. These figures are not as large as the root mean square errors, but the low values result from the low percentage errors at high feed flow rates. The percent errors between the model and experimental data at s = 15 and n = 1 are −5.63, −1.55, −1.22, and −1.21% at feed flow rates of 5, 10, 20, and 30 L/min, respectively. Therefore, we cannot accurately predict the performance of the membrane at low feed flow rates using the permeance data of pure gases. In conclusion, the fluctuation in RMSPE for the CO2 concentration in the residue was severe when pure gas permeances were used in the simulation, and the CO recovery could not be predicted with low error at low feed flow rates, indicating that pure gas permeances are inappropriate as permeance data for simulating the performance of a glassy polymer membrane.
In general, the permeability of gas increases as the pressure increases due to the increase in the diffusivity through the membrane [56,57,58,59]. The permeances of CO2, CO, and N2 in the mixed-gas state, which were calculated from the experimental data using the model at the optimal finite element numbers (s = 15, n = 1) in Table 3, increased marginally and linearly with the feed pressure in the range of 4–10 bars, as shown in Figure 6. However, the permeance of H2 increases exponentially. Any regression equations obtained from the data in Figure 6 can be incorporated into the model to predict the performance of the membrane over a wide pressure range. The simulated permeances at 8 bars of CO, N2, CO2, and H2 were 5.3, 3.2, and 173.7 GPU, respectively. These results are 6.3%, 3.9%, and 7.2% larger than the permeance of each pure gas, respectively. On the other hand, the permeance (78.3 GPU) of CO2 in LDG was decreased 4.5% compared to that (84.0 GPU) in pure gas.
Figure 7 shows the simulated CO recoveries and CO2 concentrations in the residue as a function of feed flow rates at the four feed pressures. Simulations were conducted using the optimal finite element.
For the one-stage membrane process. The feed pressure was investigated below 10 bars due to the high energy cost of gas compression. We observed that the CO recovery and CO2 concentration in the residue moved in the same direction at a feed pressure of 6–10 bars. The CO recovery and CO2 concentration become lower as the feed pressure becomes higher. Therefore, there is no condition that can simultaneously satisfy a high CO recovery of >95% and a low CO2 concentration of less than 1% at feed pressures below 10 bars. The maximum CO recovery was approximately 61% when the condition of satisfying CO2 concentration was less than 1% at a feed pressure of eight bars and a feed flow rate of 5 L/min.
The permeation rate through the membrane can be increased by drawing a vacuum on the permeate side because the permeation driving force increases by reducing the permeate side partial pressure. Simulations were carried out at the optimal finite element numbers for the vacuum operation of the permeate side in the range of 0.001–0.1 bars at a feed pressure of eight bars to determine how much the performance would improve. Compared to the atmospheric operation, the CO recovery decreased by 5.4–0.8% and the CO2 concentration in the residue decreased by 0.9–4.2% at feed flow rates of 5–30 L/min and the permeate side pressure of 0.1 bar. However, at lower vacuum pressures, they did not further decrease, as shown in Figure 8. Considering the energy cost of vacuum operation, the optimal pressure is 0.1 bar if the permeate side is operated in the vacuum mode. In the LDG separation process to recover CO of chemical grade, the most important figure of merit is CO purity. In the membrane separation step of the membrane-PSA hybrid process, the desired CO2 concentration level in the residue was below 1%. The maximum CO recovery satisfying this figure of merit is approximately 78% at a feed pressure of eight bar and a feed flow rate of 10 L/min when the permeate side is operated in the vacuum mode. Compared with the atmospheric operation of the permeate side, the CO recovery and the LDG separation rate increased from 61% to 78% and 5 to 10 L/min, respectively, by vacuuming of the permeate side. Figure 9 shows the simulation results at a permeate pressure of 0.1 bar with varying feed pressures. The CO2 concentration in the residue decreased with increasing feed pressure. The maximum feed flow rate satisfying a CO2 concentration of less than 1% is 12.6 L/min at a feed pressure of 10 bars. Compared to the simulation results at a feed pressure of 8 bars, the LDG separation rate increased by 2.6 L/min. However, the CO recovery decreased when the feed pressure was increased to 10 bars and the feed pressure at the permeate side pressure is 0.1 bar.
One of our research objectives for the membrane separation of LDG is to recover CO gas with less than 1% CO2 with a recovery of >95%. However, the objective cannot be attained in the one-stage membrane process, even if the feed pressure can be raised to 50 bars, indicating that a multi-stage process is required to achieve this goal. At present, a two-stage membrane separation process is under development.
In an asymmetric hollow fiber membrane, the shell side feed mode can raise the feed gas pressure to a higher level than in the bore side feed operation. However, with a short length of the membrane module, the performance is lower than that of the bore side feed operation due to the non-ideal gas flow pattern formed from entering and exiting the membrane module in the vertical direction. In addition, the dead volume between the tube sheets of the hollow fibers and gas line fittings affects the performance of a short-length membrane. As previously mentioned, gases can flow in a cross-flow pattern in an asymmetric hollow fiber membrane because of the porous support. Based on this flow, asymmetric models have been developed in the literature, although there are many cases where symmetric models still fit better [27,60]. Based on the simulation results, it appears that cross-flow and concentration polarization do not occur in the PS membrane module used herein. However, the shell side feed mode was simulated using the cross-flow model developed in Section 2.2 under the assumption that the flow occurs to determine the extent to which the cross-flow affects the performance. Figure 10 presented the comparison of the CO2 concentrations in the residue in the shell side and bore side feed modes simulated under the same conditions (i.e., feed pressure of 8 bars, permeate pressure of 1 bar, s = 15, and n = 1). Although the CO2 concentration in the shell side feed mode was higher than that in the bore side feed mode, the difference was less than 1% in the feed flow range of 5–30 L/min. In the case of CO recovery, the difference was small enough that it can be practically ignored. The gas follow patterns in the used membrane do not substantially affect the major performance measures for LDG separation.
One of the advantages of the finite element model is its ability to simulate various types of membrane modules and flow patterns (i.e., flat, hollow fiber, and spiral wound membranes, and mixed- and plug-flow membrane modules). One flat membrane module in which the gases are perfectly mixed can be simulated with the values of s and n as one. On the other hand, the plug-flow mode is realized when the finite element number is very large. The spiral wound membranes can be simulated using the cross-flow model described in Section 2.2.
Figure 11 shows the simulation results based on the finite element number of the shell side at a feed pressure of eight bars, permeate pressure of one bar, and n = 1. The CO2 concentration in the residue decreases as the finite elements increases from 1, which corresponds to a flat membrane, to 200. However, the rate of decline slows significantly after s = 15, and there are no noticeable changes from s = 100. The maximum difference in the CO2 concentrations between s = 1 (flat membrane) and s = 200 (plug flow membrane) was 5.1% at a feed rate of 10 L/min. The CO2 concentrations at the optimal finite number (s = 15) of the LDG separation are depicted with a blue line in Figure 11. The length per finite element is 1.6 cm and the active length of the membrane module is 24 cm. The finite element model assumes that the perfect mixing of gases occurs inside the zone of 1.6-cm length. Although perfect mixing occurs at every length of 1.6 cm, it is apparent that the CO2 concentrations are close to those of the plug flow membrane (s = 200). The differences are only 0.03–0.3% in the feed flow rate range of 5–40 L/min. Meanwhile, the recovery of the less permeable CO gas was insensitive to the finite element number. The values at s = 1 and 200 were almost the same.
The numbers of model equations required to obtain the data in Figure 11 based on the finite element number are presented in Table A2 (Appendix A). The CPU run times taken to solve these equations, F 2 and the root-mean-square errors calculated using Equation (20) are listed in the table. The CPU run times were obtained using the Fortran command “call cpu_time()” on an Intel Core™ 3.3 GHz personal computer. The solution method for the model equations using the NEQNF subroutine is powerful and effective. At the best finite element number, it took only 0.02 s within a root mean square error of 8.4 × 10−19 to solve 110 equations.

4. Conclusions

Predicting the performance of a membrane process separating a multicomponent gas mixture containing CO2, which is a glassy polymer-plasticizing gas, using the permeance data of pure gases or by approximation methods is not a good approach. The permeance data of mixed gas is more accurate when measured under real operating conditions of the membrane. In this study, a general finite element model whose element numbers can be adjusted freely in both bore and shell sides was developed to simulate the permeances of LDG mixed gas and the performance of the PS membrane module. Moreover, the permeances of the LDG component gases were calculated using the model at various finite element numbers and the experimental performance data (component gas concentrations and total flow rates at the bore and shell side outlets) obtained at one feed flow rate. At the optimal finite element numbers (s = 15, n = 1), the simulated permeances at eight bars of CO, N2, CO2, and H2 were 5.3, 3.2, and 173.7 GPU, respectively, when the finite element number was 15. These results are 6.3%, 3.9%, and 7.2% larger than the permeance of each pure gas, respectively. On the other hand, the permeance (78.3 GPU) of CO2 in LDG was decreased 4.5% compared to that (84.0 GPU) in pure gas.
The differences between the pure and mixed gas permeances were due to the plasticizing effect of CO2 on the glassy polymer polysulfone. The permeance difference was the most remarkable in the case of H2, which has the smallest molecule size among the gases. There are many industrial gases in which high concentrations of H2 and CO2 are present. The representative gases are synthetic gas and many industrial byproducts such as COG, Finex gas, and petrochemical byproducts. Moreover, the pure gas permeance data could not be used when we designed a membrane process for H2 separation from these mixed gases. In this case, our model can provide more precise information on the permeance.
Then, we can obtain the performance data of the PS membrane at arbitrary feed flow rates using the same model used to simulate the permeances. Among the permeance data, the permeances simulated at a finite number of 15 on both the bore and shell sides of the membrane gave the best results in predicting the membrane performance. The differences in the CO recoveries and CO concentrations at the residue outlet between the simulated and experimental results were in the range of −1.8 to 2.0% at four feed flow rates of 5, 10, 20, and 30 L/min. The differences increased to approximately −7.3–6.9% under the same conditions when pure gas permeances were used for the simulation instead of the mixed gas permeance mentioned previously.
In this study, the NEQNF subroutine was introduced to solve the large number of nonlinear stiff equations of the model. The newly adopted solution method allows a quick calculation of the equations with an extremely small error. The mixed-flow and plug-flow membrane modules can be simulated within a few seconds with a root mean square error of less than 1.0 × 10−16. To our knowledge, it takes approximately one to two days to accurately simulate 2- and 3-stage membrane processes with a differential model using the trial-and-error shooting conventional solution method. Therefore, the model and solution method could be effectively exploited in multistage and large-scale membrane processes, where their mixed gas permeances can be calculated using the experimental data obtained under real operating conditions.
It is desirable to operate the LDG recovery process below 10 bars because of the high energy cost of the compressor. In the 1-stage PS membrane module with an active area of 1.0 m2, the achievable CO recovery with a CO2 concentration of less than 1% is in the range of 60–80% at this pressure condition. Moreover, a multistage process is required to enhance CO recovery. At present, a two-stage recirculation process and a mini-pilot plant have been designed and constructed using the developed model and solution technique, and are operated to separate and concentrate CO2 and CO from LDG.

Author Contributions

Conceptualization, J.-H.K. and J.-Y.J.; methodology, J.-H.K., J.-Y.J. and B.-R.P.; software, J.-Y.J.; validation, J.-H.K. and J.-Y.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by “Next Generation Carbon Upcycling Project” (Project No. 2017M1A2A2043446) through the National Research Foundation (NRF) funded by the Ministry of Science and ICT, Republic of Korea.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Amembrane area (m2)
Cnumber of gas components
DoOutside diameter of membrane (m)
gpermeation rate (m3/s)
hNumber of bore side finite element per finite element of shell side
LLength of membrane (m)
nratio of number of finite elements of bore side to shell side
Qpermeance (m3/m2.s.pa)
sNumber of shell side finite element
uflow rate of reside side (m3/s)
vflow rate of permeate side (m3/s)
xconcentration in residue side (volume %)
yconcentration at in permeate side (volume %)
ȳconcentration in bulk stream of permeate side (volume %)
μviscosity (m2.pa/s)
Subscript
ffeed
ii component
jj-th finite element in bore side
kk-th finite element in shell side
mmixed gas

Appendix A

Table A1. RMSPEs calculated using the permeance of pure gases.
Table A1. RMSPEs calculated using the permeance of pure gases.
snCO Recovery (%)CO Conc. (%)CO2 Conc. (%)
213.09 2.24 171.42
513.06 0.70 58.41
1013.06 0.13 16.03
1513.05 0.17 2.26
2013.05 0.27 5.43
3013.05 0.37 12.20
6013.05 0.48 18.89
Table A2. CPU run times and root-mean-square errors of the NEQNF solution method.
Table A2. CPU run times and root-mean-square errors of the NEQNF solution method.
sEquation NumberCPU Time (sec.)ΣF2Error (Root Mean Square)
1110.0163.92 × 10−345.97 × 10−18
2220.0161.51 × 10−322.62 × 10−17
151100.027.77 × 10−358.40 × 10−19
505500.171.57 × 10−385.34 × 10−21
10011001.621.59 × 10−351.20 × 10−19
200220019.483.84 × 10−351.32 × 10−19

References

  1. Global Energy Statistics 2019, CO2 Emission from Fuel Combustion. Available online: https://yearbook.enerdata.co.kr/co2-fuel-combustion/CO2-emissions-data-from-fuel-combustion (accessed on 3 September 2019).
  2. IPCC, Special Report for Global Warming of 1.5 °C. 2018. Available online: https://www.ipcc.ch/sr15/ (accessed on 28 August 2019).
  3. Monsalve-Bravo, G.M.; Smart, S.; Bhatia, S.K. Simulation of Multicomponent Gas Transport Through Mixed-Matrix Membranes. J. Membr. Sci. 2019, 577, 219–234. [Google Scholar] [CrossRef] [Green Version]
  4. Tao, J.; Wang, J.; Zhu, L.; Chen, X. Integrated Design of Multi-Stage Membrane Separation for Landfill Gas with Uncertain Feed. J. Membr. Sci. 2019, 590, 117260. [Google Scholar] [CrossRef]
  5. Haddadi, B.; Jordan, C.; Miltner, M.; Harasek, M. Membrane Modeling Using CFD: Combined Evaluation of Mass Transfer and Geometrical Influences in 1D and 3D. J. Membr. Sci. 2018, 563, 199–209. [Google Scholar] [CrossRef]
  6. Zarca, G.; Urtiaga, A.; Biegler, L.T.; Ortiz, I. An Optimization Model for Assessment of Membrane-Based Post-Combustion Gas Upcycling Into Hydrogen or Syngas. J. Membr. Sci. 2018, 563, 83–92. [Google Scholar] [CrossRef]
  7. Qadir, S.; Hussain, A.; Ahsan, M. Numerical Analysis of a Gas Separation of CH4/CO2 Using Hollow Fiber Membrane Module. Sigma J. Eng. Nat. Sci. 2018, 36, 501. [Google Scholar]
  8. McNeil, M.V.; Wilfart, F.M.; Haelssig, J.B. Modelling Hollow Fiber Membrane Modules for Anesthesia Gas Separation. Chem. Eng. Sci. 2018, 191, 479–489. [Google Scholar] [CrossRef]
  9. Daeho, K.O. Development of a Dynamic Simulation Model of a Hollow Fiber Membrane Module to Sequester CO2 from Coalbed Methane. J. Membr. Sci. 2018, 546, 258–269. [Google Scholar] [CrossRef]
  10. Lee, S.H.; Binns, M.; Lee, J.H.; Moon, J.-H.; Yeo, J.-G.; Yeo, Y.-K.; Lee, Y.; Kim, J.-K. Membrane Separation Process for CO2 Capture From Mixed Gases Using TR and XTR Hollow Fiber Membranes: Process Modeling and Experiments. J. Membr. Sci. 2017, 541, 224–234. [Google Scholar] [CrossRef]
  11. Sharifian, S.; Harasek, M.; Haddadi, B. Simulation of Membrane Gas Separation Process Using Aspen Plus® V8. 6. Chem. Prod. Process Model. 2016, 11, 67–72. [Google Scholar] [CrossRef]
  12. Binns, M.; Lee, S.; Yeo, Y.; Lee, J.H.; Moon, J.; Yeo, J.; Kim, J. Strategies for the Simulation of Multi-Component Hollow Fibre Multi-Stage Membrane Gas Separation Systems. J. Membr. Sci. 2016, 497, 458–471. [Google Scholar] [CrossRef]
  13. Shoghl, S.N.; Raisi, A.; Aroujalian, A. A Predictive Mass Transport Model for Gas Separation Using Glassy Polymer Membranes. RSC Adv. 2015, 5, 38223–38234. [Google Scholar] [CrossRef]
  14. Szwast, M.; Szwast, Z. A Mathematical Model of Membrane Gas Separation With Energy Transfer by Molecules of Gas Flowing in a Channel to Molecules Penetrating This Channel From the Adjacent Channel. Chem. Process Eng. 2015, 36, 151–169. [Google Scholar] [CrossRef] [Green Version]
  15. Ahsan, M.; Hussain, A. A Comparison of Numerical Methods Used to Solve Cross Flow Model for Multicomponent Membrane Gas Separation. World Appl. Sci. J. 2013, 22, 703. [Google Scholar]
  16. Shamsabadi, A.A.; Kargari, A.; Farshadpour, F.; Laki, S. Mathematical Modeling of CO2/CH4 Separation by Hollow Fiber Membrane Module Using Finite Difference Method. J. Membr. Sep. Technol. 2012, 1, 19. [Google Scholar]
  17. Soni, V.; Abildskov, J.; Jonsson, G.; Gani, R. A General Model for Membrane-Based Separation Processes. Comput. Chem. Eng. 2009, 33, 644–659. [Google Scholar] [CrossRef]
  18. Murad Chowdhury, M.H.; Feng, X.; Douglas, P.; Croiset, E. A New Numerical Approach for a Detailed Multicomponent Gas Separation Membrane Model and AspenPlus Simulation. Chem. Eng. Technol. 2005, 28, 773–782. [Google Scholar] [CrossRef]
  19. Wang, R.; Liu, S.L.; Lin, T.T.; Chung, T.S. Characterization of Hollow Fiber Membranes in a Permeator Using Binary Gas Mixtures. Chem. Eng. Sci. 2002, 57, 967–976. [Google Scholar] [CrossRef]
  20. Lim, S.P.; Tan, X.; Li, K. Gas/Vapour Separation Using Membranes. Chem. Eng. Sci. 2000, 55, 2641–2652. [Google Scholar] [CrossRef]
  21. Tessendorf, S.; Gani, R.P.; Michelsen, M.L. Modeling, Simulation and Optimization of Membrane-Based Gas Separation Systems. Chem. Eng. Sci. 1999, 54, 943–955. [Google Scholar] [CrossRef]
  22. Feng, X.; Ivory, J.; Rajan, V.S.V. Air Separation by Integrally Asymmetric Hollow-Fiber Membranes. AIChE J. 1999, 45, 2142–2152. [Google Scholar] [CrossRef]
  23. Li, K.; Teo, W.K. Use of Permeation and Absorption Methods for CO2 Removal in Hollow Fibre Membrane Modules. Sep. Purif. Technol. 1998, 13, 79–88. [Google Scholar] [CrossRef]
  24. Thundyil, M.J.; Koros, W.J. Mathematical Modeling of Gas Separation Permeators-For Radial Crossflow, Countercurrent, and Concurrent Hollow Fiber Membrane Modules. J. Membr. Sci. 1997, 125, 275–291. [Google Scholar] [CrossRef]
  25. Li, K.; Chen, Y.; Teo, W.K. Removal of CO2 From a Breathing Gas Mixture Using a Hollow-Fibre Permeator With Permeate Purge and Absorption. Gas Sep. Purif. 1995, 9, 93–100. [Google Scholar] [CrossRef]
  26. Giglia, S.; Bikson, B.; Perrin, J.E.; Donatelli, A.A. Mathematical and Experimental Analysis of Gas Separation by Hollow Fiber Membranes. Ind. Eng. Chem. Res. 1991, 30, 1239–1248. [Google Scholar] [CrossRef]
  27. Sidhoum, M.; SenGupta, A.; Sirkar, K.K. Asymmetric Cellulose Acetate Hollow Fibers: Studies in Gas Permeation. AIChE J. 1988, 34, 417–425. [Google Scholar] [CrossRef]
  28. Krovvidi, K.R.; Kovvali, A.S.; Vemury, S.; Khan, A.A. Approximate Solutions for Gas Permeators Separating Binary Mixtures. J. Membr. Sci. 1992, 66, 103–118. [Google Scholar] [CrossRef]
  29. Pan, C.Y. Gas Separation by High-Flux, Asymmetric Hollow-Fiber Membrane. AIChE J. 1986, 32, 2020. [Google Scholar] [CrossRef]
  30. Pan, C.Y. Gas Separation by Permeators with High-Flux Asymmetric Membranes. AIChE J. 1983, 29, 545–552. [Google Scholar] [CrossRef]
  31. Li, K.; Acharya, D.R.; Hughes, R. Mathematical Modelling of Multicomponent Membrane Permeators. J. Membr. Sci. 1990, 52, 205–219. [Google Scholar] [CrossRef]
  32. Mccandless, F.P. Iterative Solution of Multicomponent Permeator Model Equations. J. Membr. Sci. 1990, 48, 115–122. [Google Scholar] [CrossRef]
  33. Boucif, N.; Sengupta, A.; Sirkar, K.K. Hollow Fiber Gas Permeator With Countercurrent or Cocurrent Flow: Series Solutions. Ind. Eng. Chem. Fund. 1986, 25, 217–228. [Google Scholar] [CrossRef]
  34. Boucif, N.; Majumdar, S.; Sirkar, K.K. Series Solutions for a Gas Permeator With Countercurrent and Cocurrent Flow. Ind. Eng. Chem. Fund. 1984, 23, 470–480. [Google Scholar] [CrossRef]
  35. Marriott, J.I.; Sørensen, E.; Bogle, I.D.L. Detailed Mathematical Modelling of Membrane Modules. Comput. Chem. Eng. 2001, 25, 693–700. [Google Scholar] [CrossRef]
  36. Kaldis, S.P.; Kapantaidakis, G.C.; Sakellaropoulos, G.P. Simulation of Multicomponent Gas Separation in a Hollow Fiber Membrane by Orthogonal Collocation—Hydrogen Recovery From Refinery Gases. J. Membr. Sci. 2000, 173, 61–71. [Google Scholar] [CrossRef]
  37. Marriott, J.; Sørensen, E. A General Approach to Modelling Membrane Modules. Chem. Eng. Sci. 2003, 58, 4975–4990. [Google Scholar] [CrossRef]
  38. Visser, T.; Masetto, N.; Wessling, M. Materials Dependence of Mixed Gas Plasticization Behavior in Asymmetric Membranes. J. Membr. Sci. 2007, 306, 16–28. [Google Scholar] [CrossRef]
  39. Ettouney, H.M.; Majeed, O. Transport and Permeation Properties of a Ternary Gas Mixture in a Medium-Size Polysulfone Hollow Fiber Permeator. Sep. Sci. Technol. 1996, 31, 1573–1596. [Google Scholar] [CrossRef]
  40. Chernova, E.; Petukhov, D.; Boytsova, O.; Alentiev, A.; Budd, P.; Yampolskii, Y.; Eliseev, A. Enhanced Gas Separation Factors of Microporous Polymer Constrained in the Channels of Anodic Alumina Membranes. Sci. Rep. 2016, 6, 31183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Kelman, S.; Lin, H.; Sanders, E.S.; Freeman, B.D. CO2/C2H6 Separation Using Solubility Selective Membranes. J. Membr. Sci. 2007, 305, 57–68. [Google Scholar] [CrossRef]
  42. Merkel, T.; Gupta, R.P.; Turk, B.S.; Freeman, B.D. Mixed-Gas Permeation of Syngas Components in Poly(Dimethylsiloxane) and Poly(1-Trimethylsilyl-1-Propyne) at Elevated Temperatures. J. Membr. Sci. 2001, 191, 85–94. [Google Scholar] [CrossRef]
  43. Lin, H.; Zhou, M.; Ly, J.; Vu, J.; Wijmans, J.G.; Merkel, T.C.; Jin, J.; Haldeman, A.; Wagener, E.H.; Rue, D. Membrane-Based Oxygen-Enriched Combustion. Ind. Eng. Chem. Res. 2013, 52, 10820–10834. [Google Scholar] [CrossRef]
  44. Lin, H.; He, Z.; Sun, Z.; Kniep, J.; Ng, A.; Baker, R.W.; Merkel, T.C. CO2-Selective Membranes for Hydrogen Production and CO2 Capture–Part II: Techno-Economic Analysis. J. Membr. Sci. 2015, 493, 794–806. [Google Scholar] [CrossRef] [Green Version]
  45. Chenar, M.P.; Soltanieh, M.; Matsuura, T.; Tabe-Mohammadi, A.; Feng, C. Gas Permeation Properties of Commercial Polyphenylene Oxide and Cardo-Type Polyimide Hollow Fiber Membranes. Sep. Purif. Technol. 2006, 51, 359–366. [Google Scholar] [CrossRef]
  46. Katoh, T.; Tokumura, M.; Yoshikawa, H.; Kawase, Y. Dynamic Simulation of Multicomponent Gas Separation by Hollow-Fiber Membrane Module: Nonideal Mixing Flows in Permeate and Residue Sides Using the Tanks-in-Series Model. Sep. Purif. Technol. 2011, 76, 362–372. [Google Scholar] [CrossRef]
  47. Wilke, C.R. A Viscosity Equation for Gas Mixtures. J. Chem. Phys. 1950, 18, 517–519. [Google Scholar] [CrossRef]
  48. Coker, D.T.; Allen, T.; Freeman, B.D.; Fleming, G.K. Nonisothermal Model for Gas Separation Hollow-Fiber Membranes. AIChE J. 1999, 45, 1451–1468. [Google Scholar] [CrossRef]
  49. Coker, D.T.; Freeman, B.D.; Fleming, G.K. Modeling Multicomponent Gas Separation Using Hollow-Fiber Membrane Contactors. AIChE J. 1998, 44, 1289–1302. [Google Scholar] [CrossRef]
  50. Ahmad, A.L.; Adewole, J.K.; Leo, C.P.; Ismail, S.; Sultan, A.S.; Olatunji, S.O. Prediction of Plasticization Pressure of Polymeric Membranes for CO2 Removal From Natural Gas. J. Membr. Sci. 2015, 480, 39–46. [Google Scholar] [CrossRef]
  51. Ismail, A.F.; Lorna, W. Penetrant-Induced Plasticization Phenomenon in Glassy Polymers for Gas Separation Membrane. Sep. Purif. Technol. 2002, 27, 173–194. [Google Scholar] [CrossRef]
  52. Bos, A.; Pünt, I.G.M.; Wessling, M.; Strathmann, H. CO2-Induced Plasticization Phenomena in Glassy Polymers. J. Membr. Sci. 1999, 155, 67–78. [Google Scholar] [CrossRef]
  53. Wessling, M.; Schoeman, S.; van der Boomgaard, T.H.; Smolders, C.A. Plasticization of Gas Separation Membranes. Gas Sep. Purif. 1991, 5, 222–228. [Google Scholar] [CrossRef] [Green Version]
  54. Sanders, E.S. Penetrant-Induced Plasticization and Gas Permeation in Glassy Polymers. J. Membr. Sci. 1988, 37, 63–80. [Google Scholar] [CrossRef]
  55. Chiou, J.S.; Barlow, J.W.; Paul, D.R. Plasticization of Glassy Polymers by CO2. J. Appl. Polym. Sci. 1985, 30, 2633–2642. [Google Scholar] [CrossRef]
  56. Sada, E.; Kumazawa, H.; Wang, J.-S.; Koizumi, M. Separation of Carbon Dioxide by Asymmetric Hollow Fiber Membrane of Cellulose Triacetate. J. Appl. Polym. Sci. 1992, 45, 2181–2186. [Google Scholar] [CrossRef]
  57. Ahmad, F.; Lau, K.K.; Shariff, A.M.; Fong Yeong, Y. Temperature and Pressure Dependence of Membrane Permeance and Its Effect on Process Economics of Hollow Fiber Gas Separation System. J. Membr. Sci. 2013, 430, 44–55. [Google Scholar] [CrossRef]
  58. Wang, R.; Cao, C.; Chung, T. A Critical Review on Diffusivity and the Characterization of Diffusivity of 6FDA–6FpDA Polyimide Membranes for Gas Separation. J. Membr. Sci. 2002, 198, 259–271. [Google Scholar] [CrossRef]
  59. Lin, W.-H.; Chung, T. Gas Permeability, Diffusivity, Solubility, and Aging Characteristics of 6FDA-durene Polyimide Membranes. J. Membr. Sci. 2001, 186, 183–193. [Google Scholar] [CrossRef]
  60. SenGupta, A.; Sirkar, K.K. Ternary Gas Mixture Separation in Two-Membrane Permeators. AIChE J. 1987, 33, 529–539. [Google Scholar] [CrossRef]
Figure 1. Schematic of a general finite element model for asymmetric hollow fiber membrane. (a) Bore side feed and (b) shell-side feed.
Figure 1. Schematic of a general finite element model for asymmetric hollow fiber membrane. (a) Bore side feed and (b) shell-side feed.
Membranes 12 00097 g001
Figure 2. Schematic diagram of the one-stage membrane separation process for LDG separation (PI: pressure indicator; TI: temperature indicator; GC: gas chromatography; MFC: mass flow controller).
Figure 2. Schematic diagram of the one-stage membrane separation process for LDG separation (PI: pressure indicator; TI: temperature indicator; GC: gas chromatography; MFC: mass flow controller).
Membranes 12 00097 g002
Figure 3. RMSPE for the major performance indicators of the membrane according to s value at n = 1.
Figure 3. RMSPE for the major performance indicators of the membrane according to s value at n = 1.
Membranes 12 00097 g003
Figure 4. RMSPE for the major performance indicators with the variation of n at s = 15.
Figure 4. RMSPE for the major performance indicators with the variation of n at s = 15.
Membranes 12 00097 g004
Figure 5. Comparison of simulated and experimentally measured concentrations of the components in the permeate.
Figure 5. Comparison of simulated and experimentally measured concentrations of the components in the permeate.
Membranes 12 00097 g005
Figure 6. Permeances of the mixed gas components with varying feed pressures (1 GPU = 7.5 × 10−12 m3/m2.s.Pa).
Figure 6. Permeances of the mixed gas components with varying feed pressures (1 GPU = 7.5 × 10−12 m3/m2.s.Pa).
Membranes 12 00097 g006
Figure 7. Simulated CO recoveries and CO2 concentrations in the residue based on the feed flow rate and feed pressure at a permeate pressure 1 bar.
Figure 7. Simulated CO recoveries and CO2 concentrations in the residue based on the feed flow rate and feed pressure at a permeate pressure 1 bar.
Membranes 12 00097 g007
Figure 8. Permeate side-vacuuming effect on the CO recovery and CO2 concentrations in the residue.
Figure 8. Permeate side-vacuuming effect on the CO recovery and CO2 concentrations in the residue.
Membranes 12 00097 g008
Figure 9. CO recovery and CO2 concentration in the residue with the variation of feed flow rate and feed pressure at the permeate side pressure of 0.1 bar.
Figure 9. CO recovery and CO2 concentration in the residue with the variation of feed flow rate and feed pressure at the permeate side pressure of 0.1 bar.
Membranes 12 00097 g009
Figure 10. Comparison of the CO2 concentrations in the residue in the bore side and shell side feed modes.
Figure 10. Comparison of the CO2 concentrations in the residue in the bore side and shell side feed modes.
Membranes 12 00097 g010
Figure 11. Effect of the shell side finite element number on the CO2 concentration in the residue and CO recovery.
Figure 11. Effect of the shell side finite element number on the CO2 concentration in the residue and CO recovery.
Membranes 12 00097 g011
Table 1. Specifications of the PS membrane module.
Table 1. Specifications of the PS membrane module.
Active Length
(m)
ID
(μm)
OD
(μm)
No. of FiberActive Area
(m2)
0.2420036040001.0
Table 2. Gas separation results at a feed flow rate of 10 L/min, feed pressure 8 bara, and permeate pressure 1 bar.
Table 2. Gas separation results at a feed flow rate of 10 L/min, feed pressure 8 bara, and permeate pressure 1 bar.
Flow Rate (L/min)Concentration (%)
COCO2N2H2
Residue6.91774.574.8920.260.28
Permeate3.06540.1947.256.566.00
Table 3. Permeances and selectivities of the LDG mixed and pure gases.
Table 3. Permeances and selectivities of the LDG mixed and pure gases.
Number of Finite ElementPermeance (×10−10 m3/m2·s·Pa)Selectivity
ShellBorenCOCO2N2H2CO2/COCO2/N2CO/N2
101010.3964 6.1372 0.2391 15.0102 15.66 25.96 1.66
151510.3978 5.8752 0.2402 13.0289 14.77 24.46 1.66
202010.3985 5.7216 0.2408 12.2034 14.36 23.76 1.65
303010.39925.57570.241411.466013.9723.091.65
153020.3989 5.6345 0.2412 11.8931 14.13 23.36 1.65
154530.3993 5.5574 0.2416 11.5443 14.30 23.00 1.65
156040.3994 5.5194 0.2418 11.3750 13.82 22.83 1.65
Pure gas0.37426.15220.231012.147716.4326.231.62
Table 4. Experimental results at feed rates of 5, 20, and 30 L/min.
Table 4. Experimental results at feed rates of 5, 20, and 30 L/min.
Feed Flow Rate (L/min)StreamOutlet Flow Rate (L/min)Concentration (%)
COCO2N2H2
5.0 Residue2.50276.470 0.921 22.580 0.029
Permeate2.34152.454 34.112 8.929 4.506
20.0 Residue16.65870.923 9.964 18.321 0.792
Permeate3.70433.178 54.390 5.197 7.234
30.0 Residue25.83769.102 12.123 17.680 1.096
Permeate3.83130.350 56.962 4.726 7.962
Table 5. Simulation errors with respect to the experimental data at s = 15 and n = 1 for the major performance measures of the membrane.
Table 5. Simulation errors with respect to the experimental data at s = 15 and n = 1 for the major performance measures of the membrane.
Feed Flow Rate (L/min) CO Recovery (%)CO Conc. in Residue (%)CO2 Conc. in Residue (%)
Exper.Simul.Err. (%)Exper.Simul.Err. (%)Exper.Simul.Err. (%)
559.7860.98+2.0076.4776.23−0.310.920.94+2.17
1080.5780.75+0.2174.5774.53−0.054.895.07+3.57
2092.2290.56−1.8070.9270.97+0.079.969.89−0.70
3092.9993.77+0.8369.169.11+0.0112.1212.15+0.24
Root mean square
percent error (%)
1.42 0.16 2.12
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jeon, J.-Y.; Park, B.-R.; Kim, J.-H. Numerical Simulation and Optimization of 4-Component LDG Separation in the Steelmaking Industry Using Polysulfone Hollow Fiber Membranes. Membranes 2022, 12, 97. https://doi.org/10.3390/membranes12010097

AMA Style

Jeon J-Y, Park B-R, Kim J-H. Numerical Simulation and Optimization of 4-Component LDG Separation in the Steelmaking Industry Using Polysulfone Hollow Fiber Membranes. Membranes. 2022; 12(1):97. https://doi.org/10.3390/membranes12010097

Chicago/Turabian Style

Jeon, Jong-Yeol, Bo-Ryoung Park, and Jeong-Hoon Kim. 2022. "Numerical Simulation and Optimization of 4-Component LDG Separation in the Steelmaking Industry Using Polysulfone Hollow Fiber Membranes" Membranes 12, no. 1: 97. https://doi.org/10.3390/membranes12010097

APA Style

Jeon, J. -Y., Park, B. -R., & Kim, J. -H. (2022). Numerical Simulation and Optimization of 4-Component LDG Separation in the Steelmaking Industry Using Polysulfone Hollow Fiber Membranes. Membranes, 12(1), 97. https://doi.org/10.3390/membranes12010097

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