Next Article in Journal
A FEM Structural Analysis of a Francis Turbine Blade Parametrized Using Piecewise Bernstein Polynomials
Previous Article in Journal
Simulation of Multi-Phase Flow to Test the Effectiveness of the Casting Yard Aspiration System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Theoretical Investigation of Supercritical Processes: Kinetics of Phase Transitions in Binary “2-Propanol—CO2” System

Department of Chemical and Pharmaceutical Engineering, D. I. Mendeleev University of Chemical Technology of Russia, 125047 Moscow, Russia
*
Author to whom correspondence should be addressed.
Computation 2023, 11(7), 122; https://doi.org/10.3390/computation11070122
Submission received: 29 May 2023 / Revised: 19 June 2023 / Accepted: 19 June 2023 / Published: 22 June 2023
(This article belongs to the Section Computational Chemistry)

Abstract

:
Studies of phase transition kinetics are important for such supercritical processes as supercritical drying, adsorption, micronization, etc. In supercritical technologies, “organic solvent—CO2” systems are often formed, the properties of which strongly depend on the system parameters. In this article, the kinetic curves of phase transitions in the “2-propanol—CO2” system were investigated experimentally and theoretically. Experimental studies were carried out in a 250 mL high-pressure apparatus at temperatures of 313 and 333 K and pressures of 6.3 and 7.8 MPa with and without the addition of alginate porous gel. Theoretical studies were carried out using the mass transfer equation, the Peng-Robinson equation of state, and the Van der Waals mixing rules, with Python being used for the calculations. The mass transfer coefficients and equilibrium concentrations of CO2 in the liquid phase were determined using the BFGS optimization method.

1. Introduction

Currently, supercritical technologies are a promising emerging field that finds applications in a wide variety of industries [1,2,3]. Carbon dioxide is mainly used as a supercritical fluid in supercritical processes (SC-CO2) [4,5] since it is non-toxic, non-flammable, cheap, and readily available. Moreover, it has low critical parameters (Tc = 31 °C, Pc = 7.3 MPa) [6]. If an organic solvent is used in these processes, a binary “organic solvent—CO2” system is formed. To gain a deeper understanding of mass transfer processes under sub- and supercritical conditions, it is necessary to study the phase behavior and kinetics of phase transitions in such systems. In many cases, it is also important to consider a system in the presence of porous material. New relationships for mass transfer help to intensify the supercritical extraction process by using phase behavior and kinetics data of the binary “extracted substance—CO2” system [7,8]. Similar data plays an important role in optimizing and scaling up different supercritical processes. In the micronization by the supercritical antisolvent (SAS) method, the “organic solvent—CO2” system is formed. The phase behavior of this system affects the atomization process and, consequently, the final properties of the resulting material [9,10]. Additionally, data about the kinetics of phase transitions and phase behavior are important in the supercritical drying process, which is the main stage in the production of highly porous materials such as aerogels [11,12]. To produce aerogels, the organic solvent is removed from the gel by supercritical drying without its structure collapsing or shrinking. This is possible due to the supercritical state of the “organic solvent—CO2” system formed during the process. This system is homogeneous at process parameters; thus, no phase boundary occurs inside the gel, and there are no capillary forces in the gel’s porous structure. The specified condition determines the choice of an acceptable organic solvent. The most commonly used solvents are methanol, ethanol, 2-propanol, acetone, hexane, etc. [13,14,15,16]. The supercritical drying process can be divided into four stages: during pressurization, there is a pressure increase and the phase transitions of the binary “organic solvent—CO2” system. The next steps are the displacement of the organic solvent from the free volume of the apparatus and the further diffusional exchange of the organic solvent inside the gel with SC-CO2. The final step is isothermal depressurization.
The pressurization step is an important part of the supercritical drying process, since as the pressure in the system increases, the volume of the liquid phase increases due to the dissolution of CO2 in an organic solvent; otherwise, spillage may occur. This phenomenon affects the transport of organic solvent from the pores of the gels, as a significant part of the solvent leaves during this step. A detailed investigation of the pressurization step of the supercritical drying process was introduced in [17]. The article shows the influence of the parameters and properties of the gel on the kinetics of spillage. These data confirm the importance of studying the kinetics of phase transitions for supercritical processes. The behavior of the “organic solvent—CO2” system at subcritical conditions depends on the solubility of CO2 in the organic solvent. For example, as the temperature rises to the critical point, the solubility in organic polar solvents decreases in accordance with Le Chatelier’s principle [18,19]. Thus, the dissolution process is exothermic [20]. However, with increasing pressure, the solubility of CO2 in the liquid phase increases [21]. When the pressure is increased to critical parameters of the system, the distinction between liquid and gas phases disappears, and CO2 becomes indefinitely mixed with the solvent.
Different binary systems of “organic solvent—CO2” have similar phase behavior. For example, articles [22,23,24,25,26] provide data on phase equilibrium at various parameters of the binary system “2-propanol—CO2,” which is often used in supercritical processes [17,27]. As the pressure in the system increases, the composition of the gas phase almost does not change, in contrast to the composition of the liquid phase. This is due to the increase in solubility of CO2 in 2-propanol with pressure. There is a noticeable increase in liquid phase volume with pressure and corresponding mole fraction of CO2 increasing [17]. Figure 1 demonstrates such an increase for the “2-propanol—CO2” system as a function of pressure. The curves were calculated using phase equilibrium data and the Peng-Robinson equation of state (PR-EoS) with Van der Waals mixture rules [28].
From the presented data (Figure 1), it can be seen that the curves asymptotically approach the critical pressures of the mixture at different temperatures. At these points, the difference between the two phases disappears, and the system becomes homogeneous.
These data are determined by the thermodynamic state of the system and do not contain any information about the rate of the processes—about the kinetics of phase transitions. However, it is important to note that the kinetic parameters determine the course of many supercritical processes and can be used for their optimization.
When a porous body, such as a gel, is added to the binary system, the phase behavior and phase transition kinetics of the “2-propanol—CO2” system differ from the system without a gel. The gel has small-diameter pores in which diffusion is hindered since the pore diameter is commensurate with the mean free path of a molecule. Evaluation of the effect of a porous body on the rate of phase transitions is relevant for the supercritical drying process [29,30], for solvent exchange at high pressure [31], etc.
For the study of the kinetics of phase transitions and mass transfer in binary systems, various equations of state are used, such as PR-EoS [28,32], Soave-Redlich-Kwong EoS (SRK-EoS) [33,34,35], Valderrama-Patel-Teja (VPT-EoS) [36,37], and mixing rules for multicomponent systems—Van der Waals’ [28,38,39], Huron-Vidal’s [40,41], and Sandler-Wong’s [42,43]. Often, the Peng-Robinson equation with Van der Waals mixing rules is used to study systems in the sub- and supercritical state since this equation is not very complex and allows for calculations with sufficient accuracy.
In this paper, the kinetics of phase transitions of the “2-propanol—CO2” system in the presence of a model porous material—alginate gel—and without it at various process parameters were studied. A series of experiments were carried out at temperatures of 313 and 333 K and pressures of 6.30 and 7.80 MPa. The effect of an isobaric increase of temperature and an isothermal increase of pressure on the system was studied. Mass transfer coefficients were calculated using an approach based on the mass transfer equation in combination with the PR-EoS. The obtained results are necessary for the development of supercritical processes, such as supercritical drying [44], solvent exchange under pressure [31], and supercritical adsorption [45].

