Optimization in the Absorption and Desorption of CO 2 Using Sodium Glycinate Solution

: This study used sodium glycinate as an absorbent to absorb CO 2 in the bubble column scrubber under constant pH and temperature environments to obtain the operating range, CO 2 loading, and mass transfer coefﬁcient. For efﬁcient experimentation, the Taguchi method is used for the experimental design. The process parameters are the pH, gas ﬂow rate ( Q g ), liquid temperature ( T ), and absorbent concentration ( C L ). The effects of the parameters on the absorption efﬁciency, absorption rate, overall mass transfer coefﬁcient, gas–liquid molar ﬂow rate ratio, CO 2 loading, and absorption factor are to be explored. The optimum operating conditions and the order of parameter importance are obtained using the signal/noise (S/N) ratio analysis, and the optimum operating conditions are further veriﬁed. The veriﬁcation of the optimum values was also carried out. The order of parameter importance is pH > C L > Q g > T . Evidence in the 13 CNMR (Carbon 13 Nuclear Magnetic Resonance) spectra shows that the pH value has an effect on the solution composition, which affects both the absorption efﬁciency and mass transfer coefﬁcient. There are 18 experiments for regeneration, where the operating temperature is 100–120 ◦ C. The heat of regeneration was measured according to the thermodynamic data. The CO 2 loading, the overall mass transfer, and the heats of regeneration correlation are also discussed in this work. Finally, an operating policy for the CO 2 absorption process was conﬁrmed. and circulator are off, the mass of the sodium glycinate after the is measured, and samples are taken. The samples are tested using an elements analyzer (EA) (Thermo, Flash EA1112, UK).


Introduction
Global warming is a critical issue in the world. Because of this, the Paris Agreement declared that the global mean temperature increment must not exceed 2 • C by 2100, and aims to reduce the temperature to 1.5 • C. This shows that the reduction of greenhouse gas is more pressing [1][2][3]. Fossil fuels are the most significant factors in the warming process, and the energy, petrochemical, steel, transportation, and cement industries are the major contributors. Electricity is mainly derived from thermal power plants, and the thermal power plants burn fossil fuels, where CO 2 is the main source of the greenhouse effect.
Nowadays, there are three common systems for capturing CO 2 , including post-combustion, pre-combustion, and oxyfuel-combustion. In a post-combustion system, the CO 2 concentration is 3-15% [4], which is mostly removed by chemical absorption. The common solvents for chemical solvent absorption processes include amines, NaOH, aqueous ammonia, ionic liquids, hot potassium carbonate, and amino acid salt. Among these solvents, the amine is studied a lot. Monoethanolamine (MEA) is universally used as an absorbent by over a thousand commercial plants, because it has a high absorption rate, high alkalinity, regenerability, and low price, but the it also has low loading, deterioration, high corrosively, and a high regenerated energy [3]. In order to improve the absorption and stripping efficiency, it is imperative to find new, more efficient solvents, such as In Equations (1)- (5), it can be found that the pH of the solution has an effect on the absorption and the distribution of the species. Increasing the pH means decreasing the [H + ], which favors a forward reaction in Equation (1), and hence favors CO 2 absorption. On the other hand, it favors the formation of HCO 3 − , − OOC − R − NHCOO − , and − OOC − R − NH 2 . Aqueous amino acid salt, which is similar to that of aqueous amino acid zwitterion, will be protonated or deprotonated by changing the pH, as shown in reaction [18], as follows: It can be found that a low pH (pH < pK 1 ) amino acid salt represents positive charges, HOOC − R − NH + 3 . In the range of pK 1 and pK 2 , the medium pH range becomes the zwitterion form, − OOC − R − NH + 3 . At a low pH value, the reaction of amino acid with CO 2 is less. Thus, knowing how to control the pH, temperature, and flow rate in order to obtain the desired absorption rate, efficiency, and mass transfer coefficient becomes significant. This can be achieved by adding an amino acid salt solution into the scrubber, which can maintain the desired pH to a certain value, where CO 2 can be effectively captured.
In addition, the selection of the scrubber is also significant, as it affects the performance of the CO 2 absorption. Here, the bubble-column has been used, as it shows a superior performance, such as a high absorption, high mass-transfer coefficient, simple structure, and easy operation, compared with other scrubbers [2]. However, some parameters, such as the gas-flow rate, liquid-flow rate, solution temperature, and the pH of the solution have been found to be key factors affecting the absorption of CO 2 . Because of the importance of mass transfer in a bubble-column, some papers have studied the effects of the process parameters on the mass transfer coefficients [31][32][33][34][35][36]. The reported values were in the range of 0.01-0.888, depending on the conditions. A comparison of their mass transfer coefficients is listed in Table 1. However, there are no available data in the literature regarding the absorption of CO 2 using a sodium glycinate solution. In addition, the effect of pH on the solution and operating range can be effectively explored using a bubble-column scrubber. Therefore, a study of the absorption and desorption of CO 2 using a sodium glycinate solution becomes significant.

