A Combined Experimental and Computational Study on the Effect of the Reactor Configuration and Operational Procedures on Formation, Growth and Dissociation of Carbon Dioxide Hydrate

1. Introduction

Clathrate hydrates (also called gas hydrates) are nonstoichiometric crystalline ice-like substances, consisting of a hydrogen-bonded lattice formed by water molecules (host) and entrapped gas molecules (guest gas), which are not bonded with the lattice. One of the advantages of gas hydrates is the drastic reduction of the volume of the enclathrated gas compared with that in the gaseous phase. Under standard temperature and pressure, the volume of gas stored in 1 m³ of CO₂ hydrate is about 120-160 m³ [1]. The large uptake of CO₂ in clathrate hydrates renders them a promising option for either storing, transportation or even disposing and sequestration.

CO₂ hydrates have been investigated with regard to their formation conditions and potential applications in all large-scale fields of the Carbon Capture and Storage (CCS) system, including (i) the flow assurance and safe CO₂ pipeline transportations [2,3,4,5], (ii) CO₂ storage and separations [6,7,8], including Post-Combustion Capture [9,10], and Pre-Combustion Capture [11,12,13] applications, and (iii) CO₂ sequestration [14,15,16], in deep sea sediments and other subsurface storage sites to mitigate the global warming, as well as in natural gas hydrate reservoirs to facilitate CH₄ extraction and recovery, while at the same time to prevent the release of a greenhouse gas to the atmosphere [17,18,19]. Another application of gas hydrate technologies related to CO₂ capture and separation pertains to CO₂-rich natural gas upgrade [20], and biogas upgrade via CO₂-CH₄ separation that leads to biomethane production [21,22].

Gas hydrate crystals are formed when a hydrate former and water are mixed and the temperature and pressure conditions of the system reach the (H-L_w-G) equilibrium curve where the gas, liquid and hydrate phase co-exist. Conversely, gas hydrates decompose when the (P, T) conditions diverge from the hydrate-water rich liquid-gas (H-L_w-G) phase equilibrium [23]. Hydrate dissociation is an endothermic process in which the hydrogen bonds between water molecules of the hydrate lattice and the van der Waals interaction between the guest and water molecules of the lattice break to release water and gas. Several experimental methodologies have been used to study the phase behavior and equilibrium of the CO₂-H₂O system including the formation and dissociation conditions of CO₂ hydrates [23,24,25,26,27]. These conditions are determined either by visual observation using high-pressure optical cells or by monitoring the total pressure of the system as a function of temperature that allows for constructing the (P-T) phase diagrams.

Three static methods have been applied to study the gas-liquid-hydrate phase equilibrium and gas hydrate formation kinetics, in particular the Isothermal Pressure-Search, the Isobaric Temperature-Search, and the Isochoric Pressure-Search methods. The Isothermal Pressure-Search method [27,28,29] allows to determine the (P-T) phase diagram (phase envelope) of the system, and is based on conducting a cooling/heating cycle at constant volume. In this method the cell volume is kept constant and the temperature undergoes a cyclic alteration procedure consisted of a cooling step, an isothermal step, and a slow incremental heating step whereas the pressure is monitored constantly. Hydrate formation is initiated under sufficient subcooling, ΔT_sub, which is defined as the difference between the experimental temperature where gas hydrate nucleation takes place, T_exp, and the equilibrium temperature, T_eq, which corresponds to the theoretical (H-L_w-G) equilibrium conditions (T_eq, P_eq), i.e. ΔT_sub, = T_eq – T_exp, where T_exp < T_eq [30,31,32,33,34]. Hydrate crystals growth is detected from a sharp pressure drop at a constant temperature. After completion of hydrates formation, the hydrate crystals decompose through step-wise heating. At each isothermal step of the heating process, enough time is provided to the system to establish (H-L_w-G) equilibrium and thus, to derive a (P_eq, T_eq) curve that passes through consecutive equilibrium points. In the corresponding (P-T) phase diagram, the point where the heating curve joints the cooling curve and the slope of the heating curve changes abruptly, is considered as the point of complete hydrate dissociation [28]. The dissociation point is often confirmed by visual observation.

An interesting issue related to the experimental techniques that have been applied to study the (P, T) formation conditions and kinetics of growth and dissociation of CO₂ hydrates, pertains to the effects of the design and technical characteristics of the gas hydrate autoclave and experimental configuration [35,36,37]. The design, geometrical and technical features of the utilized experimental setup as well as the operational modes can affect the characteristics of mass and heat transfer. It has been reported that application of batch stirring can improve the mass transfer and heat transfer performances of the methane hydrate formation process [36]. Usage of a stirrer can enhance the mass transfer driving force, which – unless excessive agitation is applied - can promote gas hydrate nucleation and growth, thereby increasing hydrate growth rate. Englezos et al. [38] developed a kinetic model of hydrate formation by assuming that the hydrate formation process is similar to the crystallization process composed of two successive steps with different mass transfer resistances. In the first step, the dissolved guest molecules diffuse from the bulk of the liquid into the liquid water-hydrate interface through the laminar diffusion layer around a particle. The second step consists of adsorption of the guest molecules at the liquid water-hydrate interface, followed by their incorporation and stabilization of the structured framework of the water /guest gas cages. When the agitation is conducted at a high stirring rate, the mass transfer resistance and heat-transfer resistance around the crystal seed, become negligible [38]. Experimental formation data obtained in a semi-batch stirred tank reactor have implied that at a stirring rate of 400 rpm, the mass transfer resistance around the crystal seeds is eliminated, and is no longer affected by a further increase in the rpm [38]. A stirring rate of 400 rpm was chosen as the most suitable one to avoid extensive rippling on the gas-liquid interface. This stirring rate was also able for continuous removal of the hydrate seeds from the gas-liquid interface to the bulk liquid.

On the other hand, during growth of gas hydrate crystals, an amount of latent heat of formation is released, which has to be effectively transferred to the cooling bath in order to control the temperature. The heat transfer rate between the hydrate former-liquid system and the cooling medium is affected by the design and technical features of the employed PVT set-up, the stirring mode (mechanical agitation or magnetic stirring), the design of the autoclave, the geometry and dimensions of the stirrer, the stirring rate, the contact between the autoclave and the coolant (i.e. direct contact by immersion of the autoclave into a cooling bath or cooling of the cell via circulation of the coolant through a coil), the coolant circulation rate, and the applied cooling/heating rate determine the heat transfer rate. The mass and heat transfer can affect the phase transitions (and thus, the obtained phase envelopes) in the hydrate former-liquid system, gas dissolution rate, supersaturation of the aqueous solution, and the hydrate nucleation process, which in turn determines the induction period for nucleation. Insufficient heat transfer rates can reduce subcooling, which is a driving force for hydrate formation [36].

