3.1. Characterization
Amine-functionalized SBA-15 was synthesized to improve the reaction heat between CO2 and the solvent. Because the template is mesoporous silica, the pore structure can be stably maintained under various conditions of temperature and pressure. In contrast, metal organic frameworks (MOFs) or covalent organic frameworks (COFs) form Armstrong-scale pores through a lattice structure. However, MOF and COF templates are unsuitable for amine functionalization, considering that the theoretical diameter of the amine groups is 0.935 nm. Moreover, they are vulnerable to moisture and cannot remain dispersion in aqueous amine solutions. An FE-SEM image of TEPA@SBA-15 is shown in Fig. 4a. The diameter is distributed in the range of 110–250 nm, the length is 500–750 nm, and the aspect ratio is large owing to pore structure growth during the self-assembly process. The structure of TEPA@SBA-15 dispersed in EtOH was studied using the TEM image as presented in Fig. 4b. The pore dimeters are distributed in the range of 3.7–5.8 nm. The cross-section is hexagonal, and the dimensions of SBA-15 observed by TEM and SEM are similar.
The bonding groups on the surface of SBA-15 were analyzed using XPS, as shown in Figure S2. It is mainly composed of silicon and oxygen, corresponding to amorphous silica; nitrogen is not detected. As shown in Figure S3, which presents the XPS spectrum of TEPA@SBA-15, the component ratio of silicon and oxygen is similar to that of SBA-15, and a large amount of nitrogen is observed with functionalized TEPA. Specifically, TEPA@SBA-15 contains 6.7% nitrogen, as shown in Table S1. The nitrogen content is proportional to the TEPA functionalization of NH2-(CH2)2-[NH-(CH2)2]n-NH2, and it is significantly correlated with the CO2 reaction kinetics. Figure S4a shows the N2 adsorption curve of mesoporous silica. The surface area of SBA-15 reaches 619 m2/g, indicating the distributed presence of mesopores. As amine groups are synthesized in the pores, the surface of TEPA@SBA-15 contracts to 188 m2/g. As shown in Figure S4b, the average pore size of SBA-15 is 5.68 nm. TEPA@SBA-15 has a pore size of 5.54 nm, in addition to a significantly reduced peak value.
Figure S5 shows the reaction heat versus the CO2 capture amount. The initial adsorption heat of SBA-15 was 37.1 kJ/mol-CO2, which is consistent with the values reported in related literature19. This corresponds to physisorption, as the heat tends to decrease continuously in terms of the CO2 capture amount. In the case of TEPA@SBA-15, the adsorption heat increases with CO2 capture and presents a peak of 78.5 kJ/mol-CO2. Gas adsorption is initially induced at the region with high surface energy. Following that, the reaction heat increases due to the chemisorption of amine-CO2. The CO2 capture capacity of TEPA@SBA-15 is improved by approximately 70% compared with that of SBA-15, as shown in Fig. 4c.
Table 1 compares TEPA@SBA-15 with the CO2 adsorbents reported in related literature20, 21, 22, 23, 24, 25, 26, 27, 28. Although certain adsorbents have a larger CO2 adsorption capacity than TEPA@SBA-1521, 22, 23, 24, 25, 26, 27, 28, the adsorption heat of TEPA@SBA-15 is significantly higher under similar conditions. Furthermore, considering that a large amount of thermal energy is harvested from a small amount of CO2, the electric energy production can be improved. The CO2 adsorption capacity of TEPA@SBA-15 is not an important parameter as the CO2 capture capacity of TEPA@SBA/MEA + PZ is dominantly affected by the base fluid, while the adsorbents promote the exothermic reaction heat. The desorption of TEPA@SBA-15 can be performed using heat sources above 100 oC, which is identical to base fluids of MEA and PZ.
Figure 4d shows the specific reaction heats of the solvents. The reaction heat of MEA obtained in the experiment was 81.9 kJ/mol-CO2, which deviates an error margin of 2.5% from the theoretical value of 84.0 kJ/mol-CO2. The theoretical CO2 capture capacity of MEA is 8.2 mmol-CO2/g-solvent (0.5 mol-CO2/mol-solvent) considering the bicarbonate reaction. However, the capture amount of MEA is 6.05 mmol/g in a fixed bed reactor, as presented in Fig. 4e. This can be attributed to the fact that CO2 selectivity of the solvent is inversely proportional to the CO2 concentration of flue gas. MEA + PZ induces a chain reaction between the primary and secondary amines, which increases the CO2 capture capacity and reaction heat. In addition, TEPA@SBA-15 adsorbs CO2 chemically on its surface. Finally, TEPA@SBA/MEA + PZ generates a reaction heat that is 52% higher than that generated by the primary amine. However, the CO2 adsorption capacity of TEPA@SBA-15 is smaller than that of the base fluid and does not affect the capture capacity of the entire solvent.
Figure 4f shows the thermal resistance with respect to the reactor shape. The overall thermal resistance of a fixed bed reactor can be estimated by the forced convective heat transfer between the gas (gas mixture) and liquid (solvent). The bubble volume and relative velocity are obtained depending on the porous nozzle diameter, and the heat transfer coefficient is evaluated based on the Nusselt number correlation. In the case of the membrane reactor, the natural convective heat transfer coefficient is obtained by considering the length scale and thermal conductivity of the solvent. In addition, the forced convective heat transfer coefficient is estimated for the gas channel. A comparison between the two reactor shapes show that the overall thermal resistance of the membrane reactor is 3.5 times higher than that of the fixed bed reactor. Specifically, CO2 membrane (thermal barrier) is extremely thin (8 µm), resulting in negligible thermal conduction. However, separating the gas channel and absorbent reservoir significantly increases the thermal resistance. Figure 5a presents the TGA of the CO2 membrane, and the values of 0.1, 0.3, 0.5 represent the weight fraction of the impregnated TEPA. Weight loss occurs at 70–72 oC regardless of the composition, because EtOH was included for uniform dispersion of TEPA. Evaporation of EtOH is independent of thermal stability of the membrane, and it can sufficiently maintain a stable structure at 20–150 oC, which is the operating range of CO2 absorption/regeneration. In Fig. 5b, the membrane can be manufactured in the form of a thin film. CO2 molecules can easily permeate the thin membrane by the solubility difference of amine-phase due to TEPA-impregnated structure. Moreover, the cross-linked resin structure of PEGDMA functions as a thermal barrier between gas and liquid phases by maintaining its solid state. CO2 selectivity and thermal energy harvesting performance can be controlled by varying the thickness of the membrane layer.
3.2. Point-source CO2 capture and thermal energy harvesting
The CO2 capture and thermal energy harvesting performance were evaluated for a CO2 membrane reactor (TEPA@PEGDMA) using a functionalized solvent (TEPA@SBA/MEA + PZ). Figure 6 shows the parameter variations over time during the CO2 capture process. As the inlet pressure was fixed at 3 bar in all cases, CO2 penetrates the membrane with a partial pressure difference and induces a CO2-amine reaction. Figure 6a shows that the initial pressure of the solvent reservoir was 1.75 bar, and the driving force is 7.5 kN, which is sufficient for the permeation of CO2. Here, the driving force is calculated by considering the pressure difference and surface area of the membrane. The steady section is defined as the working region (it is marked in Fig. 6(a)). The amount of CO2 captured in the working region is calculated as the working capacity, as presented in Fig. 6b. As the solvent begins to saturate, the pressure in the reservoir increases and the absorption rate decreases. As the inlet pressure increases to a value greater than 3 bar, the CO2 capture kinetics and capacity remain unaffected. This is because the permeation amount of CO2 does not depend on the partial pressure difference, and the CO2 selectivity of the functionalized solvent is decided by the temperature. However, it must be noted that a stable working region can only be obtained when the minimum inlet pressure is maintained.
Figure 6c presents the temperature variation during the point-source CO2 capture process, for which the measurement locations are marked. The temperatures at the exterior of the reactor and solvent reservoir are similar; the percentage of error between the two values is smaller than 0.2%, indicating that lumped system modeling is reasonable. Fluctuation of the solvent temperature occurs as natural convective motion is promoted by the temperature rise. The average temperature of the membrane surface is 4.5% higher than that of the solvent. In addition, the temperature variation in the longitudinal direction is negligible as CO2 gas is uniformly permeated by the same driving force, and the natural convective motion due to density difference maintains thermal equilibrium within the absorbent reservoir. Furthermore, the inlet gas temperature is maintained at 20 oC. In contrast, the temperature of the outlet gas gradually increases, as a result of the heat loss owing to mass transfer. For MEA in the membrane reactor, the maximum temperature reaches 54.2 oC at an inlet gas flow rate of 5.0 10− 5 kg/s. This is 28.4% higher than the maximum temperature of the MEA (42.2 oC) in the fixed bed reactor.
Figure 7a summarizes the working capacity of the point-source CO2 capture process. Although the working capacity of the membrane reactor is lower than that of the fixed bed reactor (CO2 capture capacity of materials in Fig. 4e), the performance trends are similar in terms of materials. Specifically, TEPA@SBA/MEA + PZ has a higher working capacity than MEA. The CO2 working capacity remains unaffected by the mass flow rate of the inlet gas. Furthermore, the large fluctuation in terms of the time step in the transient analysis can be attributed to the pressure instability caused by membrane penetration. However, the error in reproducibility of the performance remains within 5%, as demonstrated by the results of repeated experiments.
Figure 7b shows the variation in the maximum temperature during the CO2 capture process. The maximum temperature increases in proportion to the mass flow rate of the flue gas because the heat loss decreases as the saturation time is shortened. The forced convective heat transfer increases with the flow rate of the inlet gas, and the thermal resistance of the gas channel decreases. However, the effect of the gas channel is small because the thermal resistance of the solvent reservoir (natural convective heat transfer) is the dominant factor impacting the overall thermal resistance. As the flow rate increases, the heat loss over time is not significantly affected, and the reaction time is shortened to maximize the harvesting of thermal energy. The maximum temperature of TEPA@SBA/MEA + PZ reaches 67 oC at a gas flow rate of 8.33 10− 5 kg/s; the corresponding temperature rise in the fixed bed reactor is 47.5 oC. Figure 7c shows the TEHD calculated based on the temperature variation; the TEHD of TEPA@SBA/MEA + PZ is 15% higher than that of MEA. The TEHD of TEPA@SBA/MEA + PZ in the fixed bed reactor is approximately 200 kJ/kg. In the case of the membrane reactor, the TEHD is 445 kJ/kg at an inlet mass flow rate of 8.33 10− 5 kg/s, which indicates that heat loss was reduced by 245 kJ/kg. As shown in Fig. 7d, the CO2 capture performance of MEA and TEPA@SBA/MEA + PZ remains constant even after repeated absorption/desorption cycles.
3.3. Theoretical performance analysis
The numerical analysis was performed to analyze the working mechanisms of CO2 capture and thermal energy harvesting process in the membrane reactor of TEPA@PEGDMA. Furthermore, if the CO2 capture-driven energy harvesting process can be predicted, the system can be optimized through the parametric analysis. In Fig. 8a, the results of the temperature variation of MEA + PZ deviate from the experimental results by an error of 0.3%. As the CO2 capture kinetics in the working region are steady, the performance is predicted by considering the material characteristics. However, after that, the numerical value deviates from the experimental results. Thus, the data was collected only up to the temperature peak, and the subsequent section were not considered.
Figure 8b shows that TEHD increases in proportion to the mass flow rate as heat loss decreases. The effect of the mass flow rate can be discussed based on the heat transfer mechanisms. The Nusselt number (Nu) correlation for gas channel is given by Eq. (10).
$$\text{N}\text{u}= 2+1.8·{\text{R}\text{e}}^{0.5}·{\text{P}\text{r}}^{0.33}$$
10
The Prandtl number (Pr; defined as the ratio of momentum diffusivity to thermal diffusivity) for flue gas is 0.72 that is similar to ambient air. The Reynolds number (Re; defined as the ratio of inertial force to viscous force) of the gas channel is smaller than 1,000, and corresponds to laminar flow. When the inlet gas flow rate increases from 5 to 8.33 10− 5 kg/s, the heat transfer coefficient changes from 17.9 to 20.5 W/m2 K (thermal resistance changes from 0.12 to 0.095 K/W). However, the overall thermal resistance decreases only by 3.7%, and the heat loss over time is not significantly affected. In addition, the saturation time for CO2 capture process is reduced by 66.6% that corresponds to a heat loss saving of 61.3 kJ/kg.
To summarize, as the inlet mass flow rate increases, a larger amount of thermal energy can be harvested. In contrast, the CO2 concentration at the outlet increases considering that the driving force of CO2 passing through the membrane is constant. For instance, the inlet gas requires a flow rate of 8.33 10− 5 kg/s or less to maintain the outlet concentration at 3 vol%, which is the design condition of the conventional chemical absorption process. An optimal design is required based on the target outlet CO2 concentration.
3.4. Electric energy conversion
Thermal energy was converted to electricity using TE devices composed of Bi2Te3. Figure 9a shows the thermopower for a closed circuit by connecting resistors (the optimal resistors vary in terms of the working temperature range). The temperature of the hot side was varied, while that of the cold side was kept almost constant at 20 oC. The thermopower in working temperature range of 20–70 oC is presented in Fig. 9b, and optimal performance was obtained at a resistance of 5 Ω. Under these conditions, the system generated 0.12–0.94 W of thermopower, which is 44% higher than that generated by a resistance of 1 Ω for the same temperature difference.
Figure 9c shows the thermopower curve for a heat generation system that simulates the thermal energy harvesting process in point-source CO2 capture. As the solvent temperature fluctuates severely during the CO2 capture process, a constant heat generation system is considered to reduce the uncertainty. Thermopower is generated when the temperature difference is greater than 5 oC. The TEPA@SBA/MEA + PZ mixture generated a maximum power of 1.29 W. The average efficiency of thermoelectric conversion is 1.69%, although it varies in terms of materials. Figure 9d shows that the electric energy harvesting density (EEHD) of TEPA@SBA/MEA + PZ is 7.5 kJ/kg, i.e., TEPA@SBA/MEA + PZ can generate 7.5 J/g of electricity while capturing 6.6 mmol/g (0.29 g/g) of CO2. When only thermal energy harvesting is considered instead of thermo-electric conversion, TEHD is 445 kJ/kg-solvent. Furthermore, CO2 regeneration energy of amine absorbent is 4.2 kJ/g-CO2. Considering that 0.29 kg of CO2 is captured per 1 kg of solvent (working capacity), the energy consumption for CO2 regeneration is 1,218 kJ/kg-solvent. Then, the TEHD can cover 36% of thermal energy consumption in the post-combustion CO2 capture process.
3.5. Industrial applications
The industrial applications of point-source CO2 capture-driven electric energy harvesting system were considered. The data were obtained from thermal power plants units 1–7 of Korea Western Power Co., Ltd., located in Taean, South Korea. The CO2 emission of the reference system is 0.8 ton/MWh. Figure 10a shows that the fixed bed reactor using MEA as a working fluid can reduce CO2 emission to 0.18 ton/MWh. The mass transfer area of the membrane reactor is lower than that of the fixed bed reactor, and the outlet CO2 concentration is marginally higher. The CO2 capture capacity of MEA + PZ is higher than that of MEA owing to the effect of the secondary amine. However, the outlet CO2 concentration is also higher owing to the low reaction rate. The CO2 outlet concentration of TEPA@SBA/MEA + PZ is 7.4% lower than that of MEA + PZ. Because amine-functionalized silica is dispersed in the fluid, the CO2 working capacity remains unaffected, and the capture kinetics and selectivity are improved. Note that, the emission per unit energy generation is estimated based on the outlet concentration value of lab-scale results.
Figure 10b shows the contribution of electric energy harvesting to the industrial systems. The fixed bed reactor with MEA does not produce electricity, and the energy consumption for CO2 desorption reaches 4.2 kJ/g-CO2. In other words, from the perspective of energy index, 0.723 MWh of thermal energy is consumed to produce 1 MWh of electricity, which considerably reduces the efficiency of energy production. In contrast, a membrane reactor using MEA as the working fluid can harvest 23 MJ of electricity per 1 MWh of electricity production. In addition, the working fluid also significantly impacts the performance. For instance, TEPA@SBA/MEA + PZ improves the specific reaction heat by 52% compared with MEA and harvests 12.2% more electricity. This corresponds to 0.6% of the total electricity production in a power plant. Furthermore, in typical CO2 capture applications reported in recent studies3, the energy consumption of 4.2 kJ/g is only partially reduced (i.e., 10–20%). However, the system proposed in this study produces electricity, which is highly feasible and significantly increases the energy efficiency.
In a membrane reactor (TEPA@PEGDMA), TEPA@SBA/MEA + PZ absorbs 6.6 mmol/g of CO2 and simultaneously produce 7.5 kJ/kg of electricity. However, the feasibility of the present system is examined for industrial applications by considering additional energy consumption compared to a fixed bed reactor. The flue gas discharged from power plants is at ambient pressure (1 bar). In contrast, the inlet flue gas of the present system must be maintained at 3 bar to permeate the membrane. The work required for gas compression was obtained by applying the isotherm ideal gas equation to \(\int PdV\) where P is pressure and V is volume. When the inlet flow rate is 8.33 10− 5 kg/s, the gas compression work corresponds to 0.87 kJ. The temperature of cold side of the TE devices was maintained by forced convection of air. The work required to operate the fan until the temperature peak is 0.35 kJ (the heat transfer coefficient of air to maintain the temperature of the cold side was designed, and the air velocity was obtained from Nusselt number correlation. Based on the generated kinetic energy of air, isentropic work for operating fan was estimated). The net EEHD reaches 6.28 kJ/kg considering the additional work for present system compared to the fixed bed reactor. Note that, the reported thermo-electric conversion efficiency is 4% in the working temperature range (25–70 oC)29, and the EEHD can be improved by more than 2 times