Methodology
In order to reach their purpose, a mass transfer model accompanying a two-film model, thermodynamic data, and an energy duty model were used. The framework used in here is proposed in Figure 1. This study is divided into three parts, as follows: the first part uses sodium glycinate solution to capture CO 2 . In order to explore the effects of the process variables on the absorption efficiency (E), absorption rate (R A ), overall mass transfer coefficient (K G a), gas-liquid molar flow rate ratio (γ), absorption factor (ϕ), and CO 2 loading (α) in a steady state condition, a bubble-column scrubber was adopted to study the capture of CO 2 under a constant pH and constant solution temperature environments. The second part uses the Taguchi experimental design and Taguchi analysis to find the optimum operating conditions and the order of importance, which are validated further. In addition, two optimum operating conditions are used for the uncontrolled experiments. The third part uses the scrubbed liquids collected from optimum conditions for the regeneration experiments, so as to evaluate the energy duty of the sodium glycinate solution. Finally, the correlation equations were obtained for further discussion. The third part uses the scrubbed liquids collected from optimum conditions for the regeneration experiments, so as to evaluate the energy duty of the sodium glycinate solution. Finally, the correlation equations were obtained for further discussion.

Determination of RA and KGa
In a bubble-column scrubber, a gas mixture containing CO2 (A) and N2 (B) flowing into a bubble column from the bottom comes into continuous contact with the liquid solvent flowing into the column from the top. The gas mixture and liquid simultaneously come into contact within the column counter. If a binary A-B system is assumed, then RA and KGa can be obtained using the mass balance with two-film models at a steady-state condition [2]. Therefore, RA can be written as follows: where FA1 is the CO2 gas molar flow rate, which can be evaluated when Qg (gas flow rate), PA1 (partial pressure of CO2 at inlet), and T (liquid temperature in the buble column) are available. In addition, VL is the volume of liquid in the scrubber, yA1 is the CO2 gas molar fraction at the inlet, and yA2 is the CO2 gas molar fraction at the outlet. On the other hand, KGa becomes the following:

Determination of R A and K G a
In a bubble-column scrubber, a gas mixture containing CO 2 (A) and N 2 (B) flowing into a bubble column from the bottom comes into continuous contact with the liquid solvent flowing into the column from the top. The gas mixture and liquid simultaneously come into contact within the column counter. If a binary A-B system is assumed, then R A and K G a can be obtained using the mass balance with two-film models at a steady-state condition [2]. Therefore, R A can be written as follows: where F A1 is the CO 2 gas molar flow rate, which can be evaluated when Q g (gas flow rate), P A1 (partial pressure of CO 2 at inlet), and T (liquid temperature in the buble column) are available. In addition, V L is the volume of liquid in the scrubber, y A1 is the CO 2 gas molar fraction at the inlet, and y A2 is the CO 2 gas molar fraction at the outlet. On the other hand, K G a becomes the following: Considering the ideal gas law for the inlet and outlet gases, F A1 /F A2 could be replaced with (P 1 /P 2 ) (T 2 /T 1 ) (y A1 /y A2 ).

Determination of Gas-Liquid Flow Rate Ratio
Assuming that the input gas obeys the ideal gas law, the total gas flow rate can be calculated by Equation (9), as follows: where Q g is the total gas volumetric flow rate and P 1 is the total pressure at the inlet. On the other hand, the liquid flow rate can be estimated by Equation (10), as follows: where Q L is the total liquid volumetric flow rate determined at a steady state condition and C L is the concentration of sodium glycinate solution. Once F G and F L have been determined, γ can be obtained, as follows:

Determination of Scrubbing Factor
The scrubbing factor has been defined in a previous study [2], as shown below: where V is the volume of the scrubber. The scrubbing factor can be evaluated when γ, E, V, and y 1 are available.