In what concerns the gas hydrate dissociation process, this essentially involves multiphase heat and mass transfer coupled with intrinsic decomposition kinetics and gas and water flow behavior [28]. The hydrate dissociation rate depends on the intrinsic dissociation rate, heat transfer rate, and mass transfer rate [39]. A number of theoretical models has been developed to predict the intrinsic decomposition rate of gas hydrates [40,41,42,43]. On the basis of comparisons with experimental data on gas hydrate dissociation, Hong et al.[44] and Davis et al.[45] developed theoretical models for heat transfer-limited hydrate dissociation showing that heat transfer plays the dominant role in dissociation, rather than intrinsic kinetics.

Therefore, elucidation of the effects of the technical characteristics of the hydrate formation autoclave and PVT configuration, on the phase behavior of CO₂-H₂O systems and (P, T) conditions for CO₂ hydrate growth, can enable (i) CO₂ hydrates growth at a higher temperature (i.e. shorter metastability region / lower subcooling) and (ii) less experimental time due to the shortening of the induction period. Suitable process modes and conditions can enhance the hydrate nucleation kinetics and facilitate the two aforementioned effects. As a result, less energy- and time-consuming processes can be developed without any usage of kinetic promoters. Also, the potential to form CO₂ hydrate slurries with thermal stability up to 4-5^oC in a reproducible manner, by applying processes that are based on efficient operational modes and suitable PVT configurations, is of practical importance.

The data derived from gas hydrate formation experiments can be compared with suitable Equation of States that calculate the incipient (P, T) conditions for CO₂ hydrate formation (i.e. the conditions at which the first hydrate crystals emerge in the aqueous phase), which coincide with the (P, T) conditions of complete thermal hydrate dissociation. On the other hand, it is interesting to investigate the potential of simple calculation approaches to predict the solubility of CO₂ hydrate even under hydrate formation conditions, by elaborating experimental equilibrium data, and evaluate their prediction accuracy with accurate yet complex thermodynamic models or gas hydrate-related calculation models incorporated in commercial simulation software [46,47,48,49,50]. By calculating the CO₂ solubility and CO₂ fugacity in the gas phase at each equilibrium temperature (using an accurate EoS), the CO₂ mole fractions and partitioning of CO₂ among the gas, hydrate, and aqueous phase can be approximated.

The objective of this work was to investigate the effect of (i) the applied experimental procedures, and (ii) the characteristics of two different PVT configurations on the (P, T) conditions of formation and dissociation of CO₂ hydrates, and on the respective phase envelopes. The experimental procedures involved (i) the continuous temperature cycling mode, and (ii) the incremental (step-wise) temperature cycling mode, in the framework of the Isochoric Pressure Search method. Both experimental modes proceeded through consecutive stages of initial equilibration, cooling, isothermal cooling, and heating of the studied CO₂-H₂O systems. In addition, two different high pressure PVT configurations were used, of which one encompassed a stirred tank reactor and the other incorporated an autoclave of constant volume with magnetic agitation. A qualitative interpretation of the observations has been provided with the aim to gain an insight into the effect of the design and technical features of the two experimental setups on the consecutive steps of hydrates nucleation, crystal growth and thermal decomposition, as well as the form and stability of the hydrate phase. The experimental findings were complemented by two calculation approaches, one of which was a thermodynamic simulation model and the other was a developed semi-empirical calculation method, with the aim to resolve the obtained phase envelopes including CO₂ dissolution and solubility evolution, over the entire range of the applied (P, T) conditions.

2. Materials and Methods

2.1. Chemicals

CO₂ gas with a minimum listed purity of 99.998% (SOL) was used in this work. Deionized water was used after being subjected to a thorough degassing procedure.

2.2. Experimental Setups

The phase behavior and gas hydrate formation conditions of the CO₂-H₂O system have been studied using two different experimental setups. The first high pressure PVT apparatus consists of two basic parts, namely a reactor cell (autoclave) inside an air bath-thermostated cabinet, and a gas preparation manifold. A detailed description of the PVT apparatus has been provided in the work of Kastanidis et al. [51]. The second PVT configuration includes a high pressure stirred tank reactor immersed in a thermostated bath. The reactor shares the same gas preparation manifold with the first PVT apparatus, the description of which is given below.

2.3. Gas Preparation Manifold

The gas preparation manifold is the system in which CO₂ or any other hydrate former can be stored and delivered to the PVT cell under high pressures. It includes (i) a reservoir vessel having a volume of 24 L and maximum operating pressure of 12 bar, four toggle valves, and one pressure transducer (WICA D-10) for pressures up to 20 bar (with a measurement precision of ±0.02 bar), and (ii) a pneumatic gas booster (JULY GB25) bearing a pilot air valve (Swagelok), which is used to receive gas from the reservoir vessel, and deliver it to the PVT cell at several target pressures. The gas booster has a pressure multiplying ratio of 20 in relation to the pressure of the pilot air. Evacuation of the reservoir vessel down to a pressure of 10^-8 bar can be attained by means of a turbomolecular pump (Leybold, Turbovac 361 C), assisted by a rotary pump (Leybold, D16B), that is connected to one of the toggle valves. The three additional toggle valves are used for supplying single pure gases into the reservoir vessel, and are connected to gas cylinders equipped with pressure regulators.

2.4. Description of the First PVT Apparatus

2.4.1. High Pressure PVT cell

The high-pressure PVT cell, which has an internal volume of approximately 325 ml, is constructed from a special alloy (Alloy 20). The PVT cell consists of a cup and a container that can be easily separated from each other. The cup is firmly mounted in an air-bath cabinet. Two feed-through connectors, each one bearing a K-type thermocouple with a measurement precision of ±0.1^oC, are mounted on top of the cup. The thermocouples monitor the temperature of the gas and liquid phases inside the container. On top of the cup there are also two additional feed-through connectors which allow small quantities to be removed from the gas and liquid phase inside the cell via the respective gas and liquid sampling valves and loops. Also, the top of the cup bears an optical window where an external camera (Supereyes™) is mounted. The camera video-monitors the gas/liquid interface during the experiments, allowing to optically invigilate the entire gas hydrate formation and dissociation processes. The bottom of the PVT cell is in contact with a magnetic agitator with controllable rotation speed, which causes the rotation of a magnetic bar to facilitate the mixing of the components of the liquid phase.

2.4.2. Thermostated Air Cabinet

The PVT cell is housed in an insulated cabinet with air circulation. The temperature in the cabinet is precisely controlled by a PID temperature controller (Jumo dTRON 316) and can be kept constant with an accuracy of ±0.1^oC. The insulated cabinet interior also accommodates a high-pressure transducer (WICA D-10) for pressures in the range 0-200 bar with a measuring accuracy of ±0.2 bar, and a pressure transducer for the low-pressure range 0-10 bar (WICA D-10) with a measuring accuracy of ±0.01 bar. In addition, the cabinet contains the sampling loops for the gas and liquid phase with precisely defined volumes, and a pressure relief cell for the liquid sampling. All high pressure (24 MPa) diaphragm valves (Swagelok) are accessible from the front panel of the cabinet. An overview picture of all the equipment located inside the insulated cabinet can be found in the work of Kastanidis et al. [51].