2. Materials and Methods

2.1. Chemicals

The following chemicals were used: Alginic acid sodium salt (AlgNa, Algin, Alginic acid sodium salt from brown algae) (CAS 9005-38-3), Sigma-Aldrich, Saint Louis, MO, USA; CaCl2·2H2O (CAS 10035-04-8), purity ≥ 99.5%, Sigma-Aldrich; 2-propanol, purity ≥ 99.5%, RusHim, Moscow, Russia; CO2, purity ≥ 99.8%, Itc-Pex, Moscow, Russia; deionized water (H2O), prepared by laboratory water treatment equipment (Arium 61215, Sartorius) Göttingen, Germany. All chemicals were used directly without further purification.

2.2. Preparation of Alginate Gels

In this work, it was decided to use an alginate-based gel as a porous material, which is formed during the production of organic aerogels [46,47]. The production of the gel is carried out by using the dripping method in accordance with [44]. To prepare a 1 wt% solution of AlgNa in water, the required amount of AlgNa powder is poured into water and mixed on a magnetic stirrer. The resulting mixture is continuously stirred for 24 h to completely dissolve the AlgNa. Spherical particles are formed by dripping the 1 wt% solution of AlgNa drop-by-drop through a needle into the crosslinking agent (5 wt% solution of CaCl2) with constant stirring. The resulting gel particles are kept in a CaCl2 solution for 24 h, so that the chemical reactions proceed in full.
To study the kinetics of phase transitions, it is necessary to carry out solvent exchange in the gels with 2-propanol. To do this, the gel samples are transferred from the crosslinking agent solution to a 10 wt% 2-propanol solution, then to 30 wt%, 50 wt%, 70 wt%, and 90 wt% and twice to 100 wt% solutions. The gels are kept in each solution for 2–3 h with stirring. The exchange process takes a total of 14–21 h.

2.3. Experimental Study of the Kinetics of Phase Transitions

The study of the kinetics of phase transitions was performed using an experimental setup for supercritical processes [48]. The schematic diagram of the setup is shown in Figure 2. In order to organize the supply of 2-propanol to the apparatus setup, it was upgraded [31]. An additional alcohol supply line was added to the setup, which consists of a dosing container 9, a membrane dosing pump 10 (Lewa, LDB1, Leonberg, Germany), a pipeline, and a shut-off valve. This line is connected to one of the upper pipes of the apparatus. The setup also has a Coriolis flow meter (Bronkhorst, mini Cori-flow M13, Ruurlo, The Netherlands) on the carbon dioxide supply line, which monitors the mass flow of the gas and determines its total mass. In addition, a video recording system is added to the setup to record the behavior of the system during phase transitions in the device with sight glasses 5.
The experimental study of the kinetics of phase transitions involves measuring the change in liquid phase volume when establishing equilibrium in the “2-propanol—CO2” system. A 250 mL high-pressure apparatus (5) is heated to a specific temperature using the temperature control system (6). If the kinetics are being studied in the presence of a porous body, alginate gel is loaded into the apparatus and then sealed. Next, a predetermined amount of 2-propanol is fed into the apparatus from container 9 using a diaphragm pump 10 at a flow rate of 50–100 mL/min. The apparatus is then thermostated for 20 min before starting the video recording system to detect the position of the phase boundary in the apparatus. The apparatus is then fed with CO2 from tank 1 with a flow rate of no more than 1000 g/h through condenser 2 and heating element 4, after which pressurization time does not exceed 2 min. The carbon dioxide flow rate is regulated using a metering valve, while the mass flow rate and total mass of CO2 are determined by the Coriolis mass flowmeter and control system, respectively. The video recording of the system continues for 30–90 min before draining the 2-propanol from the apparatus through separator 7 to collect the liquid phase. Finally, the gas phase flow rate is measured using rotameter 8.
This study investigates the “2-propanol—CO2” system under different conditions, including 313 K at 6.3 MPa and 7.8 MPa and 333 K at 6.3 MPa and 7.8 MPa.
The experiments are conducted in two stages, with repeated experiments at specified process parameters for each stage except for the “2-propanol—CO2” system in the presence of gel at parameters 333 K and 7.8 MPa.

2.4. Image Analysis

The process of the “2-propanol—CO2” system phase transitions was studied using a video recording system. The resulting video was converted into photographic images of the system with an interval of 2 s, which were used to determine the volume of the liquid phase at a certain point of time to obtain the kinetics curves of the phase transitions of the system. For this purpose, the height of the liquid phase h and the diameter of the apparatus d were determined for each image (Figure 3). The height of the liquid phase volume in the images and the errors of its measurements were received using the method presented in Appendix A. The error in determination of the height of the liquid phase volume for all repeated experiments does not exceed 5%.
To determine the actual height of the liquid phase H, we took into account the ratio of the scale of the apparatus in the image and in the real setup (the diameter of the apparatus in the image d and the real diameter of the apparatus D):
H = D · h d
By using Equation (1) to calculate the real height of the liquid phase level, we can determine the volume of the liquid phase V at a certain point in time:
V = H D 2 · H · D H 2 + D 2 4 · π 2 + a r c t g H D 2 H · D H 2 · L
Next, the change of the liquid phase volume in relation to the liquid phase initial volume ΔV/V0 is calculated:
Δ V V 0 = V V 0 V 0 · 100 %
Based on the obtained data on the dependence of the liquid phase volume change in time, we construct kinetics curves of the “2-propanol—CO2” system phase transitions under various conditions.

2.5. Estimation of Mass Transfer Coefficients

