Effect of Pressure on the Selectivity of Supercritical CO₂ Extraction During the Fractionation of a Fatty Acid Ethyl Ester Mixture: Numerical Simulation and Experiment

1. Introduction

One of the main sources of atmospheric pollution in large cities is motor transport running on fossil fuels. The exhaust gases from such engines contain components hazardous to human health: carbon monoxide (CO), which disrupts tissue respiration; nitrogen oxides (NOx) and sulfur oxides (SOx), causing respiratory diseases; as well as carcinogenic compounds such as polycyclic aromatic hydrocarbons (e.g., benzo(a)pyrene) and soot [1,2,3,4,5]. In light of increasingly stringent environmental regulations and global climate goals, the transition to alternative, cleaner fuels has become highly relevant.

Biodiesel fuel – a mixture primarily composed of methyl or fatty acid ethyl esters (FAEE), obtained mainly through the transesterification of vegetable oils (rapeseed, soybean, palm) or animal fats with lower alcohols (methanol, ethanol) – represents a key renewable energy source [6,7,8,9,10]. Unlike petroleum-based diesel, biodiesel offers several significant environmental advantages: its combustion leads to a substantial reduction in CO emissions, unburnt hydrocarbons, soot, and almost completely eliminates SOx emissions [11,12,13]. Furthermore, biodiesel is biodegradable and exhibits excellent lubricating properties. The growing demand for environmentally friendly fuel is confirmed by market dynamics: according to a report by the International Energy Agency, published on January 29, 2026, at the India Energy Week exhibition, biofuel production could reach 500-700 million liters per year by 2030 under the main and accelerated scenarios, respectively [14]. Despite these environmental benefits, the widespread commercial adoption of biodiesel is hindered by two groups of problems: technological (high production costs, complexity of removing catalysts and glycerol) and qualitative (instability of properties depending on the feedstock). The latter includes the issue of a high content of high-viscosity saturated fatty acid esters, such as palmitate and stearate, which is characteristic, for example, of biodiesel derived from shea butter or palm oil [1,6,15]. These saturated esters possess a high cloud point and pour point, as well as increased viscosity, unlike esters of unsaturated fatty acids (ethyl oleate, etc.). This directly deteriorates the low-temperature properties of the fuel (cold flow properties), promotes deposit formation in the fuel system, and can lead to non-compliance with stringent quality standards such as EN 14214 (ν=3.5–5.0 mm²/s) or ASTM D6751 (ν=1.9–6.0 mm²/s) [16,17].

The ever-increasing demand for biodiesel necessitates the improvement of technologies to reduce its cost and enhance its quality. A promising direction is the intensification of the transesterification reaction itself, including the use of supercritical fluid media (e.g., supercritical ethanol), which simplifies the process, minimizes catalyst usage, and improves kinetics [18,19,20,21,22]. However, an equally important stage is the subsequent purification and fractionation of the resulting ester mixture to isolate the target high-quality fraction. In this context, supercritical fluid extraction (SFE) based on carbon dioxide (scCO₂) stands out as an environmentally friendly, energy-efficient, and highly selective purification technology that does not require organic solvents [23,24,25,26]. The process relies on the unique properties of scCO₂, which in its supercritical state (Tc = 31.1 °C, Pc = 7.38 MPa) combines liquid-like high density with low viscosity and gas-like high diffusivity [27,28,29]. These properties can be tuned by varying pressure and temperature, making scCO₂ an ideal agent for the selective separation of multicomponent mixtures, such as biodiesel, for the purpose of removing high-viscosity saturated fractions [30,31,32,33,34,35].

Despite its obvious advantages, the design and optimization of industrial SFE plants for biodiesel fractionation remains a challenging task. This is related to a fundamental problem at the intersection of fluid physics and heat and mass transfer: an insufficient understanding of the interrelationship between the hydrodynamics of the supercritical flow, heat transfer, and the kinetics of selective mass transfer within an extraction column.