The temperature of the air-bathed cabinet interior as well as of the gas and liquid phases inside the reactor cell are displayed by three digital indicators/controllers (Jumo) located on the front panel. The panel also contains digital readouts of the three pressure transducers attached to the reactor cell, the pressure relief cell for the liquid phase, and the primary vessel in the gas preparation manifold. A thermal bath with a water/ethylene-glycol mixture as coolant, and controlled by an automatic temperature controller (type ECO Silver RE 630S, Lauda), was used to thermostatize the PVT cell with a maximum deviation of ±0.02^oC over the temperature range of our measurements. The coolant from the bath circulates through a thermally insulated cooling coil which passes though the cabinet interior and is firmly mounted around the PVT cell. Both the coil and the cell are firmly enclosed by an insulating jacket. The coil can be easily removed from the PVT cell whenever its container needs to be removed for cleaning.

2.4.3. Monitoring and Data Acquisition

The readings of the two thermocouples for the gas and liquid phases inside the reactor cell, as well as those of the pressure transducers mounted on the reactor and relief cells are monitored, saved and acquired in a data acquisition unit using a virtual instrument application (VI) developed in Labview environment. The communication is performed via an RS232 interface.

2.5. Description of the PVT Configuration Incorporating a Stirred Reactor

2.5.1. Stirred Reactor

A high-pressure stirred tank reactor has been used for the experimental investigation of CO₂ hydrates formation and (H-L_w-G) equilibrium (Figure 1). The configuration includes a Parr 4561 bench top high-pressure stirred reactor made of T316 stainless steel, equipped with a magnetic drive for internal stirring with controllable speed of the stirrer motor (A1120HC drive model), a Parr stirring control unit (model 4565M), a refrigeration unit, a single piston sampling cylinder, and a high-pressure automated water pump. The stirred reactor encompasses a bomb cylinder (cell) with inner free volume of 100 ml and I.D. 2.1 inches, a gas inlet valve, a liquid sampling valve, a gas release valve, a mounted J-type thermocouple (with precision of ±0.3°C over the range 0°C - +200°C) a part of which reaches to the cell interior, a pressure gauge and a burst disc. In all experimental runs, the incoming gas was introduced through the gas inlet valve which is connected to a tube that reaches below the surface of the liquid and near the bottom of the cell.

The pressure of the gas phase inside the cell is monitored by using an electronic pressure transducer (type PXM4202-6KGI, Omega Engineering Inc.) with measurement precision of 0.1 bar, operation range of 0-6000 psig, and analog signal output 4-20 mA, which is mounted on the outlet of the gas release valve. The same automated thermal bath which is connected to the first large PVT apparatus, is also used for thermostating the stirred reactor of the second PVT setup. In this case the small reactor cell is immersed in the bath in direct contact with the coolant.

2.5.2. Data Monitoring and Acquisition System

The temperature of the refrigeration bath and the reactor cell interior are monitored by using a digital thermometer (type P655 instrument, Dostmann electronic) encompassing an A/D converter, that receives the analog signals from the two J-type thermocouples, one mounted inside the reactor cell (as mentioned in Section 2.5.1), and one partially immersed into the coolant in the refrigeration bath. The signals from the P655 thermometer’s digital output are transmitted to a RS232 interface in a PC. The analog readings from the pressure transducer are transmitted to a BNC-2110 A/D card with 15 BNC connectors for analog and digital I/O signal connections that transmits the digitalized signals to the PC through a PCI card. A VI application in the employed NI Labview (32-bit) software was used for data logging, storage, and acquisition.

2.6. Water Supply Configuration

A high-pressure continuous flow, pulse free pump (VP-6K Ambient pump model, Vindum Engineering Inc.) has been used to deliver water into the reactor cell of both PVT setups, with a precisely controllable flow rate and pressure, and accurately monitor the gas phase pressure over the entire course of the experiments via internal pressure transmitters with accuracy +/- 0.1% of full scale. The water amounts were introduced to the reactor cell through the high-pressure single piston sampling cylinder mentioned in Section 2.5.1 (type ProLight Ti-690 64 MB, Proserv), with inlet and outlet valves, which is connected to the water pump. The cylinder has total internal volume of 629 ml and maximum pressure limit of 690 barg. The side of the cylinder that is fixed to the liquid sampling valve of the reactor cell is filled with deionized that has been degassed for 5 hours under vacuum using a rotary pump and intense magnetic stirring. The pump has a maximum pressure rating of 6500 psi, minimum flow rate of 0.00002 ml/min, pressure resolution 0.1 psi, and flow rate and volume accuracy +/- 0.1% of set point. The pump allows for logging and acquisition of water flow rate and pressure data and thus, it can be used complementary to the readings from the pressure transducers mounted on the two PVT setups, which are monitored using the employed NI Labview data recording software.

2.7. Experimental Procedure

The Isochoric Pressure-Search method has been applied in both PVT configurations to determine the (P, T) conditions for the formation and dissociation of gas hydrates. This method allows the construction of the phase diagram of the studied CO₂-H₂O system including the point (P, T) of the complete thermal dissociation of gas hydrates. At first, the gas primary vessel and PVT apparatus were thoroughly evacuated using a turbomolecular pump (Leybold, Turbovac 361 C) assisted by a rotary pump (Leybold, D16B). CO₂ gas was then loaded into the gas primary vessel up to a pressure of 5 bar. The reactor cell was cleaned with methanol and then repeatedly with deionized water, followed by drying using compressed air. The cell was then tightly sealed and thoroughly evacuated along with the gas feed line to remove the atmospheric gases while pressure and temperature data were recorded. After the evacuation step, the temperature of the thermal bath was set to 23-25^oC and the empty reactor cell of the small PVT setup was immersed in the bath, whereas in the large PVT apparatus the temperature of the reactor cell was regulated by the coolant circulating through the cooling coil in firm contact with the cell. Subsequently, pure CO₂ was fed into the reactor cell and pressurized to a selected target pressure, at bath temperature, by using the gas booster described in Section 2.3. The (P, T) data were recorded every 30 s, and after remaining stable for 3 hours, their average values were considered as the feed (P, T) conditions.

