Next Article in Journal
Trimetallic Nanoalloy of NiFeCo Embedded in Phosphidated Nitrogen Doped Carbon Catalyst for Efficient Electro-Oxidation of Kraft Lignin
Previous Article in Journal
New Hybrid Nanocomposites with Catalytic Properties Obtained by In Situ Preparation of Gold Nanoparticles on Poly (Ionic Liquid)/Poly (4-Vinylpyridine) Nanofibers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Simulation of the Impact of Feed Gas Perturbation on CO2 Removal in a Polymeric Hollow Fiber Membrane

Department of Chemical & Petroleum Engineering, United Arab Emirates University, Al-Ain City P.O. Box 15551, United Arab Emirates
Polymers 2022, 14(18), 3783; https://doi.org/10.3390/polym14183783
Submission received: 15 August 2022 / Revised: 31 August 2022 / Accepted: 6 September 2022 / Published: 9 September 2022
(This article belongs to the Section Polymer Physics and Theory)

Abstract

:
A membrane contactor is a device that attains the transfer of gas/liquid or liquid/liquid mass without dispersion of one phase within another. Membrane contactor modules generally provide 30 times more surface area than can be achieved in traditional gas absorption towers and 500 times what can be obtained in liquid/liquid extraction columns. By contrast, membrane contactor design has limitations, as the presence of the membrane adds additional resistance to mass transfer compared with conventional solvent absorption systems. Increasing mass transfer in the gas and solvent phase boundary layers is necessary to reduce additional resistance. This study aims to increase the mass transfer in the gas phase layer without interfering with membrane structure by oscillating the velocity of the feed gas. Therefore, an unsteady state mathematical model was improved to consider feed gas oscillation. The model equation was solved using Comsol Multiphysics version 6.0. The simulation results reveal that the maximum CO2 removal rate was about 30% without oscillation, and at an oscillation frequency of 0.05 Hz, the CO2 percent removal was almost doubled.

1. Introduction

Power plants established on fossil fuels are the main source of CO2 emissions and hence, a primary environmental and climate change concern [1]. Reducing CO2 emissions from power plants through carbon capture is a reasonable solution. Various CO2 capture techniques have been implemented, including cryogenic separation, solvent absorption, and, lately, membrane solvent-gas contactors [2,3,4]. Membrane solvent-gas technology has the advantages of gas separation and compact sizes with resealable CO2 absorption [5,6]. By contrast, membrane separation technologies introduce mass transfer resistances in the liquid phase, gas phase, and membrane layer [7]. Membrane contactors are also used for purposes other than absorption. Membrane contactors have been used to strip dissolved gases other than CO2, such as methane from wastewater treatment anaerobic effluents, natural gas, and flue gas [8,9]. The treatment method has been used in converting chemical oxygen to biogas. As methane is the primary source of biogas, the membrane separation process is utilized to liberate dissolved methane from anaerobic effluents [10]. Potassium glycinate solvent is an effective solvent as an aqueous solution used to absorb CO2 inside a solvent-gas hollow fiber membrane contactor [11,12]. Feed gas oscillating flow was induced experimentally, and the mass transfer correlations were improved in terms of Sherwood number. The study revealed that inducing mixing with the gas phase enhanced the absorption rate in the membrane contactor [13]. Thin polymeric hollow fiber films are another type of practical hollow fiber membrane for CO2 removal [14,15].
Computational fluid dynamics (CFD) and the design of experiments were performed on a membrane distillation. Simulation results revealed that the membrane module’s length strongly impacts the flux [16]. Mathematical modeling was performed for membrane contactors for CO2 absorption in different solvents. Various solvents were used in the CO2 absorption process, such as 4-diethylamino-e-butanol (DEAB), diethanolamine (DEA), and methyl diethanolamine (MDEA) in hollow fiber modules. DEAB absorption rate was competitive with MEA aqueous solutions [17]. The gas velocity in membrane gas absorption plays a significant role in improving separation performance, and the addition of baffles improves SO2 absorption performance at the critical velocity [18].
The ammonia-based CO2 capturing process in a hollow fiber membrane contactor shows better mass transfer performance in the dense membrane phase than conventional gas/liquid absorption processes [19]. The membrane outer diameter has a substantial impact on membrane separation performance. Membrane modules with fibers finer than 1.2 mm in diameter increased mass transfer resistance and decreased the membrane contactor specific area per unit volume. A membrane with a 0.9 mm outer diameter or less was more suitable [20]. A significant parametric study for different physical solvents revealed the importance of physical solvents in absorption and regeneration processes compared with the conventional packed column technology [21]. Developing gas-liquid membrane conductors reduces capital cost and energy consumption of conventional CO2 absorbers and separation columns [22,23,24]. The membrane material, type of solvent, configuration of solvent and gas placed in the lumen or shell side, and the operating temperature plays significant roles in conquering the technical challenges of membrane wetting and may result in a decrease in the mass transfer coefficient [25,26]. Nanoparticles such as SiO2 and Al2O3, along with MEA aqueous solution, help reduce the energy consumption of CO2 stripping. The highest mass transfer resistance is in the liquid phase, and increased liquid velocity increases the desorption flux [27].
Based on the above findings, there are difficulties in increasing the mass transfer coefficient based on the membrane structure, composition, and configuration. Accordingly, the current work focused on modeling and simulating the gaseous solvent membrane contactor based on the oscillating feed gas flow rate to increase CO2 removal without affecting membrane structure and composition.