As noted in studies devoted to SFE for extracting biologically active compounds, local flow parameters, such as turbulence intensity (TI), have a decisive influence on the wall temperature distribution, density stratification, and, consequently, the efficiency of mass transfer [36,37,38]. In the context of biodiesel fractionation, this means that flow hydrodynamics (flow structure, mixing) directly determines the separation selectivity of saturated and unsaturated fatty acid esters and, ultimately, the viscosity and low-temperature properties of the final product. In the modern literature, considerable attention is paid to the modeling of SFE processes; however, the emphasis is shifted either towards purely hydrodynamic analysis of apparatus (separators, heat exchangers) [39,40] or the study of processes involving a solid phase (extraction from solid raw materials, aerogel drying) [41,42,43,44]. Studies specifically addressing liquid-liquid extraction or fractionation of ester mixtures in scCO₂ are often limited to phenomenological models or a simplified approach that does not account for the three-dimensional flow structure and its impact on local selectivity [45,46,47]. The selection of an adequate multicomponent and multiphase model in computational fluid dynamics (CFD) is particularly relevant. As shown in comparative studies, approaches such as the Eulerian Multifluid VOF (EMVOF) model can correctly describe complex phase interactions [48,49,50,51,52,53,54,55,56,57,58]. However, their validation and application specifically to the task of selective biodiesel fractionation in scCO₂ with product quality prediction remain practically unexplored in the literature.

Consequently, there exists a scientific and technical gap in the field of fundamental CFD modeling that quantitatively links the hydrodynamics of supercritical scCO₂ flow with the kinetics of selective mass transfer and the final physicochemical properties of the fractionated fatty acid ester mixture.

The aim of the present work is the development and experimental validation of a comprehensive CFD model for the fractionation process of a mixture of the most common fatty acid esters (FAEE) found in biodiesel, namely ethyl oleate and ethyl palmitate, in supercritical carbon dioxide. The model is aimed at reducing product viscosity through the selective removal of high-viscosity saturated esters. The investigation seeks to deepen the fundamental understanding of the "hydrodynamics-mass transfer-product properties" relationship in supercritical systems and will provide a practical tool for designing energy-efficient and selective technologies for renewable fuel purification

2. Materials and Methods

For the experimental implementation of the fractionation of an ethyl oleate and ethyl palmitate mixture using the scCO₂ extraction process, the following reagents were used:

– Carbon dioxide with a purity of not less than 99.0%, LLC "Snabtekhmet" (Kazan, Republic of Tatarstan, Russia);

– Ethyl oleate, pure, Sisco Research Laboratories Pvt. Ltd (Mumbai, Maharashtra, India) ($n_D^{20}$=1.450, ρ²⁰=870 кг/м³);

– Ethyl palmitate with a main substance content of not less than 97.0%, LLC TD "KHIMMED" (Moscow, Russia).

The SFE process was implemented on an experimental setup, the schematic diagram and visual appearance of which are presented in Figure 1. The methodology and description of the experiment are detailed in the authors' previous work [59].

Analysis of the obtained extract samples for the FAEE mixtures was performed using a "Chromatek-Crystal 5000" gas chromatograph (CJSC SKB "Khromatek", Yoshkar-Ola, Republic of Mari El, Russia). Samples were preliminarily dissolved in isopropanol at a 1:1 mass ratio. Heating was carried out from 40 to 210 °C, increasing the temperature at a rate of 10 °C/min. The injected sample volume was 0.5 µL.

When conducting CFD modeling in Ansys Fluent, the object of the study was the laboratory SFE setup. Figure 2 shows the extractor drawing (a), based on which the geometry was constructed with the indication of the experimental boundary conditions (b).

The computational mesh for the extractor model was developed considering the geometric features of the extractor and the requirements for modeling multiphase flows under supercritical conditions, based on the conducted literature review [39,40,41,42,43,44,48,49,50,51,52,53,54,55,56,57,58]. It consists of 418,639 nodes and 135,929 cells. The input and dimensional data, as well as the boundary conditions and assumptions, are provided in Table 1.

