Optimization of BOG Reliquefaction Process of Carbon Dioxide Considering Nitrogen Content
Article information
Abstract
As the importance of carbon capture and storage (CCS) technology increases for greenhouse gas reduction, CO2 transportation technology is also becoming crucial. Efficiently handling boil-off gas (BOG) generated during voyages is particularly important. This paper proposes a process incorporating two-stage separation and mixed refrigerants to reliquefy CO2 BOG containing nitrogen efficiently. This process was optimized based on specific power consumption (SPC) criteria and compared with the conventional single-stage separation using an ammonia refrigerant. The two-stage separation allows the removal of non-condensable gases, such as nitrogen, before expansion, and the use of mixed refrigerants enables more efficient heat exchange than ammonia refrigerants, improving the reliquefaction rate. As a result, the proposed process reduced SPC by up to 8.8 % with a nitrogen content of 5 mol% and up to 74.7 % with a nitrogen content of 15 mol%. In addition, the proposed process achieved a reliquefaction rate of over 74.2 % across all nitrogen content ranges of 5–15 mol%.
Nomenclature
Abbreviation
CCS: Carbon capture and storage
BOG: Boil-off gas
LNG: Liquefied natural gas
PT diagram: Pressure temperature diagram
EOS: Equation of state
CPA: Cubic plus association
SPC: Specific power consumption
LCO2: Liquefied CO2
BOR: Boil-off rate
PSO: Particle swarm optimization
JT: Joule Thomson
Symbols
Wtotal: Total energy consumption (kW)
ṁ: Mass flow rate (t/h)
r: Return LCO2 fraction
V: Volume (m3)
v: Molar volume (m3/kgmole)
x: Vapor mole fraction
PRComp: Pressure ratio of the compressor
P: Pressure (bar)
T: Temperature (°C)
C2: Ethane
C3: Propane
nC4: Normal butane
iC4: Iso butane
1. Introduction
With the adoption of the Paris Agreement in 2015, which aimed to limit the increase in average global temperature to below 2 °C compared to pre-industrial levels and restrict it to no more than 1.5 °C, reducing greenhouse gas emissions has become a significant global issue in response to rapid climate change caused by global warming. IRENA (2022) mentions carbon capture and storage (CCS) as one of the key technologies for achieving greenhouse gas emission reduction targets and expects it to account for approximately 20% of greenhouse gas reductions by 2050. Therefore, the need for research to effectively capture, transport, and store carbon is increasing.
Pipelines and liquefied CO2 carriers can be considered for transporting captured CO2 to storage sites efficiently. In particular, shipping was reported to be more economical for long-distance transport. In the case of shipping, liquefaction reduces the volume of CO2, allowing more mass of CO2 to be transported per unit volume. Decarre et al. (2010) reported that transport by ship is economically more advantageous than pipeline transport when the distance to the storage site exceeds 1,000 km.
The problem of boil-off gas (BOG) can occur while transporting liquefied CO2 by ship because liquefied CO2 must be stored at temperatures lower than room temperature. Although onboard storage tanks are insulated, they cannot block heat ingress completely. Therefore, liquefied CO2 evaporates due to external heat influx, and the resulting BOG increases the pressure in the storage tank, potentially causing safety accidents. If the generated BOG is vented to the atmosphere to prevent overpressure, it results in cargo loss and greenhouse gas emissions. A BOG reliquefaction process to prevent this can be crucial in minimizing greenhouse gas emissions and ensuring the safety of the CCS supply chain.
Existing studies have mainly focused on the liquefaction process of captured CO2. Alabdulkarem et al. (2012) proposed a new liquefaction system applying various refrigerants to improve the efficiency of the CO2 liquefaction process for CCS. Lee et al. (2015) modeled the process for liquefying captured CO2 and studied changes in process efficiency according to seawater temperature variations Seo et al. (2015) used ammonia, propane, and R-134a as refrigerants in the CO2 liquefaction process. They reported that ammonia showed the best efficiency. Seo et al. (2016) modeled the Linde–Hampson process and a process using external refrigerants to liquefy captured CO2 within the CCS value chain. They conducted a cost analysis according to liquefaction pressure. Deng et al. (2019) modeled a liquefaction process using ammonia refrigerant and performed economic evaluations of captured CO2 based on different scenarios. Aliyon et al. (2020) modeled and analyzed energy, exergy, and economic aspects of the Linde–Hampson process, a process using ammonia refrigerant and an absorption refrigeration process for the liquefaction of captured CO2. Nevertheless, these studies focused only on the CO2 liquefaction process and did not consider the reliquefaction process of CO2 BOG occurring during ship transport.
Some studies on CO2 BOG reliquefaction systems have been conducted but have limitations. Yoo (2017) studied a CO2 BOG reliquefaction system using LNG cold energy on an LNG-fueled CO2 carrier. Lee et al. (2017a) proposed a CO2 BOG reliquefaction process and presented optimal designs considering the seawater temperature and compressor discharge temperature. Lee et al. (2024) modeled a CO2 BOG reliquefaction system using the Linde–Hampson process, a process using ammonia refrigerant, and a subcooling process, comparing them based on the liquefaction pressure. On the other hand, these studies have the limitation of assuming CO2 BOG as pure CO2. In the CO2 capture process, various impurities can be mixed into the liquid CO2. According to Jeon et al. (2015), substances with low boiling points as impurities, such as nitrogen, rapidly vaporize within the liquid CO2 storage tank, making CO2 BOG a mixture containing more impurities than the liquid form. CO2 BOG mixed with nitrogen requires lower temperatures for liquefaction; therefore, an effective reliquefaction process capable of cooling to lower temperatures than pure CO2 is needed. This leads to increased energy consumption and costs, making research on efficient reliquefaction methods for impurity-containing CO2 BOG an important topic in liquefied CO2 transport.
Some studies examined reliquefying impurity-containing CO2 BOG, but they had limitations. Chu et al. (2012) reported that BOG generated during the transport of captured CO2 from thermal power plants contains impurities, such as nitrogen and oxygen, and proposed an ethylene–propane cascade cycle as a reliquefaction process. Lee et al. (2017b) proposed and compared four processes using CO2 as a refrigerant to reliquefy CO2 BOG containing 36 mol% nitrogen. Using the cascade cycle increased the process complexity because separate compression trains were required for each refrigerant. Cooling below −65.6 °C was necessary to achieve a reliquefaction rate of over 98% for impurity-containing BOG, which poses a risk of dry ice formation. In addition, using CO2 as a refrigerant required compressing CO2 to high pressures exceeding 100 bar compared to external refrigerants. To achieve a CO2 reliquefaction rate of 90%, they proposed a process applying a distillation column to separate nitrogen and CO2, but this increased energy consumption by approximately 1.5 fold compared to other processes.
The use of mixed refrigerants in the liquefaction of mixtures has also been studied widely. Existing low-temperature liquefaction processes, such as natural gas liquefaction and LNG BOG reliquefaction, use liquefaction processes using mixed refrigerants, reporting higher efficiency compared to single-refrigerant liquefaction processes. Using mixed refrigerants in natural gas liquefaction can increase efficiency compared to nitrogen single-refrigerant cycles or cascade cycles by reducing the temperature difference between hot and cold streams in the heat exchanger (Moein et al., 2015; Zhang et al., 2020). According to Lee et al. (2023), using mixed refrigerants in LNG BOG reliquefaction reduced power consumption by 30% compared to using single nitrogen refrigerants. Xu et al. (2013) reported that the composition of mixed refrigerants plays a significant role in the heat exchange efficiency of the process of natural gas liquefaction. On the other hand, such processes using mixed refrigerants have mainly been applied to natural gas liquefaction, and studies optimizing the composition of mixed refrigerants for impurity-containing CO2 BOG mixtures are difficult to find.
This paper proposes an optimal process design that enables efficient BOG reliquefaction even under temperature constraints to avoid dry ice formation by incorporating a mixed refrigerant liquefaction process into the reliquefaction of impurity-containing CO2 BOG mixtures. The CO2 BOG reliquefaction process using ammonia single refrigerant proposed by Seo et al. (2015) was compared. In the reliquefaction process, mixed refrigerants were applied, and a process incorporating a vapor–liquid separator before expansion was proposed to improve the reliquefaction rate.
2. Background
2.1 Characteristics of Carbon Dioxide (CO2) BOG
Fig.1 presents a pressure-temperature (PT) diagram generated using the commercial software Multiflash V7.0 for CO2 containing nitrogen impurities ranging from 0 to 15 mol%. The equation of state (EOS) used to plot the graph was CPA-Infochem, which has been reported to predict the phase behavior of CO2 mixtures accurately (Tsivintzelis and Kontogeorgis, 2015).
Pure CO2 can be liquefied and stored at various pressures between its triple point (5.1 bar, −56.6 °C) and critical point (73.8 bar, 31.1 °C). Seo et al. (2016) reported that storage pressures above 25 bar are uneconomical when evaluating the cost-effectiveness at various pressures within the CCS value chain, including CO2 liquefaction systems, storage tanks, and CO2 carriers. The most commonly mentioned storage pressures between the triple point and 25 bar are two conditions: pressures near the triple point of CO2 at 6.5–7.0 bar (Chu et al., 2012; Lu et al., 2023; Yoo, 2017) or at 15 bar (Jackson and Brodal, 2019; Lee et al., 2024).
In addition, the dew point temperature decreases as nitrogen impurities increase, requiring lower temperatures for liquefaction. On the other hand, dry ice can form at low temperatures near the triple point. Therefore, in the design of the liquefaction process, care must be taken to ensure that the operating temperature does not fall below the gas–solid phase transition temperature. Consequently, it is necessary to consider this when designing the reliquefaction process and develop a design that can achieve a low dew point temperature while avoiding the range where dry ice forms.
2.2 Carbon Dioxide Liquefaction Process
Fig. 2 shows the single-refrigerant liquefaction processes commonly used in CO2 liquefaction. CO2 is compressed to high pressure through a CO2 compressor and then cooled to low temperatures via heat exchange with a refrigerant in a heat exchanger (Fig. 2(a)). Subsequently, it undergoes partial liquefaction by expanding to the tank pressure through a Joule–Thomson valve. The gas and liquid are separated in a vapor–liquid separator; the gas returns to the front end of the compressor after heat exchange, and the liquid flows into the CO2 storage tank. Fig. 2(b) presents a process where low-temperature refrigerant is produced through a refrigeration cycle applying three-stage expansion, and liquefied CO2 is produced by heat exchange with CO2. Although the process becomes more complex than Fig. 2(a) due to multi-stage expansion in the refrigeration cycle, energy consumption can be reduced (Chen and Morosuk, 2021). On the other hand, these conventional reliquefaction processes are designed assuming the liquefaction of pure CO2. When liquefying CO2 containing impurities such as nitrogen, a lower liquefaction temperature is required compared to pure CO2, so it is necessary to design a separate process.
2.3 Efficiency Indices of Carbon Dioxide (CO2) BOG Reliquefaction Processes
There are various indices to evaluate the performance of CO2 BOG reliquefaction processes, with Specific Power Consumption (SPC) being widely used (Lee et al., 2017b; Lee et al., 2024). The SPC refers to the total energy consumed to produce one ton (t) of liquefied CO2 and was calculated as shown in Eq. (1)
Wtotal is the total power consumption (kW) of the compressors within the reliquefaction system, and ṁreliq is the amount of reliquefied CO2 (t/h). In this case, the unit of SPC is kWh/tCO2. The lower the total power consumption relative to the amount of liquefied CO2, i.e., a lower SPC value indicates a better reliquefaction performance.
The CO2 reliquefaction rate (Return LCO2 fraction, r) is used to assess the proportion of liquefied CO2 resulting from the reliquefaction process and was calculated using Eq. (2) (Lee et al., 2017b).
2.4 Particle Swarm Optimization (PSO)
PSO is an optimization method developed by Kennedy and Eberhart (1995) that does not require differentiation and can solve non-convex problems. The PSO algorithm starts with randomly initialized particles that update their positions and velocities to find the optimal position. Each particle individually tracks its optimal solution and leverages the discoveries of other particles, allowing the entire swarm to converge quickly on the optimal position. Through this approach, various design variables can be optimized simultaneously, effectively solving nonlinear and multidimensional problems in liquefaction processes (Khan and Lee, 2013). This PSO method has been used in many optimization studies of liquefaction processes (Khan and Lee, 2013; Lee et al., 2024).
3. Process Simulation Conditions
3.1 Carbon Dioxide Storage Conditions
Among the storage pressures for CO2, the most frequently mentioned are 6.5 bar and 15 bar. Considering the density of liquefied CO2 in saturated conditions, the density was 1154.6 kg/m3 and 1031.7 kg/m3 at 6.5 bar and 15 bar, respectively. This means that at 6.5 bar, the amount of liquefied CO2 that can be stored per unit storage volume increased by approximately 12% compared to 15 bar (Jeon et al., 2015). For this reason, the storage pressure of the CO2 tank was assumed to be 6.5 bar in this paper.
The unit capacity of the CO2 storage tank applied in this study was set to 100,000 m3 following the assumption of Chu et al. (2012), and it was assumed that the tank is filled up to an initial liquid level of 95% and operated. The boil-off rate (BOR), which is the amount of BOG generated during one day, is typically evaluated at 0.12%–0.2% (Chu et al., 2012; Yoo, 2017). In this study, a BOR of 0.13% was assumed, resulting in a BOG generation amount of 6,000 kg/h. The pressure of the generated BOG was assumed to be the same as the storage pressure, 6.5 bar.
The saturation temperature of CO2 at 6.5 bar was −50.3 °C. However, in the case of large-scale storage tanks for liquefied cryogenic fluids, heated gas accumulates at the top due to external heat ingress, existing as superheated gas. Therefore, it was assumed that the temperature of the BOG increased to −28.5 °C when entering the reliquefaction facility (Chu et al., 2012).
3.2 Carbon Dioxide BOG Composition
The CO2 loaded onto the ship exists as both boil-off gas and liquid CO2 after loading, and the volumetric ratio of the gas space in the storage tank is called the ullage. The ullage can be estimated based on the volumetric ratio of liquid and gas using Eq. (3).
3.3 Carbon Dioxide BOG Reliquefaction Process
Fig. 4 presents process flow diagrams of the CO2 BOG reliquefaction processes considered in this study. In Cases 1, 2, and 3, the CO2 BOG is first compressed to high pressure and then cooled to low temperatures through heat exchange with a refrigerant in a heat exchanger. Subsequently, it undergoes partial liquefaction by expanding to the tank pressure via a valve, and it is separated into gas and liquid phases in a vapor–liquid separator. The non-condensable gases are then vented to the atmosphere, and the liquefied CO2 is returned to the tank.
For Cases 1 and 2, ammonia refrigerant applied widely in CO2 liquefaction process studies was used as the refrigerant. Case 3 used a mixed refrigerant as the refrigerant, selecting ethane (C2), propane (C3), normal butane (nC4), and isobutane (iC4) as the components. The boiling point of pure CO2 at 6.5 bar is −50.3 °C, and propane, with a boiling point of approximately −42.1 °C at atmospheric pressure, is sufficient to liquefy it. On the other hand, a lower cooling temperature may be necessary for liquefaction for mixtures containing substances with lower boiling points, such as nitrogen. Therefore, ethane with a boiling point of approximately −88.5 °C at atmospheric pressure was included. Furthermore, normal butane and isobutane, with boiling points of approximately 0.5 °C and −11.7 °C, respectively, were included to improve the heat exchange efficiency in the high-temperature section of the heat exchanger.
Cases 4, 5, and 6 are process diagrams where an additional separator was installed in Cases 1, 2, and 3, respectively, to enhance the separation efficiency. Installing the additional separator serves to pre-release nitrogen, which vaporizes more readily than CO2, sending high-purity CO2 to the valve.
The Peng–Robinson equation of state was used to calculate the thermodynamic properties according to the pressure of liquid CO2 (LCO2). Each process was simulated using the commercial software Aspen HYSYS V11 under the conditions specified in Table 2 to ensure an equivalent comparison.
3.4 Optimization Method for the Carbon Dioxide BOG Reliquefaction Process
Optimization was performed using the Particle Swarm Optimization (PSO) algorithm to find the optimal operating conditions for each reliquefaction process and nitrogen content case. The Specific Power Consumption (SPC)―the energy consumption per unit production of liquefied CO2― was set as the objective function to be minimized to achieve the optimal performance of the CO2 reliquefaction process.
While satisfying the conditions in Table 2, optimization constraints were established to ensure that the temperature within the process remained above −55 °C even after the Joule–Thomson expansion of CO2 to prevent liquid ingress and dry ice formation. If the constraints are violated, a penalty value is assigned to prevent the solution from being selected as the optimal value. The minimum SPC of the reliquefaction process was determined by considering key process variables, such as temperature, pressure, and refrigerant composition, as optimization variables for the reliquefaction process. These are summarized in Table 3.
4. Results
Among the optimization results of the reliquefaction systems, Table A1, A2 and Table 4 lists the cases where the nitrogen content in the CO2 BOG was 5 mol% and 15 mol%. According to Table 4, when the nitrogen content in the CO2 BOG was varied from 5–15 mol%, Cases 1, 2, and 3 showed a decrease in the total power consumption of the compressors as the compression ratio of the second compressor and the refrigerant flow rate decreased. In addition, the amount of reliquefied CO2 decreased from 4,030, 4,080, and 4,296 kg/h at 5 mol% nitrogen to 570 kg/h at 15 mol% nitrogen. This is because the amount of cooling in the heat exchanger decreases, and the purge amount increases after Joule–Thomson expansion.
Fig. 5 shows the vapor fraction values of the stream after the Joule–Thomson expansion in each reliquefaction process. The vapor fraction after expansion increased as the nitrogen content in CO2 increased, reaching up to 0.91 in Cases 1, 2, and 3. This was attributed to the minimum temperature constraint in the CO2 reliquefaction process. If more nitrogen with a lower saturation temperature than CO2 (−175.7 °C at 6.5 bar) is included, the phase equilibrium of the mixture shifts, resulting in the need for a lower temperature to produce the same amount of liquid under the same pressure conditions. In the BOG reliquefaction process, however, the lower temperature limit must be constrained to prevent dry ice formation. In this study, the formation of liquid after expansion was limited because the lower limit temperature is restricted to −55 °C.
Therefore, in Cases 1, 2, and 3, as the nitrogen content increases, most of the CO2 is purged due to the temperature constraint, resulting in decreases in the compression ratio of the second compressor, the refrigerant flow rate, and the amount of reliquefied CO2. Consequently, in Cases 1, 2, and 3, when the nitrogen content varies from 5–15 mol%, the reliquefaction rate decreases from 69.4–74.0% to 10.6% (Fig. 6). Accordingly, although the total power consumption of the compressors decreased, the SPC increased exponentially from 87.5, 84.3, and 80.9 kWh/t to 380.0, 372.0, and 330.17 kWh/t, respectively (Fig. 6).
On the other hand, in Cases 4, 5, and 6, when the nitrogen content in the CO2 BOG was varied from 5 mol% to 15 mol%, the SPC increased from 86.2, 82.9, and 79.8 kWh/t to 112.6, 109.3, and 96.2 kWh/t, respectively. On the other hand, the increase was relatively low, only approximately 20–30% (Fig. 6). This is because, in Cases 4, 5, and 6, the application of two-stage vapor–liquid separation keeps the vapor fraction after expansion low, at approximately 0.11–0.22, even though the total power consumption of the compressors increases to meet the lower saturation temperatures required as the nitrogen content rises (Fig. 5). This allows for an increased amount of reliquefied CO2 compared to Cases 1, 2, and 3.
By applying two-stage vapor–liquid separation, lighter nitrogen gas can be separated from CO2 before expansion, allowing for Joule–Thomson expansion of CO2 with relatively higher purity. This makes it possible to increase the amount of reliquefied CO2 even at the limited temperature of −55 °C compared to Cases 1, 2, and 3. Consequently, a reliquefaction rate of over 63.0% was achieved in Cases 4 and 5, even when the nitrogen content was increased to 15 mol%. In Case 6, which uses mixed refrigerants, a reliquefaction rate of 74.2% can be achieved under the same conditions.
Reliquefaction is possible when using mixed refrigerants because of the more efficient heat exchange. Fig. 7 shows the thermal composite curves of Cases 4 and 6 when the nitrogen content is 15 mol%. The temperature difference between the hot and cold composite curves is, on average, 18 °C with a maximum of 38.5 °C in Case 4, whereas in Case 6, it is, on average, 11 °C with a maximum of 32. 2 °C. This suggests that more efficient heat exchange is occurring in Case 6 using mixed refrigerants. This efficiency is achieved because mixed refrigerants composed of substances with different boiling points can utilize the latent heat of each refrigerant in the BOG reliquefaction process.
As a result, the proposed two-stage vapor–liquid separation process increased the amount of liquefied CO2 compared to the single vapor–liquid separation process. When the nitrogen content changed from 5 mol% to 15 mol%, the reliquefaction rate of Case 1, which uses ammonia refrigerant in a single vapor–liquid separation process, decreased sharply from 69.4% to 10.6%. In contrast, Case 6, which applies mixed refrigerants and a two-stage vapor–liquid separation process, achieved a reliquefaction rate of over 74.2% in all cases.
Maintaining the reliquefaction rate also led to significant variations in the energy consumption required for the liquefaction process. When the nitrogen content was 5 mol%, the SPC of Case 1 was 87.5 kWh/t, whereas the SPC was 79.8 kWh/t in Case 6, showing an 8.8% reduction in SPC. As the nitrogen content increased, the effect of the two-stage separation process became more pronounced. When the nitrogen content was 15 mol%, the SPC of Case 6 was 96.2 kWh/t, showing a 74.7% reduction compared to the SPC of 380.0 kWh/t in Case 1.
5. Conclusion
Optimization was performed to minimize the specific power consumption (SPC) for cases where the nitrogen content in CO2 BOG varies from 5.0 to 15.0 mol%, considering the use of ammonia as a refrigerant (Cases 1 and 2) and the use of mixed refrigerants (Case 3). In the reliquefaction processes applying single-stage vapor–liquid separation (Cases 1, 2, and 3), the amount of reliquefied CO2 decreased sharply as the impurity content increased, leading to a significant drop in efficiency. This is because higher amounts of nitrogen, which has a lower saturation temperature, shift the phase equilibrium, requiring lower temperatures for reliquefaction. On the other hand, a temperature limit of −55 °C was imposed to avoid the possibility of dry ice formation, restricting liquid formation.
Therefore, cases were also optimized and compared by applying a two-stage vapor–liquid separation process (Cases 4, 5, and 6) to address the issue of decreasing reliquefaction rates with increasing nitrogen content. In the reliquefaction processes with an additional vapor–liquid separator before expansion (Cases 4, 5, and 6), efficient reliquefaction was possible because lighter nitrogen was removed before the Joule–Thomson expansion compared to CO2. As a result, the processes with an added vapor–liquid separator (Cases 4, 5, and 6) could achieve an SPC reduction of 1.6–74.7% compared to Cases 1, 2, and 3.
Among these, Case 6 allowed for efficient heat exchange using mixed refrigerants, making it possible to expect an additional SPC reduction of 7.4–14.6% compared to Cases 4 and 5. Consequently, in the process of applying the vapor–liquid separator and mixed refrigerants (Case 6), an SPC reduction of 8.8–74.7% could be expected compared to the existing process design. Only Case 6 showed a reliquefaction rate of over 74.2% across all nitrogen content ranges from 5% to 15%.
When impurities, such as nitrogen, increased in the CO2 BOG, the reliquefaction rate decreased gradually, and the amount of non-condensable gas vented to the atmosphere increased. This leads to increased greenhouse gas emissions and causes economic losses due to cargo loss from the perspective of CO2 carriers. Therefore, future studies will be needed to design reliquefaction processes considering economic aspects during CO2 transportation.
Notes
Youngsub Lim serves as an editorial board member of the Journal of Ocean Engineering and Technology but played no role in deciding the publication of this article. No potential conflicts of interest relevant to this study are reported.
This research was supported by Korea Institute of Marine Science & Technology Promotion(KIMST) funded by the Ministry of Oceans and Fisheries(2520000243), and the Korea Research Institute of Ships and Ocean engineering, grant from Endowment Project of “Technology Development of Onboard Carbon Capture and Storage System and Pilot Test” funded by Ministry of Oceans and Fisheries (PES5110). The Research Institute of Marine Systems Engineering and Institute of Engineering Research at Seoul National University provided research facilities.