The absorption of carbon dioxide is an important process in many practical applications such as reduction of greenhouse gases, separation and purification processes in the chemical and petroleum industries, and capture of radioactive isotopes in the nuclear fuel cycle The goal of this research is to develop a dynamic model to simulate CO2 absorption by using different alkanolamines as absorption solvents. The model is based upon transient mass and energy balances for the chemical species commonly present in CO2 gas-liquid absorption. A computer code has been written to implement the proposed model. Simulation results are discussed. The reported model simulates well the response to dynamic changes in input conditions. The proposed model can be used to optimize and control the separation of carbon-14 in the form of CO2 in the nuclear industry.

1. Theoretical Considerations

1.1. Introduction

The off-gas stream produced from reprocessing used nuclear fuel containing a host of radioactive gases including tritium, iodine, methyl iodide, krypton, xenon, and carbon dioxide. These gases are hazardous to human health, and some are of considerable economic value. Currently, there is interest in the efficient capture of these gases for recovery or storage purposes [1].

Presently, postcombustion separation of carbon dioxide from flue-gas streams is investigated using membrane separation, low-temperature distillation, physical adsorption by activated carbon and zeolites, chemical absorption by mineral oxides to produce the corresponding carbonates, and chemical absorption by suitable solvents [2]. Chemical absorption by solvents is the most common form of the CO2 capture process. Amines, and mixtures containing amines, are the most used solvents as these compounds react quickly with CO2, maintaining high mass transfer-driving forces [2]. Industrially important amines are monoethanolamine (MEA), diethanolamine (DEA), diisopropanolamine (DIPA), and the tertiary amine N-methyldiethanolamine (MDEA). Most industrial processes are operated with aqueous amine solutions, but solvents consisting of a mixture of water and a nonaqueous solvent, for example, sulfolane in the Shell-Sulfinol process [3], are also frequently used [4]. Blending of different amines is considered to be attractive because in this way the high absorption capacity of tertiary amines can be combined with the high absorption rates of primary or secondary amines [5]. Blends are also more flexible than singular amines because the relative concentration of the amines can be varied [6]. Other combinations have been considered, such as MEA and 2-amino-2-methyl-1-propanol (AMP). The latter amine is an attractive alternative since it offers higher absorption capacity and lower regeneration energy. Blending MEA with AMP is considered to combine all favorable characteristics of both amines and overcome the unfavorable characteristics [7]. Some studies have been carried out considering the influence of the solvent mixed with the alkanolamines. For example, Usubharatana and Tontiwachwuthikul [8] studied the kinetics of CO2 capture using methanol mixed into solutions of MEA.

The goal of this work is to develop a theoretical model for dynamic simulations of CO2 absorption using a generic blend of amines and/or other solvents. The proposed model can be used to simulate chemical absorption using a wide range of CO2 input stream concentrations. Through computational implementation, the derived model can be used for design, optimization, and control applications.

1.2. Model Organization

In this work, we will follow an approach similar to the one used by Greer et al. [9] in the dynamic simulation of the absorption/desorption of carbon dioxide from monoethanolamine (MEA). A dynamic model of the absorption process will be developed using a generic mixture of solvents to increase the absorption rate.

Mass transfer rates are calculated using the thin film model with chemical reactions confined to the liquid thin film. Henry’s law is used for modeling the vapor-phase equilibrium of CO2. Henry’s constant values for the different gas-phase species are calculated using fugacity ratios obtained by the Peng–Robinson equation of state (EOS). A kinetic scheme is proposed to represent the complex chemical reactions between CO2 and the solvents. Liquid and vapor energy balances are developed to calculate the liquid and vapor temperature, respectively. A schematic of the column is shown in Figure 1. The model of the absorption tower is developed by taking a small slice of the tower of height dz. The height of the packing goes from z = 0 to z = H, where H is the packing height.

1.3. Kinetic Studies

The most important part of the proposed model is the simulation of the complex chemical reactions involving the solvents and CO2. In the case of primary and secondary amines, the reaction mechanism is well understood as was originally proposed by Caplow [10] and reintroduced by Danckwerts [11]. CO2 reacts with the amine through a two-step process. The first step proceeds through the formation of a zwitterion intermediate:

This step is slow and considered to be the rate-controlling step; it is followed by a very fast removal of a proton by a base:

It is important to notice that all bases present in the liquid phase will participate in reaction (2); therefore, Bi represents a generic base participating in the reaction with the zwitterion. In this mechanism, the overall forward reaction rate equation can be derived using the quasi-steady-state assumption for the zwitterion intermediate [12]:where k1, k−1, and kbi are the respective specific constants and represents the summation over the reaction rates of all the bases reacting with the zwitterion present in the solution. In dilute aqueous solutions, for example, the amine, OH, and water act as a base, while in nonaqueous solvents, only the amine can be considered a base [12].

In the case of high amine concentration in the solvent, Equation (3) is simplified to

For aqueous MEA solutions, the overall reaction rate is of second order, and is of first order with respect to the amine. This finding indicates that the deprotonation of the zwitterion by the bases present in the solution is very fast compared to the reverse reaction. Therefore, Equation (4) is considered the main reaction in the absorption of CO2 in high-concentration aqueous alkanolamine solutions. In the case of low-concentration solutions, more complex rate equations hold [2, 13]. Alvarez-Fuster et al. [14] and Sada et al. [15] showed that changes in the solvent lead also to changes in the order of reaction.

A different mechanism applies to the reaction of CO2 with ternary amines. According to Littel et al. [16], the reaction of CO2 with tertiary amines can be described satisfactorily using the base-catalysis reaction mechanism proposed by Donaldson and Nguyen [17]:

This mechanism is essentially a base-catalyzed hydration of CO2; thus, tertiary amines cannot react directly with CO2. This finding was confirmed by Versteeg and Van Swaaij [18], who studied the absorption of CO2 into a solution of MDEA and ethanol without water.

In all the cases discussed above, the following CO2 reactions are also present:

Reaction (6) is very slow and can be neglected in most circumstances. Reaction (7), however, is fast and can enhance mass transfer even when the concentration of the hydroxyl ion is low [19].

1.4. Reaction Scheme
1.4.1. Introduction

The components in the gas phase are RNH2, R1R2R3N, CO2, H2O, N2, and O2, while the species considered in the liquid phase are , RNCOO, RNH2, R1R2R3N, R1R2R3NH+, CO2, H2O, N2, O2, , OH, and H3O+. The ionic species remain in the liquid phase, while the others are transferred from one phase to the other according to the scheme shown in Figure 2.

Only gas-phase mass transfer resistance is considered for liquid solvents (H2O, R1R2R3N, and RNH2), while liquid-only mass transfer resistance for the gas species has been considered: CO2, N2, and O2 ([2], among others). Under certain conditions, however, both resistances have to be considered for the reacting species CO2. In our computer code, we consider these situations by using a global mass transfer coefficient given bywhere is the overall mass transfer coefficient based upon liquid-phase concentrations; is the liquid-phase mass transfer coefficient; is the gas-phase mass transfer coefficient; is the concentration-based Henry’s constant, in which and are the liquid- and gas-phase i-species concentrations, respectively; and E is the enhancement factor defined in the subsection below. Typical values of these coefficients are given in Table 1. In Table 2, all the chemical species considered in the model are included.

1.4.2. Reactions

Mandal et al. [20] and Benamor and Aroua [21] proposed the following set of reactions occurring in the aqueous primary amine solution. A gas-liquid equilibrium equation and six chemical equilibria equations are introduced to describe the chemistry involved in CO2 absorption:

Following Bosch et al. [6] who studied blends of alkanol amines, we add reaction (5) catalyzed by ternary amines and reaction (16) to account for the decomposition of the ternary ammonium ion:

Reactions (5), (9), and (11) are the main reactions in CO2 absorption by amines. The other reactions are required to complete the total reaction mechanism [6]. The CO2 loading rate, defined as the ratio of CO2 to alkanolamines, determines the relative weight of all reactions. For CO2 loading rates below 0.5, Equation (9) is the main reaction. In the case of loading rates above 0.5, reaction (11) predominates, while reaction (5) will be important only for large amounts of ternary amine in the feed. In this work, we are interested in small loading rates and significant amounts of ternary amine present; therefore, the reactions with the amines will be the most important ones. A list of all reactions used in the model is included in Appendix. Equation (15) is considered by using the corresponding Henry’s law constant for CO2.