AISI 321 stainless steel was specified as the extractor material, and borosilicate glass was specified as the material for the packing. The scCO₂ supply flow rate was 2 mL/min. The mixture of ethyl oleate and ethyl palmitate (25.28/74.72 wt.%) was used as the target component. The choice of a 2 mL/min flow rate and this volume ratio was based on previous experimental extraction studies. The selection of pressures of 11 and 14 MPa is based on the phase diagrams of the binary mixtures CO₂-ethyl palmitate and CO₂-ethyl oleate (Figure 3), where at these pressures, type I-II phase behavior according to the Williams classification [60] is formed, with unbroken critical curves and the presence of critical points at these pressures and t = 40 °C.

The following main equations and models were used in the CFD modeling:

– Continuity equation (mass balance) for each phase:

$$S_{i }={\displaystyle \frac{\partial \left({\alpha }_{i }{\rho }_i\right)}{\partial t}}+\nabla ·\left({\alpha }_{i }{\rho }_i{u}_i\right)$$

(1)

where: αᵢ – volume fraction of phase i, ρᵢ – density of phase i; uᵢ – velocity of phase i; Sᵢ – mass transfer source term.

– Navier-Stokes equation (momentum balance) for the multiphase system.

$${\displaystyle \frac{\partial \left({\alpha }_{i }{\rho }_i{u}_i\right)}{\partial t}}+\nabla ·\left({\alpha }_{i }{\rho }_i{u}_i\right)=-{\alpha }_{i }\nabla P+\nabla {\tau }_{i }+{\alpha }_{i }{\rho }_{i}g+F_{drag }$$

(2)

where: P – pressure; τᵢ – viscous stress tensor; F_drag – interphase drag force;

– Species transport equation (concentration of the FAEE mixture in scCO₂).

$${\displaystyle \frac{\partial \left({\alpha }_{{co}_2 }{\rho }_{{co}_2}{Y}_{FAEE}\right)}{\partial t}}+\nabla ·\left({\alpha }_{{co}_2 }{\rho }_{{co}_2}{u}_{{co}_2}{Y}_{FAEE}\right)=\nabla ·\left(D_{Eff.}\nabla Y_{FAEE}\right)+R_{Ext.}$$

(3)

where: Y_FAEE – mass fraction of the FAEE mixture, D_Eff – effective diffusion coefficient, R_Ext – extraction rate.

Furthermore, interphase interaction models were considered, such as the drag force between scCO₂ and the FAEE mixture:

– Schiller-Naumann model:

$$F={\displaystyle \frac{3}{4}}{C}_D{\displaystyle \frac{{\alpha }_{{co}_2 }{\alpha }_{FAEE}{\rho }_{{co}_2}}{d_p}}·\left|u_{{co}_2} -u_{FAEE}\right|·\left(u_{{co}_2} -u_{FAEE}\right)$$

(4)

where: С_D – drag coefficient, dependent on the Reynolds number.

– Mass transfer model – LDF (Linear Driving Force):

$$R_{Ext.}=k_m\left(Y_{FAEE,eq}-Y_{FAEE}\right)$$

(5)

where: ${\mathrm{k}}_{\mathrm{m}}$ – mass transfer coefficient, which for the conditions of this work can be assumed to be in the range from 4∙10⁻³ to 7∙10⁻² s⁻¹; Y_FAEE,eq – equilibrium concentration of the FAEE mixture in scCO₂.

To describe the properties of scCO₂ and the FAEE mixture, the Peng-Robinson equation of state was used::

$$P={\displaystyle \frac{RT}{V_m-b}}-{\displaystyle \frac{a\left(T\right)}{V_m^2+2b{V}_m-b^2}}$$

(6)

where: a, b – parameters depending on the critical properties of the components, V – molar volume.