To determine the mass transfer coefficients across the phase interface, a mathematical description of the kinetics of the phase transitions is proposed. To do this, the change in the liquid phase volume due to mass transfer across the interface is calculated. At the initial time, the liquid phase is pure 2-propanol, and the gas phase is pure CO2. At the same time, it is assumed that the system is in isobaric and isothermal conditions. The total mass transfer through the phase boundary is calculated using the mass transfer equation:
d M = K · x * x · d F · d τ
where dM—the mass transfer rate, kg; K—the mass transfer coefficient, kg/(m2·s); (x* − x)—the driving force of the process; dF—the effective mass transfer area, m2; and —the time step, s.
The gas and liquid phase densities are calculated using the PR-EOS with Van der Waals mixing rules [28] for multicomponent systems. The equation can be represented as a polynomic for further solution:
Z 3 1 B · Z 2 + A 2 B 3 B 2 · Z A B B 2 B 3 = 0
where
A = a · P R 2 · T 2 ;   B = b · P R · T ;   Z = V m o l · P R · T
a = 0.4572 · R · T c 2 P c · 1 + m · 1 T T c 2 ; b = 0.0778 · R · T c P c
m = 0.3746 + 1.54226 · ω 0.2699 · ω 2
where P—pressure, MPa; R—gas constant, J/(mol·K); T—temperature, K; Vmol—molar volume, m3/mol; Tc—critical temperature, K; Pc—critical pressure, MPa; and ω—acentric factor.
The problem is then simplified by solving the cubic equation using Vieta’s trigonometric formula [49] to find the Z-coefficient of compressibility.
The equation is used to calculate the phase boundary F, through which mass transfer takes place:
F = L · C
where L—the length of the phase boundary, m; C—the width of the phase boundary, m.
Figure 4 shows a schematic representation of the length of the phase boundary, equal to the length of the apparatus, the width of the phase boundary, and the phase boundary surface.
The width of the boundary C is the chord of the circle of the apparatus, which is determined by geometric calculations in terms of the radius of the apparatus R and the height of the liquid level in the apparatus h:
C = 4 · R 2 R h 2
The height of the phase boundary is calculated based on the recalculation of the volumes of the liquid and gaseous phases. The process continues until the equilibrium mole fraction is established in both the liquid and gas phases. The calculated change of the liquid phase volume is determined by Equation (3) using the calculated volume of the liquid phase at a certain point of time.
The mass transfer coefficient and equilibrium mole fraction of CO2 in the liquid phase (2-propanol) are determined by minimizing the error between the calculated and experimental phase transition kinetics curves. The error is calculated using the following equations:
E r r o r   i = Δ V V 0 c a l c Δ V V 0 e x p Δ V V 0 e x p   · 100 %
E r r o r = E r r o r   i N
where Error (i)—error of the calculated data at a certain point of time, Error—average error of the calculated data for all experimental points, and N—the number of experimental points.
The minimization procedure was performed using the Broyden–Fletcher–Goldfarb–Shanno algorithm (BFGS). The algorithm was realized using the Python programming language and the SciPy library. A flowchart of the algorithm for determining the mass transfer coefficient and equilibrium mole fraction of CO2 in the liquid phase is presented in Section 3.2. Calculated curves of the phase transition kinetics.

3. Results and Discussion

3.1. Experimental Curves of Phase Transition Kinetics

Phase transition kinetics were studied at various external parameters: pressure 6.30 and 7.80 MPa and temperature 313 and 333 K. The specified parameters were selected to study the system in a heterogeneous region. It is important to note that it is complicated to reach target parameters, but they can be easily measured with high accuracy. Therefore, the values of the measured parameters for each experiment will be presented below. During the study of the “2-propanol—CO2” system, 50 mL of 2-propanol was preloaded into the apparatus. When studying the “2-propanol—CO2” system in the presence of gel, 80 mL of 2-propanol and 40 mL of gel were preloaded into the apparatus. The average particle diameter of the gel that is loaded into the apparatus was 2.7 mm with a standard deviation of 0.14.
After the research, the gels were dried using supercritical drying. The porosity of the obtained aerogels was 95%, the specific surface area was 442 m2/g, the average pore diameter was 22 nm, and the pore volume was 2.19 cm2/g.
Using the approach presented above, a set of phase transition kinetics curves was obtained. Table 1 shows the parameters of the two-component “2-propanol—CO2” system phase transition kinetics.
Table 2 shows the parameters of the two-component “2-propanol—CO2” system phase transition kinetics in the presence of alginate gel.
Experiment 4–2 was unsuccessful; therefore, we have excluded the data related to it from the article.
Furthermore, Figure 5 shows the apparatus at various points of time, in which the “2-propanol—CO2” system liquid phase volume changes in the absence of alginate gel at a temperature of 312.68 K and a pressure of 7.98 MPa (Table 1, experiment 2–2). Figure 6 shows a binary system in the presence of the gel at a temperature of 313.19 K and a pressure of 7.86 MPa (Table 2, experiment 2–2).
To obtain the phase transitions kinetics curves, at least 40 photographic images of the “2-propanol—CO2” system liquid phase were obtained in the apparatus at various points of time. As a result of the processing of photographic images (see Section 2.4), the change of the liquid phase volume was determined, and the kinetic curves of the “2-propanol—CO2” system phase transitions were obtained without and with the presence of the alginate gel (Figure 7 and Figure 8).
The kinetic curves of the phase transitions show that the change in the process parameters (temperature and pressure) significantly affects the change of the liquid phase volume (Figure 7). With an isobaric increase in temperature, the solubility of CO2 in the solvent decreases, which is confirmed by experimental studies. With an isothermal increase in pressure, the solubility of CO2 increases. It is also important to note that the values of the maximum change of the liquid phase volume (when the liquid phase volume stops changing and reaches a constant level) are not the same in different experiments with the same target parameters. The measured temperature and pressure of the experiments are different from each other due to the complexity of reaching the process parameters. The maximum volume is observed at approximately 7.80 MPa and a temperature of 313 K (Figure 7 and Figure 8b), since these parameters are close to the critical temperature and pressure of the binary mixture.
The “2-propanol—CO2” system without gel reaches the equilibrium state in less than 30 min after the start of the process (Figure 7), in contrast to the system in the presence of the alginate gel (Figure 8). This is due to the fact that the gel is a porous material; therefore, mass transfer slows down because of the diffusion inside pores.

3.2. Calculated Curves of the Phase Transition Kinetics

To study the effect of a porous body on the kinetics of the “2-propanol—CO2” system phase transitions, the mass transfer coefficients of the binary system and the equilibrium mole fraction of CO2 in the liquid phase (Section 2.5) were determined for various process parameters.
To implement this algorithm for calculating the mass transfer coefficient and the equilibrium mole fraction of CO2 in the liquid phase, the following data must be entered:
  • dVexp—a list consisting of all experimental data about changes in the liquid phase volume (see Section 2.4);
  • texp—a list consisting of the times at which the experimental data of the liquid phase volume change were determined;
  • Tcr,i/j, Pcr,i/j—critical temperature and pressure of substances, i—CO2, j—2-propanol;
  • T, P—experimental process temperature and pressure;
  • D, L—geometric dimensions of the apparatus, diameter and length, respectively;
  • V0—the liquid phase volume at the initial time is equal to the volume of pure 2-propanol;
  • x—the mole fraction of CO2, at the initial time is zero;
  • F—phase contact surface area;
  • mi/j—mass of the substance, i—CO2, j—2-propanol (const);
  • t—process time;
  • dt—time step; and
  • tend—the time of the end of the process, which is equal to the time of the experimental study of the kinetics of phase transitions.