Subsequently, a certain volume of deionized and degassed water is loaded into the autoclave through the single piston cylinder which is connected to the high-pressure water pump. Upon completion of water addition to the cell, the contained liquid is stirred either by using a magnetic agitator in the large PVT configuration or via the magnetic drive bearing a four-blades impeller when the small PVT setup is used. After charging the required amount of water to the autoclave and setting the stirring rate, the pressure began to decrease (at constant temperature), as a result of gas dissolution in water, until absorption equilibrium was reached, as indicated by stabilization of the pressure. This way of supplying the gas and water components makes it possible to calculate the exact amount of gas inside the cell by determining the compressibility factor, Z, at the steady feed conditions (P, T) of the hydrate former (whose volume is known), using the Peng-Robinson (P-R) or other Equation of State. In addition, this method allows to correctly determine the H₂O:CO₂ mole ratio in the cell based on the supplied amounts of CO₂ and water.

Accordingly, each CO₂-H₂O binary system was slowly cooled down to 13.7-13.8^oC, and kept under isothermal conditions to establish absorption equilibrium. Then, the system was further cooled to ~0^oC or -1^oC to study CO₂ dissolution and form gas hydrates by applying either an incremental (step-wise, SW) mode (~1^oC temperature drop per reduction step) or a continuous cooling (CC) mode. The system was kept at the lower temperature limit for 24 or 48 hours to study hydration and (H-L_w-G) equilibrium. Subsequently, a step-wise or continuous heating stage was conducted to thermally decompose the formed hydrate phase and close the phase envelope. The duration of the individual isothermal steps in the SW cycling mode was not predetermined but varied depending on the time required for the system to reach absorption (at the cooling stage) or desorption equilibrium (on heating) or gas-liquid-hydrate equilibrium.

The applied H₂O:CO₂ mole ratios (denoted hereafter as MR) ranged between 20.48-21.91. The loading conditions of the two components, experimental parameters, and employed PVT configuration for all experimental runs, are listed in Table 1.

3. Results

3.1. Experimental Pressure vs. Temperature Phase Diagrams

3.1.1. Experimental Runs Performed in Continuous Temperature Cycling Mode

The (P, T) phase diagrams derived from the experiments conducted in the continuous temperature cycling mode are presented in Figure 2. The steep and drastic pressure drop steps during the cooling stage, which imply the extensive formation of CO₂ hydrates, were affected by the cooling rate and the presence of perforated plates on the rotation axis of the magnetic drive, in what concerns (i) the onset temperature of each abrupt pressure reduction step, and (ii) the pressure reduction profile that led to the final (P_eq, T_eq) conditions of the (H-L_w-G) equilibrium. The pressure reduction reached to the same equilibrium pressure (12.07, 12.13, and 12.58 bara) at the lowest applied temperature limit of 0.02, 0.07, and 0.24^oC, (i.e. the same (P_eq, T_eq) equilibrium point, as is also shown in Figure S1). The equilibrium pressures were close to the P_eq values at similar temperatures which have been referred in the literature (P_eq = 12.7 bara at 0.35^oC [52]; P_eq = 13.8 bara at 0.7^oC [23]).

It is worthy to mention that the faster applied cooling rate (1^oC/2h) significantly reduced the subcooling and induction period of the CO₂ hydrate nucleation stage, thereby accelerating the initiation of the subsequent stage of hydrate crystals growth. The latter stage is related to the extensive absorption and transfer of CO₂ from the gas phase to the hydrate crystals, which accounts for the observed drastic pressure reduction.

In addition, a similar (P-T) phase diagram was obtained from a repetition experiment, when upon completion of the experiment with MR_20.21, the same CO₂-H₂O system was initially equilibrated at room temperature for 48 hours and was then subjected to the same continuous temperature cycling process, with the same cooling-heating rate of 1^oC/2h. The (P-T) phase diagrams of the two experiments are compared in Figure S2. The small difference in the initial pressure between these two runs is due to CO₂ absorption during the initial equilibration step in the repetition experiment. As it can be observed, the usage of perforated plates facilitated the formation and growth of hydrate crystals since the surface area of these plates provided an additional large number of nucleation sites and thus, it served as kinetic promoter by reducing the induction time of the nucleation phase which leads to hydrates crystal growth.

Regarding the heating curves, in all experiments the pressure increased as soon as the heating stage commenced and kept increasing smoothly up to the point where the heating curves joint the respective cooling curves. The evolution of pressure with raising temperature was common for the three runs. The observed pressure evolution might be indicative of a heat-sensitive hydrate phase, that consists of isolated hydrate crystals dispersed within the bulk aqueous liquid, the external layers of which are gradually decomposed with increasing temperature. This dissociation releases CO₂ molecules into the aqueous liquid in a continuous but not abrupt manner, which are subsequently desorbed into the gas phase. Furthermore, the temperature at which the heating curves (i.e. hydrate thermal dissociation) and cooling curves joint, ranged between 8.15-8.45^oC.

3.1.2. Experimental Runs Performed in Incremental Temperature Cycling Mode

The pressure and temperature data acquired during the cooling, isothermal and heating stage, from the runs performed with the studied CO₂-H₂O systems via the SW-mode, have been used to identify the (P, T) conditions to reach (H-L_w-G) equilibrium, and construct the phase envelope of the respective systems. The (P, T) phase diagrams are presented in Figure 3, whereas Figure 4 depicts the experimental (P, T) results at equilibrium conditions (either (L_w-G) or (H-L_w-G) equilibrium). The respective phase envelopes are formed from the equilibrium data acquired at the cooling and heating stage. The SW-mode experiments were conducted using feed H₂O:CO₂ mole ratios within the range of 20.48-21.91.

Regarding the stirred tank reactor, with the exception of one experiment performed with MR_21.37, in the other three runs no hydrate formation has taken place at all during the entire cooling period, as observed in Figure 3a and Figure 4a. In these three runs, the CO₂-H₂O system at the cooling stage did not undergo any transition from a two-phase (L_w-G) to a three-phase (H-L_w-G) system, even though (P, T) conditions were applied at which this phase transition should have taken place. In opposite, regarding the runs performed in the large PVT apparatus, application of similar H₂O:CO₂ feed mole ratios allowed for CO₂ hydrate formation, as observed in Figure 4b. Formation and growth of gas hydrates was deduced by the steep and drastic pressure reduction at the cooling stage. CO₂ hydrate crystals exhibit substantially higher capacity for CO₂ absorption compared to liquid water.

The situation where only the gas and aqueous phase coexist under hydrate formation conditions, corresponds to a metastable equilibrium between the gas phase and a super-saturated aqueous liquid, which contains a large population of dispersed hydrate nuclei. This metastability is associated with a supercooled liquid where the hydrate nuclei aggregation/clustering process that normally leads to a necessary critical size for crystal growth, involves a relatively long induction period, because of prohibited clustering of hydrate nuclei. As a consequence, the subsequent step of hydrate crystal growth cannot commence during the successive cooling steps. Hence, at the cooling stage, the supersaturation of the aqueous phase was inadequate to rapture the metastability and force the gas-liquid system to undergo the phase transition that would lead to growth of stable hydrate crystals.