Subsequently, the solution methods in Ansys Fluent were selected. The choice of solvers was the Pressure-Based Solver (preferred for incompressible/weakly compressible flows) and Phase Coupled SIMPLE (for multiphase problems). Discretization schemes: spatial (second-order upwind for momentum; QUICK for the species transport equation) and temporal (First Order Implicit for stability with larger time steps). The flow regime in this system is laminar; however, in this problem setup, a strictly laminar model led to deteriorated solution stability. Therefore, to account for sub-grid mixing and improve stability, as well as to achieve solution convergence, the "k-ε Realizable" model with "Scalable Wall Functions" was applied to account for near-wall effects.

Within the framework of numerical modeling in ANSYS Fluent, it was necessary to replace the standard library material, used by default to represent the organic phase (diesel vapor), with user-defined materials corresponding to the actual composition of the extracted mixture – ethyl oleate and ethyl palmitate. These compounds are absent from the standard Fluent database, so a user-defined material library and a set of UDF functions were implemented to ensure reproducible property definition without manual data entry in the interface and without being tied to a specific project.

For integration into the calculation environment, a user-defined material database file in the Scheme language (faee_materials.scm) was created, including definitions for the individual components – ethyl oleate and ethyl palmitate – as well as a separate mixture material faee_25_75, corresponding to the mass ratio of 25.28/74.72 wt.%. This approach allows loading the created materials using Fluent's standard tools and using them in the calculation setup on par with standard materials. During the import stage, reference property values at a fixed temperature are set in the library, ensuring correct initialization of the calculation and ease of unit control.

To account for the temperature dependence of properties in the operating temperature range of –20…100 °C, a compiled UDF module in C was developed. This module is connected to Fluent as a compiled UDF and assigned to the relevant materials via the Materials menu for the following quantities: density, viscosity, heat capacity, and thermal conductivity. The UDF uses analytical approximations based on tabulated data of temperature property profiles, ensuring stable calculation of properties in each computational cell during the iterative solution of the momentum and energy equations. Additionally, a temperature range limit was introduced in the code to prevent incorrect extrapolation beyond the specified interval.

The properties of the faee_25_75 mixture are formed directly in the UDF based on the properties of the individual components using the following mixing rules: heat capacity and thermal conductivity were calculated as mass-weighted averages; density was calculated via specific volume; viscosity was determined using the logarithmic mixing rule.

This approach allows modeling the mixture as a single working fluid without including multicomponent transport, simplifying the problem setup and reducing computational costs while maintaining correct temperature-dependent properties.

After importing the created materials into the calculation project, ethyl oleate and ethyl palmitate (and their mixture faee_25_75) were used as components of the organic phase instead of the surrogate component "diesel vapor". This allowed the calculation to be linked to the actual composition of the extracted product and eliminated the uncertainty associated with the conditional properties of a generalized material. To describe the extraction process, the model included interfacial mass transfer between CO₂ and the organic phase (CO₂ dissolution/absorption), with the composition of the organic phase specified via the created mixture. Thanks to the implemented approach, the calculation setup ensures the correct participation of ethyl oleate and ethyl palmitate (and their binary mixture) in the mass and energy transfer processes during supercritical CO₂ extraction.

It should also be noted that the developed property model has limitations: it is intended for single-phase liquid flow. For the ethyl oleate/ethyl palmitate system at reduced temperatures, crystallization and the appearance of a solid phase are possible; therefore, calculations in the lower part of the range (below the solidification onset temperature) require careful interpretation. If it is necessary to simulate the phase transition in Fluent, it is recommended to use the standard "Solidification & Melting" (enthalpy-porosity) model with specified solidus and liquidus temperatures, or to restrict the calculation to the temperature range corresponding to a stable liquid phase.