Thermal Data Estimated
Heat duty of solvent regeneration includes three parts, as follows: Q(GJt −1 ) = Q abs (GJ/t) + Q sen (GJ/t) + Q sol (GJ/t) They can be determined when the thermal data are available. The thermal data used here includes the heat capacity, heat of absorption, and latent heat, which have been reported in the literature, as shown in Table 2. In here, ∆H abs is heat of absorption, which is related to loading α and temperature [22,29]. In addition, C p is the heat capacity of the solvent, ∆T is the temperature difference, m sol is the mass of regeneration solvent, ∆m CO 2 is the mass loss of CO 2 after stripping. Here, C p can be adopted, as shown in the literature [37], and can be correlated with the weight fraction and temperature. Finally, ∆H vap is the heat of evaporation, . Q c is the cooling rate, ∆m 1 is the solvent loss during stripping, and t is the evaporation time. By collecting the thermal data reported in the literature [27], ∆H vap can be expressed as a function of the weight fraction (ω). The three thermal data equations were shown below: ∆H abs (kJ/mol) = 5.8048 exp( 583.84 T(K) )α 0.5228 (14) C P (kJ/kg · K) = 3.1702 + 9.067 × 10 −4 T (K) (15) ∆H vap (kJ/mol · K) = 4.244ω 2 − 4.937ω + 2.023 (16) Figure 1 shows the structure diagram of the Taguchi experimental design for the absorption and solvent regeneration studies. Three parts have been explored in the absorption study: the first is the Taguchi experiment, the second is the verification of the optimum conditions, and the third is the uncontrolled study. Finally, the solvent regeneration study is explored, using the solvents collected in the verification experiments as test solvents.

Experimental Design for Absorption Study
The first part of this study is CO 2 absorption. The process parameters used in here are the pH of the solution, gas flow rate (Q g ), liquid temperature (T), and solvent concentration (C L ), and each factor has four levels, as shown in Table 3. Theoretically, there should be 4 4 = 256 experiments. Because of this, a Taguchi experimental design [16] was adopted, and the number of experiments was reduced to L 16 (4 4 ) = 16; thus, we could save on the time and experimental cost. The signal/noise (S/N) ratio is calculated according to the larger-the-better value and smaller-the-better values from the data obtained in a steady state condition, expressed as Equations (17) and (18), respectively. The sequence of importance and optimum conditions are found using the S/N analysis. Where E, R A , K G a, and ϕ take the larger-the-better value, and γ takes the larger-the-better value and the smaller-the-better value, searching for the gas and liquid flow rate ratio range. Table 3 shows the factors and levels of the absorption experiment. Table 4 shows the orthogonal table for the absorption experiment. and where S/N is the objective function, n is the number of data, and z i is the number of measures. Table 3. Parameters and levels in the absorption study. Q g -gas flow rate; C L -solvent concentration.  Table 4. Orthogonal array in the absorption study.

Experimental Design for Solvent Regeneration
The second part of this study is the scrubbed solvent regeneration; this experimental procedure is similar to that reported in the literature [9,38]. The test solvents that were selected are from the optimum conditions obtained in here. A total of six solvents were adopted each for the three level temperatures (100 • C, 110 • C, and 120 • C).

Absorption of CO 2
The equipment for the absorption experiment is shown in Figure 2, including the bubble column, cooling circulator, tubing pump, CO 2 detector, heater, mass flow controller, and pressure gauge. The sodium glycinate solution absorbs the CO 2 in the bubble column, under a constant pH environment, then the simulated flue gas concentration (15%) and temperature (50 • C) from the coal-fired power plant are pumped into the column. First, the absorbent concentration for the experiment is prepared-the CO 2 detector (Guardian Plus, D600, Hartford, CT, USA), pressure indicator (SR-321), pH controller (Suntex, PC-310, New Taipei City, Taiwan), and heating controller (5020 Data Acquisition System) are switched on; the mass flow controller (N 2 and CO 2 ) (Bronkhorst, F-201CB DMFC) is adjusted to the correct flow for the experiment; and the CO 2 concentration is changed to 15%. The heater is switched on, and the temperature is controlled at 50 • C. The cooling circulator (Deng Yng, D-620, New Taipei City, Taiwan) is switched on and the temperature for the experiment is set. When the gas inlet temperature is fixed at 50 • C, the CO 2 concentration is fixed at 15%, and the corrected pH electrode is plugged in the bubble column for operation. The initial CO 2 concentration, pH, liquid temperature, gas temperature, gas inlet pressure, and absorbent feed volume are recorded. volume in the column is recorded. The tubing pump controller (EASY-LOAD, 7518-00) and pH electrode controller are switched on for timing. As the pH decreases after the solution absorbs CO2, the pH controller injects the sodium glycinate solution into the column automatically, so that the pH is controlled. The experimental data are recorded once every 5 min, and 10 mL of solution is taken as the sample every 10 min. When the CO2 concentration measured by the CO2 detector reaches a steady state, the experiment can be ended.

