Onboard CO2 Capture Process Design using Rigorous Rate-based Model

Article information

J. Ocean Eng. Technol. 2022;36(3):168-180
Publication date (electronic) : 2022 June 30
doi : https://doi.org/10.26748/KSOE.2022.006
1Graduate student, Department of Naval Architecture and Ocean Engineering, Seoul National University, Korea
2Professor, Department of Naval Architecture and Ocean Engineering, Seoul National University, Korea
Corresponding author Yutaek Seo: +82-2-880-7329, yutaek.seo@snu.ac.kr
Received 2022 March 15; Revised 2022 June 5; Accepted 2022 June 8.

Abstract

The IMO has decided to proceed with the early introduction of EEDI Phase 3, a CO2 emission regulation to prevent global warming. Measures to reduce CO2 emissions for ships that can be applied immediately are required to achieve CO2 reduction. We set six different CO2 emission scenarios according to the type of ship and fuel, and designed a monoethanolamine-based CO2 capture process for ships using a rate-based model of Aspen Plus v10. The simulation model using Aspen Plus was validated using pilot plant operation data. A ship inevitably tilts during operation, and the performance of a tilted column decreases as its height increases. When configuring the conventional CO2 capture process, we considered that the required column heights were so high that performance degradation was unavoidable when the process was implemented on a ship. We applied a parallel column concept to lower the column height and to enable easy installation and operation on a ship. Simulations of the parallel column confirmed that the required column height was lowered to less than 3 TEU (7.8 m).

1. Introduction

Efforts to respond to climate change are spreading to all industries around the world. In 2015, the Paris Agreement was adopted, aiming toward worldwide efforts to keep the global average temperature rise below 2 °C above the pre-industrial level and further limit the future temperature rise to below 1.5 °C (Bodansky, 2016). To implement this, the International Maritime Organization (IMO) has established and implemented regulations to reduce greenhouse gas (GHG) emissions from ships. According to the third IMO GHG study, CO2 emitted from ships worldwide in 2012 accounted for 2.2% of the total CO2 emissions (IMO, 2014). This exceeds the CO2 emissions of Germany, Canada, and Korea (Olivier et al., 2017). According to the fourth IMO GHG study, CO2 emitted from ships worldwide in 2018 accounted for 2.89% of the total CO2 emissions, showing an increasing trend of the proportion of CO2 emissions from ships. To respond to this, the 2019 Marine Environment Protection Committee (MEPC) 74 determined the introduction time of the Energy Efficiency Design Index (EEDI) Phase 3. The EEDI is an operational efficiency indicator for ships and refers to the CO2 emission when 1 ton of a ship operates for 1 sea mile (1.852 km). The EEDI Phase 3 requires a reduction of 30% or more of CO2 emissions from ships from 2025 to 2030 compared with 2008, and Phase 4, which will be applied after 2030, requires a CO2 reduction of more than 40%. However, the MEPC 75 in 2020 proposed to reinforce emission regulations. Accordingly, part of the EEDI Phase 3, which was originally scheduled to be introduced in 2025, was moved forward to 2022.

Various methods are being devised to reduce GHGs emitted from ships to achieve the IMO’s CO2 emission reduction strategy. According to DNV-GL (2017), the main CO2 reduction methods that have been attempted so far are classified into the following four categories: liquefied natural gas (LNG) with the use of alternative fuels e.g., hydrogen, increased energy efficiency, the reduction of navigation speed, and carbon pricing. Among alternative fuels, LNG is a representative fuel, and extensive reviews have been reported on liquefied petroleum gas, biodiesel, bio methanol, liquefied biogas, hydrogen, and nuclear power. Methods to increase energy efficiency include the development of a new hull form, recycling waste heat, engine overhaul, the development of a hybrid engine, and main engine air lubrication. Combining all these methods can achieve fuel savings of 21% to 37% per ship (Kristensen, 2012). The energy efficiency improvement of ships will continue to rise gradually until 2050 considering improvement measures, such as hull form improvement, the optimization of ship speed and operation, propulsion system, and low/zero-carbon fuels. However, from an operational point of view, the energy efficiency improvement is expected to peak in 2035, and carbon reduction by alternative fuels will dominate afterward (DNV-GL, 2018). Excluding alternative fuels, it is estimated that 20%–30% of current CO2 emissions can be reduced by currently applicable technical and operational measures. In the future, carbon reduction by alternative fuels, such as hydrogen, ammonia, and biodiesel, should be promoted. However, developing and applying related technologies are currently challenging tasks, and the corresponding infrastructure is also insufficient. Therefore, as the effective date of the EEDI Phase 3 has been partially advanced to 2022, an onboard CO2 capture technology that can be applied immediately is required to achieve the target CO2 emission reduction.

Several researchers have studied onboard CO2 capture technology in various ways. Zhou and Wang (2014) proposed a method for capturing and fixing CO2 as calcium carbonate using calcium hydroxide solution and sodium hydroxide. This method was applied to a bulk carrier with an 18,660 kW engine, and the effectiveness and economic feasibility were evaluated. Luo and Wang (2017) simulated the monoethanolamine (MEA)-based post-combustion CO2 capture process and the CO2 storage liquefaction process for a cargo ship with a 17 MW engine. They observed that the carbon reduction rate could only reach 73% when conventional marine energy systems were integrated with the CO2 capture process owing to the limited heat and electricity supply to the CO2 capture process. They also observed that the cost of CO2 capture more than doubled when an additional gas turbine was installed to achieve a carbon reduction rate of 90%. Feenstra et al. (2019) simulated the CO2 capture process based on MEA and piperazine (PZ) using Aspen Plus for 1,280 kW and 3,000 kW class marine engines. Furthermore, they calculated the capital expenditure and operating expenditure required to capture CO2 from ship exhaust gas through the process and suggested the addition of a CO2 capture process and a CO2 storage tank to the existing cargo ship design. Lee et al. (2021) proposed a new EEDI estimation method considering the CO2 capture process and applied it to a 53,200 DWT class ship. They simulated the N-methyldiethanolamine- and PZ-based CO2 capture process and the liquefaction process of the captured CO2 using Aspen Plus and considered a design that placed a liquefied CO2 storage tank on a ship. The calculation results confirmed that the carbon capture ratio required in the CO2 capture process was higher than the actual EEDI reduction rate.

2. Process Model Framework

2.1 Rate-based Model

In this study, the CO2 capture process based on an MEA solution was simulated using Aspen Plus v10, a commercial process simulator. Moreover, simulation was performed using a rate-based model for better accuracy. The conventional equilibrium model commonly used for distillation column simulation assumes that gas and liquid phases reach complete equilibrium at each stage and adjusts the performance of the distillation column by introducing an efficiency correction factor in each phase. However, such perfect gas-phase and liquid-phase equilibrium states are rare in actual processes. In contrast, the rate-based model assumes that there are several layers of thin film at the gas–liquid interface according to the film theory as shown in Fig. 1. The Maxwell–Stefan is calculated for this thin film to actualize the realistic heat and mass transfer process at the gas–liquid interface (Al-Baghli, 2001). Through this process, the rate-based model can more closely simulate the chemical reaction in the actual tray column or packed column. When simulating processes with active chemical reactions, such as the CO2 capture process using amines, the rate-based model shows higher reproducibility than the equilibrium model (Zhang and Chen, 2013).

Fig. 1.

Schematic of the behaviors of the liquid and gas phases of the inner stage of the absorber through the film theory

2.2 Thermodynamic Model

Chemical reactions in the liquid phase must be considered to simulate the chemical equilibrium between gas and liquid accurately. The CO2 capture process using amines shows a nonideal behavior owing to its own chemical reaction and ions participating in the reaction. The electrolyte nonrandom two-liquid Redlich–Kwong state model was used to simulate the activity coefficient, Gibbs energy, enthalpy, and entropy of the liquid phase. Furthermore, the fugacity coefficient of the weather (Agbonghae et al., 2014) was simulated using the perturbed-chain statistical associating fluid theory model. The CO2 absorption reaction using MEA is expressed by the equilibrium Eqs. (1)(5) below.

(1) 2H2OH3O++OH
(2) 2H2O+CO2HCO3+H3O+
(3) H2O+HCO3CO32+H3O+
(4) MEAH++H2OMEA+H3O+
(5) MEACOO+H2OMEA+HCO3

2.3 Reaction Kinetic Model

In Aspen Plus, the reaction rate rj of a chemical reaction is expressed using a power law as follows:

(6) rj=kj0exp(jR[1T1298.15]i=1Naiαij)
where rj is the reaction rate of the chemical reaction, kj0 is the pre-exponential factor, j is the activation energy, R is the gas constant, T is the absolute temperature of the system, ai is the activity coefficient, and αij is the reaction order.

The chemical reaction of carbamate and bicarbonate formed by the CO2 absorption reaction using MEA are shown in Eqs. (7)(10).

(7) MEA+CO2+H2OMEACOO+H3O+
(8) MEACOO+H3O+MEA+CO2+H2O
(9) CO2+OHHCO3
(10) HCO3CO2+OH

Zhang et al. (2009) converted the molarity-based speed constant to activation-based speed constant through the experimental data of Hikita et al. (1979) and Pinsent et al. (1956). Therefore, the calculation method of the equilibrium constant is based on “mole gamma.” This result is summarized in Table 1.

Reaction rates and pre-exponential factors for the absorption of CO2 using MEA

3. Validation of Process Simulation

3.1 Target Plant for Simulation Validation

The CO2 capture process was simulated based on pilot plant operational data (Stec et al., 2015) to validate the process simulation prior to simulating the onboard CO2 capture process. The pilot plant to be verified implements a post-combustion carbon capture process based on a 30 wt% MEA solution. Table 2 lists the physical quantities and composition of the acid gas, Table 3 lists the physical quantities of the aqueous amine solution, and Table 4 lists the detailed operating characteristics of the absorber and stripper. The acid gas contains 13.5 mol% of CO2, and the flow rate can be changed in the range of 200–400 kg/h. The absorber is a packed column with a diameter of 0.33 m and a height of 5.1 m filled with Sulzer Mellapak 500Y and 750Y. The stripper is a packed column with a diameter of 0.28 m and a height of 4.3 m filled with Sulzer Mellapak 750Y. The pilot plant was operated by applying various process improvement methods. The process simulation was validated using the amine process operation data of the most popular standard method.

Properties of the acid gas for the pilot plant

Properties of lean solvent for pilot plant

Specifications of columns

3.2 Result of Simulation Validation

Various correction coefficients, correlation coefficients, and correlation methods of the rate-based model should be carefully selected and adjusted to construct a realistic process model that simulates real chemical reactions well. Table 5 summarizes the main tunable parameters and correlation method of the Aspen Plus rate-based model used to simulate pilot plant operation data. According to Zhang et al. (2009), when the Onda correlation method (Onda et al., 1968) is used, there is a possibility of underestimating the interfacial area. Therefore, the Bravo correlation method (Bravo et al., 1985) was used instead of the Onda correlation method which is generally used as the interfacial area method and mass transfer coefficient method. (Agbonghae et al., 2014). The Stichlmair method (Stichlmair et al., 1989) was used for liquid holdup, and the heat transfer coefficient was used by the Chilton and Colburn method (Chilton and Colburn, 1934). According to Zhang et al. (2009), the prediction accuracy of the countercurrent flow model is the highest, but this model requires a large number of calculations and sometimes shows unstable calculation results. Therefore, the Vplug flow model, which produces stable results, was used.

Design parameters of the columns

Table 6 lists the simulation results based on the configuration of the Aspen Plus model described above in comparison with the pilot data. Fig. 2 compares the temperature profile inside the actual pilot plant absorber with that of the simulation. The temperature decrease at the fourth stage of the pilot plant absorber could not be simulated. However, it can be seen that the maximum temperature bulge of 7th to 9th stages formed by the CO2 absorption reaction of the amine aqueous solution and the overall temperature trend were well simulated. From the above results, it can be confirmed that the simulation method predict the actual data from pilot-plant very well.

Operating conditions of columns

Fig. 2.

Absorber temperature profile of the pilot plant

4. Simulation Result of Onboard Carbon Capture Process

4.1 Selection of Onboard CO2 Emission Scenarios

The EEDI is an efficiency indicator determined as follows. First, the CO2 emissions from the main engine, the CO2 from the auxiliary

(11) EEDI=(j=1nfj)(i=1nMEPME(i)×CFME(i)×SFCME(i))+(PAE×CFAE×SFCAE)+PTI+EFFfi×fc×Capacity×fw×Vref
(12) EEDI(modified)=(j=1nfj)(i=1nMEPME(i)×(CFME(i)×fCO2(i))×SFCME(i))+(PAE×CFAE×SFCAE)+PTI+EFFf×Capacity×Vref
engines, and the additional CO2 from power take-in (PTI) are added together. Then, this value minus the CO2 reduction by energy abatement devices, such as waste heat recovery system, is divided by the size and speed of the ship (Eq. (11)). P denotes the output of the engine (kW), and ME and AE denote the main and auxiliary engines, respectively. SFC refers to the fuel consumption rate (g/kWh) of an engine, and CF refers to the tonnage of CO2 emitted from the ship when 1 ton of specific fuel is consumed. Capacity refers to deadweight tonnage, and Vref refers to the standard speed of the ship (knot). Finally, f is the efficiency index of the ship by ship type. Detailed figures and explanations related to each item in the can be found in IMO (2018). According to the current EEDI calculation method, even if CO2 is directly captured and removed from the ship exhaust gas, this figure is not reflected in the EEDI of the ship. Lee et al. (2021) proposed an improved EEDI calculation method as expressed in Eqs. (12) and (13), which reflects the amount of removed CO2 in EEDI when CO2 contained in exhaust gas emitted from a ship to the atmosphere is removed using the onboard carbon capture process. In this method, the carbon capture process can be reflected in the EEDI by adding Eq. (13) for the CO2 reduction rate to the term denoting the CO2 emission from the main engine. mcapture refers to the mass of CO2 removed, and mexhaust refers to the mass of total CO2 included in the exhaust gas. In this study, the EEDI of a ship was calculated using the EEDI calculation method proposed by Lee et al. (2021). All the EEDI calculations below have been performed according to the calculation method of Lee et al. (2021).
(13) fCO2=1mcapturemexhaust

The CO2 emitted from ships can be estimated using the following equation:

(14) E=Pj×LFj×EFj
where E is the hourly CO2 emissions (kg/h), Pj is the average maximum output of the engine by ship type (kW), LFj is the load factor of the engine by ship type, and EFj is the emission factor (g/kWh) of the engine according to the exhaust gas composition. Tables 7, 8, and 9 show the average maximum output (kW) of diesel engines by ship type, the average load factor of diesel engines by ship type, and the emission factor according to fuel, respectively. Table 10 summarizes the composition of exhaust gas discharged from MAN B&W’s ME-GI (main engine electronic control gas injection) engine.

Average power of diesel engines by ship type (kW) (U.S. Environmental Protection Agency, 2009)

Average load factor of diesel engines by ship type (U.S. Environmental protection Agency, 2009; ENTEC, 2007)

CO2 emission factor of four-stroke engine (g/kWh) (Kristensen, 2012)

Flue gas compositions of MAN B&W ME-GI engine (Kristensen, 2012)

The IMO presented the EEDI reference line as an exponential function as shown in Table 11 according to the ship type and tonnage, and specified the target reduction rate compared with the EEDI reference line according to the ship type and tonnage. Table 12 summarizes the EEDI Phase 3 target reduction rates by major ship type.

EEDI reference line (MEPC 215(63); IMO, 2012)

Required EEDI for phase 3 (MEPC 75/18; IMO, 2020)

According to the third IMG GHG study, three ship types, i.e., bulk carrier, tanker, and container ship, account for 60% of the total onboard carbon emissions. Therefore, a total of six onboard emission scenarios were selected for the case of using diesel and LNG as fuel for each of the three ship types: bulk carrier, tanker, and container ship. Scenarios 1, 2, and 3 were selected for bulk carriers, container ships, and tankers using diesel as fuel, respectively. In addition, scenarios 4, 5, and 6 were selected for bulk carriers, container ships, and tankers using LNG as fuel, respectively. Table 13 lists the EEDI values for each of the six scenarios based on the above data related to CO2 emission from ships. The tonnage of the ship used in the EEDI calculation was interpolated using the tonnage-engine output data of the Environmental Protection Agency (2009). For the Vref of the ship, the speed of the ship was set to 0.75 MCR (maximum continuous rated) proposed by Notteboom and Carriou (2009). Furthermore, the EEDI reference according to the tonnage and ship type for each scenario was calculated. Finally, the target EEDI required for the EEDI Phase 3 was calculated, and the required CO2 reduction rate to achieve this goal is listed in Table 13.

EEDI by ship emission scenario

Then, based on the above data on ship CO2 emissions, the CO2 emissions and total exhaust gas flow rate of each scenario were calculated and summarized in Table 14. As listed in Table 13, the CO2 reduction rate required in the scenario using diesel as fuel is 40%–50%. Hence, the target CO2 reduction rate in the diesel scenario was set to 50%. Similarly, the CO2 reduction rate was set to 30% because the CO2 reduction rate required in the scenario where LNG is used as fuel is 20%–30%. The target CO2 reduction rates for each scenario are summarized in Table 14.

CO2 emissions and exhaust gas flow rate by scenarios

4.2 Simulation of Single Packed Column Process

Fig. 3 shows a flowchart of the simplest and most traditional type of CO2 capture process using an aqueous amine solution. The process consists of two packed columns or tray columns and a heat exchanger. As the aqueous amine solution (lean amine) passes through the absorber, it absorbs CO2 from the acid gas. The aqueous amine solution that has absorbed CO2 (rich amine) is introduced into the stripper through the heat exchanger. The reboiler of the stripper performs high-temperature distillation to separate CO2, which is then discharged to the top of the stripper. The regenerated aqueous amine solution passes through the heat exchanger and returns to the absorber.

Fig. 3.

Conventional CO2 capture process based on aqueous amine solution

For six onboard CO2 emission scenarios, an onboard CO2 capture process that removes 70% to 80% of CO2 trapped in acid gas using a 30 wt% MEA solution was simulated based on the various coefficients and methods of the rate-based model used for simulation validation in the previous section. The removal rate of the CO2 capture process using the MEA solution is generally set to 90%. However, when the removal rate is decreased, the amount of fluid flowing into the absorber increases. Hence, the diameter of the absorber increases, whereas the height of the absorber decreases. A removal rate lower than 90% was set to design an absorber with a lower height, and the reason will be described below. The absorber is packed with the filler material of Mellapak 250Y, and both diesel and LNG fuel usage scenarios are packed to 7.8 m. The exhaust gas cooled to 40 °C flows into the bottom of the absorber, and the MEA solution cooled to 45 °C flows into the top of the absorber. The stripper is packed with Mellapak 250Y to 4.1 m in the diesel scenario and 4.0 m in the LNG scenario. The MEA solution is heated to 109.5 °C through the reboiler at the bottom of the stripper. For both absorber and stripper, the maximum flooding rate was set to 75%. Detailed operation information of each column is summarized in Table 15, and the process operation flow rates and target CO2 removals are summarized in Table 16. The composition of flue gas CO2 of diesel fuel is 4.74 mol%, and the composition of flue gas CO2 of LNG fuel is 3.52 mol% (Table 10). When diesel is used as fuel, the flow rate of flue gas and CO2 emission are higher than those when LNG is used as fuel.

Design specifications for single column

Operating conditions for a single column

The lean loading (mol CO2/mol MEA) of the MEA solution flowing into the absorber should be considered carefully because it is one of the major factors influencing the overall process, including the flow rate of the MEA solution, the size of the absorber and stripper, and the amount of heat consumed in the reboiler of the stripper. Fig. 4 shows the effect of lean loading on the flow rate of the MEA solution required and the amount of heat consumed in the reboiler. The lowest energy consumption is 3.80 GJ/t CO2 for the lean loading of 0.19 mol CO2/mol MEA. Therefore, the lean loading of the MEA solution simulating the onboard CO2 capture process was fixed at 0.19 mol CO2/mol MEA, which requires the lowest amount of energy. At a lean loading below 0.19 mol, the required flow rate of the MEA solution decreases, whereas the CO2 absorption performance and speed increase; however, the amount of heat consumed to regenerate the MEA solution increases exponentially. This is because more energy is required to achieve a lower level of CO2 loading when rich MEA that has absorbed CO2 is regenerated in the stripper. If the lean loading increases, the required flow rate of the MEA solution increases. This in turn increases sensible heat, the overall size of the equipment, and the risk of column flooding.

Fig. 4.

Effect of lean loading on the MEA solution flow rate (red line) and reboiler duty (blue line)

Table 17 lists the size results of simulating the CO2 capture process according to the aforementioned six onboard CO2 emission scenarios. The size simulation of the device was performed using the Aspen process economic analyzer installed as an add-on of Aspen Plus. The flow rate of exhaust gas varied considerably according to each scenario, and the column diameter changed significantly as a result. In both the diesel and LNG scenarios, the absorber height was the same at approximately 12.2 m. The height of the stripper was approximately 8.5 m for the diesel scenario and approximately 8.4 m for the LNG scenario. This is because the lean solvent loading was the same in both scenarios, but the partial pressure of CO2 of the aqueous amine solution flowing into the stripper was different. Furthermore, the overall column height was significantly increased by additional devices installed in the columns compared with the height of the filler in each case.

Simulation result of single column

The onboard CO2 capture process assumed in this study aims to demonstrate the possibility of implementing a simple CO2 capture process without extensive remodeling for ships beyond EEDI Phase 3 among existing ships. However, the calculated column height exceeded 8 m in all scenarios, and the height of the absorber in particular reached 12.2 m. Assuming that all the above scenarios are for small- and medium-sized ships, there is a possibility that such a high column height may exceed the height of the engine room or the height of the stack. This indicates that, to implement the CO2 capture process, extensive repair and maintenance of the ship are required, such as making a hole in the ship’s deck or expanding the stack. Moreover, the hull is inevitably inclined owing to six-degree-of-freedom motions, such as roll, pitch, and yaw, which is unavoidable for the ships floating on the water during operation. Consequently, as the absorber installed on the ship is tilted, it is likely that the CO2 absorption performance will decrease (Di et al., 2018). According to Son et al. (2017), maldistribution of liquid inside a tilted column increases significantly as the packed column height increases. This interferes with smooth process operation and decreases column performance. Therefore, in this study, a parallel column CO2 capture process with two absorbers and two strippers was newly simulated with the objective of implementing a CO2 capture process in which the height does not exceed 3 TEU (TEU indicates the size of a standard container, and the height of 3 TEU is 7.8 m).

4.3 Simulation of Parallel Packed Column Process

Absorber intercooling is an improved absorber process operation method widely used to increase the efficiency of the on-shore CO2 capture process using amines. This method forms a temperature imbalance inside the absorber because the CO2 absorption reaction by amines is an exothermic reaction. In general, the highest temperature in the absorber is generated in the middle part of the absorber. The absorber intercooling process decreases the overall temperature of the fluid inside the absorber by installing a cooler in this highest-temperature part to provide a cooling circulating flow inside it. CO2 absorption reaction using amines occurs more actively in an environment with a low temperature owing to the solubility of gas. Hence, the installation of an intercooler increases the CO2 absorption capacity of the aqueous amine solution. Therefore, the intercooling process can have positive effects, such as reductions in the flow rate of the aqueous amine solution required for CO2 capture, the overall size of the equipment, the energy consumption of the reboiler, and process operation cost.

A parallel absorber using an intercooling device was devised to lower the column height required for the onboard CO2 capture process based on the basic concept of the absorber intercooling process. As shown in Fig. 5(a), the regenerated aqueous amine solution flows into the upper part of the first absorber, and the gas from which CO2 is removed is discharged through the upper part of the first absorber. The aqueous amine solution that has slightly absorbed CO2 comes out from the bottom of the first absorber, passes through the cooler, and flows into the top of the second absorber at a low temperature. The exhaust gas flows into the bottom of the second absorber, and the aqueous amine solution that has completely absorbed CO2 is discharged from the bottom of the second absorber.

Fig. 5.

Process flow diagram of CO2 capture process using parallel columns

Fig. 6(a) shows the temperature gradient inside the absorber in a single absorber process. A relatively low temperature is maintained at the top and bottom of the absorber because the inflow of the low-temperature MEA solution and exhaust gas continues. However, the temperature of the fluid inside the absorber rises owing to the CO2 absorption reaction of MEA, which is an exothermic reaction, and the highest temperature of 72.5 °C can be observed near stage 5. Fig. 6(b) shows the temperature gradient inside the absorber in the parallel absorber process. Stages 1 to 10 represent the first absorber, and stages 11 to 20 represent the second absorber. A low temperature is maintained in stages 1 and 20, where the MEA solution and exhaust gas are introduced, as well as in stages 10 and 11, which pass through the cooler. Hence, the maximum temperature of the entire absorber remains at approximately 61 °C. The absorption performance of the MEA solution also increased, which was able to reduce the total MEA solution flow rate by 12%.

Fig. 6.

Temperature profiles of absorbers

In Figs. 6 (a) and (b), the temperature of the incoming MEA solution is 40 °C in both cases. However, the temperatures of the single absorber in stage 1 are 52.8 °C (liquid) and 60.8 °C (gas), and the temperatures of the parallel absorber in stage 1 are 45.2 °C (liquid) and 49.6 °C (gas). This is because stage 5, the point at which the reaction occurs most actively in the absorber, is close to stage 1, the inlet of the MEA solution; therefore, the temperature of stage 1 changes relatively significantly according to the maximum internal temperature of the absorber. As the maximum temperatures inside the single absorber and the parallel absorber differ by approximately 10 °C, the temperature of stage 1, which is close to temperature bulge, also shows a difference of approximately 5–10 °C. The basic concept of stripper interheating is the same as that of absorber intercooling. A heater is installed in the middle of the stripper to raise the overall temperature of the stripper. The amine solution that absorbs CO2 is heated in the stripper to separate the absorbed CO2. The regeneration efficiency increases if a higher temperature is maintained up to the top of the stripper. This can be expected to reduce the energy consumption in the reboiler and the size of the stripper. Based on this concept, the rich MEA solution is introduced at the top of the first stripper and the separated CO2 is discharged as shown in Fig. 5(b). The MEA solution heated by the reboiler at the bottom of the first stripper flows into the top of the second stripper. At the bottom of the second stripper, the MEA solution after CO2 separation flows to the absorber through the heat exchanger.

Fig. 7 shows the energy consumption required to achieve the same lean CO2 loading according to the ratio of energy consumed by the reboiler of each stripper. The x-axis represents the ratio of energy consumption of the reboiler of the second stripper to the total energy consumption. The y-axis represents the total energy consumption. The amounts of energy consumed by the first and second reboilers, respectively, are shown in different colors. A relatively high energy consumption can be observed at 0%–10% when the reboiler of the second stripper is used less or is not used. This indicates that inefficient energy input is required to heat the bottom of the second stripper sufficiently only with the reboiler of the first stripper. The overall energy consumption decreases as the proportion of energy consumption of the reboiler of the second stripper gradually increases, reaching a minimum of 3.76 GJ/t CO2 at 60%. Beyond the proportion of 60%, the total energy consumption increases again. This appears to be because a relatively large amount of energy is input to heat the top of the first stripper sufficiently with the reboiler of the second stripper. Therefore, in the parallel absorber process, the ratio of energy consumption of the reboiler of the first stripper and the reboiler of the second stripper was simulated as 4:6. Tables 18 and 19 list the detailed operational information of the parallel packed column process.

Fig. 7.

Total duty change according to reboiler duty ratio

Design specifications for parallel absorbers

Design specifications for parallel strippers

Table 20 lists the process simulation size results for the six scenarios of the CO2 capture process for parallel packed columns. The diameters of both the absorber and stripper changed according to the flow rate of each scenario. When compared with the single absorber process, the flow rate of the MEA solution decreased, and the overall diameter of the column decreased as a result. The absorber heights in both scenarios were approximately 7.6 m, which were much lower than the height of the single packed column, i.e., 12.2 m. Thus, the column height was simulated to be lower than the intended height of 3 TEU (7.8 m). The height of the stripper was simulated as approximately 6.5 m in the diesel scenario and approximately 6.4 m in the LNG scenario, which were lower than 8.5 m and 8.4 m for the single packed column, respectively, and also lower than the target height of 3 TEU.

Simulation results of parallel columns

4. Conclusions

A technology to reduce CO2 emissions that can be immediately applied is required owing to the expedited implementation of the EEDI Phase 3. Thus, in this study, an MEA-based onboard CO2 capture process was simulated and the required column size was examined. Before the simulation, the average CO2 emission according to the ship type was calculated, and the amount of CO2 reduction required for the implementation of the EEDI Phase 3 was calculated. Six scenarios were defined according to ship type and fuel, and the MEA-based onboard CO2 capture process required for these scenarios was simulated. The process simulation was conducted using the rate-based model of Aspen Plus v10, a commercial process simulator. The operational data of the pilot plant were used to verify the accuracy and reproducibility of the process simulation. Various correction factors and correlation methods of the rate-based model were carefully adjusted and selected based on the operation data of the pilot plant. The results of process simulation were validated by comparing them with the operation data. The MEA-based CO2 capture process required for the six onboard CO2 emission scenarios selected above was simulated through the adjusted rate-based model. The onboard CO2 capture process was designed to include the most widely used basic types of absorber, stripper, and heat exchanger. Furthermore, the CO2 reduction rate of 70%–80% was targeted based on lean loading with the minimum energy consumption in the reboiler. The result of process simulation confirmed that an absorber and a stripper with a diameter of approximately 0.7 m–3.6 m and a height of 8.4 m–12.2 m were required depending on the scenario.

A lower column height in the onboard CO2 capture process is more advantageous owing to the characteristics of ships. The column height of up to 12.2 m in the onboard CO2 capture process simulated above was too high for implementation on a ship. Therefore, to lower the column height, an onboard CO2 capture process using parallel packed columns was newly designed and simulated. Both the absorber and stripper were designed to have two parallel columns, and each column was connected by applying the intercooling and interheating concepts. Consequently, the onboard CO2 capture process could be configured with columns having a diameter of approximately 0.7 m–3.5 m and a height of 6.4 m–7.6 m. Therefore, immediate onboard application would be possible because the column height is less than 3 TEU or 7.8 m. However, the CO2 capture process with parallel packed columns increases the space required for installation, as the required number of columns increases compared with that in the conventional process. This is a disadvantage, as ships have several spatial constraints. Therefore, the size of the engine room, the location of the stack, and the arrangement of the deck structure should be carefully considered to apply the parallel column process to an actual ship.

The CO2 capture process using amines requires a considerable amount of heat to regenerate the aqueous amine solution owing to the nature of the process. It is assumed that the thermal energy required for the onboard amine CO2 capture process is supplied through the waste heat of the engine and exhaust gas. However, an auxiliary engine is required to generate additional power if sufficient waste heat is not supplied to regenerate all the aqueous amine solution. This additional power generation results in higher CO2 emissions than before. Therefore, the target CO2 removal to achieve the EEDI Phase 3 will be higher than that initially calculated. This aspect requires a close review because it leads to additional increases in load and size in the capture process. Therefore, further studies are required to optimize the energy flow of the process and to design the process considering the additional energy required by the implementation of the onboard CO2 capture process.

Notes

No potential conflict of interest relevant to this article was reported.

This research was supported by a grant (RS-2022-00143644) funded by Ministry of Land, Infrastructure and Transport of Korean government.

References

Agbonghae EO, Hughes KJ, Ingham DB, Ma L, Pourkashanian M. 2014;Optimal Process Design of Commercial-Scale Amine-Based CO2 Capture Plants. Industrial & Engineering Chemistry Research 53(38):14815–14829. https://doi.org/10.1021/ie5023767 .
Al-Baghli NA, Pruess SA, Yesavage VF, Selim MS. 2001;A Rate-based Model for the Design of Gas Absorbers for the Removal of CO2 and H2S Using Aqueous Solutions of MEA and DEA. Fluid Phase Equilibria 185(1–2):31–43. https://doi.org/10.1016/S0378-3812(01)00454-X .
Bodansky D. 2016;The Legal Character of the Paris Agreement. Review of European, Comparative & International Environmental Law 25(2):142–150. https://doi.org/10.1111/reel.12154 .
Bravo JL, Rocha JA, Fair JR. 1985;Mass Transfer in Gauze Packings. Hydrocarbon Processing (International ed) 64(1):91–95.
Chilton TH, Colburn AP. 1934;Mass Transfer (Absorption) Coefficients Prediction from Data on Heat Transfer and Fluid Friction. Industrial & Engineering Chemistry 26(11):1183–1187. https://doi.org/10.1021/ie50299a012 .
DNV GL. 2017. Low Carbon Shipping Towards 2050.
DNV GL. 2018. Emission Pathway 2015–2030.
Di XN, Wang WH, Chen SJ, Huang Y. 2018;Effect of Tilt on Mass Transfer and Hydrodynamic Performance in a Packing Column. Chemical Engineering and Processing - Process Intensification 123:89–99. https://doi.org/10.1016/j.cep.2017.10.028 .
ENTEC. 2007. Ship Emissions Inventory – Mediterranean Sea, Final Report Entec Limited. London:
Environmental Protection Agency (EPA). 2009. Proposal to Designate an Emission Control Area for Nitrogen Oxides, Sulfur Oxides and Particulate Matter USA: Environmental Protection Agency.
Feenstra M, Monteiro J, Akker JT, Abu-Zahra MRM, Gilling E, Goetheer E. 2019;Ship-based Carbon Capture Onboard of Diesel or LNG-fuelled Ships. International Journal of Greenhouse Gas Control 85:1–10. https://doi.org/10.1016/j.ijggc.2019.03.008 .
Hikita H, Asai S, Katsu Y, Ikuno S. 1979;Absorption of Carbon Dioxide into Aqueous Monoethanolamine Solutions. AIChE Journal 25(5):793–800. https://doi.org/10.1002/aic.690250507 .
IMO. 2012. Guidelines for Calculation of Reference Lines for Use with the Energy Efficiency Design Index (EEDI). Resolution MEPC 215(63)International Maritime Organization. London:
IMO. 2014. Greenhouse Gas Study 2014, Executive Summary and Final Report International Maritime Organization. London:
IMO. 2018. Guidelines on the Method of Calculation of the Attained Energy Efficiency Index (EEDI) for New Ships. Resolution MEPC 308(73)International Maritime Organization. London:
IMO. 2020. Fourth IMO GHG Study 2020 – Final Report. MEPC 75/18 International Maritime Organization. London:
Kristensen HO. 2012;Energy Demand and Exhaust Gas Emissions of Marine Engines. Clean Shipping Currents 1(6):18–26.
Lee S, Yoo S, Park H, Ahn J, Chang D. 2021;Novel Methodology for EEDI Calculation Considering Onboard Carbon Capture and Storage System. International Journal of Greenhouse Gas Control 105:103241. https://doi.org/10.1016/j.ijggc.2020.103241 .
Luo X, Wang M. 2017;Study of Solvent-based Carbon Capture for Cargo Ships through Process Modelling and Simulation. Applied Energy 195:402–413. https://doi.org/10.1016/j.apenergy.2017.03.027 .
Notteboom T, Carriou P. 2009. Fuel Surcharge Practices of Container Shipping Lines: Is It About Cost Recovery or Revenue Mking? Proceedings of the International Association of Maritime Economists (IAME) Conference Copenhagen, Denmark:
Olivier JGJ, Schure KM, Peters JAHW. 2017;Trends in Gloval CO2 and Total Greenhous Gas Emissions. PBL Netherlands Environmental Assessment Agency
Onda K, Takeuchi H, Okumoto Y. 1968;Mass Transfer Coefficients between Gas and Liquid Phases in Packed Columns. Journal of Chemical Engineering of Japan 1(1):56–62. https://doi.org/10.1252/jcej.1.56 .
Pinsent BRW, Pearson L, Roughton FJW. 1956;The Kinetics of Combination of Carbon Dioxide with Hydroxide Ions. Transactions of the Faraday Society 52:1512–1520. https://doi.org/10.1039/TF9565201512 .
Son Y, Kim G, Lee S, Kim H, Min K, Lee KS. 2017;Experimental Investigation of Liquid Distribution in a Packed Column with Structured Packing under Permanent Tilt and Roll Motions Using Electrical Resistance Tomography. Chemical Engineering Science 166:168–180. https://doi.org/10.1016/j.ces.2017.03.044 .
Stec M, Tatarczuk A, Więcław-Solny L, Krótki A, Ściążko M, Tokarski S. 2015;Pilot Plant Results for Advanced CO2 Capture Process Using Amine Scrubbing at the Jaworzno II Power Plant in Poland. Fuel 151:50–56. https://doi.org/10.1016/j.fuel.2015.01.014 .
Stichlmair J, Bravo JL, Fair JR. 1989;General Model for Prediction of Pressure Drop and Capacity of Countercurrent Gas/Liquid Packed Columns. Gas Separation & Purification 3(1):19–28. https://doi.org/10.1016/0950-4214(89)80016-7 .
Zhang Y, Chen CC. 2013;Modeling CO2 Absorption and Desorption by Aqueous Monoethanolamine Solution with Aspen Rate-based Model. Energy Procedia 37:1584–1596. https://doi.org/10.1016/j.egypro.2013.06.034 .
Zhang Y, Chen H, Chen CC, Plaza JM, Dugas R, Rochelle GT. 2009;Rate-Based Process Modeling Study of CO2 Capture with Aqueous Monoethanolamine Solution. Industrial & Engineering Chemistry Research 48(20):9233–9246. https://doi.org/10.1021/ie900068k .
Zhou P, Wang H. 2014;Carbon Capture and Storage—Solidification and Storage of Carbon Dioxide Captured on Ships. Ocean Engineering 91:172–180. https://doi.org/10.1016/j.oceaneng.2014.09.006 .

Article information Continued

Fig. 1.

Schematic of the behaviors of the liquid and gas phases of the inner stage of the absorber through the film theory

Fig. 2.

Absorber temperature profile of the pilot plant

Fig. 3.

Conventional CO2 capture process based on aqueous amine solution

Fig. 4.

Effect of lean loading on the MEA solution flow rate (red line) and reboiler duty (blue line)

Fig. 5.

Process flow diagram of CO2 capture process using parallel columns

Fig. 6.

Temperature profiles of absorbers

Fig. 7.

Total duty change according to reboiler duty ratio

Table 1.

Reaction rates and pre-exponential factors for the absorption of CO2 using MEA

Related species Reaction Direction kj0 (kmol/m3s) j (kJ/mol)
MEACOO Forward 3.02×1014 41.20
Reverse (absorber, 303–353K) 5.52×1023 69.05
Reverse (stripper, 363–393K) 6.56×1027 95.24
HCO3 Forward 1.33×1017 55.38
Reverse 6.63×1016 107.24

Table 2.

Properties of the acid gas for the pilot plant

Acid gas properties

Pressure (bar) 1.05
(kPa) 206.3
Temperature (°C) 45
Flow rate (kg/h) 289

Compositions (mol fraction)

CO2 0.135
H2O 0.055
N2 0.7
O2 0.11

Table 3.

Properties of lean solvent for pilot plant

Properties of lean MEA solvent
Pressure (bar) 2.0
(kPa) 301.3
Temperature (°C) 40
Flow rate (kg/h) 1358.3
MEA concentration (wt%) 30
CO2 loading (mol CO2/mol MEA) 0.36

Table 4.

Specifications of columns

Absorber Stripper
Diameter (mm) 330 280
Height (mm) 8400 4300
Segment number 20 20
Packing type Sulzer Mellapak Sulzer Mellapak
500Y, 750Y 750Y
Reboiler temperature (°C) - 108

Table 5.

Design parameters of the columns

Absorber Stripper
Global tuning factors Reaction condition factor 0.9 0.9
Film discretization ratio 5 5
Flow model Vplug Vplug
Interfacial area 1.4 1.1
Film resistance Discretize film Discretize film
Liquid phase Number of discretization points 5 5
Liquid holdup Stichlmair89 Stichlmair89
Vapor phase Film resistance Consider film Consider film
Mass transfer coefficient Brf-85 Brf-85
Correlation methods Heat transfer coefficient Chilton and Colburn Chilton and Colburn
Interfacial area Brf-85 Brf-85

Table 6.

Operating conditions of columns

Pilot plant Aspen Plus model
Absorber Diameter (mm) 330 330
Height (mm) 8400 8400
Stripper Diameter (mm) 280 280
Height (mm) 4300 4300
Lean solvent loading mol CO2/mol MEA 0.36 0.36
Rich solvent loading mol CO2/mol MEA 0.50 0.503
CO2 removal rate % 84 86.1
Stripper reboiler duty GJ/t CO2 3.98 4.05

Stripper reboiler temperature °C 108 108

Table 7.

Average power of diesel engines by ship type (kW) (U.S. Environmental Protection Agency, 2009)

Bulk carrier Container ship Passenger ship General cargo RORO Tanker Reefer
8,000 30,900 39,600 9,300 11,000 9,400 9,600

Table 9.

CO2 emission factor of four-stroke engine (g/kWh) (Kristensen, 2012)

Fuels MDO (Marine Diesel Oil) LNG (Liquified Natural Gas)
Emission factor (g/kwh) 609 426

Table 10.

Flue gas compositions of MAN B&W ME-GI engine (Kristensen, 2012)

Composition Unit Diesel mode (MDO) Gas mode (LNG)
N2 mol% 80.21 81.83
O2 mol% 15.05 14.65
CO2 mol% 4.74 3.52
CO ppm 51.0 55.0
HC ppmC 46.5 143.5
NOx ppm 1002 1044

Table 11.

EEDI reference line (MEPC 215(63); IMO, 2012)

Reference line = a×bc

Ship type a b c
Bulk carrier 961.79 Deadweight tonnage (DWT) 0.477
Gas carrier 1120.00 0.456
Tanker 1218.80 0.488
Container ship 174.22 0.201
General cargo 107.48 0.216
Reefer 227.01 0.244
RORO 1405.15 0.498

Table 12.

Required EEDI for phase 3 (MEPC 75/18; IMO, 2020)

Ship type Bulk carrier Tanker Container ship
DWT reduction rate (%) 20,000– 10,000–20,000 20,000– 4,000–20,000 200,000– 120,000–200,000 80,000–120,000 40,000–80,000 15,000–40,000 10,000–15,000
30 0–30 30 0–30 50 45 40 35 30 15–30

Table 13.

EEDI by ship emission scenario

Case No. Fuel Ship type Capacity (DWT) Vref EEDI EEDI reference EEDI phase 3 Required CO2 reduction rate
1 Diesel (MDO) Bulk carrier 45688 14.5 6.250 5.759 4.031 (−30%) 43.9%
2 Container ship 47213 24 21.59 20.03 13.02 (−35%) 45.2%
3 Tanker 63750 12 6.323 5.512 3.859 (−30%) 40.2%

4 Gas (LNG) Bulk carrier 45688 14.5 5.029 5.759 4.031 (−30%) 28.5%
5 Container ship 47213 24 17.37 20.03 13.02 (−35%) 22.5%
6 Tanker 63750 12 5.088 5.512 3.859 (−30%) 27.2%

Table 14.

CO2 emissions and exhaust gas flow rate by scenarios

Case No. CO2 emissions (kg/h) Exhaust gas flow rate (kg/h) Target CO2 reduction rate Target CO2 reduction (kg/h)
1 3564.0 51447.4 50% 1782.0
2 15054.5 211963.3 7527.3
3 4293.4 60450.7 2146.7

4 2940.0 55341.2 30% 882.0
5 12112.8 228005.6 3633.8
6 3454.5 65025.9 1036.4

Table 15.

Design specifications for single column

Absorber Stripper
Packing type Mellapak 250Y Mellapak 250Y
Pressure 1.1 barg (210.3 kPa) 0.2 barg (121.3 kPa) (min.)
0.8 barg (181.3 kPa) (max.)
Temperature 42.7 °C (min.) 85 °C (condenser)
73.3 °C (max.) 109.5 °C (reboiler)
Lean solvent loading 0.19 -
Rich solvent loading 0.324 (Diesel) -
0.310 (LNG) -
CO2 removal rate 78.6% (Diesel) -
70.1% (LNG) -
Packing height 7.8 m 4.1 m (Diesel)
4.0 m (LNG)
L/G ratio 2 (Diesel) -
1.4 (LNG) -
Reboiler duty - 3.80 GJ/t CO2

Table 16.

Operating conditions for a single column

Case No. Target CO2 reduction (kg/h) CO2 removal rate in absorber Process inlet CO2 flow rate (kg/h) Process inlet exhaust gas flow rate (kg/h)
1 1782.0 2267.2 32727.8
2 7527.3 78.6% 9576.7 138244.8
3 2146.7 2731.2 39425.8

4 882.0 1258.2 23683.9
5 3633.8 70.1% 5183.7 97576.4
6 1036.4 1478.5 27829.9

Table 17.

Simulation result of single column

Case No. Absorber diameter (m) Absorber height (m) Stripper diameter (m) Stripper height (m)
1 1.7760 1.4340
2 3.6342 12.192 2.9473 8.5344
3 1.9408 1.5740

4 0.9604 0.7482
5 1.9880 12.192 1.5188 8.3820
6 1.0562 0.8190

Table 18.

Design specifications for parallel absorbers

Absorber 1 Absorber 2
Packing type Mellapak 250Y Mellapak 250Y
Pressure 1.05 barg (206.3 kPa) 0.97 barg (198.3 kPa)
Temperature 45.2 °C (min) 48.0 °C (min)
61.0 °C (max) 53.1 °C (max)
Lean solvent loading 0.19
Rich solvent loading 0.346 (Diesel)
0.331 (LNG)
CO2 removal rate 79.1% (Diesel)
70.6% (LNG)
Packing height 4.0 m
L/G ratio 1.7 (Diesel)
1.3 (LNG)

Table 19.

Design specifications for parallel strippers

Stripper 1 Stripper 2
Packing type Mellapak 250Y Mellapak 250Y
Pressure 0.2 barg (121.3 kPa) (min) 0.7 barg (171.3 kPa) (min)
0.68 barg (169.3 kPa) (max) 0.9 barg (191.3. kPa) (max)
Temperature 99.7 °C (min) 118.2 °C (min)
118.6 °C (max) 121 °C (max)
Packing height 2.3 m (Diesel)
2.2 m (LNG)
Reboiler duty 3.76 GJ/t CO2
3.55 GJ/t CO2

Table 20.

Simulation results of parallel columns

Case No. Absorber diameter (m) Absorber height (m) Stripper diameter (m) Stripper height (m)
1 1.7052 1.3452
2 3.4878 7.5820 2.7648 6.5360
3 1.8626 1.4764

4 0.9394 0.7019
5 1.9066 7.5820 1.4247 6.4196
6 1.0182 0.7608