The solution was initialized from the inlet, and a fixed number of iterations (1500) was set to guarantee full convergence. Convergence was monitored by the scaled residuals for the continuity, x-, y-, z-velocity, energy, k, and ε equations. Analysis of the convergence history showed that a stable (steady-state) solution was achieved at the 723rd iteration. By this point, the residuals of all equations had fallen below the specified criteria. Simultaneously, the integral parameters (mass flow rate and average temperature) reached a plateau, and their relative change over the subsequent 50 iterations did not exceed 0.1%, indicating complete stabilization of the solution. Further iterations, up to the 1500th, did not lead to any significant changes, confirming the achievement of a convergent solution.

3. Results and Discussion

The conditions and results of the experimental scCO₂ extraction separation of the binary mixture composed of ethyl oleate and ethyl palmitate are presented in Figure 4 and Table 2 as averaged values after triplicate experiments.

The authors established that the selectivity of the extraction separation of the ethyl oleate/ethyl palmitate mixture in supercritical CO₂ can be effectively regulated by varying the pressure and process time. The results presented in Table 2 show that increasing the pressure and the duration of the process leads to an increase in the content of ethyl palmitate in the extract relative to ethyl oleate. This is explained by the fact that, according to the phase diagrams (Figure 3), moving away from the critical point increases the density of the CO₂-ethyl palmitate binary mixture, which is directly related to the increased solubility of this component and, consequently, its yield. Thus, the obtained data have practical significance for the development of technologies aimed at reducing the viscosity of biodiesel fuel.

Statistical processing of the results presented in Table 2 was performed using Student's t-test. The significance of differences between the concentration values of ethyl oleate and ethyl palmitate under different extraction conditions was assessed at a significance level of p < 0.05. It was found that increasing the pressure from 11 to 14 MPa at a fixed process time leads to a statistically significant increase in the content of ethyl palmitate in the extract. This confirms the reproducibility of the experiment and the reliability of the identified patterns.

Based on the conducted experiments, the results of numerical modeling of hydrodynamics and heat transfer in the SFE extractor were analyzed for two operating pressures: 11 and 14 MPa, and t = 40 °C. The analysis focused on the fields of key parameters determining the fractionation process of the FAEE mixture: pressure, scCO₂ flow velocity, temperature, and the distribution of species mass fractions.

The distribution of static pressure in the apparatus is a key factor determining the density and, consequently, the solvent power of scCO₂. Figure 5 shows the pressure fields for the two operating modes.

The main pressure drop (hydraulic resistance) occurs in the packing zone, where the flow interacts with the surface of the spheres ( Figure 5 ). At an inlet pressure of 14 MPa, the absolute pressure values throughout the volume are naturally higher than at 11 MPa. However, more importantly, the pressure gradient in the packing zone at 14 MPa is 15-20% higher. This indicates a change in the flow regime and a more uniform penetration of scCO₂ through the packing layer at higher pressure, which is favorable for phase contact.

The flow structure directly influences the contact time and the efficiency of mass transfer between scCO₂ and the liquid FAEE mixture. Figure 6 presents the scCO₂ flow velocity fields.

Analysis of the velocity fields reveals a significant influence of pressure on hydrodynamics. At 11 MPa, the formation of a pronounced high-velocity channel is observed in the upper-central part of the apparatus, accompanied by the emergence of extensive zones with low velocities near the walls and in the inter-sphere space of the packing. This pattern indicates non-uniform flow distribution and the risk of stagnant zones where mass transfer would be limited.

When the pressure is increased to 14 MPa, the velocity profile becomes significantly more uniform. The average velocity in the inter-sphere channels increases, and the scatter of values across the apparatus cross-section decreases. This is associated with the increase in scCO₂ density and, consequently, a decrease in its kinematic viscosity, which promotes more efficient flow distribution. For the fractionation process, a uniform velocity field is critically important as it ensures homogeneous contact of the solvent with the entire mass of the fractionated liquid, prevents premature "breakthrough" of unreacted zones, and facilitates achieving uniform distribution.