Regeneration Test
The equipment for the regeneration experiment is shown in Figure 3, including a cooling circulator, oil bath, magnetic stirrer heating plate, three-neck round flask, and a ball condenser tube. The apparatus is similar to that reported in the literature [9,37]. Firstly, 0.05 kg of solution after absorption is prepared, and the ball condenser tube, three-neck round flask, and cooling circulator are assembled. The experimental heating temperature and cooling water tank temperature (5 °C) are set. When the heating temperature and cooling circulator temperature are stable, the 0.05 kg solution is poured in and the magnetic stirrer is switched on. The experimental time is 60 min, and the temperature change is recorded once every five minutes. When the experiment is finished, the heating controller and cooling circulator are switched off, the mass of the sodium glycinate solution after the experiment is measured, and samples are taken. The samples are tested using an elements analyzer (EA) (Thermo, Flash EA1112, Cambridge, UK). Subsequently, the solvent is poured into the column until it overflows, and the initial solution volume in the column is recorded. The tubing pump controller (EASY-LOAD, 7518-00) and pH electrode controller are switched on for timing. As the pH decreases after the solution absorbs CO 2 , the pH controller injects the sodium glycinate solution into the column automatically, so that the pH is controlled. The experimental data are recorded once every 5 min, and 10 mL of solution is taken as the sample every 10 min. When the CO 2 concentration measured by the CO 2 detector reaches a steady state, the experiment can be ended.

Regeneration Test
The equipment for the regeneration experiment is shown in Figure 3, including a cooling circulator, oil bath, magnetic stirrer heating plate, three-neck round flask, and a ball condenser tube. The apparatus is similar to that reported in the literature [9,37]. Firstly, 0.05 kg of solution after absorption is prepared, and the ball condenser tube, three-neck round flask, and cooling circulator are assembled. The experimental heating temperature and cooling water tank temperature (5 • C) are set. When the heating temperature and cooling circulator temperature are stable, the 0.05 kg solution is poured in and the magnetic stirrer is switched on. The experimental time is 60 min, and the temperature change is recorded once every five minutes. When the experiment is finished, the heating controller and cooling circulator are switched off, the mass of the sodium glycinate solution after the experiment is measured, and samples are taken. The samples are tested using an elements analyzer (EA) (Thermo, Flash EA1112, Cambridge, UK).

Results and Discussions
The calculated data for E, RA, KGa, γ, and φ are listed in Table 5. All of the data can be analyzed to search for the optimum conditions and the importance of parameters by using the S/N ratio.

Effects of pH and Temperature on the Solution Composition
The effect of pH and T on the composition of solution in the scrubbed solution can be determined by 13 CNMR spectra, as shown in Figure 4a-d for No. 1, 5, 9, and 13, respectively. In order for a comparison, the data was adopted using the same gas flow rate, at the rate of 3 L/min, but different pHs and temperatures of the solutions we used. In Figure 4, the major peaks were presented at 42, 45, 164, 174, and 179 ppm, which shows 3 / CO HCO in 162 ppm (C) [37,38]. The peak at 174 ppm shifted to 179 ppm, and disappeared when the pH value increased to 11, indicating the disappearance

Results and Discussions
The calculated data for E, R A , K G a, γ, and ϕ are listed in Table 5. All of the data can be analyzed to search for the optimum conditions and the importance of parameters by using the S/N ratio.

Effects of pH and Temperature on the Solution Composition
The effect of pH and T on the composition of solution in the scrubbed solution can be determined by 13 CNMR spectra, as shown in Figure 4a-d for No. 1, 5, 9, and 13, respectively. In order for a comparison, the data was adopted using the same gas flow rate, at the rate of 3 L/min, but different pHs and temperatures of the solutions we used. In Figure 4, the major peaks were presented at 42, 45, 164, 174, and 179 ppm, which shows GLY/GLH + at 42 ppm and 174 ppm (A , B ), GLYH + CO − 2 in 45 ppm and 179 ppm (A, B), and the HCO − 3 /CO 2− 3 in 162 ppm (C) [37,38]. The peak at 174 ppm shifted to 179 ppm, and disappeared when the pH value increased to 11, indicating the disappearance of carbamate, and hence an increase in the absorption of CO 2 . The effect of pH on the absorption of CO 2 can be explained in E and K G a, as shown in Table 4; the values of E and K G a are 55. 33

