Research Article  Open Access
Carbon Dioxide Absorption Modeling for OffGas Treatment in the Nuclear Fuel Cycle
Abstract
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 CO_{2} 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 CO_{2} gasliquid 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 carbon14 in the form of CO_{2} in the nuclear industry.
1. Theoretical Considerations
1.1. Introduction
The offgas 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 fluegas streams is investigated using membrane separation, lowtemperature 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 CO_{2} capture process. Amines, and mixtures containing amines, are the most used solvents as these compounds react quickly with CO_{2}, maintaining high mass transferdriving forces [2]. Industrially important amines are monoethanolamine (MEA), diethanolamine (DEA), diisopropanolamine (DIPA), and the tertiary amine Nmethyldiethanolamine (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 ShellSulfinol 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 2amino2methyl1propanol (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 CO_{2} capture using methanol mixed into solutions of MEA.
The goal of this work is to develop a theoretical model for dynamic simulations of CO_{2} 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 CO_{2} 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 vaporphase equilibrium of CO_{2}. Henry’s constant values for the different gasphase 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 CO_{2} 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 CO_{2}. 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]. CO_{2} reacts with the amine through a twostep process. The first step proceeds through the formation of a zwitterion intermediate:
This step is slow and considered to be the ratecontrolling 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, B_{i} 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 quasisteadystate assumption for the zwitterion intermediate [12]:where k_{1}, k_{−1}, and k_{bi} 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 CO_{2} in highconcentration aqueous alkanolamine solutions. In the case of lowconcentration solutions, more complex rate equations hold [2, 13]. AlvarezFuster 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 CO_{2} with ternary amines. According to Littel et al. [16], the reaction of CO_{2} with tertiary amines can be described satisfactorily using the basecatalysis reaction mechanism proposed by Donaldson and Nguyen [17]:
This mechanism is essentially a basecatalyzed hydration of CO_{2}; thus, tertiary amines cannot react directly with CO_{2}. This finding was confirmed by Versteeg and Van Swaaij [18], who studied the absorption of CO_{2} into a solution of MDEA and ethanol without water.
In all the cases discussed above, the following CO_{2} 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 RNH_{2}, R_{1}R_{2}R_{3}N, CO_{2}, H_{2}O, N_{2}, and O_{2}, while the species considered in the liquid phase are , RNCOO^{−}, RNH_{2}, R_{1}R_{2}R_{3}N, R_{1}R_{2}R_{3}NH^{+}, CO_{2}, H_{2}O, N_{2}, O_{2}, , OH^{−}, and H_{3}O^{+}. 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 gasphase mass transfer resistance is considered for liquid solvents (H_{2}O, R_{1}R_{2}R_{3}N, and RNH_{2}), while liquidonly mass transfer resistance for the gas species has been considered: CO_{2}, N_{2}, and O_{2} ([2], among others). Under certain conditions, however, both resistances have to be considered for the reacting species CO_{2}. 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 liquidphase concentrations; is the liquidphase mass transfer coefficient; is the gasphase mass transfer coefficient; is the concentrationbased Henry’s constant, in which and are the liquid and gasphase ispecies 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 gasliquid equilibrium equation and six chemical equilibria equations are introduced to describe the chemistry involved in CO_{2} 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 CO_{2} absorption by amines. The other reactions are required to complete the total reaction mechanism [6]. The CO_{2} loading rate, defined as the ratio of CO_{2} to alkanolamines, determines the relative weight of all reactions. For CO_{2} 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 CO_{2}.
1.4.3. Kinetic Parameters
Kinetic data were collected using a general primarysecondary amine that follows the zwitterion mechanism represented by Caplow [10] and Danckwerts [11]. The forward reaction (9) is thought to occur through a twostep mechanism. Initially, a CO_{2} 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 k_{1f} is calculated using the method in [22]:
The bicarbonate ion formation, reaction (11), is the most important reaction for CO_{2}/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 CO_{2}/alkanolamines above 0.5, this becomes the dominant reaction for CO_{2} 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 CO_{2} 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 CO_{2} molar flow term () is given as follows [2, 9, 24]:where a_{w} 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 k_{if} are the forward specific reaction rate constants. The influence of the different reactions on the total rate of CO_{2} 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 gasliquid 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 C_{NaOH} in Equation (23) to calculate the enhancement factor when all RNH_{2} is consumed, and reaction (3) is the dominant CO_{2} 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 liquidphase 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 CO_{2} 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 (RNH_{2} and R_{1}R_{2}R_{3}N) 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/molCO_{2} and 20 kJ/molCO_{2}, respectively. The first value was taken from Greer [2] and the second from Pinsent et al. [28].
A twoequation 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, u_{g} is the superficial velocity inside the gas phase, and is the heat released by the chemical reaction. The CO_{2} 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 RNH_{2}, R_{1}R_{2}R_{3}N, and H_{2}O components when the resistance is assumed to be in the gasliquid film [2]. A similar expression for CO_{2}, O_{2}, and N_{2} 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; Z_{g} is the gasphase 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 CO_{2} absorption using primary/binary amine solvents [30]. An explicit finitedifference 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/m^{3}). The residence time of the liquid phase (t_{r} = H/u_{l}) was used to define a dimensionless time (t = time/t_{r}). 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 CO_{2} axial concentration profiles against those values reported by Greer [2] for dimensionless time (t) equal to 0.328 (2000 s) and steadystate 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 CO_{2} 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 CO_{2} 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 CO_{2} concentration is depicted in Figure 4, where typical axial concentration profiles are presented. At dimensionless time equal to 0, a gas mixture containing CO_{2} is injected. The CO_{2} 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 CO_{2} 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 CO_{2} following reaction (9). Two amine ions, the ammoniumlike 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 leftside 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 CO_{2} moves upwards in the gas phase, until a maximum is reached at the bottom, coinciding with the highest CO_{2} concentration at the bottom of the column (Z = 0). The concentration values of the ammoniumlike 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 CO_{2} concentrations: (i) a step which permanently increases the CO_{2} 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 steadystate 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 CO_{2} input concentration. The results are shown in Figure 8 where the CO_{2} gasphase 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 CO_{2} concentration until the new steadystate 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 CO_{2} input concentration is set to its original value (1.7 mol/m^{3}), and the CO_{2} concentration starts decreasing throughout the column. This behavior is more noticeable at the bottom of the column, as the extra amount of CO_{2} 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 CO_{2} 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 steadystate curve and the CO_{2} concentration profile at a particular time represents the remaining extra CO_{2} 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 CO_{2} input concentrations.
3. Conclusions
A model for the dynamic behavior of reactive CO_{2} 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 CO_{2} gasliquid 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 CO_{2} 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 CO_{2} 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.
Appendix
A. Overall CO_{2} Mass Transfer Coefficient
The computer code developed to implement our model calculates an overall mass transfer coefficient based upon liquidphase concentrations () and compares its value with the liquidlayer 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 pseudosteady 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 (R_{geni}) is calculated by a molar balance using
Nomenclature
a_{T}:  Mass transfer specific area (m^{2}·m^{−3}) 
a_{w}:  Packing specific area (m^{2}·m^{−3}) 
B_{i}:  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}) 
D_{i}:  Diffusivity of component i (m^{2}·s^{−1}) 
E:  Enhancement factor (dim.) 
E_{a}:  Activation energy (J·mol^{−1}) 
:  Fugacities of component i in phase j (dim. or Pa) 
H^{cc}:  Concentrationbased Henry’s constant (dim.) 
Ha:  Hatta number (dim.) 
k_{bi}:  Specific reaction rate constant of base i (m^{3}·mol^{−1}·s^{−1}) 
k_{i}:  Specific reaction rate constant (m^{3}·mol·s^{−1} or s^{−1}) 
:  Preexponential factor reaction i (m^{3}·mol·s^{−1} or s^{−1}) 
k_{if}:  Forward reaction specific rate constant (m^{3}·mol·s^{−1} or s^{−1}) 
k_{ir}:  Reverse reaction specific rate constant (m^{3}·mol·s^{−1} or s^{−1}) 
:  Phase j mass transfer coefficient of species i (m/s) 
K_{i}:  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}) 
r_{i}:  Reaction rate for species i (mol·m^{−3}·s^{−1}) 
R:  Gas constant (J·mol^{−1}·K^{−1}) 
t:  Time (s) 
T_{i}:  Temperature of the i phase (K) 
u_{i}:  Superficial velocity of the i phase (m·s^{−1}) 
U_{T}:  Global heat transfer coefficient (J·m^{−2}·K^{−1}·s^{−1}) 
z:  Column height (dim.) 
Z_{g}:  Compressibility factor (dim.) 
ΔH_{R}:  Heat released by chemical reaction (J·mol^{−1}) 
:  Fugacity coefficients for component i in phase j (dim.) 
Dim.:  Dimensionless. 
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.
Acknowledgments
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 # NFE1203822), 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 UTBattelle, LLC, under Contract DEAC0500OR22725 with the U.S. Department of Energy.
References
 A. Ladshaw, S. Yiacoumi, R. Lin, Y. Nan, L. L. Tavlarides, and C. Tsouris, “A mechanistic modeling framework for gas phase adsorption kinetics and fixedbed transport,” AIChE Journal, vol. 63, no. 11, pp. 5029–5043, 2017. View at: Publisher Site  Google Scholar
 T. Greer, “Modeling and simulation of post combustion CO_{2} capturing,” Telemark University College, Faculty of Technology, Porsgrunn, Norway, 2008, M.Sc. thesis. View at: Google Scholar
 A. L. Kohl and R. B. Nielsen, Gas Purification, Gulf Publishing, Houston, TX, USA, 5th edition, 1997.
 G. F. Versteeg and W. P. M. van Swaaij, “Solubility and diffusivity of acid gases (CO_{2}, N_{2}O) in aqueous alkanolamine solutions,” Journal of Chemical & Engineering Data, vol. 33, pp. 29–34, 1988. View at: Publisher Site  Google Scholar
 D. P. Hagewiesche, S. S. Ashour, H. A. AlGhawas, and O. C. Sandall, “Absorption of carbon dioxide into aqueous blends of monoethanolamine and methyldiethanolamine,” Chemical Engineering Science, vol. 50, no. 7, pp. 1071–1079, 1995. View at: Publisher Site  Google Scholar
 H. Bosch, G. F. Versteeg, and W. P. M. van Swaaij, “Gasliquid mass transfer with parallel reversible reactionsIII. Absorption of CO_{2} into solutions of blends of amines,” Chemical Engineering Science, vol. 44, no. 11, pp. 2745–2750, 1989. View at: Publisher Site  Google Scholar
 R. Sakwattanapong, A. Aroonwilas, and A. Veawab, “Reaction rate of CO_{2} in aqueous MEAAMP solution: experiment and modeling,” Energia Procedia, vol. 1, no. 1, pp. 217–224, 2009. View at: Publisher Site  Google Scholar
 P. Usubharatana and P. Tontiwachwuthikul, “Enhancement factor and kinetics of CO_{2} capture by MEAmethanol hybrid solvents,” Energia Procedia, vol. 1, no. 1, pp. 95–102, 2009. View at: Publisher Site  Google Scholar
 T. Greer, A. Bedelbayev, J. M. Igreja, J. F. Pereira Gomes, and B. Lie, “A dynamic model for the deabsorption of carbon dioxide from monoethanolamine solution,” in Proceedings of SIMS, Amsterdam, Netherlands, October 2008. View at: Google Scholar
 M. Caplow, “Kinetics of carbamate formation and breakdown,” Journal of the American Chemical Society, vol. 90, no. 24, pp. 6795–6803, 1968. View at: Publisher Site  Google Scholar
 P. V. Danckwerts, “The reactions of CO_{2} with ethanolamines,” Chemical Engineering Science, vol. 34, no. 4, pp. 443–446, 1979. View at: Publisher Site  Google Scholar
 G. F. Versteeg and W. P. M. van Swaaij, “On the kinetics between CO_{2} and alkanolamines both in aqueous and nonaqueous solutionsI. Primary and secondary amines,” Chemical Engineering Science, vol. 43, no. 3, pp. 573–585, 1988. View at: Publisher Site  Google Scholar
 R. J. Littel, W. P. M. van Swaaij, and G. F. Versteeg, “Kinetics of CO_{2} with primary and secondary amines in aqueous solutionI. Zwitterion formation deprotonation kinetics for DEA and DIPA in aqueous blends of alkanolamines,” Chemical Engineering Science, vol. 47, no. 8, pp. 2027–2035, 1992. View at: Publisher Site  Google Scholar
 C. AlvarezFuster, N. Midoux, A. Laurent, and J. C. Charpentier, “Chemical kinetics of the reaction of carbon dioxide with amines in pseudo mn^{th} order conditions in aqueous and organic solutions,” Chemical Engineering Science, vol. 35, no. 8, pp. 1717–1723, 1980. View at: Publisher Site  Google Scholar
 E. Sada, H. Kumazawa, Y. Osawa, M. Matsuura, and Z. Q. Han, “Reaction kinetics of carbon dioxide with amines in nonaqueous solvents,” Chemical Engineering Journal, vol. 33, no. 2, pp. 87–95, 1986. View at: Publisher Site  Google Scholar
 R. J. Littel, W. P. M. van Swaaij, and G. F. Versteeg, “Kinetics of carbon dioxide with tertiary amines in aqueous solution,” AIChE Journal, vol. 36, no. 11, pp. 1633–1640, 1990. View at: Publisher Site  Google Scholar
 T. L. Donaldson and Y. N. Nguyen, “Carbon dioxide reaction and transportation in aqueous amine membranes,” Industrial & Engineering Chemistry Fundamentals, vol. 19, no. 3, pp. 260–266, 1980. View at: Publisher Site  Google Scholar
 G. F. Versteeg and W. P. M. van Swaaij, “On the kinetics between CO_{2} and alkanolamines both in aqueous and nonaqueous solutionsII. Ternary amines,” Chemical Engineering Science, vol. 43, no. 3, p. 587, 1988. View at: Publisher Site  Google Scholar
 H. Bosch, G. F. Versteeg, and W. P. M. van Swaaij, “Gasliquid mass transfer with parallel reversible reactionsI. Absorption of CO_{2} into solutions of sterically hindered amines,” Chemical Engineering Science, vol. 44, no. 11, pp. 2723–2734, 1989. View at: Publisher Site  Google Scholar
 B. P. Mandal, M. Guhab, A. K. Biswasb, and S. S. Bandyopadhyaya, “Removal of carbon dioxide by absorption in mixed amines: modelling of absorption in aqueous MDEAMEA and AMPMEA solutions,” Chemical Engineering Science, vol. 56, no. 2122, pp. 6217–6224, 2001. View at: Publisher Site  Google Scholar
 A. Benamor and M. K. Aroua, “Modeling of CO_{2} solubility and carbamate concentration in DEA, MDEA and their mixtures using the DeshmukhMather model,” Fluid Phase Equilibria, vol. 231, no. 2, pp. 150–162, 2005. View at: Publisher Site  Google Scholar
 A. Jamal, A. Meisen, and C. Jim Lim, “Kinetics of carbon dioxide absorption and desorption in aqueous alkanolamine solutions using a novel hemispherical contactor—II: experimental results and parameter estimation,” Chemical Engineering Science, vol. 61, no. 19, pp. 6590–6603, 2006. View at: Publisher Site  Google Scholar
 S. Freguia and G. T. Rochelle, “Modeling of CO_{2} capture by aqueous monoethanolamine,” AIChE Journal, vol. 49, no. 7, pp. 1676–1686, 2003. View at: Publisher Site  Google Scholar
 A. Bedelbayev, T. Greer, and B. Lie, “Model based control of absorption tower for carbon dioxide capturing,” in Proceedings of SIMS, pp. 1420–1440, Amsterdam, Netherlands, February 2008. View at: Google Scholar
 N. Haimour, A. Bidarian, and O. Sandall, “Kinetics of the reaction between carbon dioxide and methyldiethanolamine,” Chemical Engineering Science, vol. 42, no. 6, pp. 1393–1398, 1987. View at: Publisher Site  Google Scholar
 J. E. Critchfield, “CO_{2} absorption/desorption in methyldiethanolamine solutions promoted with MEA and D2EA. Mass transfer and reaction kinetics,” University of Texas, Austin, TX, USA, 1988, Ph. D. thesis. View at: Google Scholar
 H. A. Rangwala, B. R. Morrell, A. E. Mather, and F. D. Otto, “Absorption of CO_{2} into aqueous tertiary amine/MEA solutions,” in Proceedings of AIChE National Meeting, vol. 8, New Orleans, LA, USA, March 1988. View at: Google Scholar
 B. R. W. Pinsent, L. Pearson, and F. J. W. Roughton, “The kinetics of combination of carbon dioxide with hydroxide ions,” Journal of the Chemical Society, Faraday Transactions, vol. 52, pp. 1512–1520, 1956. View at: Publisher Site  Google Scholar
 Y. Liu, L. Zhang, and S. Watanasiri, “Representing vapourliquid equilibrium for an aqueous MEACO_{2} system using the electrolyte nonrandomtwoliquid model,” Industrial & Engineering Chemistry Research, vol. 38, no. 5, pp. 2080–2090, 1999. View at: Publisher Site  Google Scholar
 L. L. Tavlarides, R. Lin, Y. Nan et al., Sorption Modeling and Verification for OffGas Treatment, NEUP, Washington, DC, USA, 2015, https://neup.inl.gov/.
 R. H. Perry and D. W. Green, Chemical Engineers Handbook, McGrawHill, New York, NY, USA, 7th edition, 1999.
 G. F. Versteeg, J. V. Holst, P. P. Politiek, and J. P. Niederer, “CO_{2} capture from flue gas using amino acid salt solutions,” June 2006, http://www.co2cato.nl/doc.php?lid=317. View at: Google Scholar
 C. Gòmez, D. O. Borio, and N. S. Schbib, “Simulation of an industrial packed column for reactive absorption of CO_{2},” Latin American Applied Research, vol. 33, pp. 201–205, 2003. View at: Google Scholar
 E. P. Van Elk, M. C. Knaap, and G. F. Versteeg, “Application of the penetration theory for gasliquid mass transfer without liquid bulk differences with systems with a bulk,” IChemE Symposium Series, vol. 152, 2006. View at: Google Scholar
 J. Gabitto and C. Tsouris, “Dynamic simulation of CO_{2} absorption in a packing tower,” in Proceedings of AIChE Annual Meeting, vol. 11, Atlanta, GA, USA, November 2014. View at: Google Scholar
 A. Lawal, M. Wang, P. Stephenson, and H. Yeubg, “Dynamic modeling of CO_{2}absorption for post combustion capture in coalfired power plants,” Fuel, vol. 88, no. 12, pp. 2455–2462, 2009. View at: Publisher Site  Google Scholar
 J. H. Richardson, J. H. Harker, and J. R. Backhurst, Chemical Engineering, vol. 2, ButterworthHeinemann, Publishers, Oxford, UK, 5th edition, 2008, Particle Technology and Separation Processes.
 R. Billet and M. Schultes, “Prediction of mass transfer columns with dumped and arranged packing’s: updated summary of the calculation method of billet and schultes,” Chemical Engineering Research and Design, vol. 77, no. 6, pp. 498–504, 1999. View at: Publisher Site  Google Scholar
 A. Aboudheir, P. Tontiwachwuthikula, A. Chakmab, and R. Idema, “Kinetics of the reactive absorption of carbon dioxide in high CO_{2}loaded, concentrated aqueous monoethanolamine solutions,” Chemical Engineering Science, vol. 58, no. 2324, pp. 5195–5210, 2003. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Jorge Gabitto and Costas Tsouris. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.