To verify the assumption of isothermality of the process and assess the influence of operating parameters on the temperature field inside the extractor, temperature profiles along the apparatus height were analyzed. Figure 7 shows the temperature distributions along the central axis and near the extractor walls for the two investigated pressures – 11 and 14 MPa. Since the modeling results for pressures of 11 and 14 MPa showed coinciding temperature profiles, Figure 7 presents generalized dependencies for both pressures. Line 1 corresponds to the temperature on the central axis, line 2 – near the inner surfaces of the extractor walls.

As can be seen from the graph, the temperature at the wall remains practically constant along the entire height and equals the set value, which fully corresponds to the thermostatic boundary condition. Temperature deviations on the central axis may arise due to the compressibility of the working fluid or due to internal (viscous) friction, but with thermostated walls, they are quickly suppressed by heat exchange with the wall. The temperature profiles obtained for pressures of 11 and 14 MPa coincide – the differences between them do not exceed 0.1 °C, which is within the error of the numerical solution. This indicates that, within the framework of the adopted model, the temperature field is determined solely by the boundary conditions and does not depend on the operating pressure in the investigated range. The absence of significant temperature gradients confirms the correctness of the isothermal process assumption and allows, in the further analysis of hydrodynamics and mass transfer, to consider pressure changes as the main factor influencing density, viscosity, and, consequently, the efficiency of fractionation.

The analysis of temperature fields confirmed that the process proceeds under isothermal conditions, and the influence of temperature gradients on the properties of the supercritical fluid is negligible. This allows proceeding to consider the factor directly determining the extraction efficiency – density. It is density, dependent on pressure, that directly affects the solvent power of the fluid and the selectivity of extracting mixture components. Figure 8 shows the density distribution fields in the extractor for the 11 and 14 MPa modes, allowing a visual assessment of the thermodynamic conditions at each point of the apparatus.

Analysis of the density fields allows the following key observations. Throughout the working volume, the density at 14 MPa systematically exceeds the values at 11 MPa, which fully corresponds to thermodynamic expectations. Along the flow direction (from bottom to top), a smooth increase in density is observed, caused by the gradual saturation of CO₂ with esters and the minor influence of hydraulic resistance. The density gradient is more pronounced at 14 MPa, which is associated with the greater compressibility of the medium in this parameter range. In the near-wall zones, the density is somewhat lower than on the central axis, which is explained by the lower intensity of mass transfer near the walls. The difference between the density values at the center and the wall does not exceed 15–20 kg/m³ for both pressures, indicating a relatively uniform density field across the apparatus cross-section and being a favorable factor for the uniform course of the extraction process. The obtained density fields unequivocally show that increasing the operating pressure from 11 to 14 MPa leads to a significant increase in density at all points of the apparatus. Since the solvent power of scCO₂ with respect to lipophilic compounds, such as FAEE, has a direct correlation with density, it can be stated that the 14 MPa mode provides a higher extraction potential. This conclusion, based on the visualization of density fields, creates a solid thermodynamic basis for interpreting the possible increase in yield and selectivity in this mode.

After establishing the patterns of distribution of parameters determining the thermodynamic state of the system, the next stage was the analysis of mass transfer – the process responsible for the selective extraction of target components. For this purpose, separate CFD calculations of the transport of ethyl oleate and ethyl palmitate in supercritical CO₂ at pressures of 11 and 14 MPa were performed. Figure 9 shows the obtained profiles of the mass fraction distribution of ethyl oleate and ethyl palmitate along the extractor height. For each pressure, dependencies are presented reflecting the change in concentration of these components along the central axis and averaged values for the near-wall region (left and right walls). Since the mass fraction of carbon dioxide at each point is determined as the complement to unity, it is not shown in the figure. The presented dependencies provide a quantitative representation of the process selectivity and allow linking the observed changes in composition with the hydrodynamic and thermodynamic conditions discussed earlier. Averaging the values for both walls allows assessing the transverse non-uniformity of component distribution across the apparatus cross-section. The mass fraction values in Figure 9 are given relative to the zero level at the extractor inlet, which allows a clear visualization of the accumulation dynamics of each component along the apparatus height.