1.4.3. Kinetic Parameters

Kinetic data were collected using a general primary-secondary amine that follows the zwitterion mechanism represented by Caplow [10] and Danckwerts [11]. The forward reaction (9) is thought to occur through a two-step mechanism. Initially, a CO2 molecule and an alkanolamine molecule form a zwitterion intermediate which, in a second step, reacts with another alkanolamine molecule. The second step is much faster than the first step; hence, the first step is rate limiting and of second order. The reaction rate is given by

The specific forward rate constant k1f is calculated using the method in [22]:

The bicarbonate ion formation, reaction (11), is the most important reaction for CO2/alkanolamine loadings above 0.5. It is of second order, given by

The forward rate for the formation of bicarbonate is significantly fast, but the overall rate is usually quite small due to the low concentration of OH ions used. At loadings of CO2/alkanolamines above 0.5, this becomes the dominant reaction for CO2 removal. The forward rate is calculated from [23]:

In order to complete the rate calculations, we collected literature data corresponding to the equilibrium rate constants of reactions (9)–(14) from Bedelbayiev et al. [24] and Greer [2]. In order to deal with reactions (5) and (16) for ternary alkanol amines, we followed many investigators including Haimour et al. [25], Critchfield [26], Littel et al. [16], and Rangwala et al. [27] who fit the rate constant of the reaction as a function of temperature using

All values of kinetic parameters used are listed in Table 3. A full description of all the rate equations and calculation of generation terms is also included in Appendix.

1.4.4. Enhancement Factor

The carbon dioxide absorption is accompanied by strong chemical reactions. Therefore, the calculation of the CO2 flux term requires the use of an enhancement factor (E) to account for the enhanced mass transfer. The enhancement factor is defined as the mass transfer rate under reactive absorption divided by the mass transfer rate under nonreactive absorption conditions [31]. The CO2 molar flow term () is given as follows [2, 9, 24]:where aw is the interphase area per unit volume. The enhancement factor (E) is a function of the Hatta number (Ha) defined aswhere is the diffusion coefficient and kif are the forward specific reaction rate constants. The influence of the different reactions on the total rate of CO2 absorption is considered by the enhancement factor E. The Hatta number is the ratio of the rate of homogeneous reaction relative to the rate of gas dissolution. Ha is also a measure of the amount of dissolved gas that reacts inside the diffusion film near the gas-liquid interface compared to the amount that reaches the bulk of the solution without reacting. When Ha = 0, we have purely physical absorption. The higher the value of the Hatta number, the stronger the effect of chemical reactions on mass transfer. In the case of Ha > 2, the enhancement factor E is directly equal to Ha [31].

Versteeg et al. [32] recommended a different definition of the Hatta number including instead of CNaOH in Equation (23) to calculate the enhancement factor when all RNH2 is consumed, and reaction (3) is the dominant CO2 removal reaction. Other authors have used both formulations with similar results [24, 33, 34]. In this work, however, we preferred to use the general definition of the Hatta number, Equation (23).

1.5. Mass Transfer Model
1.5.1. Mass Balances

The mass balance of component i in the liquid phase was calculated using the following equation [2, 35]where is the superficial liquid-phase velocity, represents moles of species i generated/consumed by interphase reaction per unit volume, and is the mass flow of component i from the liquid phase into the gas phase. In the model presented in this work, the effect of reaction on the CO2 absorption is considered through the use of the enhancement factor. In the cases of ionic species, there is no interphase mass transfer; therefore, for these reactions, Equation (24) becomes

The generic amines (RNH2 and R1R2R3N) are the only chemical species for which we have to calculate a generation term plus an interfacial mass transfer term ( and ).

1.6. Energy Balances

The reactions given by Equations (8) and (11) are highly exothermic; therefore, an energy balance has to be solved in order to consider temperature changes. The heat of reactions for Equations (9) and (11) used was 65 kJ/molCO2 and 20 kJ/molCO2, respectively. The first value was taken from Greer [2] and the second from Pinsent et al. [28].

