A simulation study of CO<sub>2</sub> gas adsorption with bottom ash adsorbent from palm oil mill waste using computational fluid dynamic (CFD)

Novi Sylvia; Aden Syahrullah Tarigan; Rozanna Dewi; Yunardi Yunardi; Yazid Bindar; Mutia Reza; Novi Sylvia; Aden Syahrullah Tarigan; Rozanna Dewi; Yunardi Yunardi; Yazid Bindar; Mutia Reza

doi:10.3934/environsci.2024022

AIMS Environmental Science

2024, Volume 11, Issue 3: 444-456. doi: 10.3934/environsci.2024022

Previous Article Next Article

Research article

A simulation study of CO₂ gas adsorption with bottom ash adsorbent from palm oil mill waste using computational fluid dynamic (CFD)

1.
Department of Chemical Engineering, Malikussaleh University, Lhokseumawe, 24351, Indonesia
2.
Department of Chemical Engineering, Syiah Kuala University, Banda Aceh, 23111, Indonesia
3.
Department of Chemical Engineering, Faculty of Industrial Technology, Institut Teknologi Bandung, Jl. Ganesha 10 Bandung 40132, Indonesia
4.
Department of Chemical Engineering, Institut Teknologi Kalimantan, Indonesia

Received: 04 March 2024 Revised: 01 April 2024 Accepted: 05 June 2024 Published: 11 June 2024

Biogas is a cost-effective, efficient, environmentally friendly, and renewable alternative energy source. While biogas contains CH₄, it also contains impurities in the form of 27–45% CO₂ gas. Therefore, it is necessary to purify biogas by removing CO₂ gas as it can reduce the calorific value of CH₄, the main component of biogas. The process of methane purification can be achieved through adsorption. Bottom ash, derived from palm oil mill waste, shows great potential for methane purification by effectively adsorbing CO₂. This research simulated the methane purification process using the computational fluid dynamics (CFD) method with the student version of the ANSYS R20 software. The study utilized an adsorbent made from bottom ash obtained from palm oil mill waste. The main objective was to investigate the performance of bottom ash as an adsorbent for removing CO₂ gas in a continuous gas flow within an adsorption column. The study involved varying the column bed height (4 cm, 8 cm, 12 cm) and gas flow rate (10 L/min, 15 L/min, 20 L/min). The results showed that the highest efficiency in removing CO₂ gas was 84.53% with a bed height of 12 cm and a flow rate of 10 L/min, while the lowest efficiency was 47.87% with a bed height of 4 cm and a flow rate of 20 L/min. Furthermore, the highest adsorption capacity for CO₂ gas was 1.64 mg/g with a bed height of 12 cm and a flow rate of 10 L/min, while the lowest capacity was 0.93 mg/g with a bed height of 4 cm and a flow rate of 20 L/min. The linearization of adsorption isotherm data indicated that the CO₂ gas adsorption process using bottom ash adsorbent followed the Langmuir model.

Keywords:

Citation: Novi Sylvia, Aden Syahrullah Tarigan, Rozanna Dewi, Yunardi Yunardi, Yazid Bindar, Mutia Reza. A simulation study of CO2 gas adsorption with bottom ash adsorbent from palm oil mill waste using computational fluid dynamic (CFD)[J]. AIMS Environmental Science, 2024, 11(3): 444-456. doi: 10.3934/environsci.2024022

Related Papers:

[1]	Lukluatus Syavika, Anugrah Ricky Wijaya, Alif Alfarisyi Syah . Synthesis and development of Fe₃O₄/SiO₂/CaCO₃ nanocomposite adsorbents for ammonia adsorption in the shrimp pond waste. AIMS Environmental Science, 2024, 11(6): 883-899. doi: 10.3934/environsci.2024044
[2]	Novi Sylvia, Husni Husin, Abrar Muslim, Yunardi, Aden Syahrullah, Hary Purnomo, Rozanna Dewi, Yazid Bindar . Design and performance of a cyclone separator integrated with a bottom ash bed for the removal of fine particulate matter in a palm oil mill: A simulation study. AIMS Environmental Science, 2023, 10(3): 341-355. doi: 10.3934/environsci.2023020
[3]	Premrudee Kanchanapiya, Supachai Songngam, Thanapol Tantisattayakul . The Adsorption of Perfluorooctanoic Acid on Coconut Shell Activated Carbons. AIMS Environmental Science, 2022, 9(2): 128-139. doi: 10.3934/environsci.2022010
[4]	Alif Alfarisyi Syah, Anugrah Ricky Wijaya, Irma Kartika Kusumaningrum . Preparation and development of amidoxime-modified Fe₃O₄/SiO₂ core-shell magnetic microspheres for enhancing U(VI) adsorption efficiency from seawater. AIMS Environmental Science, 2024, 11(1): 21-37. doi: 10.3934/environsci.2024002
[5]	Meriatna, Zulmiardi, Lukman Hakim, Faisal, Suryati, Mizwa Widiarman . Adsorbent performance of nipa (nypafruticans) frond in methylene blue dye degradation: Response surface methodology optimization. AIMS Environmental Science, 2024, 11(1): 38-56. doi: 10.3934/environsci.2024003
[6]	Nguyen Trinh Trong, Phu Huynh Le Tan, Dat Nguyen Ngoc, Ba Le Huy, Dat Tran Thanh, Nam Thai Van . Optimizing the synthesis conditions of aerogels based on cellulose fiber extracted from rambutan peel using response surface methodology. AIMS Environmental Science, 2024, 11(4): 576-592. doi: 10.3934/environsci.2024028
[7]	Isaiah Kiprono Mutai, Henry Kirimi Kiriamiti, Milton M M'Arimi, Robert Kimutai Tewo . Fluoride removal from water using Al(OH)$ _3 $-surface modified diatomite mixed with brick: optimization, isotherm and kinetic studies. AIMS Environmental Science, 2024, 11(3): 360-379. doi: 10.3934/environsci.2024017
[8]	Winai Meesang, Erawan Baothong, Aphichat Srichat, Sawai Mattapha, Wiwat Kaensa, Pathomsorn Juthakanok, Wipaporn Kitisriworaphan, Kanda Saosoong . Effectiveness of the genus Riccia (Marchantiophyta: Ricciaceae) as a biofilter for particulate matter adsorption from air pollution. AIMS Environmental Science, 2023, 10(1): 157-177. doi: 10.3934/environsci.2023009
[9]	Zainuddin Ginting, Adi Setiawan, Khairul Anshar, Ishak Ibrahin, Jalaluddin, Indri Riski Hasanah, Fitriyaningsih, Zetta Fazira . Optimization of the production of grade A liquid smoke from Patchouli solid residue using a distillation-adsorption process. AIMS Environmental Science, 2025, 12(1): 106-119. doi: 10.3934/environsci.2025005
[10]	Nasrul Arahman, Cut Meurah Rosnelly, Sri Mulyati, Wafiq Alni Dzulhijjah, Nur Halimah, Rinal Dia'ul Haikal, Syahril Siddiq, Sharfina Maulidayanti, Muhammad Aziz, Mathias Ulbricht . Investigation of microplastics in community well water in Banda Aceh, Indonesia: a separation technique using polyethersulfone-poloxamer membrane. AIMS Environmental Science, 2025, 12(1): 53-71. doi: 10.3934/environsci.2025003

Abstract

Mathematics subject classification:
Magnitude	Unit	Description	Magnitude	Unit	Description
q_e	mg/g	Sorption capacity	ρ	Kg m^-3	Fluid density
C₀	g/L	Initial concentration	α	m^-1	Viscous resistance coefficient
C_e	g/L	Final concentration	t	sec	Time
m	g	Mass	u	m sec^-1	Gas velocity (x-direction)
V	L	Volume	v	m sec^-1	Gas velocity (y-direction)
C₂	m	Inertial resistance Coefficient	w	m sec^-1	Gas velocity (z-direction)
µ	Ns m^-2	Fluid viscosity	ɛ		Bed voidage fraction (Porosity)

1. Introduction

The rapid development of technology and population growth have led to a significant increase in energy demand. However, due to the limited availability of fossil energy sources, this can result in an energy crisis and unmet energy needs. An indication of this energy crisis is the recent scarcity of oil, such as kerosene, petrol, and diesel. This scarcity is a result of the high and ever-increasing fuel demand each year. As a consequence, the price of crude oil worldwide continues to rise annually. Resolving this conflict requires serious attention. Therefore, there is a need for alternative renewable energy sources that are cost-effective, efficient, and environmentally friendly ^[1].

One such alternative energy source that meets these criteria is biogas. Biogas is not only rich in methane gas, but also has a simple and environmentally friendly production process. It is produced through the fermentation of organic waste, such as garbage, food scraps, animal feces, and waste from the food industry. Methane gas (CH₄), carbon dioxide gas (CO₂), oxygen gas (O₂), hydrogen sulfide gas (H₂S), hydrogen gas (H₂), and other gases are the components of biogas. Among these elements, methane gas (CH₄) and carbon dioxide gas (CO₂) play a crucial role in determining the quality of biogas. A high content of CH₄ results in a high calorific value of the biogas, while a high content of CO₂ leads to a low calorific value. Therefore, in order to increase the calorific value of biogas, the CO₂ content must be removed ^[2,3].

Various methods have been employed to reduce CO₂ emissions, including physical and chemical absorption, cryogenics, membrane separation, microalgae bio-fixation, and adsorption. Among these methods, the sorption method that utilizes amine compounds is commonly used to separate the gas component from the gas mixture by passing it through a liquid. However, this separation process has its drawbacks, such as amine mixture degradation due to impurities in the gas phase at low temperatures (below 50℃), equipment corrosion, high energy requirements for regeneration, degradation of solvent quality due to impregnation with other products, and amine reduction through evaporation. Among these methods, the adsorption method stands out as the most favorable due to its economic viability, ability to remove organic matter, and absence of toxic side effects. Furthermore, the adsorption method shows promise as an alternative for CO₂ sequestration, as it is relatively inexpensive, simpler, and does not produce liquid waste ^[4,5,6].

Removing CO₂ gas levels can be achieved using biosorbents made from boiler bottom ash. Palm oil mill boiler bottom ash is one type of adsorbent that can be used for gas purification because it contains carbon (C), oxygen (O), magnesium (Mg), aluminium (Al), silica (Si), phosphorus (P), sulfur (S), potassium (K), calcium (Ca), and iron (Fe). However, the primary components are carbon (C), oxygen (O), and silica (Si). It has been discovered that bottom ash contains the silica oxide (SiO₂) compound, which has demonstrated excellent activity in the CO₂ adsorption process ^[7,8,9,10].

Studying the adsorption process of CO₂ is of great interest, both experimentally and through simulation. Computational fluid dynamics (CFD) is a simulation method that is known for its cost-effectiveness and ability to validate experimental data. CFD is a field within fluid mechanics that employs numerical methods and algorithms to solve and analyze problems related to fluid flow, heat transfer, and chemical reactions through computer simulations ^[7,11].

Building upon the aforementioned information, this research involves adsorption utilizing a fixed-bed column system with adsorbents created from palm oil mill boiler waste, specifically bottom ash. Computational fluid dynamics (CFD) with Ansys 2020 R2 commercial software was employed in this study. CO₂ gas was introduced into the fixed-bed column, with a predetermined CO₂ composition. The height of the adsorbent bed and the incoming gas flow rate were varied. Figure 1 illustrates the research scheme. The simulation results were analyzed, and the impact of separation factors under different conditions was calculated. The objective of this study was to analyze the influence of bed height and the incoming gas flow rate on sorption capacity, sorption efficiency, and adsorption isotherms ^[8].

Figure 1. The mechanism of fixed-bed column simulation in ANSYS.

DownLoad: Full-Size Img PowerPoint

2. Research and method

2.1. Model geometry

The 3-D model geometry for the simulation in this study was developed with the help of design modeling using the application integrated in ANSYS Fluent. This was used to model a fixed-bed column to analyse the dynamic behavior of the adsorption process. For this study, the length of the cylinder was fixed at a 30-cm height and 6.4-cm diameter. The 3-D design for the fixed-bed column is illustrated in Figure 2. The boundary conditions of the study are shown in Figure 3.

Figure 2. The fixed-bed column geometry.

DownLoad: Full-Size Img PowerPoint

Figure 3. The boundary condition geometry.

DownLoad: Full-Size Img PowerPoint

2.1.1. Mesh generation

Before conducting simulations on the fixed-bed column geometry, it is crucial to determine the mesh value as it significantly impacts the generated cells. The assessment of mesh quality is carried out to identify the appropriate mesh size that balances computation time and result accuracy ^[12]. The accuracy value can be seen in Figure 4. As for this study, when viewed from the quality of the mesh through the scewness mess metrics spectrum, the quality value is obtained in the range of 0.00049633–0.355 with a total shell of 600,000. The mesh method used in this study is multizone with a prism mesh type. From the value range obtained before entering the post-processing stage, the mesh geometry value prepared is classified as very good. For a view of the fixed-bed column geometry at the mesh stage, see Figure 5.

Figure 4. The mesh value of the fixed-bed column geometry.

DownLoad: Full-Size Img PowerPoint

Figure 5. The validity of the mesh quality.

DownLoad: Full-Size Img PowerPoint

2.2. Mathematical model

The dynamic process of an adsorption column can be modeled and characterized using an integrated CFD model. To formulate a general adsorption model to fit the adsorption mechanism, the following assumptions are made ^[6]:

- Heat transfer within the bed is negligible;

- Competitive adsorption between CO₂ and N₂ is assumed with a mixtures model to represent mass transfer into the bed pellets;

- The mass transfer coefficient is equated to the external fluid film resistance and macropore diffusion;

- Porosity is uniform and constant;

- The equation used in the adsorption column is the continuity equation.

The generalized 3-D continuity equation for unsteady-state fluid flow is shown in Eq 1 ^[11,13,14]:

$\frac{\partial \rho }{\partial t}+\frac{\partial \rho u}{\partial x}+\frac{\partial \rho v}{\partial y}+\frac{\partial \rho w}{\partial z} = 0$

(1)

The Navier-Stokes equations for fluid flow through porous media where heat sources S_ix, S_iy, and S_iz are added to Eqs 2–4 to model flow resistance in 3-D dimensions ^[12,15] are:

Navier-Stokes equation in the x-direction

$\begin{array}{l} \frac{{\partial (\rho u)}}{{\partial t}} + \frac{{\partial (\rho uu)}}{{\partial x}} + \frac{{\partial (\rho vu)}}{{\partial y}} + \frac{{\partial (\rho wu)}}{{\partial z}} = - \frac{{\partial P}}{{\partial x}} + 2\frac{\partial }{{\partial x}}(\mu \frac{{\partial u}}{{\partial x}}) \hfill \\ - \frac{2}{3}\frac{\partial }{{\partial x}}\left[ {\mu \left( {\frac{{\partial u}}{{\partial x}} + \frac{{\partial v}}{{\partial y}} + \frac{{\partial w}}{{\partial z}}} \right)} \right] + \frac{\partial }{{\partial y}}\left[ {\mu \left( {\frac{{\partial u}}{{\partial y}} + \frac{{\partial v}}{{\partial x}}} \right)} \right] + \frac{\partial }{{\partial z}}\left[ {\mu \left( {\frac{{\partial u}}{{\partial z}} + \frac{{\partial w}}{{\partial x}}} \right)} \right] + {S_{ix}} \hfill \\ \end{array}$

(2)

Navier-Stokes equation in the y-direction

$\begin{array}{l} \frac{{\partial (\rho v)}}{{\partial t}} + \frac{{\partial (\rho vv)}}{{\partial y}} + \frac{{\partial (\rho uv)}}{{\partial x}} + \frac{{\partial (\rho wv)}}{{\partial z}} = \frac{{\partial P}}{{\partial y}} + 2\frac{\partial }{{\partial y}}\left( {\mu \frac{{\partial v}}{{\partial y}}} \right) \hfill \\ - 2\frac{\partial }{{\partial y}}\left[ {\mu \left( {\frac{{\partial u}}{{\partial x}} + \frac{{\partial v}}{{\partial y}} + \frac{{\partial w}}{{\partial z}}} \right)} \right] + \frac{\partial }{{\partial x}}\left[ {\mu \left( {\frac{{\partial v}}{{\partial x}} + \frac{{\partial u}}{{\partial y}}} \right)} \right] + \frac{\partial }{{\partial z}}\left[ {\mu \left( {\frac{{\partial v}}{{\partial z}} + \frac{{\partial w}}{{\partial y}}} \right)} \right] + {S_{iy}} \hfill \\ \end{array}$

(3)

Navier-Stokes equation in the z-direction

$\begin{array}{l} \frac{{\partial (\rho w)}}{{\partial t}} + \frac{{\partial (\rho ww)}}{{\partial z}} + \frac{{\partial (\rho uw)}}{{\partial x}} + \frac{{\partial (\rho vw)}}{{\partial y}} = \frac{{\partial P}}{{\partial z}} + 2\frac{\partial }{{\partial z}}\left( {\mu \frac{{\partial w}}{{\partial z}}} \right) \hfill \\ - \frac{2}{3}\frac{\partial }{{\partial z}}\left[ {\mu \left( {\frac{{\partial u}}{{\partial x}} + \frac{{\partial v}}{{\partial y}} + \frac{{\partial w}}{{\partial z}}} \right)} \right] + \frac{\partial }{{\partial x}}\left[ {\mu \left( {\frac{{\partial w}}{{\partial x}} + \frac{{\partial u}}{{\partial z}}} \right)} \right] + \frac{\partial }{{\partial y}}\left[ {\mu \left( {\frac{{\partial w}}{{\partial y}} + \frac{{\partial v}}{{\partial z}}} \right)} \right] + {S_{iz}} \hfill \\ \end{array}$

(4)

The S_i porous media momentum source can calculate the pressure gradient in the packed bed and generate a pressure drop proportional to the fluid velocity, as shown in Eq 5 ^[12]:

${S_i} = \frac{\mu }{\alpha }{\mu _i} + {C_2}\left( {\frac{1}{2}\rho {u_i}\left| {\left. {{u_i}} \right|} \right.} \right)$

(5)

where

${C_2} = \frac{{1.75}}{{{D_p}}}\frac{{(1 - \varepsilon )}}{{{\varepsilon ^3}}}$

(6)

$\alpha = \frac{{D{{_p^2}^{}}}}{{150}}\frac{{{\varepsilon ^3}}}{{{{(1 - \varepsilon )}^2}}}$

(7)

where C₂ and α are the inertial resitance and viscous resitance, which are estimated by Eqs 6 and 7. We can calculate the percent absorption and absorption capacity (q_e) using the following Eq 8 ^[7]:

${q_e} = \frac{{({C_0} - {C_e})V}}{m}$

(8)

$removal\;\% = \frac{{({C_0} - {C_e})}}{{{C_0}}}$

(9)

where q_e is the amount of heavy metal adsorbed (mg/g), C₀ and C_e are the initial and final concentrations of gas, respectively, V is the volume of solution (L), and m is the mass of adsorbent (g).

2.3. Boundary conditions

The adsorbent material used in this study is bottom ash, with a particle size of 46 µm. The gas employed is a mixture of CO₂ and N₂, with variable flow rates when passed through the fixed-bed column. The experimental procedure is outlined in Figure 1. The mole fraction of CO₂ and N₂ in the gas mixture is 95% and 5%, respectively. The system operates at a pressure of 1 atm and a temperature of 25 ℃. Detailed specifications for the gas flow rate and bed height are provided in Tables 1 and 2.

Table 1. Simulation data of the adsorption process.

Operating Conditions	Value
Inlet gas flow rate (L/min)	10, 15, 20
Bed height (cm)	4, 8, 12
Inlet gas composition	95 % CO₂; 5 % N₂
Column height (cm)	30
Porosity, ɛ	0.88

| Show Table

DownLoad: CSV

Table 2. The spatial discretization schemes.

Variabel	Discrete scheme
Pressure-velocity coupling	Coupled
Pressure	Presto
Density	Quick
Momentum	Second order upwind
Species	Second order upwind
Energy	Second order upwind

| Show Table

DownLoad: CSV

3. Results

3.1. The effect of bed height and flow rate on CO₂ removal efficiency

The impact of flow rate on CO₂ removal efficiency is illustrated in Figure 6(a). The findings demonstrate that a higher flow rate corresponds to a lower CO₂ removal efficiency. This is attributed to the fact that a higher CO₂ gas flow rate reduces the contact time between CO₂ and the adsorbent, thereby leading to a decrease in the percentage of CO₂ removal efficiency ^[16,17,18]. Similar findings were also observed in a study conducted by Dupre et al. ^[19], wherein the flow rate significantly influenced the results of CO₂ gas absorption. Decreasing the flow rate extends the breakthrough time. By employing a low flow rate, an extended breakthrough time can be achieved, which in turn enhances the contact time and mass transfer of CO₂ gas or the diffusion process on the adsorbent ^[11].

Figure 6. The effect of flow rate (a) and bed height (b) on CO₂ removal efficiency.

DownLoad: Full-Size Img PowerPoint

Based on Figure 6(b), it is evident that as the height of the adsorbent bed in the adsorption column increases, the efficiency of CO₂ removal also increases. This correlation can be attributed to the surface area and active sites of the adsorbent, where taller beds provide a larger surface area and greater adsorption capacity. Additionally, the highest bed contains a larger quantity of adsorbent particles, thereby facilitating enhanced contact between the solid and CO₂ gas as it flows through the column ^[14]. A higher height of the bed can increase the contact time or delay the breakthrough time of the adsorbent, which means the bed is able to be operated for a longer period of time without changing the type of adsorbent being used. Instead, for a lower bed, the breakthrough time will increase or approach faster. Therefore, the performance decreases ^[4]. The results of this study can be seen in Figures 6 and 7, where the highest CO₂ removal efficiency was obtained in a bed with a height of 12 cm by 84.53% with a flow rate of 10 L/min, while the lowest CO₂ removal efficiency was obtained 47.87% in a bed with a height of 4 cm with a flow rate of 20 L/min ^[16,17].

Figure 7. The effect of flow rate and bed height on CO₂ gas absorption capacity.

DownLoad: Full-Size Img PowerPoint

3.2. The effect of flow rate and bed height on CO₂ adsorption capacity

Based on Figure 7 (a), (b), it is evident that the flow rate has an impact on the adsorption capacity of CO₂ gas. Specifically, when the flow rate in the adsorption column is higher, the absorption capacity of CO₂ gas decreases. This can be attributed to the fact that increasing the flow rate reduces the contact time and mass transfer between the adsorbent and CO₂. Additionally, a higher flow rate reduces the amount of time the adsorbate (CO₂) spends within the adsorbent, leading to a lower adsorption capacity for CO₂ gas. Conversely, decreasing the flow rate of the incoming gas enhances the contact time and mass transfer between the adsorbent and CO₂. In other words, a lower flow rate within the adsorption column allows for a longer residence time of CO₂ molecules, leading to a higher amount of CO₂ being adsorbed and more time for the adsorbate (CO₂) to be absorbed by the adsorbent. Ultimately, this results in an increased CO₂ adsorption capacity ^[1]. The highest adsorption capacity was obtained at a flow rate of 10 L/min with a bed height of 12 cm with an adsorption capacity of 1.64 mg/g, while the lowest adsorption capacity was at a rate of 20 L/min with a bed height of 4 cm with an absorption capacity of 0.93 mg/g ^{[3,6,18,19,20,21]}.

3.3. The contour of CO₂ removal by the simulation result

The simulated fixed-bed column was made with a cylindrical geometry, with the flow entering evenly from the bottom and exiting at the top during the adsorption process ^[7,17]. Figure 8 shows the schematic geometry of the fixed-bed column, shaft zone, and gas mixture inflow and outflow. A two-dimensional (2-D) axisymmetric CFD model was created to view the contours of concentration and velocity in the fixed-bed column, which can be seen in Figure 9 for concentration and Figure 10 for velocity.

Figure 8. The skema fixed-bed column.

DownLoad: Full-Size Img PowerPoint

Figure 9. The contours of CO₂ concentration distribution with flow rate and bed height: (a) 10 L/m at 4 cm, (b) 15 L/m at 4 cm, (c) 20 L/m at 4 cm, (d) 10 L/m at 8 cm, (e) 15 L/m at 8 cm, (f) 20 L/m at 8 cm, (g) 10 L/m at 12 cm, (h) 15 L/m at 12 cm, and (i) 20 L/m at 12 cm.

DownLoad: Full-Size Img PowerPoint

Figure 10. CO₂ velocity contours with flow rate and bed height: (a) 10 L/m at 4 cm, (b) 10 L/m at 8 cm, (c) 10 L/m at 12 cm, (d) 15 L/m at 4 cm, (e) 15 L/m at 8 cm, (f)15 L/m at 12 cm, (g) 20 L/m at 4 cm, (h) 20 L/m at 8 cm, and (i) 20 L/m at 12 cm.

DownLoad: Full-Size Img PowerPoint

Based on Figure 9, it is evident that the contour results show a significant decrease in CO₂ gas concentration when passing through the bottom ash adsorbent. Figure 10 illustrates the changes in concentration that occur in the porous zone. In this zone, a portion is visualized as an adsorbent bed that is responsible for absorbing CO₂ gas. The decrease in concentration is a result of the CO₂ gas being absorbed through the contact between CO₂ and bottom ash. This absorption process takes place in the porous zone area, which has a surface area of 0.93 m2/g and a pore volume of 0.12 cm3/g ^[3,10,19,22]. Bottom ash has 3 types of pores, namely, macropore (55 nm and > 55 nm), mesopore (7 nm and 8 nm), and micropore (<1 nm) ^{[4,6,8,9,16,17,20]}. Gas adsorption capability is dependent on the texture of the adsorbent. CO₂ gas adsorption requires adsorbents with micropore size (<1 nm) because the gas has a small molecular size of < 4 nm. The decrease in CO₂ gas concentration in the porous zone also occurs because it is dependent on Vanderwaals interactions related to the polarity of adsorbate and adsorbent molecules ^{[5,7,9,10,13,15,18,20,21]}.

Based on the contours results in Figure 10, it can be seen that the higher the bed, the lower the CO₂ velocity. The decrease in velocity in the bed section is because the greater the gas velocity, the greater the pressure drop that occurs ^{[9,11,12,13,15,16,17,21,23]}.

3.4. Adsorption isotherms

According to Table 3 of the Langmuir isotherm, the R² value is 0.995, the qm value is 0.567, and the Langmuir constant is -2.694. In contrast, the Freundlich isotherm yields an R² value of 0.965. Based on the obtained y-value, the n-value is -1.19, and the Freundlich constant is 11.74. By comparing the determination values of R in the Langmuir isotherm and the Freundlich isotherm, this study suggests that the bonding of CO₂ with the bottom ash adsorbent follows the Langmuir model. This indicates that the absorption process occurs in a monolayer, as supported by the coefficient of determination R being closest to 1 ^{[1,2,4,7,12,13,16,18,21]}.

Table 3. The adsorption isotherm model at 1 atm pressure and 25 ℃ temperature.

Type	Equation Nonliner		Linear
Langmuir		$q=\frac{q_m \times k_L \times C_e}{1+k_L \times C_e}$	$\frac{C_e}{q}=\frac{1}{q_m \times k_L}+\frac{C_e}{q_m}$
q	0.567
k	-2.694
R	0.995
Freundlich		$q=k \times C_e{ }^{\frac{1}{n}}$	$\log q=\log k+\frac{1}{2} \log C_{e}$
k	11.74
n	-1.19
R	0.965

| Show Table

DownLoad: CSV

3. Conclusions

The study found that increasing the height of the bed resulted in higher efficiency and capacity for removing CO2 gas. However, as the flow rate increased, the efficiency and capacity decreased. The highest CO2 gas removal efficiency was 84.53% with a bed height of 12 cm and a flow rate of 10 L/min, while the lowest was 47.87% with a bed height of 4 cm and a flow rate of 20 L/min. Similarly, the highest CO2 gas adsorption capacity was 1.64 mg/g with a bed height of 12 cm and a flow rate of 10 L/min, and the lowest was 0.93 mg/g with a bed height of 4 cm and a flow rate of 20 L/min. The results also showed that the adsorption process using bottom ash as an adsorbent following the Langmuir model, indicating that it occurred in a monolayer manner. Bottom ash, a waste product from palm oil mill boilers, can be used as an alternative adsorbent for CO2 adsorption. In future research, a fixed-bed column system and bottom ash from palm oil mills could be used for CO2 adsorption from biogas.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

The authors are grateful to Malikussaleh university of chemical engineering simulation and process laboratory team

Conflict of interest

The authors have no conflicts of interest to declare.

References

[1]	Li H, Tan Y, Ditaranto M, et al. (2017) Capturing CO₂ from biogas plants. Energy procedia 114: 6030-6035. https://doi.org/10.1016/j.egypro.2017.03.1738 doi: 10.1016/j.egypro.2017.03.1738
[2]	Li M, Zhu Z, Zhou M, et al. (2021) Removal of CO₂ from biogas by membrane contactor using PTFE hollow fibers with smaller diameter. J Membr Sci 627: 119232. https://doi.org/10.1016/j.memsci.2021.119232 doi: 10.1016/j.memsci.2021.119232
[3]	Oyedeji OA (2019) Understanding and modeling the formation of syngas contaminants during biomass gasification, University of Tennessee
[4]	Cho YM, Yang YM, Park DS, et al. (2013) Study on CO₂ adsorption on LiOH-modified Al₂O₃. Appl Mech Mater 284: 342-346. https://doi.org/10.4028/www.scientific.net/AMM.284-287.342. doi: 10.4028/www.scientific.net/AMM.284-287.342
[5]	Ochedi FO, Liu Y, Adewuyi YG (2020) State-of-the-art review on capture of CO₂ using adsorbents prepared from waste materials. Process Saf Environ 139: 1-25 https://doi.org/10.1016/j.psep.2020.03.036. doi: 10.1016/j.psep.2020.03.036
[6]	Chen Q, Rosner F, Rao A, et al. (2019) Simulation of elevated temperature solid sorbent CO₂ capture for pre-combustion applications using computational fluid dynamics. Appl energy 237: 314-325. https://doi.org/10.1016/j.apenergy.2019.01.042 doi: 10.1016/j.apenergy.2019.01.042
[7]	Pham TD, Xiong R, Sandler SI, et al. (2014) Experimental and computational studies on the adsorption of CO₂ and N₂ on pure silica zeolites. Microporous Mesoporous Mater 185: 157-166. https://doi.org/10.1016/j.micromeso.2013.10.030 doi: 10.1016/j.micromeso.2013.10.030
[8]	Sylvia N, Fitriani F, Dewi R, et al. (2021) Characterization of bottom ash waste adsorbent from palm oil plant boiler burning process to adsorb carbon dioxide in a fixed bed column. Indones J Chem 21: 1454-1462. https://doi.org/10.22146/ijc.66509 doi: 10.22146/ijc.66509
[9]	Sethupathi S, Chung LL, Younas M, et al. (2017) Gaseous pollutant removal using solid wastes adsorbents, " in : Air, Gas, and Water Pollution Control Using Industrial and Agricultural Solid Wastes Adsorbents: CRC Press, 3-20.
[10]	Promraksa A, Rakmak N (2020) Biochar production from palm oil mill residues and application of the biochar to adsorb carbon dioxide. Heliyon 6: e04019. https://doi.org/10.1016/j.heliyon.2020.e04019 doi: 10.1016/j.heliyon.2020.e04019
[11]	Sylvia N, Mutia R, Dewi R, et al. (2019) A computational fluid dynamic comparative study on CO₂ adsorption performance using activated carbon and zeolite in a fixed bed reactor. IOP Conf Ser Mater Sci Eng 536: 012042 https://iopscience.iop.org/article/10.1088/1757-899X/536/1/012042 doi: 10.1088/1757-899X/536/1/012042
[12]	Nouh S, Lau KK, Shariff AM (2010) Modeling and simulation of fixed bed adsorption column using integrated CFD approach. J Appl Sci 24: 3229-3235. 10.3923/jas.2010.3229.3235 doi: 10.3923/jas.2010.3229.3235
[13]	Li S, Deng S, Zhao L, et al. (2018) Mathematical modeling and numerical investigation of carbon capture by adsorption: literature review and case study. Appl Energy 221: 437-449. 10.1016/j.apenergy.2018.03.093 doi: 10.1016/j.apenergy.2018.03.093
[14]	Mesfer MKA, Danish M, Khan MI, et al. (2020) Continuous fixed bed CO₂ adsorption: breakthrough, column efficiency, mass transfer zone. Processes 8: 1233. DOI: 10.3390/pr8101233 doi: 10.3390/pr8101233
[15]	Sylvia N, Dewi R, Aprilia B, et al. (2023) Simultaneous fine particulate matter separation and CO₂ adsorption in a cyclone separator with a fixed bed bottom ash from a palm oil mill boiler: a simulation study. Math Model Eng Probl 10: 597-604. https://doi.org/10.18280/mmep.100229 doi: 10.18280/mmep.100229
[16]	Akpasi SO, Isa YM (2022) Effect of operating variables on CO₂ adsorption capacity of activated carbon, kaolinite, and activated carbon–Kaolinite composite adsorbent. Water-Energy Nexus 5: 21-28. https://doi.org/10.1016/j.wen.2022.08.001 doi: 10.1016/j.wen.2022.08.001
[17]	Shafeeyan MS, Daud WM, Shamiri A, et al. (2015) Modeling of carbon dioxide adsorption onto ammonia-modified activated carbon: kinetic analysis and breakthrough behavior. Energy Fuels 29: 6565-6577. https://doi.org/10.1021/acs.energyfuels.5b00653. doi: 10.1021/acs.energyfuels.5b00653
[18]	Kan X, Chen X, Shen Y, et al. (2019) Box-Behnken design based CO₂ co-gasification of horticultural waste and sewage sludge with addition of ash from waste as catalyst. Appl Energy 242: 1549-1561. https://doi.org/10.1016/j.apenergy.2019.03.176 doi: 10.1016/j.apenergy.2019.03.176
[19]	Dupre KR, Vyas A, Goldfarb JL, et al. (2019) Investigation of computational upscaling of adsorption of SO₂ and CO₂ in fixed bed columns. Adsorption 25: 773-782. https://doi.org/10.1007/s10450-019-00050-4 doi: 10.1007/s10450-019-00050-4
[20]	Umenweke GC, Afolabi IC, Epelle EI, et al. (2022) Machine learning methods for modeling conventional and hydrothermal gasification of waste biomass: A review. Bioresour Technol Rep 17: 100976. https://doi.org/10.1016/j.biteb.2022.100976 doi: 10.1016/j.biteb.2022.100976
[21]	Yang H, Li W, Liu J, et al. (2019) Polyethylenimine-impregnated resins: The effect of support structures on selective adsorption for CO₂ from simulated biogas. Chem Eng J 35: 822-829. https://doi.org/10.1016/j.cej.2018.08.220 doi: 10.1016/j.cej.2018.08.220
[22]	Gorbounov M, Petrovic B, Ozmen S, et al. (2023) Activated carbon derived from Biomass combustion bottom ash as solid sorbent for CO₂ adsorption. Chem Eng Res Des 194: 325-343. https://doi.org/10.1016/j.cherd.2023.04.057 doi: 10.1016/j.cherd.2023.04.057
[23]	Liang F, Wei S, Wu L, et al. (2023) Improving the value of CO₂ and biogas slurry in agricultural applications: A rice cultivation case. Environ Exp Bot 208: 105233. https://doi.org/10.1016/j.envexpbot.2023.105233 doi: 10.1016/j.envexpbot.2023.105233