Analysis of the component mass fraction distribution profiles ( Figure 9 ) allows drawing the following conclusions: at a pressure of 11 MPa, the mass fraction of ethyl palmitate at the extractor outlet (Y = 380–400 mm) reaches 0.69, while the fraction of ethyl oleate is 0.51. This indicates the preferential extraction of the saturated ester – ethyl palmitate – under these conditions. Increasing the pressure to 14 MPa leads to an even more pronounced increase in palmitate concentration: its mass fraction rises to 0.76, whereas the fraction of oleate increases only to 0.58. Thus, throughout the entire apparatus height, the content of ethyl palmitate in the extract consistently exceeds that of ethyl oleate, and this difference amplifies with increasing pressure.

The identified pattern is explained by the difference in the solubility of the components in supercritical CO₂. According to the phase diagrams ( Figure 3 ), ethyl palmitate exhibits higher solubility in the investigated parameter range, which determines its preferential extraction. Increasing the pressure further increases the density of CO₂ and, consequently, its solvent power, which particularly affects the more soluble component – ethyl palmitate. Thus, both the calculated and experimental data confirm that the scCO₂ extraction process allows for the efficient separation of the binary ester mixture, enriching the extract with the high-viscosity saturated component – ethyl palmitate – and the raffinate with the fraction enriched in unsaturated ethyl oleate having reduced viscosity. This directly corresponds to the aim of the work – reducing the viscosity of biodiesel fuel.

4. Conclusion

This study addressed the pressing issue of improving biodiesel quality by reducing its viscosity through the selective removal of high-viscosity saturated fatty acid esters from a mixture with unsaturated esters. Fractionation using supercritical carbon dioxide was proposed as a promising method, the efficiency of which directly depends on the hydrodynamic and thermodynamic conditions within the extractor. A three-dimensional CFD model of the supercritical extraction process was developed and experimentally verified for a laboratory-scale extractor with a volume of 0.092 L using the ANSYS Fluent software environment. Since the target components – ethyl oleate and ethyl palmitate – are absent from the standard material database, a custom property library and compiled User-Defined Function (UDF) routines were implemented for the first time. These enable the reproducible specification of temperature-dependent density, viscosity, heat capacity, and thermal conductivity for both the individual substances and their binary mixture using mixing rules, ensuring a physically correct description of the working fluid's properties.

A comparative CFD modeling study was conducted for two operating pressures – 11 and 14 MPa – at a constant temperature of 40 °C. Analysis of the pressure, velocity, temperature, and density fields, as well as the distributions of component mass fractions, revealed a comprehensive positive effect of increasing pressure on the process. Increasing the pressure to 14 MPa led to a 15–20% increase in density throughout the extractor volume. Considering the direct correlation between density and solvent power, this creates a fundamentally higher extraction potential. Concurrently, a homogenization of the flow velocity field was observed, along with a reduction in hydraulic resistance within the packing layer and a decrease in the extent of stagnant zones. These changes promote more uniform and intensive contact between phases, which is critically important for efficient mass transfer.

Both the calculated and experimental data confirm that ethyl palmitate is extracted preferentially over ethyl oleate, and this selectivity intensifies with increasing pressure. The mass fraction of ethyl palmitate at the extractor outlet increases from 0.69 at 11 MPa to 0.76 at 14 MPa, while the fraction of ethyl oleate only rises from 0.51 to 0.58. Experimental concentrations at 14 MPa and an extraction time of 30 min reach 94.3 wt.% for ethyl palmitate and 5.7 wt.% for ethyl oleate, showing good agreement with the calculated profiles. The mass fraction of the target ester mixture in the extract increases by 9.5 wt.% upon raising the pressure, providing quantitative confirmation of process intensification.