Taguchi Analysis
Using an S/N ratio analysis, using data E as an example, the larger-the-better value can be determined using Equation (17). The calculated values are shown in Table 6. The results show that the importance of the parameters is B (Q g ) > A (pH) > D (C L ) > C (T), while the optimum condition is A4B1C2D1. Other analysis data, such as R A , K G a, γ, and ϕ, can be determined using the same procedure. All of the results are collected in Table 7, showing the different results in an optimum condition and the importance of the parameters. In order to understand the importance of the parameters of the whole system, the ratios of 4, 3, 2, 1, for the sequence 1, 2, 3, and 4, respectively are set. It was found that the pH has highest rations (3.67) and T has lowest rations (1.17). Because of this, the importance of the parameters of the whole system became A > D > B > C. Figure 5 shows the importance of the parameters in a radar chart for whole system and for the individual data.  Table 7. Optimum conditions and importance of parameters in the absorption process, according to the Taguchi analysis.

Taguchi Analysis
Using an S/N ratio analysis, using data E as an example, the larger-the-better value can be determined using Equation (17). The calculated values are shown in Table 6. The results show that the importance of the parameters is B (Qg) > A (pH) > D (CL) > C (T), while the optimum condition is A4B1C2D1. Other analysis data, such as RA, KGa, γ, and φ, can be determined using the same procedure. All of the results are collected in Table 7, showing the different results in an optimum condition and the importance of the parameters. In order to understand the importance of the parameters of the whole system, the ratios of 4, 3, 2, 1, for the sequence 1, 2, 3, and 4, respectively are set. It was found that the pH has highest rations (3.67) and T has lowest rations (1.17). Because of this, the importance of the parameters of the whole system became A > D > B > C. Figure 5 shows the importance of the parameters in a radar chart for whole system and for the individual data.  Table 7. Optimum conditions and importance of parameters in the absorption process, according to the Taguchi analysis.

Verification of Optimum Conditions
Verification tests of the optimum conditions were carried out using the same procedure as in the Taguchi experiment. All of the data evaluated are listed in Table 8. All of the conditions that satisfied the optimum values are shown in red. This demonstrates the feasibility and reliability of the Taguchi experimental design.

Verification of Optimum Conditions
Verification tests of the optimum conditions were carried out using the same procedure as in the Taguchi experiment. All of the data evaluated are listed in Table 8. All of the conditions that satisfied the optimum values are shown in red. This demonstrates the feasibility and reliability of the Taguchi experimental design.

Uncontrolled Experimental Results
In order to better understand the absorption process in a pilot scale, uncontrolled experiments were carried out. Two runs were conducted here. One, for the condition of optimum E was selected (No. 17), and the other was the condition of optimum K G a (No. 19). However, pH and T were uncontrolled for both. Figure 6 shows the elapsed times in the pH and T for No. 23, while the controlled experiment (No. 17) is also presented in this figure as blue points for comparison. It was found that the change in pH was nearly the same for both of the systems, while the temperature in No. 23 was lower than in No. 17. However, the system for No. 23 reached a steady state after 40 min, which was similar to that of No. 17. Table 9 shows the data for both No. 17

Uncontrolled Experimental Results
In order to better understand the absorption process in a pilot scale, uncontrolled experiments were carried out. Two runs were conducted here. One, for the condition of optimum E was selected (No. 17), and the other was the condition of optimum KGa (No. 19). However, pH and T were uncontrolled for both. Figure 6 shows the elapsed times in the pH and T for No. 23, while the controlled experiment (No. 17) is also presented in this figure as blue points for comparison. It was found that the change in pH was nearly the same for both of the systems, while the temperature in No. 23 was lower than in No. 17. However, the system for No. 23 reached a steady state after 40 min, which was similar to that of No. 17. Table 9 shows the data for both No. 17 and No. 23. It can be seen that the values evaluated for the controlled experiments (No. 17 and No. 19) were higher than those obtained in the uncontrolled experiments (No. 23 and No. 24). The reason for this was that the gasliquid contact time for both were different. The contact time for the former was longer than that for the latter.

Loading of CO2 at the Final Stage
The loading of CO2 (α) can be determined using TOC analysis (Teledyne Tekmer, Phoenix 8000, Mason, OH, USA), and the results are shown in the last column of Table 4. The obtained values are in the range of 0.1081-0.6061 mol-CO2/mol-amine, depending on the operating conditions. It can be found that the value is related to the pH, T, and flow rate ratio (γ). The results show that the lower the pH, the higher the α value; however, it was also found that the higher the γ, the higher the α.