2. Mathematical Model

We built a two-dimensional unsteady state mathematical model to describe carbon dioxide absorption from methane in potassium glycinate (PG) aqueous solution using a hollow fiber membrane absorption process under the influence of an oscillating gas flow. The model was solved in two scenarios; the first method is to validate the model predictions in which the gas is fed into the tube lumen side (Figure 1a). In the second scenario, the gas is fed to the shell side (Figure 1b). The model equations listed below describe the most common scenario taking place in a membrane contactor. The potassium glycinate solvent is fed to the lumen at z = L, and the gas stream enters the shell side currently at z = 0 (Figure 1b). Specific assumptions were considered, such as isothermal operation, laminar flow, and Henry’s law, which is applicable at the liquid-gas interface and non-wetting mode (gas-filled membrane pores).

2.1. Hollow Fiber Lumen (Solvent Flow)

Considering the case where gas is fed into the shell side and liquid solvent to the lumen side, the unsteady steady-state gas transport in the lumen solvent side i is CO 2 ,   CH 4 .
C i , t t = D i , t [ r ( 1 r C i , t r ) + 2 C i , t z 2 ] + r i , t + v z , t C i , t z
where c i , t is the molar concentration of CO2 in the liquid solvent flowing in the membrane tube side and D i is the diffusion coefficient of CO2 in the tube side.
The aqueous solution velocity ( v z , t ):
v z , t = 2 Q t n π r 1 2 ( 1 ( r r 1 ) 2 )
where Q t is the liquid solvent feed rate, n is the number of fibers, and r1 is the inner radius of the tubes.
The boundary conditions:
at z = 0 , C C O 2 , t = 0 (the CO2 concentration in the fresh solvent is zero);
at z = L , C C O 2 , t z = 0 (the convective flux is zero at the exit of the tube);
at r = 0 , C C O 2 , t   r = 0 (axial symmetry assumption);
at r = r 1 , C C O 2 , t = m   C C O 2 , m (at the gas-liquid interface, the solubility of CO2 in PG).
The forward reaction rate is expressed as follows [28].
r C O 2 = 2.42 × 10 16 exp ( 8544 T ) exp ( 0.44 C P g ) C P g C C O 2
where T(K) is the liquid temperature, C p g , are the concentrations of PG, and C C O 2 is the concentration of CO2.

2.2. Membrane Skin

In this section, the model equations describe the unsteady steady-state gas diffusion across the membrane walls where the convective flux is neglected, and only diffusion is considered:
C i , m t = D i , m [ r ( 1 r C i , m r ) + 2 C i , m z 2 ]
The boundary conditions of the membrane skin:
at z = 0 , C i , m z = 0 (at the inlet of the module, the membrane walls are solid, and diffusion flux is assumed insignificant);
at z = L , C i , m z = 0 (at the exit of the membrane module, the membrane walls are solid, and diffusion flux is ignored);
at r = r 1 , D i , m C i , m r = D i , t C i , t r (at the liquid-liquid interface, the fluxes are assumed equal);
at r = r 2 , C i , m = C i , s (at the membrane-gas interface, the membrane resistance is neglected, and the concentration of the gas in the shell side is assumed to be equal to that in the membrane walls).

2.3. The Shell of the Module (Gas Stream)

