1. Introduction
It is now widely recognized that climate change, standing as the most pressing challenge of our time, is driven primarily by anthropogenic activities. Among the unequivocal evidence there are rising global temperatures, melting polar ice caps, ocean acidification phenomena, and increasingly erratic weather patterns [
1]. In the pursuit of mitigating these impacts, the scientific community has tirelessly explored various strategies to reduce carbon dioxide (CO
2) emissions such as [
1,
2]:
Improvement of fossil fuel-based energy efficiency;
Enhancement of nuclear and renewable energy as well as increasing use of biofuel-based energy;
Development of environmental engineering works such as afforestation and reforestation.
On the other hand, during this period of ecological transition and contextual energy crisis, the belief of making our current economies independent from fossil fuels seems a utopia as the above mitigation strategies alone may fall short in achieving the necessary emission reductions [
3,
4]. Carbon capture and storage (CCS) would serve as a critical bridge, allowing industries with substantial point-source CO
2 emissions, such as power generation and heavy manufacturing, to significantly curtail their carbon footprint.
After the capture process, CO
2 is usually compressed and transported via pipelines to the storage or utilization facilities. In geographical regions without extensive pipeline infrastructure, CO
2 could be transported in the liquid phase via ships. CO
2 could then be injected and stored underground in geological formations such as depleted oil and gas wells, saline aquifers, and deep coal seams. CO
2 could also react with certain minerals to form stable carbonates, providing a long-term storage solution. Moreover, CO
2 utilization (CCU) involves converting the captured CO
2 into valuable products, such as bulk chemicals, plastics, and synthetic fuels. This strategy, although energy and/or hydrogen intensive, offers advantages in offsetting the cost of carbon capture by creating revenue streams [
2,
3,
4,
5].
The rising emphasis on achieving a net-zero society by mid-century, as articulated in international agreements and embraced by nations and corporations alike, has propelled CCS into the forefront of climate change solutions. Policies and financial incentives aimed at fostering CCS deployment further underscore its integral role in achieving the decarbonization goals [
3,
4,
5]. The key role of CCUS in reducing the CO
2 emissions can be justified by a number of reasons:
So far, diverse CCS technologies have been proposed, however the quest for an efficient, scalable and low-cost solution remains ongoing. Although the deployment of CCS projects has seen a notable progress, the amount of CO
2 captured and sequestered is overall relatively modest. Up to now, only a few Mt of CO
2 have been captured globally, with various projects demonstrating the feasibility of CCS in different industrial settings. However, meeting the scale required for a meaningful climate impact would require a remarkable upscaling of CCS units. The International Energy Agency (IEA) has anticipated a substantial increase in the capture and storage of CO
2, reaching Gt-scale levels by mid-century [
5,
6]. According to the IEA report, the Sustainable Development Scenario will meet the energy-related goals by fully eliminating CO
2 emissions by 2070 [
5,
6]. Butnar et al. [
11] examined different scenarios for the European decarbonization, demonstrating that Europe needs a large-scale CCS industry to achieve the future goals. In the 1.5°C scenario, the average CO
2 captured by CCS would be in the range of 230–430 Mt
CO2 y
–1 by 2030 and 930–1200 Mt
CO2 y
–1 by 2050. Koelbl et al. [
12] compared 18 integrated assessment models showing that, although carbon capture rates will fluctuate in a wide range in the period 2020–2100, no model predicts less than 600 Gt
CO2 of cumulated captured emissions.
The selection of a specific capture technology depends on several factors such as the CO
2 concentration of in the gas stream, the target CO
2 purity and CO
2 recovery, the allowed energy requirements and additional constraints based on the specific industrial process. Three main strategies have been presented as follows. In post-combustion capture, CO
2 is captured from the flue gases arising from the combustion of fossil fuels [
13]. Pre-combustion capture involves capturing CO
2 from H
2-rich gas mixtures generated by the gasification of fossil fuels prior to the combustion [
14]. Oxy-fuel combustion involves burning fossil fuels in an atmosphere enriched with oxygen rather than air [
14]. This results in a flue gas stream predominantly composed of CO
2 and water vapor, making it easier to separate CO
2. A further strategy is offered by chemical looping combustion, which separates CO
2 by using metal oxides as oxygen carriers being circulating between two reactors [
15].
Among the main carbon capture technologies, amine-based absorption has gained prominence for its efficacy in capturing CO
2 from industrial flue gases. Examples of amine-based solvents are monoethanolamine (MEA), diethanolamine (DEA), and proprietary amines. After the absorption process, the CO
2-rich solvent is regenerated by providing heat to release the captured CO
2 [
16]. Alternative solvents such as ionic liquids have also been proposed as they offer some advantages with respect to amines, including low volatility, low regeneration temperatures and potential for tunability [
17]. In pressure swing adsorption (PSA) processes, nanoporous adsorbents, such as zeolites and metal organic frameworks, selectively adsorb CO
2 at high pressure. The adsorbent is then regenerated at low pressures, sometimes under vacuum, to release the captured CO
2 [
18,
19]. Temperature swing adsorption (TSA) processes rely on temperature variations to desorb CO
2 from the adsorbent material. Besides, membrane technology based on polymeric, metallic or ceramic materials, use selective permeable membranes to separate CO
2 from other flue gases by exploiting differences in molecular size, shape, or affinity [
20].
As the demand for more sustainable and energy-efficient solutions grows, attention has turned to cryogenic carbon capture as a promising technology in this field. Cryogenic separation involves cooling the flue gas stream to very low temperatures, causing CO
2 to condense or desublimate while the other gases remain in the gaseous phase. Cryogenics refers most closely to processes that occur at temperatures below 120 K (–153°C), such as the condensation of nitrogen and oxygen; however, the term is often used to indicate in a generic way low-temperature separations [
21]. While the idea of cryogenic separation has been explored for decades, recent advancements and successful pilot projects have reignited interest in its potential as a game-changing carbon capture technology [
21,
22]. The earliest works on cryogenic flue gas treatment were developed by the Bechtel Power Corporation on cryogenic SO
2 capture from a coal plant [
22]. In the last few decades, many patents were filed in relation to CH
4 purification. Holmes and Ryan [
23] proposed a conventional cryogenic distillation for the separation of CH
4 from CO
2 with extractive distillation [
23] using an entrainment solvent (heavy hydrocarbon such as n-butane) to control the CO
2 freezing. Amongst other technologies, the CryoCell® process was developed by Cool Energy and tested in collaboration with other industrial partners including Shell Global Solutions in order to cryogenically remove CO
2 from natural gas [
24].
In the field of post-combustion carbon capture, cryogenic separation is considered energy-intensive but highly effective. The numerous advances in the cryogenic sector have pushed research to also adopt CCS solutions to treat combustion flue gases. However, several works were related to vapor-liquid equilibrium and therefore limited to gas purifications having a high concentration of CO
2 in the feed stream [
25,
26]. For instance, Berstad et al. [
27] analysed low-temperature CO
2 capture processes from IGCC plants by condensation and phase separation, identifing feasible design configurations and determinging a specific energy requirement of around 300 kJ kg
CO2–1 with a CO
2 capture recovery of 85%.
More recently, various scientific works have identified the importance of vapor-solid equilibrium to achieve high CO
2 purities and recoveries, even from flue gases at low pressures and low CO
2 concentrations, and have thus laid the foundation for the thermodynamic understanding and engineering principles of cryogenic desublimation. For instance, some research groups have proposed cryogenic packed beds operating in parallel and undergoing different cycle steps [
28,
29]. Other researchers have developed prototypes and the first pilot plants based on CO
2 desublimation, showcasing the feasibility of implementing cryogenic carbon capture [
30]. Overall, according to the literature reviewed in the present work, cryogenic carbon capture (CCC) based on CO
2 desublimation from post-combustion flue gases offers several advantages compared to other technologies as follows:
The technology features extremely high selectivity, and thus minimizes the competitive capture of other components in the flue gases;
The separation process is driven by temperature differences and does not rely on solvents, adsorbents or membranes;
The CO2 product is generally extracted in the liquid phase, so as to avoid the downstream compression work;
Scalability potential and wide ranges of applications have been demonstrated. In addition, CCC can be integrated into existing industrial processes with minimal retrofit requirements;
CCC is characterized by low water consumption and offers a large heat integration potential. This, in turn, minimizes the disposal of wastes and enhances safety and environmental aspects.
While these advantages make CCC an intriguing option, it is also important to acknowledge that several challenges exist, including the need for very low operating temperatures, the high cost of equipment materials, and the overall feasibility of large-scale implementation. Ongoing research and developments are essential to address these challenges and optimize CCC for widespread application. In the last few years various review articles have been published discussing both modelling and experimental aspects of cryogenic desublimation methods for post-combustion carbon capture [
21,
26,
31,
32,
33].
Shen et al. [
26] focused on the cryogenic capture systems from the perspective of constructing new cryogenic capture system structures, exploring the optimal system parameters, and analyzing the challenges faced. Font-Palma et al. [
21] presented a first comprehensive review on the topic. Song et al. [
31] systematically discussed several CCC technologies, including standalone process schemes and hybrid schemes coupled with cold sources such as LNG regasification [
31]. Asgharian et al. [
32] published a review paper on process modeling of CCC, focusing on thermodynamics and energy-related aspects. Aneesh and Sam [
33] evaluated additional theoretical aspects of CCC technologies such as mass transfer mechanisms.
All the above works often focussed on a limited research area, either thermodynamic aspects or process configurations. Furthermore, these works often included process schemes that are based on vapor-liquid equilibrium and/or hybrid schemes to purify CH4 rather than capturing CO2. This means that the presented results are heterogeneous and do not represent a systematic classification of the processes based on cryogenic desublimation. In addition, desublimation-based technology is experiencing a rising interest with the first demonstration pilot plants being built and operating, which could confirm the theoretical and modelling aspects investigated in the previous papers.
Based on these considerations, a novel and comprehensive review work on this topic is necessary. This paper delves into the evolution of cryogenic carbon capture by desublimation, emphasizing both strengths and limitations of the latest process schemes. It meticulously explores the engineering principles behind desublimation, ranging from thermodynamics to mass transfer phenomena. By critically evaluating the achievements of pilot plants and commercial solutions, we aim to elucidate the transformative potential of cryogenic carbon capture, shedding light on recent breakthroughs and real-world applications. This will be achieved by discussing the main process schemes and comparing the key performance indicators of the proposed technologies.
2. Fundamentals of Cryogenic Desublimation
To predict the phase equilibria involved in cryogenic desublimation processes, cubic equations of state (EoS) including Soave–Redlich–Kwong (SRK), Peng–Robinson (PR) and Patel–Teja (PT) can be used coupled with a relatively small amount of experimental data to determine their mixing rules parameters [
34]. Several works studied theoretically and experimentally CO
2 desublimation in the oil and gas industry. De Guido et al. [
35] demonstrated the reliability of their thermodynamic modelling in representing solid-liquid-vapor equilibrium (SLVE) with a focus on obtaining CO
2 solubility predictions in hydrocarbon-rich mixtures (notably for natural gas purification). The same authors also described the solid-vapor equilibrium (SVE) for CO
2-CH
4 and CO
2-CH
4-N
2 mixtures applying both selected cubic EoS and the Gibbs energy minimization method [
36]. A parallel field of studies was also related to the prediction of the Joule-Thomson effect, where Wang et al. [
37] developed and validated an improved 25-parameter model. Below is reported the most commonly used set of equations to describe SVE of CO
2 in a flue gas mixture as well as to predict CO
2 dew and freezing points. In particular, Equation (1) represents the isofugacity condition for the SVE, while the PR EoS shown in Equation (2) is used to calculate the fugacity coefficients of CO
2 in both solid and vapor phases. The sublimation vapor pressure of CO
2 can be determined by Equation (3), as tested by Jensen et al. [
30].
In the above equations
,
,
,
,
and
represent the mole fraction of CO
2 in the vapor phase, the fugacity coefficient of CO
2 in the vapor phase, the total pressure, the sublimation vapor pressure of CO
2, the fugacity coefficient of CO
2 at the saturated solid pressure and the molar volume of solid CO
2, respectively. It should be noted that in Equation (1) the solid phase is represented by pure CO
2. Several modifications of the above equations are reported in the literature, especially for the binary interaction parameters of the cubic EoS [
30,
38,
39]. In the works of both Jensen et al. [
30] and Pellegrini et al. [
38], SVE calculations have been implemented by minimizing the Gibbs energy of the system. Most theoretical data have also been validated through dedicated experimental campaigns [
30,
38].
In this work similar results based on SVE were achieved using Aspen Plus V10.0 software. To obtain these data the RGibbs reactor was selected. The CO2 solid phase is formed as a result of a chemical equilibrium reaction from the gas phase. The results were obtained for binary systems of CO2 and N2 considering two CO2 concentrations in the feed: 5 mol% as representative for natural gas-fired power plant flue gases and 15 mol% as representative for coal-fired power plant flue gases.
Figure 1a shows the CO
2 recovery in the solid phase against the flue gas temperature according to the SVE model at three different pressure levels and two inlet CO
2 concentrations.
Figure 1b shows the flue gas temperature required to obtain a certain CO
2 recovery in the solid phase at a given pressure starting with a flue gas having a CO
2 composition of 15%. In
Figure 1b simulation results are compared against the calculations carried out by Pellegrini et al. [
38] and Baxter et al. [
40]. In case of low CO
2 concentrations (5%) in the feed, high CO
2 recoveries are achievable only by operating at temperatures below 150 K or considering flue gases at higher pressures. However, these temperature levels could be raised to around 170 K if the flue gases contain 15% of CO
2 at relatively low pressures.
Overall, it is possible to obtain a minimum work of separation in the order of 1 MJ kg
CO2–1 at these temperatures and for inlet CO
2 concentrations of around 15%. Yurata et al. [
41] carried out process simulations involving CO₂-H₂ mixtures demonstrating the thermodynamic superiority of desublimation processes against other unit operations. Energy consumption values of 1–2 MJ kg
CO2–1 represent an enormous advantage compared to the state of the art of technology, especially considering that for many cryogenic solutions the compression/liquefaction work necessary for the CO
2 transport is minimized.
According to chemical thermodynamics, the minimum work of separation required to separate CO
2 from the other flue gases (inert) involves overcoming the thermodynamic barrier associated to the gas mixing. This corresponds to the opposite of the Gibbs energy of mixing between two products (
p1,
p2) and the feed (
f) in a generic separator, as indicated by Equation (4) [
13,
42]. In addition, the subsequent compression and/or liquefaction of CO
2 is a crucial aspect in evaluating the efficiency and feasibility of carbon capture processes [
13,
22]. A more accurate calculation can be attained by including the theoretical work of compression and/or liquefaction to deliver CO
2 in a dense phase. For instance, purified CO
2 must be usually compressed to very high pressures to overcome the pressure drops in the pipelines and to be injected in the form of supercritical fluid into reservoirs at depths of the km order of magnitude. The minimum work required to compress the purified CO
2 is calculated by Equation (5) and corresponds to a reversible isothermal compression. Considering a hydrostatic pressure of roughly 1 MPa per 100 m of depth, the battery limits for CO
2 delivery would be around 10–15 MPa. Alternatively, a liquefaction process can be considered by including the theoretical heat removal by liquefaction in a reversible refrigeration cycle, as reported in Equation (6).
In chemical thermodynamics, the minimum work of separation required to separate CO
2 from the other flue gases (inert) involves overcoming the thermodynamic barrier associated to the gas mixing. This corresponds to the opposite of the Gibbs energy of mixing between two products (
p1,
p2) and the feed (
f) in a generic separator, as indicated by Equation (4) [
13,
42]. In addition, the subsequent compression and/or liquefaction of CO
2 is a crucial aspect in evaluating the efficiency and feasibility of carbon capture processes [
13,
22]. A more accurate calculation can be attained by including the theoretical work of compression and/or liquefaction to deliver CO
2 in a dense phase. For instance, purified CO
2 must be usually compressed to very high pressures to overcome the pressure drops in the pipelines and to be injected in the form of supercritical fluid into reservoirs at depths of the km order of magnitude. The minimum work required to compress the purified CO
2 is calculated by Equation (5) and corresponds to a reversible isothermal compression. Considering a hydrostatic pressure of roughly 1 MPa per 100 m of depth, the battery limits for CO
2 delivery would be around 10–15 MPa. Alternatively, a liquefaction process can be considered by including the theoretical heat removal by liquefaction in a reversible refrigeration cycle, as reported in Equation (6).
Where the specific thermal duty for condensation
coincides to the latent heat of vaporization. The minimum work of separation is reported in
Figure 2 for an inlet CO
2 concentration of 15% as a function of the CO
2 recovery in the solid phase. The calculated value for separating 90% of CO
2 into a pure CO
2 stream and a second N
2-enriched stream is 140 kJ
e kg
CO2–1. By considering a liquefaction occurring at 50 bar, and a work of compression also at 50 bar, the resulting total minimum work of separation is around 400 kJ
e kg
CO2–1. In fact, it has been reported that optimized compression and liquefaction with external refrigeration (e.g., ammonia refrigeration cycle) consume about 100 kWh t
CO2–1 [
43]. These values can be compared to the actual work of separation, which can be expressed through the II law efficiency shown in Equation (7). Quantifying the contributions from an energy point of view allows us to use a second law efficiency and therefore to compare a theoretical separation work with the thermal or mechanical energy used in a real system. In fact, the numerator of Equation (7) represents the useful product while the denominator represents the equivalent work including the transformation of thermal energy into mechanical energy (higher quality). More refined methods are presented in the cited literature [
13,
44,
45]. The actual work of separation could then represent a better comparison with the real energy penalties obtained in the CCS projects deployed worldwide.
To account for the energy consumption in a cryogenic desublimator we can refer to the schematic shown in
Figure 3 adapted from Swanson et al. [
22], who performed theoretical calculations for energy penalties associated to CCC. In the figure, the flue gas is pre-cooled by recovering the internal cold energy, including the cold CO
2 product and the cold CO
2-lean flue gas streams. Recuperative heat exchangers are expected to pre-cool the incoming flue gas. Depending on the operating conditions, the recovery ration and the selected temperature pinch within the recuperative heat exchanger, the process may require an external refrigeration cycle. As can be seen from the figure, there are two types of heat removal strategies: i) cold energy recovery from internal process streams such as the clean gas or the recovered liquid CO
2; ii) refrigeration by expending mechanical/electrical energy, i.e. using a refrigeration cycle. According to Swanson et al. [
22], and as reported in some real process schemes, two types of refrigeration cycles could be implemented: one with a higher coefficient of performance (COP) rejecting the heat to an internal stream (e.g., the solid CO
2 is used as a low-temperature heat sink for a portion of the heat that must be rejected by the refrigeration system), and another one with a lower COP rejecting the heat at ambient temperature. Rejecting the heat to a cold stream is less costly than the rejecting the heat at ambient temperature. However, the choice of the refrigeration cycle will depend on the equipment design and the dichotomy between capital costs and operational costs due to the low temperatures involved.
Downstream of the separation, the energy costs to compress CO
2 to the required pipeline pressures is usually very limited, since CO
2 is already in the liquid phase. For instance, when pressurizing 1 kg of liquid CO
2 at 220 K from atmospheric pressure to 120 bar, the resulting enthalpy change is only 4 kJ. This leads to an overall energy consumption of around 450 kJ kg
CO2–1 with 90% CO
2 recovery, which is comparable to the total minimum works of separation shown in
Figure 2. Although the real desublimation processes described in the following sections disclosed slightly higher values of energy consumption, cryogenic processing can be considered very close to the thermodynamic optimum. For instance, in amine-based absorption systems, the energy required for solvent regeneration constitutes a significant portion of the overall work of separation and, if summed to the compression and transportation energy requirements, can total approximately 4–6 MJ kg
CO2–1 [
13,
16]. This means that, compared to the benchmark processes, cryogenic processes can push the II law efficiency of carbon capture from around 20% to more than 50%.
Desublimation can occur from the gaseous bulk phase towards a cold surface or through direct contact with a cold liquid (cryogenic fluid). From the point of view of the design of unit operations, the flue gas can be treated as a stream of non-condensable gases saturated with CO
2 and the formalism presented by classic chemical engineering textbooks can be used to obtain the basic design of the equipment for heat and mass transfer [
42]. The heat and mass balances to determine the gas temperature
TG and the gas composition
yCO2 along the axial dimension
z can be drawn according to the film theory for gas desublimation in a mixture of non-condensable gases in countercurrent contact with a cold liquid, as exhibited in
Figure 4. The figure shows the schematics of the film (considering that the liquid cannot vaporize) and the control volumes of the gas flowing from bottom to top and the liquid flowing from top to bottom of a desublimating separator considering 1-dimensional case of a direct-contact desublimating separator. The CO
2 flux occurs through the interface and the separated solid CO
2 fall down into the descending liquid phase that is not volatile (not present in the gaseous phase).
Heat and mass balances can be generalized by applying them to a control volume between the flue gas mixture from which CO
2 is separated and the cryogenic liquid. Considering a separator section of height
dz the following equations represent the heat fluxes exchanged for sensible heat (Equation (8)) and latent heat (Equation (9)) by the gas phase. Equation (10) describe the heat flux transferred to the cryogenic liquid. These equations are also valid for a heat exchanger with a wall between liquid and gas phases under the assumption of the wall having a negligible thermal resistance. By neglecting the sensible heat of solid CO
2 and assuming that the specific surface related to mass transfer coincides with the specific surface related to heat transfer, the balance of fluxes
=
can be used to calculate the interface temperature.
The mathematical modelling and the related algorithms to define the temperature – concentration profile of a desublimator are similar to those determined for vapor mixtures containing both condensable and non-condensable gases. Simplified methodologies to calculate the direct contact heat and mass transfer apparatus can be found for systems characterized by a Lewis number approximately equal to 1 (air-water humidification and dehumidification systems) [
42]. In this case (CO
2-flue gases-cryogenic liquid), the heat and material balances should be numerically integrated to determine the temperature of the phase interface as well as the temperatures of the liquid and gas phases as well as the CO
2 concentration profile. In the field of desublimation there are no complete theoretical works combining thermodynamic aspects and heat and mass transfer phenomena. A few modelling works are available for desublimation equipment although most of these works have not been fully validated against experimental data. Only Asgharian et al. [
32] provided a detailed numerical framework for the modelling of the main components constituting a cryogenic carbon capture process, including CO
2 separators, storage tanks, heat exchangers and turbomachinery.
Table 1 details the main modelling features of literature works related to CCC by desublimation.