For the runs with a H₂O:CO₂ mole ratio of 21.66 and 21.91, the pressure drop became rapid and prominent below ~3^oC. As observed in the phase envelopes of the experiments with MR_21.37, MR_21.66 and MR_21.91 (Figure 3a,b), the temperature difference between the point of complete hydrate decomposition and the onset of hydrate crystal growth (i.e. the degree of subcooling) is either moderate, as in the two latter runs (ΔΤ_sub = 2.74^oC and 4.18^oC respectively), or high, as in the experiment with MR_21.37 (ΔΤ_sub = 6.74^oC). The onset temperature of hydrate crystals growth (i.e. rapture of metastability) has been taken as the temperature of the sudden heat release (temperature jump).

Further evidence for gas hydrate formation is provided by the first derivative of pressure with respect to temperature (dP/dT) vs. temperature curves at the SW cooling stage, as shown in Figure S3. In the case of the experiments with MR_21.37 and MR_21.91, an intense maximum is observed at 0.46 and 2.17^oC respectively, which starts to emerge at 1.34 and 4.29^oC. Another important finding is the different profiles of the heating curves obtained from the two PVT configurations. This dissimilarity could arise from the different internal design and liquid agitation mode of the crystallization autoclaves, which affect the mass and heat transfer rates (lower temperature reduction rate in the large PVT autoclave compared to the cell of the stirred reactor, for the same cooling rate), as will be discussed later.

It is mentioned that the SW-mode proceeds though consecutive equilibrium steps of gas dissolution, hydrate formation, hydrate dissociation and gas desorption, which result in prolonged processes. On one hand, application of intense agitation will increase the gas dissolution rate and thus, the concentration of gas hydrate nuclei within the bulk liquid and/or on the gas/liquid interface during the cooling stage. However, it is possible that excessive agitation and shear, especially for quite long periods, destabilize and decompose the clusters of gas hydrate nuclei, thereby hindering their aggregation that would initiate the growth of hydrate crystals. As a consequence, the subcooling degree and the nucleation induction time will increase and the gas/liquid system may even be unable to undergo the phase transformation needed to establish (H-L_w-G) equilibrium. The hindrance of the hydrate nuclei clustering is possibly the reason for the lack of a CO₂ hydrate phase for the experiments with MR_20.48, MR_21.37 (repetition experiment), and MR_21.67. For the run with MR_21.37, the conditions at which (H-L_w-G) equilibrium was established were identified at (i) T = 0.46^oC and P = 12.83 bara, (ii) T = -0.25^oC and P = 12.4 bara.

3.2. Correlation of the Pressure Evolution During the Cooling Stage with Crystallization Effects

In is known that formation of CO₂ hydrates is a two-phase process that includes (i) formation of fine nuclei consisted of CO₂ molecules enclathrated within ice-like crystalline lattices (crystal seeds), which grow in size via clustering, and (ii) growth of bulky hydrate crystals [53]. The duration of the nucleation phase is defined as the induction period. During this period the aqueous phase is supersaturated with CO₂. The nucleation phase is ended by the formation of stable hydrate nuclei clusters of critical size. As soon as the growing nuclei exceed the critical size, they begin to grow spontaneously as bulky hydrate crystals, and in this way the crystal growth phase commences. The greater the subcooling / supersaturation, the greater the driving force for hydrate nucleation and the rate of hydrate crystal growth.

The effect of the subcooling degree on the extent and rate of CO₂ hydrate growth and associated pressure drop, can be assessed by the phase envelopes of the runs with MR_21.37, MR_21.66, and MR_21.91. In the first run, the evolution of pressure during the cooling stage was quite different from the case of the two latter runs. Since the three experiments were performed using the same stirring rate (400 rpm) and at nearly the same initial pressure (31.04-31.40 bar, see Table 1), the sharper pressure reduction that occurred in the first run can be associated with the greater subcooling and induction time. This ensues larger driving force for spontaneous growth of hydrate crystals and also higher crystal growth rate and thus, higher pressure reduction rate. Similar findings regarding the effect of subcooling on the hydrate growth rate have been reported in the literature, although cases with large variability in induction times within a single experimental setup have been also reported [54,55].

On the other hand, the larger induction time and subcooling can be attributed to the intense shear or turbulence imposed in the aqueous liquid by the rotation of the impeller, that could mechanically extort CO₂ molecules for the hydrate nuclei dispersed in the liquid phase, especially when prolonged cooling stages are applied. This would prevent the formation of stable hydrate nuclei clusters of critical size, until the driving force for this clustering at a sufficiently low temperature, eventually overcomes the mass transfer limitation. It is noted that for all runs performed in the stirred tank reactor by applying the SW mode, the experiment with MR_21.37 was the only case where a CO₂ hydrate phase was formed. This is indicative of the hindering effect of the excessive shear / turbulence, and prolonged process time, in the capability of CO₂-H₂O systems to undergo the phase transition that will lead to the formation of a gas hydrate phase, despite the higher heat flux.

In the two latter runs that were performed in the large PVT setup, the pressure reduction became sharp below ~3^oC (for MR_21.91) or ~4^oC (for MR_21.66). It is likely that the moderate mass transfer imposed by magnetic stirring did not hinder aggregation of hydrate nuclei dispersed in the aqueous liquid - this aggregation is an essential process for hydrate crystal growth - since no excessive shear and turbulence conditions were applied in the liquid to destabilize clustering of hydrate nuclei. This would result in the observed lowering of the subcooling degree compared to the high ΔT_sub that was observed when the stirred tank reactor was used (experiment with MR_21.37).

Besides the abrupt pressure reduction in the experiments with MR_21.66 and MR_21.91, a rapid switching in the appearance of the aqueous phase was observed, from a lucid solution to an opaque one. The surface of the stirred liquid became essentially opaque to optical observation by the mounted camera (rather suddenly) once the autoclave was cooled below nearly 3^oC, indicating the formation of a non-transparent hydrate layer. Several optical microscopy studies have revealed that during the hydrate crystals growth phase, a CO₂ hydrate film can be formed and expand laterally along the gas/liquid interface [56,57]. The lateral growth rate of this hydrate film also depends on the system’s subcooling ΔT_sub.

Prior to the sharp pressure reduction that occurred at the cooling stage for the runs with MR_21.37, MR_21.66, and MR_21.91, an instantaneous increase (jump) of the liquid’s temperature was observed, which coincided to the onset of the drastic pressure drop. Natarajan et al.[58] have reported a slight but detectable rise of the aqueous solution temperature at the turbidity point because of the heat released by the sudden formation of stable hydrate nuclei. The sudden increase in temperature can be attributed to a rapid initiation of the growth phase of large amounts of stable hydrate crystals, which is an exothermic process. It is noticed that the temperature jump was more pronounced when it took place at a lower temperature and thus, greater subcooling, possibly due to the larger population of hydrate nuclei clusters of critical size in the aqueous phase. This would cause an abrupt release of larger amounts of hydration heat during the initiation of the phase of hydrate crystals growth.