A two-equation model for the transient energy balance in the control volume depicted in Figure 1 leads to the following equations for all the components shown in the figure [9, 36]:where and are the heat capacities of component i in the mixture, is the global heat transfer coefficient, ug is the superficial velocity inside the gas phase, and is the heat released by the chemical reaction. The CO2 molar flow term () is given by Equation (22), while the enhancement factor (E) is given by the Hatta number defined by Equation (23).

1.7. Thermodynamics

The molar flow of component i from the gas phase into the liquid phase () is calculated usingwhere is the interfacial liquid equilibrium concentration. In order to evaluate , we use a formulation based upon the calculation of gas and liquid fugacity values [2]. This formulation avoids the use of an iterative procedure as the fugacity values can be calculated directly as

This expression for the diffusion molar flow is valid for the RNH2, R1R2R3N, and H2O components when the resistance is assumed to be in the gas-liquid film [2]. A similar expression for CO2, O2, and N2 can also be derived by

In Equations (29) and (30), R is the ideal gas constant; P is the pressure; T is the temperature; is the total molar concentration in the liquid phase; Zg is the gas-phase compressibility factor; , are the fugacities of component i in liquid and gas phases, respectively; and , are the fugacity coefficients for component i in the liquid and gas phases, respectively.

2. Results and Discussion

2.1. Model Validation

A FORTRAN computer code was developed based upon a previous one prepared for simulating CO2 absorption using primary/binary amine solvents [30]. An explicit finite-difference scheme was used to solve the relevant hyperbolic partial differential equations. For convenience, in our equations, all variables were made dimensionless, but the species concentrations were calculated in dimensional values (mM and mol/m3). The residence time of the liquid phase (tr = H/ul) was used to define a dimensionless time (t = time/tr). In all our simulation runs, we used geometric parameters and operating conditions taken from the literature. A list of our input data is shown in Table 4. This list includes information about the packing, the ranges of all values, and the standard set used in most of our calculations.

Validation of the computer code was achieved by comparing calculated parameters with data from the literature, especially with the work of Greer [2]. The validation process included checking against literature values [2, 5, 6, 23, 26, 39] the values of our calculated geometric parameters (), transfer coefficients (), kinetic parameters (), and equilibrium parameters (). It also included a comparison of calculated mass species concentration profiles against literature values [2, 5, 23, 26, 30].

In Figure 3, we show a comparison of calculated CO2 axial concentration profiles against those values reported by Greer [2] for dimensionless time (t) equal to 0.328 (2000 s) and steady-state condition. The basic set of parameters and operating conditions reported by Greer ([2] Table 4.2 in page 90) was used in these calculations. It can be seen in the figure that there is good agreement between our calculated data and those from Greer [2] for the CO2 axial concentration profiles; our calculated values, however, were always slightly lower than those reported by Greer [2]. Similar agreement, not shown in this work, was achieved when calculating the other chemical compounds’ axial profiles.

The accuracy of the reaction model was validated by carrying out a mass balance for the amine chemical species. A global mass balance for the amines and their reaction products at equilibrium gives

The difference between both sides in Equations (31) and (32) was used as a way of estimating the accuracy of the reaction scheme. In all our calculations, the error between the calculated concentrations was on the order of the code precision (10−6). A global mass balance for CO2 was used as another way of estimating the global accuracy of the proposed model. The biggest relative error in all our calculations was 0.01%.

2.2. Simulation Results

The time change of the CO2 concentration is depicted in Figure 4, where typical axial concentration profiles are presented. At dimensionless time equal to 0, a gas mixture containing CO2 is injected. The CO2 concentration is highest at the base of the column, Z = 0, and drops as we approach the top, Z = 1. The results in Figure 4 show that as time increases, the concentration of CO2 increases as we approach the top. We can also see that, at long times, the steady state is achieved.