Loading of CO 2 at the Final Stage
The loading of CO 2 (α) can be determined using TOC analysis (Teledyne Tekmer, Phoenix 8000, Mason, OH, USA), and the results are shown in the last column of Table 4. The obtained values are in the range of 0.1081-0.6061 mol-CO 2 /mol-amine, depending on the operating conditions. It can be found that the value is related to the pH, T, and flow rate ratio (γ). The results show that the lower the pH, the higher the α value; however, it was also found that the higher the γ, the higher the α. However, the relationship between them is hard to quantitatively determine. Therefore, a linear regression including the Taguchi experimentation is required. A total of sixteen data sets, as listed in Table 4, are adopted, and the results are shown as follows: Figure 7 shows the confidence of regression, where it is found that the measured values are close to that of the calculated values, indicating that they are reliable in regression. Therefore, the equation can be used to calculate the values obtained in the optimum conditions and uncontrolled experiments. The results are also presented in Figure 7, where a slight deviation shows that Equation (19) is predicable. In addition, the parameter importance analysis shows γ > pH > T, indicating that the liquid gas contact is significant for CO 2 loading. For the sake of caparison, two solvents, MEA [35] and aqueous ammonia [36], were used in the bubble-column scrubbers for CO 2 capture, and they are also exhibited Figure 7. It was found that most of the data falls in between the error range. In addition, the loading of the aqueous ammonia was close to that obtained in this work.
Appl. Sci. 2018, 8, x 13 of 20 However, the relationship between them is hard to quantitatively determine. Therefore, a linear regression including the Taguchi experimentation is required. A total of sixteen data sets, as listed in Table 4, are adopted, and the results are shown as follows:  Figure 7 shows the confidence of regression, where it is found that the measured values are close to that of the calculated values, indicating that they are reliable in regression. Therefore, the equation can be used to calculate the values obtained in the optimum conditions and uncontrolled experiments. The results are also presented in Figure 7, where a slight deviation shows that Equation (19) is predicable. In addition, the parameter importance analysis shows γ > pH > T, indicating that the liquid gas contact is significant for CO2 loading. For the sake of caparison, two solvents, MEA [35] and aqueous ammonia [36], were used in the bubble-column scrubbers for CO2 capture, and they are also exhibited Figure 7. It was found that most of the data falls in between the error range. In addition, the loading of the aqueous ammonia was close to that obtained in this work.

Overall Mass-Transfer Coefficient
The effects of T, pH, CL, and Qg on KGa can be found in Table 4. However, the effective correlation equation for KGa is not prone to being obtain. Here, the mass transfer coefficient can be correlated to the liquid molar flux, gas molar flux, and temperature [39,40]. In order to obtain the empirical equation, the gas and liquid flow rates are needed to translate into the molar flux, that is, Gx (mol/s·m 2 ) and Gy (mol/s·m 2 ); the values are in the range of 0.9469-2.8407 mol/s·m 2 and 0.2321-1.7936 mol/s·m 2 for Gy and Gx, respectively. Using the linear regression, a total of sixty data sets for KGa, as listed in Table 4, were used. The correlation results became the following:

Overall Mass-Transfer Coefficient
The effects of T, pH, C L , and Q g on K G a can be found in Table 4. However, the effective correlation equation for K G a is not prone to being obtain. Here, the mass transfer coefficient can be correlated to the liquid molar flux, gas molar flux, and temperature [39,40]. In order to obtain the empirical equation, the gas and liquid flow rates are needed to translate into the molar flux, that is, G x (mol/s·m 2 ) and G y (mol/s·m 2 ); the values are in the range of 0.9469-2.8407 mol/s·m 2 and 0.2321-1.7936 mol/s·m 2 for G y and G x , respectively. Using the linear regression, a total of sixty data sets for K G a, as listed in Table 4, were used. The correlation results became the following: In order to understand the confidence of regression, a plot of (K G a) cal versus (K G a) mea is shown in Figure 8, where it is found that the most data are within a 15% error and are close to the solid line, indicating that regression is reasonable. In order to verify the optimum experiments and the uncontrolled experiments, all of the data were plotted in the figure. It was found that all of the data are close to the solid line, thus, Equation (20) is reliable. In addition, from an operation range analysis, the parameter importance was found to be G x > T > G y . Because of this, controlling G x can effectively improve the K G a. This could be done by increasing the pH of the solution. In addition, the K G a values obtained in the two solvents, MEA [35] and aqueous ammonia [36], which were used in the bubble-column scrubbers for CO 2 absorption, were shown in the same figure for comparison. It was found that the reported data for the two solvents were close to that obtained in this work. This indicated that the K G a of a sodium glycinate solution is comparable to that of the MEA and aqueous ammonia solvents. The range in K G a for the three solvents was 0.04-0.4 s −1 , depending on the operating conditions. In order to understand the confidence of regression, a plot of (KGa)cal versus (KGa)mea is shown in Figure 8, where it is found that the most data are within a 15% error and are close to the solid line, indicating that regression is reasonable. In order to verify the optimum experiments and the uncontrolled experiments, all of the data were plotted in the figure. It was found that all of the data are close to the solid line, thus, Equation (20) is reliable. In addition, from an operation range analysis, the parameter importance was found to be Gx > T > Gy. Because of this, controlling Gx can effectively improve the KGa. This could be done by increasing the pH of the solution. In addition, the KGa values obtained in the two solvents, MEA [35] and aqueous ammonia [36], which were used in the bubblecolumn scrubbers for CO2 absorption, were shown in the same figure for comparison. It was found that the reported data for the two solvents were close to that obtained in this work. This indicated that the KGa of a sodium glycinate solution is comparable to that of the MEA and aqueous ammonia solvents. The range in KGa for the three solvents was 0.04-0.4 s −1 , depending on the operating conditions.

Heat of Regeneration of CO2 Loaded Sodium Glycinate Solution
The regeneration energy can be determined by the thermal data and energy balance, as shown in Equation (13), and the results are listed in Table 10 and Figure 9. It was found that the individual penalties were 0.181-0.453 GJ t −1 , 1.66-4.38 GJ t −1 , and 0.51-7.43 GJ t −1 for Qabs, Qsen, and Qsol, respectively. Therefore, the total energy required here was in the range of 3.68-10.75 GJ t −1 , which is higher than 5.0 GJ/t, comparable to that reported by Rabensteiner et al. [22]. From the energy duty distribution analysis, it can be found that the heat of absorption is below 5%, while the heat of evaporation and the sensitive heat are both close to 48%. Therefore, a reduction in Qsen and Qsol can effectively obtain a small Q.

Heat of Regeneration of CO 2 Loaded Sodium Glycinate Solution
The regeneration energy can be determined by the thermal data and energy balance, as shown in Equation (13), and the results are listed in Table 10 and Figure 9. It was found that the individual penalties were 0.181-0.453 GJ t −1 , 1.66-4.38 GJ t −1 , and 0.51-7.43 GJ t −1 for Q abs , Q sen , and Q sol , respectively. Therefore, the total energy required here was in the range of 3.68-10.75 GJ t −1 , which is higher than 5.0 GJ/t, comparable to that reported by Rabensteiner et al. [22]. From the energy duty distribution analysis, it can be found that the heat of absorption is below 5%, while the heat of evaporation and the sensitive heat are both close to 48%. Therefore, a reduction in Q sen and Q sol can effectively obtain a small Q.   However, the heat of regeneration can be correlated with the pH of the solution, regeneration temperature, and CO2 loading. A total of eighteen data was adopted for regression as shown below: It can be found that most of the data in Figure 10 shows a ± 20% error. Equation (21)   However, the heat of regeneration can be correlated with the pH of the solution, regeneration temperature, and CO 2 loading. A total of eighteen data was adopted for regression as shown below: Q(GJ/t) = 0.5110 exp(0.9343pH − 2158.9 T R (K) )α 1.002 (21) It can be found that most of the data in Figure 10 shows a ±20% error. Equation (21) can be used to predict the heat of regeneration when the loading, temperature, and pH of the solution are given. In addition, the analysis showed that the sequence of the parameter is α > pH > T R , indicating the importance of CO 2 loading. This implies that the loading control is important for reducing the heat of regeneration. In addition, in the real stripper, the efficiency could be higher, as the gas-liquid contact environment is better and the effect of pressure on the CO 2 stripping could be obviously. given. In addition, the analysis showed that the sequence of the parameter is α > pH > TR, indicating the importance of CO2 loading. This implies that the loading control is important for reducing the heat of regeneration. In addition, in the real stripper, the efficiency could be higher, as the gas-liquid contact environment is better and the effect of pressure on the CO2 stripping could be obviously.

Operating Policy
From the uncontrolled experiments, the operating pH and liquid temperature (T) are around 11 and 28 °C, respectively, after steady-state operation. The KGa and E values for No. 23 and No. 24, shown in Table 9, are 85.90 and 0.125, and 66.23 and 0.243, respectively. The different E and KGa values for the two runs are due to differences in Gx and Gy, which are 1.03 and 1 for No. 23, and 3.66 and 3 for No. 24, respectively. This indicates that at a higher mixing intensity, No. 24 has a higher KGa, but reduces its removal efficiency. From the designation viewpoint, a higher KGa has a smaller size, and hence gets more effective in the removal of CO2. If we assume the removal efficiency E = 80%, the KGa becomes 0.16 s −1 , which falls in between the two values. In this manner, the evaluated α is 0.1123, using Equation (19), and the Q is 5.08 GJ/t, using Equation (21), at TR = 373K. Therefore, a mild mixing and increase in the pH of the solution are required to give a higher KGa and a higher E. Alternatively, an agitator put inside the column could enhance the dispersion and mixing, which makes KGa and E higher.

Conclusions
This study successfully used a continuous bubble-column scrubber with a sodium glycinate solution to capture CO2. Using material balance with a two-film model, the absorption rate, overall mass-transfer coefficient, CO2 loading, and scrubbing factor could be determined at a steady-state condition. Quantitatively, the effects of the variables on the E, KGa, RA, γ, α, and φ can be explained using the empirical equations obtained by this study. The operating range for γ was found to be in the range of 0.74-11.37, which is flexible. The heat of the regeneration (Q) mechanism analysis suggests that a reduction of the evaporation heat and sensitive heat can obviously attenuate the Q in a stripper. The 13 CNMR spectra of the scrubbed solution for the glycine system showed GLY/GLYH + , GLYH + CO2 − , and HCO3 − /CO3 2− , depending on the absorption conditions. At a higher pH, it favors E and KGa, but reduces the values of γ, while α and φ show the pH on the effect of the solution chemistry. E and KGa can reach 80% and 0.16 s −1 , respectively; when γ < 1, it corresponds to a pH of

Operating Policy
From the uncontrolled experiments, the operating pH and liquid temperature (T) are around 11 and 28 • C, respectively, after steady-state operation. The K G a and E values for No. 23 and No. 24, shown in Table 9, are 85.90 and 0.125, and 66.23 and 0.243, respectively. The different E and K G a values for the two runs are due to differences in G x and G y , which are 1.03 and 1 for No. 23, and 3.66 and 3 for No. 24, respectively. This indicates that at a higher mixing intensity, No. 24 has a higher K G a, but reduces its removal efficiency. From the designation viewpoint, a higher K G a has a smaller size, and hence gets more effective in the removal of CO 2 . If we assume the removal efficiency E = 80%, the K G a becomes 0.16 s −1 , which falls in between the two values. In this manner, the evaluated α is 0.1123, using Equation (19), and the Q is 5.08 GJ/t, using Equation (21), at T R = 373K. Therefore, a mild mixing and increase in the pH of the solution are required to give a higher K G a and a higher E. Alternatively, an agitator put inside the column could enhance the dispersion and mixing, which makes K G a and E higher.

Conclusions
This study successfully used a continuous bubble-column scrubber with a sodium glycinate solution to capture CO 2 . Using material balance with a two-film model, the absorption rate, overall mass-transfer coefficient, CO 2 loading, and scrubbing factor could be determined at a steady-state condition. Quantitatively, the effects of the variables on the E, K G a, R A , γ, α, and ϕ can be explained using the empirical equations obtained by this study. The operating range for γ was found to be in the range of 0.74-11.37, which is flexible. The heat of the regeneration (Q) mechanism analysis suggests that a reduction of the evaporation heat and sensitive heat can obviously attenuate the Q in a stripper. The 13 CNMR spectra of the scrubbed solution for the glycine system showed GLY/GLYH + , GLYH + CO 2 − , and HCO 3 − /CO 3 2− , depending on the absorption conditions. At a higher pH, it favors E and K G a, but reduces the values of γ, while α and ϕ show the pH on the effect of the solution chemistry. E and K G a can reach 80% and 0.16 s −1 , respectively; when γ < 1, it corresponds to a pH of 11. In this manner, α is less than 0.12 mol/mol, which can reduce the Q value. The regression equations in K G a, α, and Q can satisfy the optimum conditions and uncontrolled experiments. Because of this, empirical equations can be used to predict the K G a, α, and Q for a pilot-scale bubble-column scrubber and heat regeneration. Finally, the order of parameter importance is pH > C L > Q g > T, which can be effectively adjusted in order to obtain the desired E, K G a, and α within the operating range in this work.

Conflicts of Interest:
The authors declare no conflicts of interest.