Apart from the effect of the utilized PVT setup on the subcooling degree and the onset of the temperature jump, the temperature profile of the steep pressure drop phase, can also be explained on the basis of the mass and heat transfer characteristics of each experimental apparatus. Although the excessive agitation in the stirred tank reactor and the direct contact of the reactor with the coolant result in enhanced heat transfer rates, that maintain constant temperature, the release of substantial amounts of hydration heat can impose a temporary increase of both temperature and pressure inside the reactor, thereby introducing a negative effect, as can be observed in Figure 3a.

3.3. Thermal Dissociation of the Hydrate Phase

The thermal dissociation stage of the hydrate phase was studied with regard to the evolution of the gas pressure with increasing temperature. In the experiments where the system established (H-L_w-G) equilibrium (i.e. those with MR_21.37, MR_21.66 and MR_21.91), the pressure exhibited two different behaviors. In the run with MR_21.37, the pressure did not change up to 0.8^oC, whereas in the other two experiments it remained either constant or even further decreased as the system was heated from the lower temperature limit up to ~4^oC, thereby demonstrating the thermal stability of the CO₂ hydrate phase up to this threshold. However, it is likely that the stability of the pressure over this temperature range is due to the growth of a thick dense hydrate layer on the gas/liquid interface rather than to a slow heat transfer and dissociation rate of large CO₂ hydrate crystals dispersed in the bulk liquid. This hydrate layer acted as a practically impermeable cap separating the gaseous from the aqueous phase, thus preventing the desorption of released CO₂ molecules from decomposed hydrate crystals below the solid layer/liquid interface towards the gas phase. Once the temperature exceeded 4^oC, the gradual decomposition and pore opening of the hydrate barrier layer, enabled the release of considerable amounts of CO₂ into the gas phase and this decomposition-desorption process continued up to 8^oC.

The prominent first derivative dP/dT peak at 4.58^oC for the experiment with MR_21.91 (7.09 bar/^oC), which implies the considerable release of CO₂ amounts into the gas phase (which started between 3.76-4.36^oC), is indicative of the dense hydrate cap assumption. At the temperature of ~8^oC the phase envelopes closed, similar to the experiment with MR_21.37. For all three experiments the “hydrate decomposition point” at which the thermal decomposition curve joins the cooling curve, was reproducible to within 0.37^oC (8.08^oC for the MR_21.37, 8.41^oC for the MR_21.66, and 8.04^oC for the MR_21.91).

3.4. Phase Transformations and Equilibrium at CO₂ Hydrate Formation Conditions

As mentioned previously, during the crystal growth phase, a rapid decline in pressure occurs owing to the extensive CO₂ absorption by the hydrate phase. The spontaneous growth of hydrate crystals leads to the formation of a thin hydrate film that propagates along the gas/liquid interface. Eventually the hydrate film covers the entire interface and then, it grows in thickness into the aqueous phase-side. It has been reported that the size of the hydrate crystals decreases and the lateral film growth rate increases with greater supercooling [54,55]. An explanation for the effect of supercooling on the lateral growth rate of the hydrate films along gas / water interfaces has been provided in literature by a film growth model involving convective transfer of the heat of hydrate formation, from the edge of the film to the surrounding phases [59].

There are two sources of CO₂ molecules that contribute to the lateral growth and thickness of a gas hydrate film; these are the gas phase and the bulk aqueous phase. Mori et al.[60] and Sugaya et al.[61] considered the hydrate film to be a thin water-permeable hydrophilic crystalline plate with numerous microperforations which are much larger than the molecular dimensions of the hydrate formers and water. In each microperforation liquid water diffuses continuously by capillary forces toward the gas phase-side surface of the hydrate layer, resulting in successive formation of hydrate crystals on this surface and therefore, in thickening of the hydrate film [62]. At the same time, the opposite surface of the hydrate film is being eroded as a result of dissociation of the hydrate crystals exposed to the bulk aqueous phase. Dissociation is enhanced with increasing rate of diffusive or convective removal of the gas molecules from the aqueous phase-side of the hydrate layer, which is inherently dependent on the flow of water relative to the film. The film thickness is inversely proportional to the mass transfer coefficient on the aqueous phase-side of the film [63].

In contrast, at low subcooling and low supersaturation levels, the population of hydrate nuclei is lower, and the lateral growth rate of the hydrate layer is controlled by the growth rate of a single crystal. As a result, the hydrate film surface becomes rougher and the lateral growth proceeds at lower rate [64,65]. The lower lateral growth rate of the hydrate layer provides the crystals in the film front with sufficient time to grow to 3D dimensions, which leads to the formation of a thicker and coarser hydrate film. The wider inter-crystalline gaps (film capillaries), in conjunction to the coarser internal texture (that provides larger hydrophilic surface area) favor the transfer of water and CO₂ molecules from the gas phase, which come into contact on the pore mouths. Eventually, if the temperature becomes quite low, the width of pores may substantially decrease, and then the hydrate film will become practically impermeable to gas permeation. In this case, the hydrate film will separate the two fluid phases from each other and as a result, the pressure decrease will cease and the system will reach a steady state, which however is not identical to the actual (H-L_w-G) equilibrium.

Indeed, if the phase diagrams of the experiments with MR_21.37, MR_20.09 (stirred tank reactor, Figure 3a and Figure 4) and MR_21.66, MR_21.91 (large PVT setup, Figure 3b) are compared, it can be deduced that the larger mass transfer and shear, established in the mechanically agitated tank reactor, prevented the formation of a mechanically and thermally stable hydrate layer. Since a dense hydrate layer could hardly appear on the gas/liquid interface under such conditions, the bulk liquid was accessible to CO₂ dissolution from the gas phase, until (H-L_w-G) equilibrium was eventually established and the pressure drop ceased, as observed in Figure 4a. This was also the reason for the larger reduction of the normalized pressure, as seen in Figure S3a. The structural instability of the CO₂ hydrate layer could account for the fact that it started decomposing gradually almost as soon as the temperature of the system was raised (MR_21.37).

On the other hand, in the case of the large autoclave, the moderate magnetic stirring resulted in slower mass transfer between the hydrate crystals and the surrounding fluid liquid. Thus, destabilization of the hydrate nuclei clusters was prevented, and enough time was provided to the individual clusters to continue enclathrating CO₂ molecules and grow in dimensions up to a critical size. Such operational conditions tended to induce a rather moderate subcooling with shortening of induction period, and produce a thicker hydrate layer with low porosity on the gas/liquid interface, which acted like an impermeable cap, separating the gas phase from the bulk liquid phase [63,64,65]