In order to study the time evolution of the concentration of the generic amine, we chose to use monoethanolamine (MEA) due to the high amount of data available for this particular chemical compound. MEA reacts with CO2 following reaction (9). Two amine ions, the ammonium-like ion and the carbamate ion , are produced by this reaction and consumed by reactions (13) and (14), respectively. The interaction between these ions and MEA is very complex and varies with time. In Figures 5 and 6, we show contour plots which describe the variations of the ionic species axial profiles with time. In Figures 5 and 6, the depth axis depicts the axial variation of the concentration from the bottom to the top, i.e., it is a view from the bottom. The horizontal axis shows the time change of the concentration profile. Low values are represented by blue colors, while high values are represented by red colors. The “back” wall represents ionic input values, while the left-side wall represents initial conditions.

The cation profile is shown in Figure 5. The figure shows that the values of the cation concentration are initially very high throughout the column and decrease as time increases. The cation concentration is equal to zero at the column top (Z = 1), passes through a maximum, and decreases as it moves to the bottom (Z = 0). As time increases, the concentration of the cation decreases at the bottom of the column. Our results show that the decrease in the cation concentration leads to an increase in both the anion concentration and the neutral solvent. This behavior is produced by the complex interaction among the different chemical reactions.

Figure 6 depicts the time evolution of the carbamate ion concentration. The anion is not present in the input solution at the top (Z = 1); its concentration increases as CO2 moves upwards in the gas phase, until a maximum is reached at the bottom, coinciding with the highest CO2 concentration at the bottom of the column (Z = 0). The concentration values of the ammonium-like ion were also higher than the carbamate concentration values in all simulations.

The time variation of the solvent MEA concentration is shown in Figure 7.

It is shown in Figure 7 that as the amount of the reaction products increases, we can expect that the amount of the free amine will decrease. An analysis of the three contour plots shows that, at short times, there is a high rate of generation of the cation by consumption of the MEA solvent. As time increases, the concentration of the cation decreases, while the concentration of the anion increases, and some increase of the solvent concentration is observed until equilibrium is achieved.

The main goal of this project is to develop a dynamic simulator for nuclear energy applications. The dynamic simulator is evaluated by how it describes sudden changes in input parameters. In order to test our model, we simulated the response to two sudden changes in input CO2 concentrations: (i) a step which permanently increases the CO2 input concentration and (ii) a discrete pulse which increases the input concentration during a fixed time; then, the input concentration is set to its original value (Figures 8 and 9). In the case of a step change in input concentration, we assumed that the column reaches the steady-state operating conditions shown in Figure 4; this is the starting time for the simulation (t = 0). At t = 0, we double the value of the CO2 input concentration. The results are shown in Figure 8 where the CO2 gas-phase concentration starts increasing in value from the bottom of the column (input) towards the top (output) until a new equilibrium is achieved at higher concentration values (t > 0.1). There is a continuous increase in CO2 concentration until the new steady-state values are reached. The area between both curves is proportional to the amount of extra mass added.

The simulation results for a discrete pulse are shown in Figure 9. We start from the final equilibrium concentration profile shown in Figure 8 (t = 0). At t = 0, the CO2 input concentration is set to its original value (1.7 mol/m3), and the CO2 concentration starts decreasing throughout the column. This behavior is more noticeable at the bottom of the column, as the extra amount of CO2 injected during the pulse travels through the column. The solid line in Figure 9 represents the final steady value to be reached, equal to the one depicted in Figure 4 for the same CO2 input concentration. At this point in time, all the extra mass injected by the pulse has been absorbed and removed from the gas phase. The area in between the steady-state curve and the CO2 concentration profile at a particular time represents the remaining extra CO2 mass injected during the pulse.

The results presented in Figures 8 and 9 show that the proposed model can be successfully used to analyze the response of the absorption column to sudden changes in CO2 input concentrations.

3. Conclusions

A model for the dynamic behavior of reactive CO2 absorption using mixtures of alkanolamine solvents has been successfully developed. The model is based upon transient mass and energy balances for several different chemical species commonly present in CO2 gas-liquid absorption. The phase equilibrium has been considered using a thermodynamic model and through the use of experimentally based Henry’s law values. Typical values for the geometric parameters of the absorber and the packing characteristics have been collected. A reaction scheme that takes into account the different reactions between CO2 and blends of amines in an alkaline environment has been proposed. A computer code has been written to implement the proposed model. The computer code has been validated by checking the values of parameters calculated and comparing results to those reported in the literature. The mass balances for CO2 have been close within a 0.01% relative error, while the alkanolamine solvents' mass balances have been closed within the computer code precision (10−6). The results have been collected, and they are logical and agree with equivalent literature results. The computer code developed in this work can describe adequately the dynamic processes occurring due to sudden changes in operating conditions. Therefore, it is a valuable tool to design, optimize, and control absorption processes in the nuclear industries.


A. Overall CO2 Mass Transfer Coefficient

The computer code developed to implement our model calculates an overall mass transfer coefficient based upon liquid-phase concentrations () and compares its value with the liquid-layer mass transfer coefficient () value to decide which one should be used at each step of the calculations. Typical results are presented in Table 2. Equation (8) was used to calculate . We can see that, in most of the calculations presented in this work, the relative error between using E and is around 10%; therefore, the overall mass transfer coefficient is used.

B. Reaction List

The 8 reactions used in the kinetic model are listed below:

C. Reaction Rates

The equations used to calculate the forward and reverse reactions areHere, is the water molar fraction in the feed solvent phase.

D. Kinetic Data

Table 3 summarizes all the rate equations used in the model. The first number in the first column is the equation number used in Appendix. The second number refers to the original equation number used in the article.

E. Simulation Parameters

In our simulation runs, we used a Montz B 200 metal structured packing [2] as typical packing. Table 4 shows typical geometric parameters and operating variables used in the simulation runs.

F. Generation Terms

Assuming the pseudo-steady state for every chemical species, we can calculate the generation terms that enter into the chemical species mass balances. In order to simplify the calculations, we determined an overall rate per reaction according to

Every generation term (Rgeni) is calculated by a molar balance using


aT:Mass transfer specific area (m2·m−3)
aw:Packing specific area (m2·m−3)
Bi:Concentration of base i participating in reaction (2) (mol·m−3)
:Interfacial liquid equilibrium concentration (mol·m−3)
:Concentration of species i in phase j (mol·m−3)
:Heat capacities of component i in phase j (J·mol−1 K−1)
Di:Diffusivity of component i (m2·s−1)
E:Enhancement factor (dim.)
Ea:Activation energy (J·mol−1)
:Fugacities of component i in phase j (dim. or Pa)
Hcc:Concentration-based Henry’s constant (dim.)
Ha:Hatta number (dim.)
kbi:Specific reaction rate constant of base i (m3·mol−1·s−1)
ki:Specific reaction rate constant (m3·mol·s−1 or s−1)
:Preexponential factor reaction i (m3·mol·s−1 or s−1)
kif:Forward reaction specific rate constant (m3·mol·s−1 or s−1)
kir:Reverse reaction specific rate constant (m3·mol·s−1 or s−1)
:Phase j mass transfer coefficient of species i (m/s)
Ki:Equilibrium constant reaction i (dim. or mol m−3)
:Overall liquid mass transfer coefficient of species i (m/s)
:Flow of component i from the liquid phase into the gas phase (mol m−3 s−1)
:Flow of component i from the gas phase into the liquid phase (mol·m−3·s−1)
:Species i generated/consumed by chemical reaction (mol·m−3·s−1)
ri:Reaction rate for species i (mol·m−3·s−1)
R:Gas constant (J·mol−1·K−1)
t:Time (s)
Ti:Temperature of the i phase (K)
ui:Superficial velocity of the i phase (m·s−1)
UT:Global heat transfer coefficient (J·m−2·K−1·s−1)
z:Column height (dim.)
Zg:Compressibility factor (dim.)
ΔHR:Heat released by chemical reaction (J·mol−1)
:Fugacity coefficients for component i in phase j (dim.)

Data Availability

The data collected during the submitted research project are stored in the authors’ servers. The data and metadata collected are stored in Excel spreadsheets and comprise less than 1 GB. They can be made public as required.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


Funding for this research provided by the Office of Nuclear Energy of the U.S. Department of Energy, under the Nuclear Energy University Program (Grant # NFE-12-03822), is gratefully acknowledged by the authors. This study was conducted at Prairie View A&M University in collaboration with the Oak Ridge National Laboratory (ORNL). ORNL is managed by UT-Battelle, LLC, under Contract DE-AC05-00OR22725 with the U.S. Department of Energy.