Reader Comments

Your name:*

Email:*
© 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Environmental Science

1.6 2.9

Metrics

Article views(1440) PDF downloads(182) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(10) / Tables(3)

AIMS Environmental Science

A simulation study of CO₂ gas adsorption with bottom ash adsorbent from palm oil mill waste using computational fluid dynamic (CFD)

Related Papers:

Abstract

1. Introduction

2. Research and method

2.1. Model geometry

2.1.1. Mesh generation

2.2. Mathematical model

2.3. Boundary conditions

3. Results

3.1. The effect of bed height and flow rate on CO₂ removal efficiency

3.2. The effect of flow rate and bed height on CO₂ adsorption capacity

3.3. The contour of CO₂ removal by the simulation result

3.4. Adsorption isotherms

3. Conclusions

Use of AI tools declaration

Acknowledgments

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

AIMS Environmental Science

A simulation study of CO2 gas adsorption with bottom ash adsorbent from palm oil mill waste using computational fluid dynamic (CFD)

Related Papers:

Abstract

1. Introduction

2. Research and method

2.1. Model geometry

2.1.1. Mesh generation

2.2. Mathematical model

2.3. Boundary conditions

3. Results

3.1. The effect of bed height and flow rate on CO2 removal efficiency

3.2. The effect of flow rate and bed height on CO2 adsorption capacity

3.3. The contour of CO2 removal by the simulation result

3.4. Adsorption isotherms

3. Conclusions

Use of AI tools declaration

Acknowledgments

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

A simulation study of CO₂ gas adsorption with bottom ash adsorbent from palm oil mill waste using computational fluid dynamic (CFD)

3.1. The effect of bed height and flow rate on CO₂ removal efficiency

3.2. The effect of flow rate and bed height on CO₂ adsorption capacity

3.3. The contour of CO₂ removal by the simulation result