As a result, the system reached a steady state that can be far from the theoretical (H-L_w-G) equilibrium. This case is different from what was observed in the run with MR_21.37, at which once the rapid consumption of CO₂ was completed and the pressure reached a constant limiting value, the system reached the (P, T) conditions for (H-L_w-G) equilibrium.

3.5. Comparison between Experimental and Computational Results on the Phase Behavior of The studied CO₂-H₂O Systems

3.5.1. Computational Results on the Incipient Conditions for CO₂ Hydrate Formation

In addition to the experimental investigation of the phase behavior and hydrate formation conditions, a commercial thermodynamic software (Multiflash v.6.1) has been used to computationally determine the incipient pressure for CO₂ hydrate formation, P_inc, and CO₂ solubility in the aqueous phase, at each equilibrium temperature, T_eq, for all experiments. The simulation model used for these calculations was the CPA Infochem Equation of State, which is available in the simulation software. Figure 5 depicts a collocation of the experimental (P_eq, T_eq) equilibrium curves for the cooling stage with the (P_inc, T_eq) curves that have been calculated by applying the CPA Infochem model and for the same equilibrium temperatures.

According to the simulation results presented in Figure 5, there is a temperature threshold at which the experimental equilibrium (P_eq, T_eq) curve for each run intersects the respective calculated (P_inc, T_eq) curve, below which the experimental pressures exceed the predicted incipient pressure at any equilibrium temperature, T_eq. This intersection point is detected at a different temperature depending on the applied H₂O:CO₂ mole ratio. For all experiments the intersection point ranges between ~7.9-8.38^oC and thus, it coincides to the “hydrate decomposition point” of the respective phase envelopes. The computational results in the case of the run with MR_21.37 showed also two more intersection points at the experimental conditions where (H-L_w-G) was established (see Section 3.1.2).

At the experimental pressures above the calculated incipient ones (at same temperatures), the studied CO₂-H₂O systems are two-phase gas hydrate-water rich liquid (H-L_w) systems, and thus, the gas phase should have been completely depleted and CO₂ would only be distributed between the hydrate phase and the aqueous liquid. The conditions for (H-L_w-G) equilibrium in CO₂-H₂O systems with H₂O:CO₂ mole ratios of 20:1, 30:1, and 40:1 have been determined by applying the CPA Infochem model, and the calculation results are shown in Figure S4. The temperature range in which (H-L_w-G) equilibrium can be attained, varies from 271.81-283.31 K for the MR_20:1 and MR_30:1, and from 271.81-280.37 K for the MR_40:1. Figure S5-S7 in Supplementary Materials presents density and volume values of the aqueous phase, as well as change of volume (%) of the aqueous phase, at the (P_eq, T_eq) equilibrium conditions of the experiments, that have been calculated using the CPA Infochem model.

3.5.2. CO₂ Solubility Calculations

CO₂ solubility in the aqueous phase at the applied (P_eq, T_eq) equilibrium conditions, has been calculated using three approaches: (i) calculation of the Henry coefficient for all equilibrium temperatures by a fitting equation, which was derived from an extended data base of experimental measurements of CO₂ solubility in water, combined with the calculation of CO₂ fugacity at absorption equilibrium using the Peng-Robinson EoS, (ii) the CPA Infochem model of the Multiflash v.6.1 software, and (iii) the simple subtraction of the CO₂ amount in the gas phase - at each phase equilibrium point (P_eq, T_eq) - from the initial CO₂ amount supplied in the autoclave at the feed conditions, using the Peng-Robinson EoS. The CO₂ solubility values calculated by the three aforementioned approaches were compared regarding their mutual deviations and with respect to their temperature dependence over the entire cooling stage (including superaturation conditions), and at the (H-L_w-G) equilibrium stage. The evolution of CO₂ solubility in the aqueous phase with decreasing temperature in the presence of hydrates was resolved, and the results of the three calculation approaches will be discussed.

3.5.2.1. Henry Coefficients at the Equilibrium Temperatures

Carroll et al.[66] performed a thorough review of the research literature on CO₂ dissolution in water and found approximately 80 experimental studies that reported dissolution equilibrium data for pressures lower than 10 bar. The authors developed a model based on the Henry's law and the experimental solubility data derived from the literature, to calculate the Henry's coefficient values, and allowed for the data to be compared and correlated. The developed model related the calculated Henry's coefficients (expressed in terms of mole fraction) as a function of absolute temperature, as follows:

$$ln{H}_{21}=-6.8346+1.2817\times {10}^4/T-3.7668\times {10}^6/T^2+2.997\times {10}^8/T^3,$$

(1)

where H₁₂ is the Henry coefficient of the system H₂O (solvent, component 1) - CO₂ (solute, component 2), in MPa/(mole fraction) and T is the absolute temperature. This correlation was found to be valid for the temperature range between 0 and 160^oC. The Henry coefficients at any experimental temperature can be calculated by applying the fitting equation.

The maximum total pressure where the experimental data from the literature were taken into account to produce the above fitting equation, was limited to about 10 bar because this model ignores the deviation from the ideal behavior of the liquid phase, which at pressures greater than 10 bar starts to become significant. At higher pressures, the solubility of CO₂ increases to the point where the activity coefficients can no longer be ignored. Also, at pressures above 10 bar, the Poynting correction (the effect of pressure on the reference volatility) becomes significant. The error resulting from neglecting the Poynting factor at a pressure of 10 bar was estimated to be about 1%. For a gas/liquid solvent system and a given temperature, the Henry coefficient expresses the solubility of the gas and is defined by the following relationship:

$$H_{21}=\underset{x_i\to 0}{\lim }{\displaystyle \frac{y_i{f}_i}{x_i}}=\underset{x_i\to 0}{\lim }{\displaystyle \frac{y_i{\widehat{\phi }}_{i}P}{x_i}}$$

(2)

where x_i: mole fraction of gas component i in the liquid phase, y_i: mole fraction of gaseous component i in the gas phase, f_i: fugacity of component i in the gas phase, φ_i: the fugacity coefficient of component i in the gas phase, and P: total pressure. The product y_i·P is the partial pressure (p_i) of component i (in this case CO₂) in the gas phase, x_i = x_CO2,aq, and y_i·f_i = φ_i·p_i. It is mentioned that in order for the calculated Henry constants to be valid, the mole fraction of the dissolved gas must be quite small (normally x_CO2,aq → 0 and φ → 1).

Subsequently, the CO₂ solubility curves which are calculated by applying the fitting equation (1) coupled with the Peng-Robinson EoS for CO₂, and the Henry coefficient equation, that is H₂₁ = φ_CO2·p_CO2 / x_CO2,aq, can be compared with simulation results for CO₂ solubility in water at various (P, T) conditions. In particular, the equilibrium (P, T) conditions of the experimental runs will be the selected parameters for comparison purposes.