The unsteady steady-state gas concentration in the shell side:
C i , s t = D i , s [ r ( 1 r C i , s r ) + 2 C i , s z 2 ] + v z , s ( C i , s z )
Here, C i , s is the molar concentration of the CO2 in the gas on the module shell side and D i s is the diffusion coefficient of CO2 in the gas phase flowing on the shell side.
The perturbation of the gas velocity on the shell side [29] is modified to consider the sine wave frequency of the gas velocity on the shell side ( v s m ), whereas v s m is the main gas velocity on the shell side before reaching the oscillating device.
v z , s = v z m   sin ( ω t ) { 1 ( r 2 r 3 ) 2 } { ( r r 3 ) 2 ( r 2 r 3 ) 2 2 l n ( r r 2 )   3 + ( r 2 r 3 ) 4 4 ( r 2 r 3 ) 2 + 4 ln ( r 2 r 3 ) }
where ω is the frequency. The appropriate boundary conditions:
z = L , C i , s = C i , 0 (the gas flows in the module counter currently; accordingly, the concentration of CO2 in the gas equals the initial concentration);
z = 0 , 2 C i , s a z 2 = 0 (the gas leaves from the entrances of the module and hence, convective flux of the exit gas stream is negligible);
r = r 2 , D i , s   C i , s a r = D i , m s   C i , m a r (at the membrane skin-shell gas interface, the diffusive flux is equal);
r = r 3 , C i , s a r = 0 (at the outer module radius, the diffusion flux is neglected).
The radius of the free surface ( r 3 ), is expressed as follows:
r 3 = r 2 ( 1 1 φ ) 0.5
The module void fraction ( φ ):
φ = R 2 n   r 2 2 R 2
R , and r 2 n are the inner radius of the module and fiber outer radius, respectively. Qt is the solvent circulation volumetric rate, and n is the number of fibers. Table 1 lists the characteristics of the hollow fiber membrane module. The experimental conditions were at atmospheric pressure and room temperature.

3. Mass Transfer Coefficient, K G

The total mass transfer coefficient, K G ( k m o l m 2 . s . k P a ) , was estimated based on the two-phase theory as the ratio of the CO2 absorption flux ( J C O 2 ) to the solute concentration gradient between gas and liquid [30]. CO2 absorption flux ( J C O 2 ) was calculated based on the following expression [31]:
J C O 2 ( m o l / m 2 . s ) = Q g , i n   y C O 2 , i n Q g , o u t   y C O 2 , o u t A
where Q g is the gas molar flow rate (mol/s) (assuming inlet and exit flow rate are the same; neglecting the effect of the amount of CO2 being absorbed on the exit gas flow rate), y C O 2 , i n and y C O 2 , o u t are the inlet and exit mole fractions of CO2, respectively, and A is the gas-–liquid interface contact area based on the internal diameter of the hollow fiber. The total mass transfer coefficient is calculated as follows [32]:
K G   ( m / s ) = J C O 2 C g , L m
where K G symbolizes the overall mass transfer coefficient in the gas phase and C g , L m is the average log mean concentration of CO2 in the bulk gas phase of the outlet and inlet concentration [33].
C g , l m = C C O 2 , i n C C O 2 , o u t ln ( C C O 2 , i n / C C O 2 , o u t )
The exit gas concentration ( C C O 2 , o u t ) is determined using the boundary line integration built using the software Comsol, version 6.0 [34].

4. Results and Discussion

4.1. Model Validation

The model predictions were validated using experimental data available from the literature [35]. The model validation was performed by comparing the overall mass transfer coefficient between the experimental data and simulation predictions. As per the experimental data used in the model validation, the gas stream was fed to the lumen side and the solvent to the shell side (Figure 1a). Figure 2 shows the effect of the feed flow rate represented by the Reynolds number on the overall mass transfer coefficient (KG). The results revealed that the mass transfer coefficient increased with an increased Reynolds number, a common trend of the membrane contactors that is attributed to the increase in gas velocity being associated with a decrease in the boundary layer thickness. This replicates the presence of supplementary mixing in the gas phase boundary layer, as the oscillating flow disrupts the thickness of the boundary layer and the flow regime around this layer [36]. The model predictions revealed that increasing the gas velocity had significant effects on the total gas mass transfer coefficients owing to the formation of a thin film for mass transfer layer at higher speeds, which resulted in increases in the mass transfer coefficient of the gas phase layer in spite of decreasing the residence time of the gas inside the membrane module at high velocities. Increase in the velocity has insignificant impact on the overall mass transfer coefficient [37].
A comparison of the gas feed rate between the tube side and shell side revealed that introducing the gas to the shell side of the membrane module enhanced the overall mass transfer coefficient compared with gas being introduced into the tube side, attributed to increases in surface contact area and gas residence time in the membrane module (Figure 3). Accordingly, the rest of the results are based on gas flow fed into the shell side and liquid solvent into the tube side (Figure 1b). This configuration matches most of the experimental work performed in the literature, where liquids are introduced to the lumen side of the membrane module for better contact area and residence time [8,33]. The increase in the mass transfer coefficient when gas is fed into the shell side is attributed to better mixing and the increase in residence time when the gas flow rate decreases. Our simulation is based on the data available in Table 1, where the gas is introduced into the shell side.

4.2. Gas Feed Perturbation

The gas feed velocity coefficient at different employed sinusoidal frequencies was investigated. As absorption occurred in a gas-liquid membrane system that was stabilized very rapidly, a time range of 60 s was investigated [35]. The frequency range was maintained while increasing the speed gradually and reaching the maximum point in 60 s. For example, at a frequency of 0.03, the gas speed gradually increased until it reached its maximum value. Frequencies (other than 0.05) cause the velocity to reach the maximum speed more than once in an oscillatory manner within the studied range. Under sinusoidal operating conditions, the velocity of the feed gas at different frequencies is described in Figure 4. The figure shows the effect of the feed gas sinusoidal frequency on CO2 concentration along the shell side of the membrane unit. The lowest exit CO2 concentration profile was at a frequency of 0.05 HZ. A smaller CO2 concentration in the exit stream indicates a higher percentage removal of CO2 from the gas stream mixture [35].
Figure 5 depicts the effect of sinewave frequencies on the percent removal of CO2 from the feed gas stream. The figure shows that the CO2 removal percentage was low (0.03 percent), and the removal rate reached the maximum at a frequency of 0.05. Further increase in the frequencies causes a decline in the percentage of CO2 removal, and that is attributed to the sinusoidal number of cycles; at a frequency of 0.2, the velocity coefficient increases within the investigated interval of 60 s around four times. The results indicated that an oscillation frequency of 0.05 achieves the highest CO2 removal rate associated with a reduction in the boundary layer thickness. This reflected the presence of supplementary mixing in the gas phase boundary layer as the oscillating flow interrupted the boundary layer thickness and flow regime around this layer [37].
The increase in membrane efficiency caused by a frequency of 0.05 shows the best removal rate as the velocity reaches a maximum and then decreases to give more room for fresher gas and increase the driving force. Further increase in frequencies leads to a greater number of cycles within the investigated period, decreasing the residence time and not giving enough time for separation [35]. The reasonable frequency of 0.05 doubled the percentage removal because the pressure wave generated by the oscillation of the feed gas produced localized increases in the partial pressure driving force for mass transfer across the membrane. The results were consistent with previous experimental observations [37].
Figure 6 shows the model predictions for the effect of a feed gas flow rate at a constant sinusoidal frequency of 0.05 on CO2 concentrations and removal flux along the dimensionless membrane length. The figure shows that, along the membrane dimensionless length, the higher the gas flow rate, the lower the CO2 removal rate, and the higher the CO2 removal flux. This is attributed to the fact that an increase in the gas flow rate usually reduces the thickness of the gas boundary layer, which is supposed to enhance the mass transfer rate and increase the CO2 removal ratio. At the same time, an increase in the feed gas flow rate reduces the residence time, which reduces the rate of mass transfer, and thus, reduces the rate of CO2 absorption along with the membrane gas-phase compartment (membrane unit envelope side). Similar results were reported in previous studies [38,39,40]. This indicates that the residence time strongly impacts the CO2 removal rate. The increase in the CO2 removal flux with the gas flow rate is credited to the increased amount of CO2 being absorbed with the gas flow rate. Accordingly, residence time strongly affects CO2 removal compared with gas mass transfer intensification.
Figure 7 indicates the influence of the gas feed rate on the total mass transfer coefficients under the presence (0.05 sinusoidal frequency) and absence of feed oscillation. The total mass transfer coefficient increased with the gas feed rate due to the amount of carbon dioxide absorbed, and the thickness of the gas boundary layer decreased with the increase in gas velocity. The total mass transfer coefficients increase with the increase in gas velocity mentioned earlier. The results are within the range of those obtained experimentally. In contrast, the effect of gas supply rate on the total mass transfer coefficient in the absence of gas supply disturbance in the sine waveform is negligible. There is a significant increase in KG compared with the non-oscillating mode of operation. The oscillating flow disturbed the boundary layer thickness and flow regime around the layer, causing further mixing in the gas phase boundary layer, and hence, enhancing the CO2 diffusion rate to the absorbing solvent.
Figure 8 compares the CO2 removal flux between the oscillatory and non-oscillatory gas feed velocity represented by the gas feed rate. The figure shows that the removal percentage increased slightly in the non-oscillating operating mode at a gas feed rate up to a flow rate of 0.25 L/min. Beyond this, the change in the removal ratio with the gas feed rate was insignificant. In contrast, under oscillating mode, a considerable increase in the mass transfer coefficient was observed, indicating the success of the feed gas oscillation mode of operation.

5. Conclusions

A two-dimensional (2D) transient mathematical model, assumed to operate under non-wetted mode conditions, was upgraded by considering feed gas oscillation conditions. The feed gas enters the membrane module at different sinusoidal waves with different frequencies. The model was utilized to study the absorption of CO2 from natural gas into a PG aqueous solution within a hollow fiber membrane unit. In addition to the perturbation of the inlet gas velocity, the model considers the radial and axial diffusion within the hollow fiber lumen, membrane walls, and module shell sides. The predicted results revealed that the feed gas oscillation frequency and flow rate of both gas and liquid streams influenced CO2 removal efficiency. We found significant improvement in the sinusoidal frequency of the feed gas stream. The optimum frequency that doubled the percentage removal of CO2 was around 0.05. The percentage removal increased from almost 30% at no oscillation to around 70% at a frequency of 0.05. The model predictions were in good agreement with experimental data.

Funding

This research was supported by the research office, United Arab Emirates University, UAE, Grant number 31N374.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article. Additional explanations or data are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nieminen, H.; Järvinen, L.; Ruuskanen, V.; Laari, A.; Koiranen, T.; Ahola, J. Insights into a Membrane Contactor Based Demonstration Unit for CO2 Capture. Sep. Purif. Technol. 2020, 231, 115951. [Google Scholar] [CrossRef]
  2. Ho, C.-D.; Chang, H.; Lin, G.-H.; Chew, T.L. Enhancing Absorption Performance of CO2 by Amine Solution through the Spiral Wired Channel in Concentric Circular Membrane Contactors. Membranes 2022, 12, 4. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, H.; Xue, K.; Cheng, C.; Gao, D.; Chen, H. Study on the Performance of CO2 Capture from Flue Gas with Ceramic Membrane Contactor. Sep. Purif. Technol. 2021, 265, 118521. [Google Scholar] [CrossRef]
  4. Ruan, H.; Wu, S.; Chen, X.; Zou, J.; Liao, J.; Cui, H.; Dong, Y.; Qiu, Y.; Shen, J. Capturing CO2 with NaOH Solution from Reject Brine via an Integrated Technology Based on Bipolar Membrane Electrodialysis and Hollow Fiber Membrane Contactor. Chem. Eng. J. 2022, 450, 138095. [Google Scholar] [CrossRef]
  5. Soroodan Miandoab, E.; Kentish, S.E.; Scholes, C.A. Non-Ideal Modelling of Polymeric Hollow-Fibre Membrane Systems: Pre-Combustion CO2 Capture Case Study. J. Memb. Sci. 2020, 595, 117470. [Google Scholar] [CrossRef]
  6. Mohammadi Saadat, M.; Norouzbahari, S.; Esmaeili, M. CO2/N2 Separation by Glycerol Aqueous Solution in a Hollow Fiber Membrane Contactor Module: CFD Simulation and Experimental Validation. Fuel 2022, 323, 124370. [Google Scholar] [CrossRef]
  7. Chen, G.; Chen, G.; Cao, F.; Zhang, R.; Gao, H.; Liang, Z. Mass Transfer Performance and Correlation for CO2 Absorption into Aqueous 3-Diethylaminopropylamine Solution in a Hollow Fiber Membrane Contactor. Chem. Eng. Process. Process Intensif. 2020, 152, 107932. [Google Scholar] [CrossRef]
  8. Ghasem, N. Modeling and Simulation of the Simultaneous Absorption/Stripping of CO2 with Potassium Glycinate Solution in Membrane Contactor. Membranes 2020, 10, 72. [Google Scholar] [CrossRef]
  9. Ghasem, N. CFD Simulation of CO2 Absorption by Water-Based TiO2 Nanoparticles in a High Pressure Stirred Vessel. Sci. Rep. 2021, 11, 1984. [Google Scholar] [CrossRef]
  10. Velasco, P.; Jegatheesan, V.; Othman, M. Modeling of Hollow Fiber Membrane Contactors (HFMCs) for the Recovery of Dissolved Methane from Anaerobic Effluents. Sep. Purif. Technol. 2022, 286, 120488. [Google Scholar] [CrossRef]
  11. Ghasem, N. Modeling and Simulation of CO2 Absorption Enhancement in Hollow-Fiber Membrane Contactors Using CNT–Water-Based Nanofluids. J. Membr. Sci. Res. 2019, 5, 295–302. [Google Scholar] [CrossRef]
  12. Yang, D.; Zhang, J.; Xie, M.; Tong, X.; Jiang, T.; Yu, W.; Qu, D. An Integrated Hollow Fiber Membrane Contactor and Chemical Precipitation to Recover N, P and K from Human Urine Wastewater. J. Environ. Chem. Eng. 2022, 10, 107844. [Google Scholar] [CrossRef]
  13. Hosseini, E.; Soroodan Miandoab, E.; Stevens, G.W.; Scholes, C.A. Absorption of CO2 from Flue Gas under Oscillating Gas Flow Conditions in Gas-Solvent Hollow Fibre Membrane Contactors. Sep. Purif. Technol. 2020, 249, 117151. [Google Scholar] [CrossRef]
  14. Liang, C.Z.; Liu, J.T.; Lai, J.-Y.; Chung, T.-S. High-Performance Multiple-Layer PIM Composite Hollow Fiber Membranes for Gas Separation. J. Memb. Sci. 2018, 563, 93–106. [Google Scholar] [CrossRef]
  15. Liang, C.Z.; Yong, W.F.; Chung, T.-S. High-Performance Composite Hollow Fiber Membrane for Flue Gas and Air Separations. J. Memb. Sci. 2017, 541, 367–377. [Google Scholar] [CrossRef]
  16. Choi, J.; Cho, H.; Choi, Y.; Lee, S. Combination of Computational Fluid Dynamics and Design of Experiments to Optimize Modules for Direct Contact Membrane Distillation. Desalination 2022, 524, 115460. [Google Scholar] [CrossRef]
  17. Saidi, M. Mathematical Modeling of CO2 Absorption into Novel Reactive DEAB Solution in Hollow Fiber Membrane Contactors; Kinetic and Mass Transfer Investigation. J. Memb. Sci. 2017, 524, 186–196. [Google Scholar] [CrossRef]
  18. Kong, X.; Xu, P.; Fu, K.; Gong, D.; Chen, X.; Qiu, M.; Fan, Y. Critical Gas Velocity of Hydrophobic Ceramic Membrane Contactors for SO2 Absorption. Chem. Eng. Sci. 2021, 231, 116327. [Google Scholar] [CrossRef]
  19. Mansourizadeh, A.; Rezaei, I.; Lau, W.J.; Seah, M.Q.; Ismail, A.F. A Review on Recent Progress in Environmental Applications of Membrane Contactor Technology. J. Environ. Chem. Eng. 2022, 10, 107631. [Google Scholar] [CrossRef]
  20. Li, M.; Zhu, Z.; Zhou, M.; Jie, X.; Wang, L.; Kang, G.; Cao, Y. Removal of CO2 from Biogas by Membrane Contactor Using PTFE Hollow Fibers with Smaller Diameter. J. Memb. Sci. 2021, 627, 119232. [Google Scholar] [CrossRef]
  21. Cesari, L.; Castel, C.; Favre, E. Membrane Contactors for Intensified Gas-Liquid Absorption Processes with Physical Solvents: A Critical Parametric Study. J. Memb. Sci. 2021, 635, 119377. [Google Scholar] [CrossRef]
  22. Houlker, S.; Davey, C.J.; Allemand, A.; Brookes, A.; Moore, A.; Vale, P.; Pidou, M.; McAdam, E.J. Reconciliation of Gas to Liquid Mass Transfer in Parallel and Transverse Flow (Cross-Flow) Hollow Fiber Membrane Contactors (HFMC) for CO2 Absorption. Sep. Sci. Technol. 2021, 56, 129–140. [Google Scholar] [CrossRef]
  23. Scholes, C.A.; Kentish, S.E.; Qader, A. Membrane Gas-Solvent Contactor Pilot Plant Trials for Post-Combustion CO2 Capture. Sep. Purif. Technol. 2020, 237, 116470. [Google Scholar] [CrossRef]
  24. Lian, S.; Song, C.; Liu, Q.; Duan, E.; Ren, H.; Kitamura, Y. Recent Advances in Ionic Liquids-Based Hybrid Processes for CO2 Capture and Utilization. J. Environ. Sci. 2021, 99, 281–295. [Google Scholar] [CrossRef] [PubMed]
  25. Kim, S.; Scholes, C.A.; Heath, D.E.; Kentish, S.E. Gas-Liquid Membrane Contactors for Carbon Dioxide Separation: A Review. Chem. Eng. J. 2021, 411, 128468. [Google Scholar] [CrossRef]
  26. Pang, H.; Chen, Z.; Gong, H.; Du, M. Fabrication of a Super Hydrophobic Polyvinylidene Fluoride–Hexadecyltrimethoxysilane Hybrid Membrane for Carbon Dioxide Absorption in a Membrane Contactor. J. Memb. Sci. 2020, 595, 117536. [Google Scholar] [CrossRef]
  27. Mohammed, H.N.; Ahmed, S.M.R.; Al-Naseri, H.; Al-Dahhan, M. Enhancement of CO2 Desorption from MEA-Based Nanofluids in Membrane Contactor: Simulation Study. Chem. Eng. Process. Process Intensif. 2021, 168, 108582. [Google Scholar] [CrossRef]
  28. Eslami, S.; Mousavi, S.M.; Danesh, S.; Banazadeh, H. Modeling and Simulation of CO2 Removal from Power Plant Flue Gas by PG Solution in a Hollow Fiber Membrane Contactor. Adv. Eng. Softw. 2011, 42, 612–620. [Google Scholar] [CrossRef]
  29. Shiravi, A.; Naderifar, A.; Raisi, A. Modeling of Carbon Dioxide Separation by Aqueous Ammonia Solution in a Hollow-Fiber Membrane Contactor. Chem. Eng. Technol. 2022, 45, 1637–1647. [Google Scholar] [CrossRef]
  30. Ghasem, N. Mass Transfer Modeling in Nanofluids: Theoretical Basics and Model Development. In Nanofluids and Mass Transfer; Elsevier: Amsterdam, The Netherlands, 2022; pp. 247–271. ISBN 9780128239964. [Google Scholar]
  31. Velasco, P.; Jegatheesan, V.; Othman, M. Effect of Long-Term Operations on the Performance of Hollow Fiber Membrane Contactor (HFMC) in Recovering Dissolved Methane from Anaerobic Effluent. Sci. Total Environ. 2022, 841, 488–496. [Google Scholar] [CrossRef]
  32. Pahnavar, N.; Keramat, F.; Azari, A.; Osfouri, S. Computational Fluid Dynamics Analysis of CO2 Absorption Intensification in an Hollow Fiber Membrane Contactor Using SiO2 and Carbon Nanotubes Nanofluids. Environ. Prog. Sustain. Energy 2022, 41, e13777. [Google Scholar] [CrossRef]
  33. Sohaib, Q.; Muhammad, A.; Younas, M.; Rezakazemi, M. Modeling Pre-Combustion CO2 Capture with Tubular Membrane Contactor Using Ionic Liquids at Elevated Temperatures. Sep. Purif. Technol. 2020, 241, 116677. [Google Scholar] [CrossRef]
  34. COMSOL. Multiphysics®, version 6.0; COMSOL AB: Stockholm, Sweden, 2022. Available online: www.comsol.com(accessed on 7 July 2022).
  35. Ghobadi, J.; Ramirez, D.; Khoramfar, S.; Jerman, R.; Crane, M.; Hobbs, K. Simultaneous Absorption of Carbon Dioxide and Nitrogen Dioxide from Simulated Flue Gas Stream Using Gas-Liquid Membrane Contacting System. Int. J. Greenh. Gas Control 2018, 77, 37–45. [Google Scholar] [CrossRef]
  36. Hosseini, E.; Soroodan Miandoab, E.; Stevens, G.; Scholes, C.A. Vibration-Induced Enhanced Mass Transfer within Membrane Contactors for Efficient CO2 Capture. Sep. Purif. Technol. 2022, 295, 121251. [Google Scholar] [CrossRef]
  37. Soroodan Miandoab, E.; Scholes, C.A. A Rigorous Membrane Gas-Solvent Contactor Model for Flowsheet Simulation of the Carbon Capture Process. Ind. Eng. Chem. Res. 2022, 61, 9381–9393. [Google Scholar] [CrossRef]
  38. Lin, Z.; Liu, Y.; Zhang, Z.; Yao, J. Preparation and Characterization of OH/SiO2-TiO2/PES Composite Hollow Fiber Membrane Using Gas-Liquid Membrane Contactor for CO2/CH4 Separation. Sep. Purif. Technol. 2022, 288, 120551. [Google Scholar] [CrossRef]
  39. Cao, Y.; Nakhjiri, A.T.; Sarkar, S.M.; Ghadiri, M. Time-Dependent Numerical Investigation of 3-Hydroxypropionic Acid Extraction Using a Microporous Membrane Contactor. Eur. Phys. J. Plus 2022, 137, 494. [Google Scholar] [CrossRef]
  40. Wu, J.; Pan, Z.; Shang, L.; Sun, X.; Sun, C.; Sun, X. Characteristics of CO2 Absorption by N,N-Dimethylethanolamine (DMEA) in Hollow Fiber Membrane Contactor. Huagong Jinzhan Chem. Ind. Eng. Prog. 2022, 41, 2132–2139. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram and the membrane contactor model domain: (a) gas fed in the lumen side, (b) gas fed in the shell side.
Figure 1. Schematic diagram and the membrane contactor model domain: (a) gas fed in the lumen side, (b) gas fed in the shell side.
Polymers 14 03783 g001
Figure 2. Model validation with the experimental data [35]. The gas phase is fed into the lumen side and the solvent into the shell side (Figure 1a).
Figure 2. Model validation with the experimental data [35]. The gas phase is fed into the lumen side and the solvent into the shell side (Figure 1a).
Polymers 14 03783 g002
Figure 3. Comparison between the oscillatory gas flows in the shell and tube sides at operating conditions (temperature of 25 °C, atmosphere pressure, liquid flow rate of 10 mL/min).
Figure 3. Comparison between the oscillatory gas flows in the shell and tube sides at operating conditions (temperature of 25 °C, atmosphere pressure, liquid flow rate of 10 mL/min).
Polymers 14 03783 g003
Figure 4. Effect of gas feed oscillation frequency on the CO2 concentration profile versus membrane dimensionless length of the gas flow on the shell side. Gas feed rate 500 mL/min, liquid feed rate 10 mL/min.
Figure 4. Effect of gas feed oscillation frequency on the CO2 concentration profile versus membrane dimensionless length of the gas flow on the shell side. Gas feed rate 500 mL/min, liquid feed rate 10 mL/min.
Polymers 14 03783 g004
Figure 5. Effect of the sinewave frequency in the feed stream on the percent removal of carbon dioxide from the natural gas stream. The liquid feed rate is 10 mL/min; the gas feed rate is 500 mL/min.
Figure 5. Effect of the sinewave frequency in the feed stream on the percent removal of carbon dioxide from the natural gas stream. The liquid feed rate is 10 mL/min; the gas feed rate is 500 mL/min.
Polymers 14 03783 g005
Figure 6. Effect of feed gas flow rate on CO2 concentrations and removal flux at the optimum frequency of 0.05.
Figure 6. Effect of feed gas flow rate on CO2 concentrations and removal flux at the optimum frequency of 0.05.
Polymers 14 03783 g006
Figure 7. Effect of gas feed rate on the CO2 overall mass transfer coefficients under the optimum oscillating frequency (0.05) and in normal conditions at a constant liquid feed rate of 10 mL/min. Gas was introduced to the shell side and solvent into the tube side.
Figure 7. Effect of gas feed rate on the CO2 overall mass transfer coefficients under the optimum oscillating frequency (0.05) and in normal conditions at a constant liquid feed rate of 10 mL/min. Gas was introduced to the shell side and solvent into the tube side.
Polymers 14 03783 g007
Figure 8. Influence of feed perturbation on the CO2 removal flux under the effect of oscillatory and non-oscillatory gas feed rate. The oscillating mode was at the optimum frequency of 0.05 HZ.
Figure 8. Influence of feed perturbation on the CO2 removal flux under the effect of oscillatory and non-oscillatory gas feed rate. The oscillating mode was at the optimum frequency of 0.05 HZ.
Polymers 14 03783 g008
Table 1. Membrane characteristics used in the mathematical model at room temperature and ambient pressure.
Table 1. Membrane characteristics used in the mathematical model at room temperature and ambient pressure.
Number of fibers20
Hollow fiber, inner radius, mm0.21
Hollow fiber, outer radius, mm0.55
Module inside radius, m 0.008
Module length, m 0.25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ghasem, N. Modeling and Simulation of the Impact of Feed Gas Perturbation on CO2 Removal in a Polymeric Hollow Fiber Membrane. Polymers 2022, 14, 3783. https://doi.org/10.3390/polym14183783

AMA Style

Ghasem N. Modeling and Simulation of the Impact of Feed Gas Perturbation on CO2 Removal in a Polymeric Hollow Fiber Membrane. Polymers. 2022; 14(18):3783. https://doi.org/10.3390/polym14183783

Chicago/Turabian Style

Ghasem, Nayef. 2022. "Modeling and Simulation of the Impact of Feed Gas Perturbation on CO2 Removal in a Polymeric Hollow Fiber Membrane" Polymers 14, no. 18: 3783. https://doi.org/10.3390/polym14183783

APA Style

Ghasem, N. (2022). Modeling and Simulation of the Impact of Feed Gas Perturbation on CO2 Removal in a Polymeric Hollow Fiber Membrane. Polymers, 14(18), 3783. https://doi.org/10.3390/polym14183783

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