The mass transfer coefficient K and the equilibrium mole fraction of CO2 in the liquid phase x* are variables and are determined using the BFGS algorithm by finding the minimum average error between the experimental and calculated data.
A flowchart for determining the mass transfer coefficient and the equilibrium mole fraction of CO2 in the liquid phase is shown in Figure 9.
Experimental and theoretical kinetic curves of the phase transitions of the system at different temperatures and pressures are shown in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17.
Table 3 shows the numerical values of the calculated data obtained for all the experimentally obtained process parameters, the equilibrium mole fraction from other researches (an interpolation procedure was used for data processing in some cases), and the average error of the calculated data of kinetic curves according to experimental ones.
Table 3 shows that the mass transfer coefficients and the equilibrium mole fraction of CO2 in the liquid phase can differ significantly since the binary system is sensitive to external parameters. The most obvious is the effect of the porous body on the rate of the process. A system with an alginate gel takes longer to reach equilibrium than a system without a gel. For the system at 313 K and 6.3 MPa, the mass transfer coefficient is reduced in 1.6 times; at 313 K and 7.8 MPa, in 5.0 times; at 333 K and 6.3 MPa, in 2.7 times; and at 333 K and 7.8 MPa, in 2.7 times. The mass transfer coefficient for the system without the gel decreases in 1.3 times with an isobaric increase in temperature by 20 degrees at 6.3 MPa and decreases in 3.1 times with an increase in temperature by 20 degrees at 7.8 MPa. The mass transfer coefficient for the system with the gel decreases in 2.2 times with an isobaric increase in temperature by 20 degrees at 6.3 MPa and decreases in 1.6 times with increase in temperature at 7.8 MPa. Moreover, the mass transfer coefficient increases for all cases with an isothermal pressure increase. It increases in 3.8 times for the system without the gel and in 1.2 times for the system with the gel at 313 K and increases in 1.6 times for both systems with and without the gel at 333 K. It can be noted that the effect of pressure is less significant in a system with the presence of a gel.
The specified behavior of the system can be directly related to the density of carbon dioxide. When changing external parameters, its density changes significantly. To demonstrate the influence of the density of carbon dioxide on the mass transfer coefficient, the data were plotted in Figure 18. The target external parameters are also marked in the picture. This applies for data without the gel and for data with the gel. It is important to note that the influence of the density of carbon dioxide is more interesting than the influence of the density of 2-propanol, as the density of 2-propanol changes insignificantly within the considered parameters.
The data clearly show that with an increase in the density, an increase in the coefficient of mass transfer occurs, that is, the rate of the process increases. This is due to the increased diffusion rate at higher densities of CO2. This is mainly true for systems without a gel. Moreover, it can be seen that the behavior of the system differs significantly at different pressures. As noted earlier, more obvious differences in the rate of the process in the system without the gel and in the system with the gel are noticeable at 7.8 MPa than at 6.3 MPa. The critical pressure of the mixture “2-propanol—CO2” at temperature 313 K is 8.29 MPa, and at temperature 333 K it is 10.35 MPa [23]. Thus, for both temperatures, the system at 7.8 MPa is closer to the critical parameters than the system at 6.3 MPa. It is obvious that under conditions near the critical parameters, a convective transport can arise. This is due to the fact that the densities of the gas and liquid phases become closer to each other, and their mutual mixing is facilitated. The presence of a porous body (alginate gel) significantly reduces the effect of convective transport. Thus, the rate of the mass transfer does not increase as noticeably with an increase in the density of CO2 and with the approach of the critical parameters for the system with alginate gel.
It can be seen that the closer the system parameters are to the critical parameters, the more difficult it becomes to maintain constant process parameters. Thus, there is a significant increase in the confidence intervals at the parameters of 313 K and 7.8 MPa.
The calculated equilibrium mole fractions of CO2 in the liquid phase are close to the literature data (Table 3). The average error of the calculated equilibrium mole fractions for all experiments is 13%. Moreover, the highest error is observed for experiments at the external target parameters 313 K and 7.8 MPa (25%), which are close to the critical parameters of the “2-propanol—CO2” system. This is due to the fact that the system in near critical conditions is less stable, and PR-EoS is worse for describing the system in a subcritical state.

4. Conclusions

Experimental studies of the kinetics of “2-propanol—CO2” system phase transitions under high pressure have been carried out. It was found that with an isobaric increase in temperature, the solubility of CO2 in 2-propanol decreases, and with an isothermal pressure increase, on the contrary, it increases. The system with alginate gel takes more than two times longer to reach a state of equilibrium. As a result of the experimental study, the kinetic curves of phase transitions were obtained using a new method with and without the presence of a porous body (alginate gel).
A new approach was proposed for calculating the mass transfer coefficient in the system under consideration. It was calculated using the Peng-Robinson equation of state with Van der Waals mixing rules. The proposed equations accurately describe the mass transfer process, and the error between the calculated and experimental data was, on average, 2%. The calculated values of the mass transfer coefficient and the analysis led to an important conclusion about the influence of external parameters on the transport rate in the heterogeneous system “2-propanol—CO2” with and without the presence of a porous body (alginate gel).
The data can be used to study the kinetic characteristics of many supercritical processes. The proposed approach can be used for optimization of processes that are carried out under conditions close to critical. This will help in the production of aerogels, and in some cases, it can be used in supercritical extraction, drying, etc.

Author Contributions

Conceptualization, A.L. and E.S.; methodology, A.L. and E.S.; software, A.L. and E.S.; validation, A.L., M.M. and E.S.; formal analysis, E.S.; investigation, E.S.; resources, M.M.; data curation, E.S.; writing—original draft preparation, A.L. and E.S.; writing—review and editing, E.S.; visualization, E.S.; supervision, A.L. and M.M.; project administration, E.S. and A.L.; funding acquisition, A.L. All authors have read and agreed to the published version of the manuscript.

Funding

Research was funded by the grant from Russian Science Foundation (Project No. 22-79-00154).

Data Availability Statement

The experimental data and code will be provided upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Supplementary Data (Liquid Phase Height and Measurement Error)

This section describes the determination of the height of the liquid phase h in the photographic image and its measurement errors. For a more accurate determination of the height, a grid of five parallel vertical lines (red lines in Figure A1) perpendicular to the diameter of the apparatus was constructed. An example of the grid is shown in Figure A1.
Figure A1. Auxiliary grid for determining the height of the liquid phase h.
Figure A1. Auxiliary grid for determining the height of the liquid phase h.
Computation 11 00122 g0a1
This grid was superimposed on the photographic images of the system in the apparatus in such a way that the phase boundary coincided with the surface shown in Figure A1. To determine the height h, the values of perpendiculars of each line that intersect with the circle that represents the wall of the apparatus (AiBi) and the distances from the center of the interface A0 to the vertical line Ai (A0Ai) are used. With the help of the obtained data, the height of the liquid phase h was recalculated according to the equation
h i = d 2 d 2 2 A 0 A i 2 + A i B i
Table A1 presents the data obtained from a single photographic image at the time τ = 20 min for the “2-propanol—CO2” system (experiment 4–2, Table 1).
Table A1. Auxiliary segment data used to determine the height of the liquid phase.
Table A1. Auxiliary segment data used to determine the height of the liquid phase.
Vertical Line №, iA0AiAiBihi
14.153.314.66
22.284.234.61
31.594.374.55
43.493.604.52
55.192.124.40
Next, the real height of the liquid phase H was calculated according to Equation (1) and its average value H ¯ at five points (n = 5). Then, the mean absolute deviation of the height values was determined using equation
σ = i = 1 n H i H ¯ n
After determining the height and the calculation error, a curve of the change of the height of the liquid phase with a vertical average deviation was constructed for each photographic image. An example of constructing the experiment curve 4–2 (Table 1) is shown in Figure A2.
Figure A2. Kinetics of the change of the “2-propanol—CO2” system liquid phase height (experiment 4—2, Table 1).
Figure A2. Kinetics of the change of the “2-propanol—CO2” system liquid phase height (experiment 4—2, Table 1).
Computation 11 00122 g0a2
When studying further the change of the liquid phase volume, the vertical average deviation is not indicated in the curves (Figure 5 and Figure 6), since they are within the limits shown in Figure A2, and the average deviation of the height of the liquid phase does not exceed 5%. The height of the liquid phase obtained in this way is used to further determine the liquid phase volume at each time using Equation (2).

