Fei Zhao1, Ling Zhu1, Zhenzhong Wang2, Yan Hou1, Jiaqing Chen1, Chunyu Wang1,3, Danyun Xu1. 1. Department of Environmental Engineering, Beijing Institute of Petrochemical Technology, Beijing 102617, China. 2. SINOPEC Research Institute of Safety Engineering, Qingdao, Shandong 266071, China. 3. College of Environmental and Energy Engineering, Beijing University of Technology, Beijing 100124, China.
Abstract
In this work, the adsorption parameters of n-butane vapor on an absorbent were tested following the fixed-bed method. According to the corresponding experiments, the maximum adsorption capacity and breakthrough time of activated carbon (AC) are 0.2674 g·g-1 and 924 min, respectively. According to the two-energy-state model formula and the classical adsorption heat formula, the values of theoretical and actual adsorption heat of AC adsorbing n-butane are 5.48 and 5.56 kJ·mol-1, respectively. The model for adsorption of n-butane by an AC fixed bed is based on the analytical solutions to the mass, momentum, and energy conservation equations. The model is built using porous media zone in ANSYS Fluent, the implementation of the model into ANSYS Fluent under user-defined functions (UDFs) is also described, the mass source term Si and energy source term S T are loaded into Fluent through UDF, and then the mass- and heat-transfer processes of AC in the absorption of n-butane are simulated. Furthermore, the predictions by ANSYS Fluent are compared with in situ experimental data, and the deviation rate of breakthrough time and temperature of six monitoring points is less than 5%. The results verify the accuracy and feasibility of computational fluid dynamics (CFD). Therefore, the model can be used to predict the engineering application of the adsorption of organic gases by various porous media.
In this work, the adsorption parameters of n-butane vapor on an absorbent were tested following the fixed-bed method. According to the corresponding experiments, the maximum adsorption capacity and breakthrough time of activated carbon (AC) are 0.2674 g·g-1 and 924 min, respectively. According to the two-energy-state model formula and the classical adsorption heat formula, the values of theoretical and actual adsorption heat of AC adsorbing n-butane are 5.48 and 5.56 kJ·mol-1, respectively. The model for adsorption of n-butane by an AC fixed bed is based on the analytical solutions to the mass, momentum, and energy conservation equations. The model is built using porous media zone in ANSYS Fluent, the implementation of the model into ANSYS Fluent under user-defined functions (UDFs) is also described, the mass source term Si and energy source term S T are loaded into Fluent through UDF, and then the mass- and heat-transfer processes of AC in the absorption of n-butane are simulated. Furthermore, the predictions by ANSYS Fluent are compared with in situ experimental data, and the deviation rate of breakthrough time and temperature of six monitoring points is less than 5%. The results verify the accuracy and feasibility of computational fluid dynamics (CFD). Therefore, the model can be used to predict the engineering application of the adsorption of organic gases by various porous media.
According to the World Health Organization (WHO), volatile organic
compounds (VOCs) are organic compounds with a Reid vapor pressure
exceeding 10.3 Pa at the normal temperature (293.15 K) and pressure
(101.325 kPa).[1] VOCs are a large group
of carbon-based chemicals that easily evaporate at room temperature,
including alkanes, aromatics, esters, alkenes, carboxylic acids, and
alcohols.[2−7] VOCs are classified as major contributors to air pollution owing
to their persistence and tendency to accumulate in the environment.[8,9]They contribute to pollution both indirectly and directly
as ozone/smog
precursors and substances toxic to the environment, respectively.[10,11]Activated carbon (AC) is one of the most effective absorbents
for
VOCs owing to its specific surface area, pore-size distribution,[8,12−16] and the chemical affinity between the adsorbents and VOC pollutants.
To further improve the adsorption efficiency, deep investigations
into modified AC materials have received more attention, with modification
processes including the rebuilding of pore structures to match the
pore size of adsorbents to the molecular diameter of gaseous VOCs.
Transport phenomena, such as heat and mass transfer in porous AC,
have an important effect on the adsorption processes and their efficiency.
An absorbent with a higher adsorption heat and lower heat-transfer
efficiency would cause an increase in the temperature of the adsorbent
bed, which is the major factor causing the adsorption performance
to degrade and affecting the regeneration efficiency. Moreover, when
the adsorbent bed temperature exceeds the flashing point of VOCs,
some serious accidents would occur, such as fire or an explosion in
the adsorption device. A computational fluid dynamics (CFD) approach
is used for modeling transport and adsorption in porous media using
Fluent. Much work has already been conducted in this area.[17−20]For example, Bahadori et al. investigated the effect of variations
in viscosity on natural convection in a saturated enclosure, predicting
the flow field and heat transfer in double-layer porous media using
the Darcy model in a CFD commercial package, but did not simulate
or analyze the concentration field.[21] Hou
et al. established the computational domain model of the onboard refueling
vapor recovery (ORVR) carbon canister and simulated the internal structure
of carbon using a porous medium model. The internal pressure fields
of the carbon canister’s adsorption and desorption processes
were also numerically simulated by Fluent,[22] while the temperature and concentration were not considered. Bai
et al. simulated fixed-bed adiabatic adsorption/desorption processes
and used the linear driving force model for the heat- and mass-transfer
rates. They also used a two-dimensional cylindrical canister and three-dimensional
automotive production canister geometry to study the adsorption/desorption
processes of carbon dioxide in a helium carrier gas on Norit B4 AC.[23] However, their simulation results differed from
the experimental results. Lin et al. presented a test process based
on the analytical simulation of the butane working capacity (BWC)
of an automotive carbon canister and based their computational model
of the fixed-bed system of a carbon canister on nonequilibrium, nonisothermal,
and nonadiabatic algorithms to simulate the real-life loading/purging
of hydrocarbon vapors from this device.[24] However, the thermal transportation and temperature fields were
not considered in their studies. Hou established the unsteady three-dimensional
mathematical model for the absorption of benzene by AC absorption,
and the effects of the inlet benzene concentration and gas velocity
on outlet benzene concentration, temperature, and pressure distribution
in the U-tube were investigated.[25,26] The fixed
bed was considered an integral solid, so the effects of bed porosity
and transfer resistance were not investigated.In this paper,
adsorption experiments, theoretical calculations
of adsorption heat, and numerical simulations were used. The objective
of this work is to investigate the mass- and heat-transfer properties
during the adsorption of n-butane on AC in the fixed
bed. The adsorption experiment was carried out on the fixed bed with
a frequently used AC. The temperature of the fixed bed and the concentration
of n-butane at the outlet were measured. Several
kinetic models were used to fit and analyze the adsorption breakthrough
curve, which was measured in the adsorption experiment. According
to the formula of adsorption heat, the theoretical adsorption heat
and the actual adsorption heat of n-butane adsorbed
on AC were calculated. Based on the experimental data, the expressions
of mass source term (S) and energy source
term (ST) for the simulation model were
obtained. The fixed-bed model was implemented following a finite volume
method in the commercial CFD package Fluent, the same initial conditions
as the adsorption experiment were setup in Fluent. The mass- and heat-transfer
processes of AC adsorbing n-butane were simulated
by loading the S and ST with user-defined functions (UDFs). Validations of the
newly introduced model predictions are demonstrated through comparisons
with experimental results.
Results and Discussion
To understand the heat and mass transfer during adsorption on the
fixed bed, a mathematical model of the adsorption of n-butane on an AC fixed bed was established based on mass, energy,
and momentum equations. The Integrated Computer Engineering and Manufacturing
(ICEM) code was used to establish the physical model and generate
the grid, and Fluent was used for simulations. Adsorption was modeled
as a flow through a porous medium, with the mass and energy source
terms defined in the UDF. The binary gas mixture flow in the porous
zone with the laminar fluid flow was simulated using the built-in
solver. The results obtained were quantitatively analyzed with a three-dimensional
numerical model that solves the nonstationary continuity, momentum,
and energy equations, which allows the spatial distributions of the
temperature and adsorbed quantities to be simulated.
Physical
Model
The geometric diagram
of the fixed bed presented in Figure consists of five parts. The inlet size was R1 = 2 mm, z1 = 15
mm; the cavity radius was R2 = 17.5 mm,
and the total length was 480 mm. The front buffer length was z2 = 20 mm, the AC particle area length was z3 = 410 mm, and the tail buffer length was z4 = 20 mm; the outlet (R3, z5) was of the same size as
the inlet. There were six temperature measuring points in the vertical
direction on the side of the fixed bed with a straight tube, which
were marked as points 1–6, respectively.
Figure 1
Geometric sketch of the
fixed bed.
Geometric sketch of the
fixed bed.The size of the physical model
was set according to the actual
size of the fixed bed. The structured and unstructured grids were
combined to divide the mesh. The number of meshes was 12 871.
The mesh of the geometric body was divided, as shown in Figure .
Figure 2
Fixed-bed model mesh-generation
diagram.
Fixed-bed model mesh-generation
diagram.
Mathematical
Model
The adsorption
of n-butane by AC is a three-dimensional, unsteady
process that follows the mass, momentum, and energy conservation laws.
Fluent was used to simulate the adsorption process, and the mathematical
models for numerical simulations were based on the mass, momentum,
and energy conservation equations. A complete mathematical model with
mass, momentum, and energy equations was established.
Mass Conservation Equation
The
mass conservation equation (eq )[27−31] can be expressed aswhere ρf is the fluid density
(kg·m–3), t is the time (s),
ε is the porosity of AC, u is the fluid local
velocity (m·s–1), and Sm is the generalized mass source term (kg·m–3·s–1), and S is
the mass source term of gas i (kg·m–3·s–1).The adsorbed phase mass conservation eqs and 3(27,32) can be expressed aswhere Y is the
mass fraction of gas i, D is the diffusion coefficient of gas i, M is the amount of the substance of gas i (kg·m–3), q is the gas concentration (kg·m–3), and ρAC is the density of AC (kg·m–3).
Momentum Conservation Equation
The momentum conservation equation (eq )[27,30−32] can be expressed
aswhere u is the
velocity of gas i (m·s–1),
μ is the viscosity (Pa·s), and SV is the momentum source term (N·m–3).The porous medium model is a momentum source term superimposed on
the momentum equation. The source term includes two parts, i.e., the
viscous and inertial loss terms, as shown in eq (37,40)where α is the Darcy
permeability coefficient
(m2), C2 is the inertial resistance
(m–1), νm is the
velocity (m·s–1), and ν is the velocity for i (x, y, or z)
(m·s–1).There are two resistance coefficients
in the porous simulation
process, i.e., viscous resistance 1/α and inertial resistance C2, which can be calculated according to Ergun
equation (eq )[30,33]where ε is the porosity of AC
(0.38), Dp is the mean particle diameter
of AC (0.00085 m), and 1/α
is the viscous resistance (m–2).
Energy Conservation Equation
The
energy conservation equation[27−29] can be expressed aswhere cp is the
specific heat capacity (J·kg–1·°C–1), T is the temperature (°C), kh is the heat-transfer coefficient of the fluid
(W·K·m–2), and ST is the generalized energy source term (W·m–3). x, y, and z are three directions of the space coordinate system and u, v, and w are the local
fluid velocities in the x, y, and z directions (m·s–1).
Adsorbed Phase Mass Balance Equation
The adsorbed phase
mass balance equation[32,34] for component i can be expressed aswhere k is the mass-transfer coefficient and q* is the gas equilibrium concentration.According to the first-order kinetics fitting equation q = 51.822 – 51.822 e–0.02386 (Figure ), it can be calculated that
Figure 14
Fitting curve of first-order
adsorption kinetics equation.
The adsorbed phase mass source
term S in eq can be expressed aswhere mC is the mass of n-butane (g), mAC is the mass
of AC (g), VAC is the volume of AC (m3), and ρAC is the density of AC (550 kg·m–3).Adsorption heat of n-butane
on AC: Qa,a = 5.56 kJ·mol–1, the energy
source term ST in eq can be expressed aswhere MC is the molar mass of n-butane
(58.12 g·mol–1).To generate a simulation
result with acceptable accuracy in comparison
to the experimental data, it is necessary to select models for the
specific situation in Fluent. When Fluent is used for calculation,
the N–S equation is set as the governing equation and the calculation
type is set to be implicit and pressure-based. The physical model
is imported into Fluent, and the species transfer model, laminar flow
model, and porous medium model were used. The inlet velocity boundary
condition was applied at the left end of the fixed bed, and a pressure
outlet boundary condition was applied at the exit. The mass source
term S and energy source
term ST were loaded into Fluent through
UDF, and the model of n-butane adsorption on AC fixed
bed was established.
Simulation Results and
Discussion
The initial parameters of the model are set according
to the experimental
conditions. The inlet condition of the model is the velocity inlet,
the boundary condition is the wall, and the outlet condition is the
outlet vent. The inlet speed was set to 0.5305 m·s–1, the inlet concentration was 0.02 (VC/V), and the initial temperature was 300 K, the viscous resistance
of the porous medium area was set to 5.5267 × 108 m–2, and the inertial resistance of the porous medium
area was set to 4.6525 × 104 m–1.It took 895 min to simulate the whole adsorption process
using Fluent. After simulations, the temperature and concentration
cloud images were obtained using the Tecplot postprocessing software.
The adsorption process was analyzed using the experimental data, simulated
results, and simulated cloud images.
Cloud
Image Analysis
In this study,
the cloud images of six points in the fixed bed were selected for
analysis.As is known, during gas-phase adsorption, a high amount
of heat is released. The isosteric and limit heat of adsorption of
AC are 17.25–21.5 and 22.5 kJ·mol–1,
respectively.[35] The heat-transfer coefficient
of AC is 0.15–0.20 W·m–1·K–1.[35−37] As this value is low, the heat of adsorption cannot
be quickly transferred to the outside environment. The heat of adsorption
accumulates internally, resulting in a local temperature increase
in the fixed bed. The temperature cloud images of the fixed bed are
shown in Figure ,
and the local temperature isoline plots of the adsorption process
are shown in Figure .
Figure 3
Temperature cloud images of the adsorption process.
Figure 4
Local temperature isoline plots of the adsorption process.
Temperature cloud images of the adsorption process.Local temperature isoline plots of the adsorption process.Figure shows that
the local heating zone of the fixed bed moves forward along the airflow.
The maximum temperature in the high-temperature area of the six monitoring
points is 306–309.5 K. As soon as AC adsorbed and saturated,
it would no longer release heat. When the mixed gas kept on passing
through the adsorption saturation region, it would act as sweeping
gas and remove the heat gathered in the fixed bed, resulting in a
temperature decrease in the saturated region. For example, Figure a shows that the
adsorption of n-butane occurs in the area of point
1, and the temperature of point 1 reaches the maximum. From Figure a,b, we can find
that when the AC in the area of point 1 is saturated, the adsorption
moves forward to the area of point 2, and the temperature of point
2 increases gradually and reaches the maximum. As the continuous sweeping
of the gas transfers the locally increased heat to the fixed bed at
the back of the unsaturated part, the temperature of point 1 decreases
to the background temperature.In Figure , the
background temperatures of point 1 and point 6 were 300.5 and 301
K, respectively. The initial temperature of the fixed bed was set
as 300 K, and the temperature of the whole bed increased approximately
1 K. From Figure ,
we can see that the maximum temperatures of points 1–6 are
7.5, 6.5, 6, 8, 9, and 9.5 K respectively. The maximum temperatures
of the first three monitoring points decreased gradually, and those
of the latter three points increased. Because under the continuous
purging of the inlet gas, the heat generated by the first half of
the fixed bed is purged to the second half, and as a result, the fixed
bed’s background temperature is increased.Figures and 4 show that the closer to the center of the fixed
bed, the higher the temperature. However, the closer to the center
of the fixed bed, the more difficult it was to transfer heat to the
outside.The n-butane concentration cloud images
of the
fixed bed are shown in Figure , and the isotope plot of the n-butane concentration
of adsorption is shown in Figure .
Figure 5
n-Butane concentration cloud images of
the adsorption
process.
Figure 6
Isotope plot of the local n-butane concentration
of the adsorption process.
n-Butane concentration cloud images of
the adsorption
process.Isotope plot of the local n-butane concentration
of the adsorption process.The region between the saturated and nonadsorbed regions was denoted
as the adsorbed belt. Figures and 6 show that the adsorbed belt
moves forward gradually because, under the sweeping action of the
airflow, when the n-butanegas was adsorbed by AC,
the AC gradually became saturated and the adsorbed belt moved forward.
The adsorbed belt was curved, with its top at the center of the fixed
bed, and the direction was consistent with the direction of the airflow.
As the inlet was located in the middle of the fixed bed, the airflow
entering the AC area was faster near the center, which causes the
AC in the middle to become saturated earlier than that on both sides,
thus forming the arc-shaped adsorbed belt.Figures and 5 show that the adsorbed belt is consistent with
the higher temperature distribution of the fixed bed. The process
of organic gas adsorption by AC was mass transfer, and the process
of organic gas adsorption to release heat was heat transfer. This
indicates that the AC adsorption process is simultaneous mass and
heat transfer and reflects the synchronization and consistency of
the fluid mass- and heat-transfer processes.
Data
Validation and Analysis
The
temperature changes of the six monitoring points of the fixed bed
of AC are shown in Figure , and the ΔTmax of each monitoring point is given in Table . The temperature curves are parabola, and
the transverse span is basically the same, which shows that the length
of the mass-transfer belt is almost the same as that in the adsorption
process, the temperature changes obviously along the airflow direction
in the bed, the maximum temperature increase area is transmitted along
the air inlet direction, and the maximum temperature increase value
on the mass-transfer belt is between 7 and 10 K. The maximum temperatures
at the simulated and experimental process points 1–6 are shown
in Figure . Simulative
data of temperature at six monitoring points of the adsorption process
are shown in Figure ; the trend of the simulation data is consistent
with that of the experimental data, and the range of the temperature
change is small. This reflects the correctness of the simulation data
for the heat-transfer process of the fixed bed.
Figure 7
Experimental data of
temperature at six monitoring points of the
adsorption process.
Table 4
Maximum Temperature Increase Values
of Six Monitoring Points
number
point 1
point 2
point 3
point 4
point 5
point 6
ΔTmax (K)
7
9.2
10
9.2
8.5
9
Figure 8
Simulative data of temperature
at six monitoring points of the
adsorption process.
Experimental data of
temperature at six monitoring points of the
adsorption process.Simulative data of temperature
at six monitoring points of the
adsorption process.The deviation rates of
the experimental and simulated temperatures
at six points are shown in Table . The deviation rate between the simulated and experimental
temperatures was between 0.162 and 1.290%, which is well below 5%,
and the deviation rate between the simulated adsorption breakthrough
time 895 min and the experimental adsorption breakthrough time 924
min was 3.139%. This demonstrates the accuracy of the simulation results
obtained using Fluent.
Table 1
Deviation Rate between
the Simulated
and Experimental Temperatures
point 1
point 2
point 3
point 4
point 5
point 6
simulated data (K)
307.5
306.5
306
308
309
309.5
experimental
data (K)
307
309.2
310
309.2
308.5
309
deviation rate (%)
0.163
0.873
1.290
0.388
0.162
0.162
Figure shows the
changes in the concentration of n-butane at the six
monitoring points and the outlet. The trend of the n-butane concentration at the seven points was consistent, reflecting
the stability of the AC adsorption process. The trend was consistent
by comparing outlet concentrations obtained in the experiments and
simulated results, which reflects the reliability of the simulation
data. The simulated breakthrough time of adsorption was approximately
895 min, and the experimental breakthrough time was approximately
924 min. The breakthrough time of the simulated process was shorter
than that of the experimental process because, in the simulation process,
the particle size of AC was fairly uniform and the distribution of
AC was regular, which results in a higher adsorption efficiency and
shorter time required to simulate adsorption saturation.
Figure 9
Comparison
of the experimental and simulated concentration change
curves.
Comparison
of the experimental and simulated concentration change
curves.
Conclusions
In this study, a multidimensional transient model for the adsorption
process of an AC fixed bed was established, and the heat and mass
transfer of the n-butane adsorbed in the fixed bed
was simulated by Fluent. According to the experimental data and simulation
results, the trends of the changes in the temperature and concentration
in fixed-bed adsorption were analyzed. From the discussions in this
paper, we can draw the following conclusions:The simulation results were compared
with the experimental results of the absorption process. The deviation
ratio between the simulated and experimental temperatures at six points
ranged from 0.162 to 1.290%. The trend of the simulated outlet concentration
was consistent with that of the experimental outlet concentration,
and the accuracy and reliability of the simulation process and results
were demonstrated. The simulation results have value as a reference
for industrial applications of AC adsorption.Temperature cloud images show that
the temperature increase at the center of the fixed bed was larger
during the whole adsorption process. n-Butane concentration
cloud images show that the adsorbed belt was curved, and the AC at
the center of the fixed bed becomes saturated more rapidly than that
on both sides. In the design of the fixed bed, the heat release can
be promoted by adjusting the arrangement of AC particles and the packing
density in the middle of the fixed bed. The simulation results can
provide a reference for the design and modification of the fixed bed.
Experimental Section
The AC was supplied by Mead Westvaco
(MWV) Co. Ltd. (Shanghai,
China). Also, its properties and parameters are shown in Table . Experimental pretreatment:
Table 2
AC Properties and Parameters
AC type
AC quality mAC (g)
valid diameter Dp (mm)
porosity ε
density
ρAC (kg·m–3)
specific heat capacity cpAC (J·kg–1·K–1)
WV-A1500
217
0. 85
0.38
550
840
Equilibrium
Adsorption Capacity Experiment
The equilibrium adsorption
capacity (qe) of AC for n-butane was determined by conducting
flow adsorption equipment (Figure ) at room temperature and 1.0 atm. After preheating
and pretreatment, approximately 0.65 g of AC was added to a U-tube
with an inside diameter of 4 mm. The total gas flow was controlled
using mass flow meters to maintain a constant flow rate of 30 mL·min–1, and the concentration of n-butane
in the intake gas ranged from 250 ppm (10–6V/V)
to 2000 ppm (10–6V/V). The U-tube mass was weighed
at 25 min intervals until the value became stable.
Figure 10
Schematic diagram of
the VOC adsorption equipment on the U fixed
bed.
Schematic diagram of
the VOC adsorption equipment on the U fixed
bed.The equilibrium adsorption capacities
(qe, gVOCs/gadsorb) of the adsorbents were measured
and calculated by eq where WA and WB are the weights of the tube after and before
adsorption, respectively, (g), and WS is
the weight of samples before adsorption (g).
Breakthrough
Curve and Temperature Experiment
The breakthrough curve and
temperature in the fixed bed were tested
in the equipment (Figure ), and a 450 mm long 304 stainless steel pipe with an inside
diameter of 35 mm was used as the adsorption column.
Figure 11
Schematic diagram of
the VOC adsorption experiment system in the
insulated fixed bed.
Schematic diagram of
the VOC adsorption experiment system in the
insulated fixed bed.After 30 min of nitrogen
purging pretreatment, a stream mixture
with 20 000 ppm (10–6V/V) of n-butane/N2 (400 mL·min–1) was passed
through the fixed bed. The breakthrough curves were obtained by recording
the concentration of n-butane at the outlet of the
adsorption column, which was continuously monitored using a gas chromatograph
(GC) equipped with a total hydrocarbon-filling column (3 × 2000
mm) and an FID detector.Six thermocouples were used to monitor
the changes in the bed temperature
during adsorption. The thermocouple positions were defined with respect
to the inlet (z = 0) and tank axis (r = 0). Thermocouple P6 was located at the entrance of
the tank (z = 45 mm, r = 0); P5 at z = 120 mm, r = 0; P4 at z = 195 mm, r = 0; P3 at z = 270 mm, r = 0; P2 at z = 345 mm, r = 0; and
P1 was located at the outlet of the tank (z = 420 mm, r = 0).
Computational
Methods
Adsorption Capacity
To trace the
adsorption capacity changes in the adsorption process, isotherms of n-butane on AC were fitted by the Langmuir model[38,39] with the following formula (eq ).where qmax is
the maximum adsorption capacity (g·g–1), b is the Langmuir constant (L·mg–1), Ce is the mass concentration after
adsorption equilibrium (mg·L–1), and qe is the equilibrium adsorption capacity (g·g–1).The adsorption capacity of n-butane on AC could be fitted well by the Langmuir model, as indicated
by the good relative coefficients (R2 =
0.9871). According to the fitting equation in Figure , qmax = 0.2674
g·g–1 and b = 0.00299 L·mg–1. Therefore, the maximum adsorption capacity qmax of AC was 0.2674 g·g–1.
Figure 12
Fitting diagram of the adsorption capacity based on the Langmuir
equation.
Fitting diagram of the adsorption capacity based on the Langmuir
equation.
Adsorption
Kinetic Parameters
Yoon–Nelson Model
The Yoon–Nelson
model[40] is a relatively simple model for
the adsorption of gases on AC. This model assumes that the diminution
rate during adsorption for every adsorbate molecule is relative to
the probability of solute breakthrough on the adsorbent. The equation
for this model is described by the following formula (eq ).where CA and C0 are the outlet and inlet concentrations of
the stream gas (ppm), respectively, τ0 is the breakthrough
time (time for CA = 0.5 C0), and k′ is the Yoon–Nelson
constant.The fitting results of the Yoon–Nelson model
are shown in Figure and Table . According
to the adsorption breakthrough curve, the total amount of n-butane adsorbed by AC is 19.135 g, i.e., nC = 0.328 mol.
Figure 13
Breakthrough curves
and mathematic models of AC samples.
Table 3
Fitting Kinetic Parameters of Adsorption
of n-Butane by AC with Different Models
Yoon–Nelson adsorption
model
first-order
dynamic
equation
temperature (K)
k′ (min–1)
τ0 (min)
R2
k1 (min–1)
R2
300
0.0635
924
0.992
0.02386
0.953
Breakthrough curves
and mathematic models of AC samples.
First-Order Kinetics
Equation Model
The first-order kinetics equation (eq ) was used to fit the
adsorption breakthrough curve,
and the results are shown in Figure and Table .Fitting curve of first-order
adsorption kinetics equation.The first-order kinetics equation[41] is
described by the following formula (eq ).where t is the time (min), qe is the equilibrium adsorption capacity (mg·g–1), q is the adsorption capacity of n-butane at time t (mg·g–1), and k1 is the first-order adsorption rate constant (min–1).According to the fitting equation of the
Yoon–Nelson adsorption
model, the adsorption breakthrough time is 924 min. It can be seen
from the correlation coefficient R2 that
the Yoon–Nelson adsorption model and first-order kinetics equation
have high goodness of fit with the adsorption breakthrough curve,
indicating that the adsorption process of n-butane
on AC is physical adsorption.
Adsorption
Heat
In the adsorption
process, the gas molecules move to the solid surface, the molecular
movement speed is greatly reduced, and the heat is released. The adsorption
heat is one of the important indicators to measure the performance
of the adsorbent.[42]The maximum temperature
increase values (ΔTmax) of the six
monitoring points in the adsorption experiment are shown in Table .When
the adsorption of n-butane by AC particles
reaches the saturation state, the heat is no longer released, and
the adsorbed n-butane molecules and free n-butane molecules in the AC are in a thermal dynamic equilibrium
state. The distribution of n-butane obeys Boltzmann
distribution, and the adsorption heat of two-energy-state model εa[43,44] formula is (eq )where k is the
Boltzman number
(1.380649 × 10–23 J·K–1), T is the temperature (K), N1 is the molecular number of free n-butane,
and N2 is the molecular number of adsorbed n-butane.The total adsorption heat Qa isUnder the standard condition,
90% of the adsorbed n-butane is in the adsorbed state
and 10% is in the free state.[43] The ambient
temperature of this adsorption experiment
is 300 K, and the ratio of the molecular number of adsorbed and free n-butane is N1/N2 = 9:1. According to the adsorption heat formula of two-energy-state
model (eq ),[32,33] the theoretical adsorption heat of AC adsorbing n-butane (Qa,) can be
calculatedwhere NA is the
Avogadro constant (6.02 × 1023).The theoretical
adsorption heat of n-butane on
AC is 5.48 kJ·mol–1.The mass of n-butane adsorbed by AC is nC = 0.328 mol, and
the maximum temperature increase is ΔTmax = 10 K (Table ). The actual adsorption heat of AC adsorbing n-butane (Qa,a) in the experiment can
be calculatedThe deviation rate between Qa, and Qa,a is 1.44%; the
actual
adsorption heat is consistent with the theoretical adsorption heat.
The results show that the adsorption of n-butane
by AC belongs to physical adsorption. The adsorption heat of n-butane
on AC is Qa,a = 5.56 kJ·mol–1.
Authors: José L Domingo; Joaquim Rovira; Lolita Vilavert; Martí Nadal; María J Figueras; Marta Schuhmacher Journal: Sci Total Environ Date: 2015-03-12 Impact factor: 7.963