Figure 6a depicts the effect of temperature on CO₂ solubility. The CO₂ mole fractions in the liquid have been calculated using the fitting equation (1) to determine the Henry constant for each temperature and the Peng-Robinson EoS to calculate the fugacity at the respective experimental equilibrium pressures. Since the fitting equation (1) has been produced by collecting a large number of CO₂ solubility data published in literature, the two calculation models (CPA Infochem and P-R EoS) can be indirectly validated with regard to their agreement with the experimental results for CO₂ absorption in water reported in literature.

3.5.2.2. Calculation of CO₂ Dissolution in Water Using the Simulation Model and Comparison with the Calculation Results Derived from the Fitting Equation-Based Approach

Figure 6b presents the CO₂ solubility, calculated by applying the CPA Infochem model, as a function of temperature, for all applied MRs. The simulation model calculates the amount (in moles) and the mole fraction of CO₂ in any phase in equilibrium. Figure S8 presents the deviations in the calculated values of CO₂ solubility in water as a function of temperature, between the CPA Infochem EoS and the H₂₁ vs. T fitting equation (1). The deviations are given as mean absolute percentage deviation (MAPD). A reasonable agreement between the fitting equation (1) and the thermodynamic model, especially at the higher equilibrium temperatures, can be deduced. A small deviation of the calculated results between the two approaches emerges at the lower temperature regime. Therefore, the simple calculation approach that is based on the fitting equation (1) can be useful for simple and rough predictions of the CO₂ solubility in CO₂-H₂O systems (and mass balances as well), even under gas hydrate formation conditions. Usage of a more complex EoS will add complexity to this approach without significantly enhancing the accuracy.

Furthermore, in the experiment with MR_21.37, the CO₂ amount (moles) that remained in the gas phase at 0.46^oC, calculated by using the CPA Infochem EoS, reached to 16.1% of its initial amount at the feed. At the lowest temperature of -0.25^oC, the respective CO₂ amount was slightly reduced to 15.65%. It is reminded that these two equilibrium temperatures correspond to (H-L_w-G) equilibrium conditions, as mentioned in the Section 3.1.2. In the experiment with MR_21.66, the remaining CO₂ amount in the gas phase at -0.32^oC, reached to 27.41% of its initial amount at the feed, even though the initial gas pressures in these two runs were almost equal (31.4 and 31.04 bara respectively, see Table 1). The difference in the remaining CO₂ amount in the gas phase between these two experiments can be attributed to the assumption of the growth of an impermeable hydrate layer on the gas/liquid interface and the subsequent establishment of a steady-state condition in the experiment with MR_21.66, which prevented the CO₂-H₂O system from reaching the actual (H-L_w-G) phase equilibrium.

3.5.2.3. Calculation of CO₂ Dissolution in Water Using the Peng-Robinson EoS

The CO₂ solubility values that have been calculated using the approaches described above (i.e. H₂₁ vs. T fitting equation, CPA Infochem model, and P-R EoS), as a function of temperature, are presented in Figure 7. The presented calculation results refer to both PVT configurations. In almost all cases, the CO₂ solubilities calculated by using the P-R EoS clearly diverge from the results obtained from the other two calculation approaches especially under the conditions at which formation of hydrates occurred.

Furthermore, the results reveal a decline in CO₂ solubility in the presence of the hydrate phase, as the temperature decreased. This reduction was more pronounced in the experiment with MR_21.37 that exhibited the greater supercooling degree (minimum x_CO2,aq ~ 0.015) and the most sharp and pronounced pressure drop. The decline of the CO₂ mole fraction in the aqueous phase is a result of enclathration of large amounts of gas at the expense of the gas dissolved in the super-saturated aqueous liquid. Decrease of CO₂ solubility with decreasing temperature over the hydrate formation region, has been reported previously by other research groups [24,67].

4. Conclusions

In this work the phase behavior, and gas hydrates formation and dissociation conditions of CO₂-H₂O systems have been experimentally investigated with respect to (i) the applied temperature cycling procedures and (iii) characteristics of the experimental PVT configuration. The operational procedures involved (i) the incremental temperature cycling mode and (ii) the continuous temperature cycling process, in the framework of the Isochoric Pressure Search method. The (P, T) conditions required to reach the (H-L_w-G) equilibrium depend on the temperature cycling mode applied as well as the design and technical characteristics of the reactor cell. The incremental temperature cycling mode imposes longer induction times and higher subcooling when the experiments are performed using the agitated tank reactor, compared to the magnetically stirred autoclave.

Qualitative considerations on the subcooling, induction time, and thermal stability of the developed hydrate phases allowed to gain an insight on the effect of the technical parameters of the two PVT setups (i.e. agitation method, internal design, mass and heat transfer coefficients, shear, perforated plates), on the growth rate thickness, and porosity of the hydrate layer on the gas/liquid interface. A thin permeable hydrate layer will allow for the system to reach actual (H-L_w-G) equilibrium whereas a thick dense layer with low porosity may act as an impermeable barrier separating the two fluid phases from each other. As a result, the system may eventually reach a steady state where the gas phase is isolated from the aqueous/hydrate phase, which is different from the actual (H-L_w-G) equilibrium.

In what concerns the continuous temperature cycling mode, the faster applied cooling rate significantly promoted the CO₂ hydrates growth. Also, the usage of perforated plates reduced the induction time of the nucleation phase and facilitated the formation and growth of hydrate crystals by providing an additional large number of nucleation sites on their surface. Moreover, the obtained heating curves joint the cooling curves at temperatures that ranged between 8.15-8.45^oC, similar to the cases of all runs performed by applying the SW operational mode which exhibited CO₂ hydrae formation, regardless of the employed PVT configuration. On the other hand, the profile of the heating curves in the CC operational mode was also similar to the heating curve of the run with M_21.37 that was performed in the same autoclave, i.e. the stirred tank reactor, indicating that usage of this type of hydrate formation set-up can only result in formation of a thin hydrate layer which is not thermally and mechanically stable, in opposite to the dense hydrate layer that is grown when the magnetically agitated PVT autoclave with moderate mass transfer is used.

The simulation results corroborated the experimental findings in what concerns the hydrate decomposition point, the metastability of the studied systems, and their ability to undergo the phase transformation that leads to CO₂ hydrates formation. CO₂ solubility at the equilibrium (P, T) conditions has been resolved using two calculation approaches: (i) calculation of the Henry coefficient for the experimental temperatures by combining a fitting equation, which was derived from an extended data base of experimental measurements of CO₂ solubility in water, with the Peng-Robinson EoS-based calculation of the CO₂ fugacity at the equilibrium (P, T) conditions, and (ii) the CPA Infochem model of a commercial simulation software. The results revealed a continuous decline in CO₂ solubility with decreasing temperature in the presence of CO₂ hydrates while the opposite trend was found in the two-phase (L_w-G) regime. Thus, the contribution of each one of the CO₂ sources (i.e. the aqueous phase and the gas phase) on CO₂ hydrates formation was assessed.