References

  1. Hegel, P.; Martín, L.; Popovich, C.; Damiani, C.; Pancaldi, S.; Pereda, S.; Leonardi, P. Biodiesel Production from Neochloris Oleoabundans by Supercritical Technology. Chem. Eng. Process. Process Intensif. 2017, 121, 232–239. [Google Scholar] [CrossRef] [Green Version]
  2. Gracia, E.; García, M.T.; Rodríguez, J.F.; de Lucas, A.; Gracia, I. Improvement of PLGA Loading and Release of Curcumin by Supercritical Technology. J. Supercrit. Fluids 2018, 141, 60–67. [Google Scholar] [CrossRef]
  3. Arango-Ruiz, Á.; Martin, Á.; Cosero, M.J.; Jiménez, C.; Londoño, J. Encapsulation of Curcumin Using Supercritical Antisolvent (SAS) Technology to Improve Its Stability and Solubility in Water. Food Chem. 2018, 258, 156–163. [Google Scholar] [CrossRef] [PubMed]
  4. Johner, J.C.F.; Hatami, T.; Meireles, M.A.A. Developing a Supercritical Fluid Extraction Method Assisted by Cold Pressing for Extraction of Pequi (Caryocar Brasiliense). J. Supercrit. Fluids 2018, 137, 34–39. [Google Scholar] [CrossRef]
  5. Kohli, R. Applications of Supercritical Carbon Dioxide for Removal of Surface Contaminants. In Developments in Surface Contamination and Cleaning: Applications of Cleaning Techniques; Elsevier: Amsterdam, The Netherlands, 2018; Volume 11, pp. 209–249. ISBN 9780128155776. [Google Scholar]
  6. Nikolai, P.; Rabiyat, B.; Aliev, A.; Abdulagatov, I. Supercritical CO2: Properties and Technological Applications—A Review. J. Therm. Sci. 2019, 28, 394–430. [Google Scholar] [CrossRef]
  7. De Souza, A.T.; Benazzi, T.L.; Grings, M.B.; Cabral, V.; Antônio da Silva, E.; Cardozo-Filho, L.; Ceva Antunes, O.A. Supercritical Extraction Process and Phase Equilibrium of Candeia (Eremanthus Erythropappus) Oil Using Supercritical Carbon Dioxide. J. Supercrit. Fluids 2008, 47, 182–187. [Google Scholar] [CrossRef]
  8. Scopel, R.; Falcão, M.A.; Lucas, A.M.; Almeida, R.N.; Gandolfi, P.H.K.; Cassel, E.; Vargas, R.M.F. Supercritical Fluid Extraction from Syzygium aromaticum Buds: Phase Equilibrium, Mathematical Modeling and Antimicrobial Activity. J. Supercrit. Fluids 2014, 92, 223–230. [Google Scholar] [CrossRef]
  9. Sierra-Pallares, J.; del Valle, J.G.; Bueno, C.M.; García, J.; Castro, F. High-Pressure Fluid Flow in Carbon Dioxide—Organic Solvent Mixing Layers. Equilibrium Models and the Prediction of Phase Transitions. J. Supercrit. Fluids 2020, 169, 105024. [Google Scholar] [CrossRef]
  10. De Marco, I.; Knauer, O.; Cice, F.; Braeuer, A.; Reverchon, E. Interactions of Phase Equilibria, Jet Fluid Dynamics and Mass Transfer during Supercritical Antisolvent Micronization: The Influence of Solvents. Chem. Eng. J. 2012, 203, 71–80. [Google Scholar] [CrossRef]
  11. Smirnova, I.; Gurikov, P. Aerogel Production: Current Status, Research Directions, and Future Opportunities. J. Supercrit. Fluids 2018, 134, 228–233. [Google Scholar] [CrossRef]
  12. Lazrag, M.; Lemaitre, C.; Castel, C.; Hannachi, A.; Barth, D. Aerogel Production by Supercritical Drying of Organogels: Experimental Study and Modelling Investigation of Drying Kinetics. J. Supercrit. Fluids 2018, 140, 394–405. [Google Scholar] [CrossRef]
  13. Chen, D.; Gao, H.; Liu, P.; Huang, P.; Huang, X. Directly Ambient Pressure Dried Robust Bridged Silsesquioxane and Methylsiloxane Aerogels: Effects of Precursors and Solvents. RSC Adv. 2019, 9, 8664–8671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rao, A.P.; Rao, A.V.; Gurav, J.L. Effect of Protic Solvents on the Physical Properties of the Ambient Pressure Dried Hydrophobic Silica Aerogels Using Sodium Silicate Precursor. J. Porous Mater. 2008, 15, 507–512. [Google Scholar] [CrossRef]
  15. Tripathi, A.; Parsons, G.N.; Khan, S.A.; Rojas, O.J. Synthesis of Organic Aerogels with Tailorable Morphology and Strength by Controlled Solvent Swelling Following Hansen Solubility. Sci. Rep. 2018, 8, 2106. [Google Scholar] [CrossRef] [Green Version]
  16. Omranpour, H.; Motahari, S. Effects of Processing Conditions on Silica Aerogel during Aging: Role of Solvent, Time and Temperature. J. Non-Cryst. Solids 2013, 379, 7–11. [Google Scholar] [CrossRef]
  17. Bueno, A.; Selmer, I.; SP, R.; Gurikov, P.; Lölsberg, W.; Weinrich, D.; Fricke, M.; Smirnova, I. First Evidence of Solvent Spillage under Subcritical Conditions in Aerogel Production. Ind. Eng. Chem. Res. 2018, 57, 8698–8707. [Google Scholar] [CrossRef]
  18. Gui, X.; Tang, Z.; Fei, W. Solubility of CO2 in Alcohols, Glycols, Ethers, and Ketones at High Pressures from (288.15 to 318.15) K. J. Chem. Eng. Data 2011, 56, 2420–2429. [Google Scholar] [CrossRef]
  19. Décultot, M.; Ledoux, A.; Fournier-Salaün, M.C.; Estel, L. Solubility of CO2 in Methanol, Ethanol, 1,2-Propanediol and Glycerol from 283.15 K to 373.15 K and up to 6.0 MPa. J. Chem. Thermodyn. 2019, 138, 67–77. [Google Scholar] [CrossRef]
  20. Liu, H.; Gao, H.; Idem, R.; Tontiwachwuthikul, P.; Liang, Z. Analysis of CO2 Solubility and Absorption Heat into 1-Dimethylamino-2-Propanol Solution. Chem. Eng. Sci. 2017, 170, 3–15. [Google Scholar] [CrossRef]
  21. Lazzaroni, M.J.; Bush, D.; Brown, J.S.; Eckert, C.A. High-Pressure Vapor−Liquid Equilbria of Some Carbon Dioxide + Organic Binary Systems. J. Chem. Eng. Data 2005, 50, 60–65. [Google Scholar] [CrossRef]
  22. Bamberger, A.; Maurer, G. High-Pressure (Vapour + Liquid) Equilibria in (Carbon Dioxide + Acetone or 2-Propanol) at Temperatures from 293 K to 333 K. J. Chem. Thermodyn. 2000, 32, 685–700. [Google Scholar] [CrossRef]
  23. Radosz, M. Vapor-Liquid Equilibrium for 2-Propanol and Carbon Dioxide. J. Chem. Eng. Data 1986, 31, 43–45. [Google Scholar] [CrossRef]
  24. Ziegler, J.W.; Chester, T.L.; Innis, D.P.; Page, S.H.; Dorsey, J.G. Supercritical Fluid Flow Injection Method for Mapping Liquid—Vapor Critical Loci of Binary Mixtures Containing CO2. In Innovations in Supercritical Fluids; ACS Symposium Series; American Chemical Society: Washington, DC, USA, 1995; Volume 608, pp. 93–110. ISBN 978-0-8412-3324-9. [Google Scholar]
  25. Elizalde-Solis, O.; Galicia-Luna, L.A.; Camacho-Camacho, L.E. High-Pressure Vapor–Liquid Equilibria for CO2+alkanol Systems and Densities of n-Dodecane and n-Tridecane. Fluid Phase Equilibria 2007, 259, 23–32. [Google Scholar] [CrossRef]
  26. Tachiwaki, T. Supercritical Fluid Drying of Aqueous 2-Propanol Suspension. Dry. Technol. 2004, 22, 325–334. [Google Scholar] [CrossRef]
  27. Egerer, K.; Brandl, F.; Golzari, N.; Beuermann, S. Recycling of Printed Circuit Boards Employing Supercritical Carbon Dioxide. Mater. Sci. Forum 2019, 959, 100–106. [Google Scholar] [CrossRef]
  28. Fateen, S.E.K.; Khalil, M.M.; Elnabawy, A.O. Semi-Empirical Correlation for Binary Interaction Parameters of the Peng-Robinson Equation of State with the van Der Waals Mixing Rules for the Prediction of High-Pressure Vapor-Liquid Equilibrium. J. Adv. Res. 2013, 4, 137–145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Muhr, A.H.; Blanshard, J.M.V. Diffusion in Gels. Polymer 1982, 23, 1012–1026. [Google Scholar] [CrossRef]
  30. Masmoudi, Y.; Rigacci, A.; Ilbizian, P.; Cauneau, F.; Achard, P. Diffusion During the Supercritical Drying of Silica Gels. Dry. Technol. 2006, 24, 1121–1125. [Google Scholar] [CrossRef]
  31. Lebedev, A.; Suslova, E.; Troyankin, A.; Lovskaya, D. Investigation of Aerogel Production Processes: Solvent Exchange under High Pressure Combined with Supercritical Drying in One Apparatus. Gels 2021, 7, 4. [Google Scholar] [CrossRef]
  32. Peng, D.-Y.; Robinson, D.B. A New Two-Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59–64. [Google Scholar] [CrossRef]
  33. Markočič, E.; Knez, Ž. Redlich-Kwong Equation of State for Modelling the Solubility of Methane in Water over a Wide Range of Pressures and Temperatures. Fluid Phase Equilibria 2016, 408, 108–114. [Google Scholar] [CrossRef]
  34. Su, C.S. Prediction of Volumetric Properties of Carbon Dioxide-Expanded Organic Solvents Using the Predictive Soave-Redlich-Kwong (PSRK) Equation of State. J. Supercrit. Fluids 2012, 72, 223–231. [Google Scholar] [CrossRef]
  35. Gregorowicz, J.; De Loos, T.W. Modelling of the Three Phase LLV Region for Ternary Hydrocarbon Mixtures with the Soave-Redlich-Kwong Equation of State. Fluid Phase Equilibria 1996, 118, 121–132. [Google Scholar] [CrossRef]
  36. Forero, G.L.A.; Velásquez, J.J. The Patel-Teja and the Peng-Robinson EoSs Performance When Soave Alpha Function Is Replaced by an Exponential Function. Fluid Phase Equilibria 2012, 332, 55–76. [Google Scholar] [CrossRef]
  37. Valderrama, J.O.; Cardona, L.F.; Rojas, R.E. Correlation of Ionic Liquid Viscosity Using Valderrama-Patel-Teja Cubic Equation of State and the Geometric Similitude Concept. Part II: Binary Mixtures of Ionic Liquids. Fluid Phase Equilibria 2019, 497, 178–194. [Google Scholar] [CrossRef]
  38. Barragán-Aroche, J.F.; Bazúa-Rueda, E.R. A Local Composition Extension of the van Der Waals Mixing Rule for PR Cubic Equation of State. Fluid Phase Equilibria 2005, 227, 97–112. [Google Scholar] [CrossRef]
  39. Grgurić, I.R.; Šerbanović, S.P.; Kijevčanin, M.L.; Tasić, A.Ž.; Djordjević, B.D. Volumetric Properties of the Ternary System Ethanol + 2-Butanone + Benzene by the van Der Waals and Twu-Coon-Bluck-Tilton Mixing Rules: Experimental Data, Correlation and Prediction. Thermochim. Acta 2004, 412, 25–31. [Google Scholar] [CrossRef]
  40. Forero, G.; Luis, A.; Velasquez, J.; Jorge, A. Modeling the Liquid-Liquid Equilibria of Polar Aprotic Solvents/Alkanes Type Mixtures Using a Modified PR EoS and the Huron-Vidal Mixing Rules. J. Mol. Liq. 2019, 292, 111380. [Google Scholar] [CrossRef]
  41. Velásquez, J.; Jorge, A.; Forero, G.; Luis, A. Calculation of Complex Phase Equilibria of DMF/Alkane and Acetonitrile/Alkane Systems Using a Modified Peng-Robinson EoS and the Huron-Vidal Mixing Rules. J. Mol. Liq. 2017, 243, 600–610. [Google Scholar] [CrossRef]
  42. López, J.A.; Cardona, C.A. Phase Equilibrium Calculations for Carbon Dioxide + N-Alkanes Binary Mixtures with the Wong-Sandler Mixing Rules. Fluid Phase Equilibria 2006, 239, 206–212. [Google Scholar] [CrossRef]
  43. Espanani, R.; Miller, A.; Jacoby, W. Prediction of Vapor-Liquid Equilibria for Mixtures of Low Boiling Point Compounds Using Wong-Sandler Mixing Rule and EOS/GE Model. Chem. Eng. Sci. 2016, 152, 343–350. [Google Scholar] [CrossRef]
  44. Menshutina, N.V.; Lovskaya, D.D.; Lebedev, A.E.; Lebedev, E.A. Production of Sodium Alginate-Based Aerogel Particles Using Supercritical Drying in Units with Different Volumes. Russ. J. Phys. Chem. B 2017, 11, 1296–1305. [Google Scholar] [CrossRef]
  45. Bolten, D.; Johannsen, M. Influence of 2-Propanol on Adsorption Equilibria of α- and δ-Tocopherol from Supercritical Carbon Dioxide on Silica Gel. J. Chem. Eng. Data 2006, 51, 2132–2137. [Google Scholar] [CrossRef]
  46. Lovskaya, D.; Menshutina, N. Alginate-Based Aerogel Particles as Drug Delivery Systems: Investigation of the Supercritical Adsorption and In Vitro Evaluations. Materials 2020, 13, 329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Gaudio, P.D.; Auriemma, G.; Mencherini, T.; Porta, G.D.; Reverchon, E.; Aquino, R.P. Design of Alginate-Based Aerogel for Nonsteroidal Anti-Inflammatory Drugs Controlled Delivery Systems Using Prilling and Supercritical-Assisted Drying. J. Pharm. Sci. 2013, 102, 185–194. [Google Scholar] [CrossRef]
  48. Lebedev, A.E.; Katalevich, A.M.; Menshutina, N.V. Modeling and Scale-up of Supercritical Fluid Processes. Part I: Supercritical Drying. J. Supercrit. Fluids 2015, 106, 122–132. [Google Scholar] [CrossRef]
  49. Lipschutz, S.; Liu, J.; Spiegel, T. Schaum’s Outline: Mathematical Handbook of Formulas and Tables, 4th ed.; McGraw-Hill Education: Sydney, Australia, 2018; ISBN 978-1-260-01053-4. [Google Scholar]
  50. Zevnik, L.; Levec, J. Hydrogen Solubility in CO2-Expanded 2-Propanol and in Propane-Expanded 2-Propanol Determined by an Acoustic Sensor. J. Supercrit. Fluids 2007, 41, 335–342. [Google Scholar] [CrossRef]
  51. Lim, J.S.; Jung, Y.G.; Yoo, K.-P. High-Pressure Vapor−Liquid Equilibria for the Binary Mixtures of Carbon Dioxide + Isopropanol (IPA). J. Chem. Eng. Data 2007, 52, 2405–2408. [Google Scholar] [CrossRef]
  52. Secuianu, C.; Feroiu, V.; Geanǎ, D. High-Pressure Vapor-Liquid Equilibria in the System Carbon Dioxide and 2-Propanol at Temperatures from 293.25 K to 323.15 K. J. Chem. Eng. Data 2003, 48, 1384–1386. [Google Scholar] [CrossRef]
  53. Khalil, W.; Coquelet, C.; Richon, D. High-Pressure Vapor−Liquid Equilibria, Liquid Densities, and Excess Molar Volumes for the Carbon Dioxide + 2-Propanol System from (308.10 to 348.00) K. J. Chem. Eng. Data 2007, 52, 2032–2040. [Google Scholar] [CrossRef]
  54. Sima, S.; Racovita, R.; Dinca, C.; Feroiu, V.; Secuianu, C. Phase Equilibria Calculations for Carbon Dioxide + 2-Propanol System. UPB Sci. Bull. Ser. B Chem. Mater. Sci. 2017, 79, 11–24. [Google Scholar]
  55. Yaginuma, R.; Nakajima, T.; Tanaka, H.; Kato, M. Densities of Carbon Dioxide + 2-Propanol at 313.15 K and Pressures to 9.8 MPa. Available online: https://pubs.acs.org/doi/pdf/10.1021/je9700028 (accessed on 2 June 2021).
  56. Hsieh, C.-M.; Windmann, T.; Vrabec, J. Vapor-Liquid Equilibria of CO2 + C1-C5 Alcohols from the Experiment and the COSMO-SAC Model. J. Chem. Eng. Data 2013, 58, 3420–3429. [Google Scholar] [CrossRef]
Figure 1. Liquid phase volume of the “2-propanol—CO2” system with increasing pressure.
Figure 1. Liquid phase volume of the “2-propanol—CO2” system with increasing pressure.
Computation 11 00122 g001
Figure 2. Schematic diagram of the equipment for conducting supercritical processes (PI—pressure probe, TC—temperature probe, TI—temperature probe).
Figure 2. Schematic diagram of the equipment for conducting supercritical processes (PI—pressure probe, TC—temperature probe, TI—temperature probe).
Computation 11 00122 g002
Figure 3. Processed image of the system for determining the volume of the liquid phase in the apparatus.
Figure 3. Processed image of the system for determining the volume of the liquid phase in the apparatus.
Computation 11 00122 g003
Figure 4. Schematic representation of the system in the apparatus.
Figure 4. Schematic representation of the system in the apparatus.
Computation 11 00122 g004
Figure 5. Kinetics of the change in the “2-propanol—CO2” system liquid phase volume without the addition of the alginate gel at a temperature of 312.68 K and a pressure of 7.98 MPa.
Figure 5. Kinetics of the change in the “2-propanol—CO2” system liquid phase volume without the addition of the alginate gel at a temperature of 312.68 K and a pressure of 7.98 MPa.
Computation 11 00122 g005
Figure 6. Kinetics of the change in the “2-propanol—CO2” system liquid phase volume without the addition of the alginate gel at a temperature of 313.19 K and a pressure of 7.86 MPa.
Figure 6. Kinetics of the change in the “2-propanol—CO2” system liquid phase volume without the addition of the alginate gel at a temperature of 313.19 K and a pressure of 7.86 MPa.
Computation 11 00122 g006
Figure 7. Phase transitions kinetics curves of a two-component “2-propanol—CO2” system with target parameters: (a) 313 K, 6.30 MPa; (b) 313 K, 7.80 MPa; (c) 333 K, 6.30 MPa, (d) 333 K, 7.80 MPa.
Figure 7. Phase transitions kinetics curves of a two-component “2-propanol—CO2” system with target parameters: (a) 313 K, 6.30 MPa; (b) 313 K, 7.80 MPa; (c) 333 K, 6.30 MPa, (d) 333 K, 7.80 MPa.
Computation 11 00122 g007aComputation 11 00122 g007b
Figure 8. Phase transitions kinetics curves of the two-component “2-propanol—CO2” system in the presence of the alginate gel at the target parameters: (a) 313 K, 6.30 MPa; (b) 313 K, 7.80 MPa; (c) 333 K, 6.30 MPa, (d) 333 K, 7.80 MPa.
Figure 8. Phase transitions kinetics curves of the two-component “2-propanol—CO2” system in the presence of the alginate gel at the target parameters: (a) 313 K, 6.30 MPa; (b) 313 K, 7.80 MPa; (c) 333 K, 6.30 MPa, (d) 333 K, 7.80 MPa.
Computation 11 00122 g008
Figure 9. Flowchart for determining the mass transfer coefficient K and the equilibrium mole fraction of CO2 in the liquid phase x.
Figure 9. Flowchart for determining the mass transfer coefficient K and the equilibrium mole fraction of CO2 in the liquid phase x.
Computation 11 00122 g009
Figure 10. Experimental and calculated data of the kinetics of phase transitions at target parameters of 313 K, 6.30 MPa: experiments conducted without an alginate gel (see Table 1).
Figure 10. Experimental and calculated data of the kinetics of phase transitions at target parameters of 313 K, 6.30 MPa: experiments conducted without an alginate gel (see Table 1).
Computation 11 00122 g010
Figure 11. Experimental and calculated data of the kinetics of phase transitions at target parameters of 313 K, 6.30 MPa: experiments conducted with an alginate gel (see Table 2).
Figure 11. Experimental and calculated data of the kinetics of phase transitions at target parameters of 313 K, 6.30 MPa: experiments conducted with an alginate gel (see Table 2).
Computation 11 00122 g011
Figure 12. Experimental and calculated data of the kinetics of phase transitions at target parameters of 313 K, 7.80 MPa: experiments conducted without an alginate gel (see Table 1).
Figure 12. Experimental and calculated data of the kinetics of phase transitions at target parameters of 313 K, 7.80 MPa: experiments conducted without an alginate gel (see Table 1).
Computation 11 00122 g012
Figure 13. Experimental and calculated data of the kinetics of phase transitions at target parameters of 313 K, 7.80 MPa: experiments conducted with an alginate gel (see Table 2).
Figure 13. Experimental and calculated data of the kinetics of phase transitions at target parameters of 313 K, 7.80 MPa: experiments conducted with an alginate gel (see Table 2).
Computation 11 00122 g013
Figure 14. Experimental and calculated data of the kinetics of phase transitions at target parameters of 333 K, 6.30 MPa: experiments conducted without an alginate gel (see Table 1).
Figure 14. Experimental and calculated data of the kinetics of phase transitions at target parameters of 333 K, 6.30 MPa: experiments conducted without an alginate gel (see Table 1).
Computation 11 00122 g014
Figure 15. Experimental and calculated data of the kinetics of phase transitions at target parameters of 333 K, 6.30 MPa: experiments conducted with an alginate gel (see Table 2).
Figure 15. Experimental and calculated data of the kinetics of phase transitions at target parameters of 333 K, 6.30 MPa: experiments conducted with an alginate gel (see Table 2).
Computation 11 00122 g015
Figure 16. Experimental and calculated data of the kinetics of phase transitions at target parameters of 333 K, 7.80 MPa: experiments conducted without an alginate gel (see Table 1).
Figure 16. Experimental and calculated data of the kinetics of phase transitions at target parameters of 333 K, 7.80 MPa: experiments conducted without an alginate gel (see Table 1).
Computation 11 00122 g016
Figure 17. Experimental and calculated data of the kinetics of phase transitions at target parameters of 333 K, 7.80 MPa: experiments conducted with an alginate gel (see Table 2).
Figure 17. Experimental and calculated data of the kinetics of phase transitions at target parameters of 333 K, 7.80 MPa: experiments conducted with an alginate gel (see Table 2).
Computation 11 00122 g017
Figure 18. Mass transfer coefficient at different densities of CO2. (ο—without gel, □—with gel).
Figure 18. Mass transfer coefficient at different densities of CO2. (ο—without gel, □—with gel).
Computation 11 00122 g018
Table 1. Parameters of the two-component “2-propanol—CO2” system kinetics of the phase transitions without gel addition.
Table 1. Parameters of the two-component “2-propanol—CO2” system kinetics of the phase transitions without gel addition.
ExperimentExternal Target ParametersMeasured External Parameters
Temperature, [K]Pressure, [MPa]Temperature, [K]Pressure, [MPa]
1–13136.30313.836.13
1–2312.996.34
1–3314.356.48
2–13137.80312.667.76
2–2312.687.98
2–3314.987.98
3–13336.30332.956.50
3–2333.386.34
3–3325.286.67
4–13337.80333.057.66
4–2333.837.96
4–3333.257.64
Table 2. Parameters of the two-component “2-propanol—CO2” system kinetics of the phase transitions in the presence of the gel.
Table 2. Parameters of the two-component “2-propanol—CO2” system kinetics of the phase transitions in the presence of the gel.
ExperimentExternal Target ParametersMeasured External Parameters
Temperature, [K]Pressure, [MPa]Temperature, [K]Pressure, [MPa]
1–13136.30313.936.36
1–2314.036.29
1–3312.806.29
2–13137.80314.487.75
2–2313.197.86
2–3314.017.78
3–13336.30333.346.32
3–2333.796.28
3–3333.526.28
4–13337.80333.727.80
4–3332.267.79
Table 3. Calculated data of the kinetics of phase transitions.
Table 3. Calculated data of the kinetics of phase transitions.
CO2 Density, kg/m3ExperimentAverage Error of the Calculated Data of Kinetic Curves (Error), %Calculated Mass Transfer Coefficient,
kg/(m2·s)
Average Value of the Mass Transfer Coefficient, kg/(m2·s)Calculated Equilibrium Mole Fraction of CO2 in the Liquid Phase, mol/molEquilibrium Mole Fraction of CO2 in the Liquid Phase from [22,23,24,50,51,52,53,54,55,56], mol/mol
Without gel
T = 313 K,
P = 6.3 MPa
1641–11.000.05770.06550.440.47
1–21.920.04460.560.52
1–31.020.08640.590.52
With gel
T = 313 K,
P = 6.3 MPa
1641–11.200.05000.03960.590.51
1–21.400.02270.560.50
1–31.650.04600.620.52
Without gel
T = 313 K,
P = 7.8 MPa
2672–10.920.10320.23610.730.88
2–23.460.28140.790.95
2–32.770.32370.810.86
With gel
T = 313 K,
P = 7.8 MPa
2672–11.270.06730.04690.650.83
2–23.230.03950.640.91
2–31.940.03380.640.85
Without gel
T = 333 K,
P = 6.3 MPa
1343–13.190.03900.04840.360.38
3–23.720.03860.400.37
3–30.790.06760.400.41
With gel
T = 333 K,
P = 6.3 MPa
1343–11.470.02240.01820.410.37
3–20.920.01640.390.37
3–32.030.01570.440.37
Without gel
T = 333 K,
P = 7.8 MPa
1844–11.560.07420.07730.450.48
4–21.180.08610.550.52
4–31.760.07160.520.48
With gel
T = 333 K,
P = 7.8 MPa
1844–12.380.03180.02880.420.51
4–32.700.02570.480.50
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Suslova, E.; Mochalova, M.; Lebedev, A. Experimental and Theoretical Investigation of Supercritical Processes: Kinetics of Phase Transitions in Binary “2-Propanol—CO2” System. Computation 2023, 11, 122. https://doi.org/10.3390/computation11070122

AMA Style

Suslova E, Mochalova M, Lebedev A. Experimental and Theoretical Investigation of Supercritical Processes: Kinetics of Phase Transitions in Binary “2-Propanol—CO2” System. Computation. 2023; 11(7):122. https://doi.org/10.3390/computation11070122

Chicago/Turabian Style

Suslova, Ekaterina, Maria Mochalova, and Artem Lebedev. 2023. "Experimental and Theoretical Investigation of Supercritical Processes: Kinetics of Phase Transitions in Binary “2-Propanol—CO2” System" Computation 11, no. 7: 122. https://doi.org/10.3390/computation11070122

APA Style

Suslova, E., Mochalova, M., & Lebedev, A. (2023). Experimental and Theoretical Investigation of Supercritical Processes: Kinetics of Phase Transitions in Binary “2-Propanol—CO2” System. Computation, 11(7), 122. https://doi.org/10.3390/computation11070122

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop