Comparative Analysis of Turbulence Models for the Description of the Operation of a Submerged Combustion Apparatus

V. A. Demin; A. V. Kostyrya; В. А. Демин; А. В. Костыря

doi:10.23947/2687-1653-2026-26-3-2266

Comparative Analysis of Turbulence Models for the Description of the Operation of a Submerged Combustion Apparatus

V. A. Demin, A. V. Kostyrya

https://doi.org/10.23947/2687-1653-2026-26-3-2266

EDN: SWMOTQ

Full Text:

PDF (Rus) PDF (Eng) HTML HTML (Rus) XML XML (Rus)

Generate QR code

Contents

Scroll to:

Abstract

Introduction. Improving the submerged combustion apparatus (SCA) requires a detailed understanding of its internal hydrodynamic behavior. There are studies on the numerical modeling of liquefied natural gas regasification units. However, these units, unlike vaporizers, do not provide high-velocity gas flow injection. Therefore, the methodology for numerical modeling of SCV (Submerged Combustion Vaporizer) used to heat concentrated brines and contaminated media has not been fully developed. The objective of this study is a comparative analysis of turbulence models for numerical simulation of aerodynamic flow paths and substantiation of recommendations for their selection. Three tasks are solved. The first is the creation of a set of numerical SCA models. The second involves numerical experiments with them. The third is the identification of their similarities and differences based on quantitative and qualitative experimental results.

Materials and Methods. The Eulerian–Eulerian approach was used to describe the motion of liquid and gas in the liquid–gas–solid particle system, while the Eulerian–Lagrangian approach was used for the solid phase. The equations of the turbulence models were sequentially incorporated into a general physicomathematical model of the system, which was closed by the empirical relations of Schiller–Naumann, Ishii–Zuber, and Rantz–Marshall. The problem was solved using the finite volume method in both steady-state and transient modes (with a time step of 0.1 s). Turbulence was simulated in homogeneous and heterogeneous settings. The computational environment was Ansys CFX. The grid cell size was 10 mm.

Results. For the eight turbulence models, the specific turbulent kinetic energy, velocity, liquid volume, and liquid temperature were calculated at the computational cell level. The maximum for the first indicator was 0.361 (k – ω), the minimum was 0.026 RNG (k – ε). For the second — 0.457 (DES) and 0.130 (k – ω), respectively. According to the third — 1.339 (k – ω); 1.186 (DES). According to the fourth — 58.7 RNG (k – ε); 28.2 (k – ω). The velocity fields were visualized for both modes. A low scatter of the deposited solid phase mass was observed (1.18–1.38 kg).

Discussion. The EARSM (Explicit Algebraic Reynolds Stresses Model) reproduces a flow structure similar to k – ε and SST (Shear Stresses Transfer), thereby confirming the reliability of their results. The k – ω model yields an unphysical result due to the lack of strict symmetry of the velocity field. The homogeneous DES (Detached Eddy Simulation) predicts maximum free-surface asymmetry. The resulting inclination angle (approximately 30°) does not correspond to the actual hydrodynamics in the SCA; therefore, the homogeneous DES is inapplicable. Three cases of greatest reliability have been established: homogeneous k – ω (flow symmetry, minimum average temperature); heterogeneous RNG k – ε (maximum average temperature); homogeneous DES (maximum surface deformation). SST is recommended for stationary modes, LES (Large Eddy Simulation) and heterogeneous SST — for dynamic processes.

Conclusion. In numerical simulations of the SCA, the calculations have proved to be conditionally stable. With a few exceptions, the turbulence models are generally equivalent. Recommendations are provided for using the models in steady-state and transient modes. In the future, the reasons for the extreme temperature and surface curvature results produced by the models will be studied in detail.

Keywords

submerged combustion, multiphase flow, submerged gas-liquid jet, turbulence

For citations:

Demin V.A., Kostyrya A.V. Comparative Analysis of Turbulence Models for the Description of the Operation of a Submerged Combustion Apparatus. Advanced Engineering Research (Rostov-on-Don). 2026;26(3):2266. https://doi.org/10.23947/2687-1653-2026-26-3-2266. EDN: SWMOTQ

Introduction. Submerged combustion apparatus (SCA) are thermal engineering devices for the direct heating and evaporation of solutions. The mixture is heated via direct contact with a jet of hot flue gases generated in the combustion chamber. The key advantage of this heating method is the absence of heat-transfer surfaces, where salt precipitation and scale deposition typically occur. During operation of the SCA, the flue gas jet contacts the liquid and breaks into bubbles, which bubble through the liquid layer and transfer heat to it. Turbulent flow occurs in the liquid due to the high-velocity inlet of the flue gas jet. Therefore, the flows in the apparatus are calculated using one of the turbulence models. Several classes of such models are known [1].

Fig. 1. Classification of turbulence models

Figure 1 shows a general classification scheme for turbulence models. Most of them belong to the RANS class (Reynolds-averaged Navier-Stokes equations). The essence of averaging is to replace the true velocity with the ensemble-averaged velocity [2]:

(1)

Here, v — true velocity; 〈〉 — average velocity; v′ — velocity fluctuation. The Navier–Stokes equation takes the form:

(2)

(3)

Equation (3) represents the Reynolds stress tensor. Equation (2) remains open until the tensor components are determined. RANS class models define the Reynolds stress tensor.

A subclass of models using the Boussinesq hypothesis [3], is based on the proposition of proportionality of the stress tensor and the strain rate tensor. In this case, the proportionality coefficient is turbulent viscosity, and its determination becomes the key to describing the liquid motion. A number of models have been proposed for calculating turbulent viscosity, differing in the number of differential equations (DEs) used. The following models are known:

– without DE (zero-order models described in [4] and [5]);

– with one DE [6];

– with two and more DE (k – ε [7], k – ω [8], k – ε – n² [9]).

An alternative turbulence modeling within the RANS framework involves a departure from the Boussinesq hypothesis [10]. From a mathematical standpoint, the Reynolds stresses represent the second-order single-point moments of the velocity fluctuations, and the transport equations for these quantities can be rigorously derived from the Navier–Stokes expressions:

(4)

Here, D_ijk — diffusion term, P_ij — generation term, F_ij — diffusion term, e_ij — dissipative term. Determination of the relationship between the terms of equation (4) is the essence of the problem of constructing Reynolds stress transfer models. A set of algebraic and differential equations can be used.

Eddy-resolving models, instead of the stress averaging procedure, implement a variable filtering procedure [11]:

(5)

Large-scale turbulence eddies are resolved directly using filtered variables. Small-scale eddies, smaller than the grid cell size, are resolved using subgrid-scale models. The most well-known of these is large-eddy simulation (LES). In the limit of a vanishing grid size, LES reduces to direct numerical simulation (DNS) of turbulence, which is currently used only for fundamental research.

The combination of RANS and eddy resolution is the basis for a class of hybrid models. Among them, we note the detached eddy simulation (DES) model [12], as well as the group of scale-adaptive simulation (SAS) models [13]. This diversity, on the one hand, allows for the flexible selection of a turbulent model suitable for a specific computational case, and on the other hand, requires careful analysis and verification of the models.

For numerical experiments, an abstract representation of the SCA is used — a cell with a submerged jet and a free surface [14]. Its diagram is shown in Figure 2.

Fig. 2. Schematic of the test cell: submerged jet and free surface

From a computational fluid dynamics (CFD) perspective, this cell is a challenging computational case. The complexity stems from five factors:

turbulent flow;
laminar flow;
local instability in the inlet region of the submerged gas-liquid jet;
dispersed phase (gas bubbles);
free deformable surface.

Recommendations for describing these factors contradict each other¹. A multi-velocity continuum is proposed for modeling bubble flows, and a single-velocity continuum for modeling the free surface. Locally unstable flows require the use of eddy-resolving turbulence models, while laminar and turbulent regions require resolution of the laminar-turbulent transition. Researchers of gas-lift [15] and bubbling [16] systems use the k – ε model, but in such systems, there is no high-speed injection of the gas-liquid jet. Thus, it is highly challenging to make a definitive choice of turbulence model for the cell under consideration based solely on previous publications. At the same time, the national economy requires reliable technology for evaporating concentrated brines. These brines are either waste from potash fertilizer production or an intermediate product from underground salt leaching.

Materials and Methods. The following turbulence models were considered for the analysis:

k – ε;
k – ε with renormalization groups (RNG k – ε);
k – ω;
transfer of shear stresses (SST);
explicit algebraic Reynolds stress model (EARSM);
large eddy simulation (LES);
detached eddy simulation (DES);
zero-order (zero equation).

These models can be considered generally applicable, as they are all included in the engineering handbook for chemical industry workers [17]. The zero-order model, obviously, is not intended for solving this class of problems. It was included in the experiment for indirect verification of the numerical model. Verification will be considered successful if the calculation using the zero-order model yields unphysical results or end with an error.

The models under consideration implement different approaches to turbulence modeling:

Reynolds stress averaging with the Boussinesq hypothesis (k – ε, k – ω, SST) and without it (EARSM);
filtering of short-wave velocity fluctuations (LES and DES).

Furthermore, these models are integrated into the following computational packages: CFX, Fluent, FlowVision, and Logos.

The equations of the turbulence models under consideration were included in a general physicomathematical model of the hydrodynamic system under study. The general model describes a multiphase flow in which the carrier phase is a liquid, and the dispersed phase consists of gas bubbles and solid particles. The Euler–Euler approach [18], was used to calculate the gas-phase flow, whereas the Euler–Lagrange approach [19] was applied for the solid phase. Interphase interaction was calculated using empirical relationships:

momentum transfer between bubbles and liquid — the Ishii–Zuber correlation [20];
momentum transfer between solid particles and liquid — the Schiller–Naumann correlation [21];
heat transfer between bubbles and liquid — the Rantz–Marshall correlation [22].

Gas bubbles with a constant diameter of 3 mm are introduced into the combustion chamber nozzle exit. Gas density is variable, which is a consequence of strong temperature variations in the gas phase along the flow line.

Bubble coalescence was not calculated. Air was chosen as the gas, as the thermal properties of flue gases were close to those of the flue gas at an excess air ratio of 1.7. The solid phase inlet is located at the hypocenter of the combustion chamber, some distance from the nozzle exit. This distance is selected so that the inlet region is located within the gas-liquid jet, in the zone of maximum thermal stress. In the actual process, the liquid evaporates precisely at this point, which consequently results in salt precipitation. An open boundary allows liquid to escape from the computational domain, and gas to penetrate into it and partially escape. In this section of the boundary, the parameters of the gas phase beyond the computational domain are efficiently set: all walls are impermeable and adiabatic, plus no-slip conditions. Additionally, for the bottom wall (Fig. 1), a condition for complete absorption of the solid phase is introduced. This simulates the settling of solid particles on the bottom of the apparatus.

The equations of the physicomathematical model were solved by the finite volume method with the computational fluid dynamics (CFD) software package Ansys CFX 2019 R3.

Description of the mathematical model. The mathematical model of the cell is based on the Eulerian-Eulerian and Eulerian-Lagrangian equations. The equations of continuity, conservation of momentum, and energy are presented below:

(6)

(7)

(8)

The viscous stress tensor from (7) and (8) is calculated using the corresponding equations of the turbulence models mentioned above.

Law of motion for Lagrangian particles:

(9)

Forces in (7) and (9):

In the presented equations, α — volume fraction; ρ — density, kg/m³; v — velocity, m/s; t — time, s; p — pressure, Pa; H — specific enthalpy, J/kg; λ — thermal conductivity, W/m∙K; T — temperature, K; Q — amount of heat transferred from phase j, J/m³; c_p — heat capacity at constant pressure, J/kg∙K; d — diameter, m; m — mass, kg. Vectors , , , represent the acceleration due to gravity, m/s², specific force of interphase interaction, N/m³, specific force of added mass, N/m³, and the force of gravity, N; τ — viscous stress tensor. Subscripts i, j indicate arbitrary phases. Subscripts g, l, s correspond to the gas, liquid, and solid phases, respectively. C_VM — coefficient of added mass, set to 0.5. C_D — drag coefficient calculated using the Schiller–Naumann correlation for solid particles or the Ishii–Zuber correlation for gas bubbles.

The boundary and initial conditions described below are added to the basic equations.

The first is the combustion chamber nozzle exit (the conditions correspond to the experimental setup [23]):

Second — open border:

Third — liquid inlet:

Fifth — initial conditions:

The free surface of the liquid is bounded by plane z = 0.8 m.

The following dimensions of the computational domain were adopted: length and width – 1.3 m, height – 1 m. The number of computational grid nodes in three dimensions was 130:130:100. A uniform tetrahedral grid was used. To verify the independence of the calculation results from the grid cell size, numerical experiments were conducted using grids of varying sizes. The results of these experiments are presented in Table 1.

Table 1

Average Specific Energy of Turbulent Pulsations Depending on the Cell Size, J/kg

Cell, mm	Energy
20	0.045
15	0.052
10	0.051

As can be seen from Table 1, the average energy of turbulent pulsations (a key parameter for this study) for grids with 15- and 10-mm cells differs by 0.001 J/kg, which is 2% of the obtained values. Therefore, both the 15- and 10-mm grids can be used for further calculations. However, the 10-mm grid demonstrated greater calculation stability, so it was chosen for further work.

Numerical experiments were conducted in steady-state and transient modes. The time step in the transient mode was 0.1 s. In both modes, the single-continuum (homogeneous) and multi-continuum (heterogeneous) formulations were alternately applied for turbulence modeling. The remaining physical quantities were computed using the heterogeneous formulation.

Research Results. The following data were obtained as a result of the computational experiments:

average specific energy of turbulent pulsations (Table 2);
average liquid velocity in the cell at the end of the calculation (Table 3);
liquid volume in the cell (Table 4);
average liquid temperature in the cell (Table 5).

In the tables, the type of data that could not be obtained during the experiment is designated as n/a. Symbol х indicates cases where the calculation was not completed due to a solver error.

Table 2

Average Specific Energy of Turbulent Pulsations, J/kg

Mode	Formulation	Model
Mode	Formulation	DES	LES	k – ε	RNG k – ε	k – ω	SST	EARSM	Zero equation
Transient	Homogeneous	0.131	n/a	х	х	х	0.211	х	х
Transient	Heterogeneous	х	х	х	х	х	0.098	х	х
Steady-state	Homogeneous	–	–	0.065	0.068	0.361	0.142	х	х
Steady-state	Heterogeneous	–	–	0.043	0.026	х	0.073	0.065	х

Table 3

Average Liquid Velocity, m/s

Mode	Formulation	Model
Mode	Formulation	DES	LES	k – ε	RNG k – ε	k – ω	SST	EARSM	Zero equation
Transient	Homogeneous	0.457	0.292	х	х	х	0.402	х	х
Transient	Heterogeneous	х	х	х	х	х	0.269	х	х
Steady-state	Homogeneous	–	–	0.236	0.244	0.130	0.240	х	х
Steady-state	Heterogeneous	–	–	0.238	0.264	х	0.278	0.164	х

Table 4

Volume of Liquid, m³

Mode	Formulation	Model
Mode	Formulation	DES	LES	k – ε	RNG k – ε	k – ω	SST	EARSM	Zero equation
Transient	Homogeneous	1.186	1.320	х	х	х	1.194	х	х
Transient	Heterogeneous	х	х	х	х	х	1.338	х	х
Steady-state	Homogeneous	–	–	1.307	1.304	1.339	1.331	х	х
Steady-state	Heterogeneous	–	–	1.297	1.299	х	1.271	1.331	х

Table 5

Average Liquid Temperature, ºС

Mode	Formulation	Model
Mode	Formulation	DES	LES	k – ε	RNG k – ε	k – ω	SST	EARSM	Zero equation
Transient	Homogeneous	37.7	42.3	х	х	х	34.4	х	х
Transient	Heterogeneous	х	х	х	х	х	33.7	х	х
Steady-state	Homogeneous	–	–	39.7	42.7	28.2	34.0	х	х
Steady-state	Heterogeneous	–	–	44.0	58.7	х	37.6	32.5	х

Figures 3 and 4 present the liquid velocity fields with the gas bubble streamlines overlaid on them.

Fig. 3. Liquid velocity field (transient mode):
a — homogeneous DES; b — homogeneous LES; c — homogeneous SST; d — heterogeneous SST

Fig. 4. Field of liquid velocities (steady-state mode):
a — SST homogeneous; b — SST heterogeneous; c — k – ε homogeneous; d — k – ε heterogeneous; e — RNG k – ε homogeneous; f — RNG k – ε heterogeneous; g — k – ω homogeneous; h — EARSM heterogeneous

Additionally, the accumulation of solids at the bottom of the cell in a transient mode was studied. The mass of solid particles touching the lower boundary of the computational domain over the entire duration of the numerical experiment (60 s) was calculated. The data are presented in Table 6.

Table 6

Mass of Settled Solid Phase, kg

Formulation	Model
Formulation	DES	LES	SST
Homogeneous	1.290	1.377	1.253
Heterogeneous	х	х	1.182

As can be seen, the mass of the settled solid phase is in the range of 1.18–1.38 kg, indicating a low scatter of values. The difference between the two SST model formulations is insignificant (6%). The homogeneous SST and DES models yielded the smallest difference in results (3%), despite the differences in flow patterns.

Discussion. A comparison of the specific turbulent kinetic energy values (Table 2) shows that they are all of the same order of magnitude, ranging from tenths to hundredths of Joules per kilogram. The homogeneous formulation yields higher fluctuation energy values than the heterogeneous one.

The average liquid velocities (Table 3) are also of the same order of magnitude — tenths of a meter per second. The spread of pulsation energies and velocities is small.

The liquid volume values (Table 4) in the computational domain are close to each other for the steady-state mode. The exception is the heterogeneous SST model, with a volume value 0.04 m³ lower than the average for the other models. In the transient mode, the DES and homogeneous SST models yield volume values of approximately 1.19 m³, while the other two yield values of approximately 1.3 m³. This indicates a low degree of conservatism for the first two models, meaning they may yield unacceptable results over long-time scales.

The distribution of mean temperature values (Table 5) shows a different trend. Thus, in the transient analysis, the difference between the two SST model formulations is minimal (0.7°C), while the difference between the SST and eddy-resolving models is several degrees. The steady-state mode is characterized by a spread in mean temperature values. The minimum (28.2°C) is given by the k – ω model, and the maximum (58.7°C) — by the heterogeneous k – ε.

It should be noted that all calculations using the zero-order model resulted in errors, as expected. This is entirely consistent with expectations and indicates indirect verification of the numerical model.

As can be seen from Tables 2–5, twelve out of twenty-four calculations (excluding calculations with the zero-order model) resulted in errors. This indicates low computational stability for the numerical model containing both a free surface and a dense bubble flow. This instability is due to the need to calculate the position of the free surface in a heterogeneous formulation. As mentioned earlier, in a homogeneous formulation in a number of numerical experiments, only turbulence is modeled. To simulate flows with a free surface, it is recommended to use a complete homogeneous formulation². However, the need to calculate bubble flow precludes the use of a homogeneous approach.

The steady-state calculation mode turned out to be more stable than the transient one. This can be explained by the accumulation of numerical errors in the transient solution [24]. This effect is especially noticeable for multiphase problems, where a larger number of equations (a multiple of the number of phases) are solved. Thus, the steady-state mode can be recommended for modeling stabilized modes, when the time dependence of variables is not significant.

Examining the velocity fields obtained in the transient mode (Fig. 3), it can be observed that the DES and LES models produce a flow structure with a greater number of eddies than the SST. The free surface of the liquid turned out to be strongly deformed in the calculations with the DES and SST models within the homogeneous formulation.

Examining the velocity fields obtained in the steady-state mode (Fig. 4), we see that none of the models produces significant free surface deformation, even the SST. The latter yields a flow pattern with a slightly more complex free surface relief. The result obtained using the k – ω model is worth noting. Here, the free surface is smoothest, and the velocity field is symmetrical, which is essentially reminiscent of the reference flow pattern in the SCA [25]. The remaining models gave similar results in both homogeneous and heterogeneous formulations (Fig. 4 а–f, h): a two-eddy flow pattern with a predominant left-to-right flow direction across the computational domain.

The liquid velocity fields correspond to regions where the volume fraction of liquid exceeds 85%. This allows assessing the state of the free surface. In all the illustrations, we see a similarly shaped region of the gas-liquid jet. Its range is virtually identical for all the models and computational modes considered. The jet diameter on the free surface varies slightly, except in cases where the free surface is highly curved (Figs. 3 a and 3 c). The liquid velocity vectors are directed tangentially to the gas phase flow lines. This is consistent with previously obtained results from studies of the hydrodynamic situation in the SCA [23].

Let us analyze the flow structure in the computational domain. The problem formulation itself is asymmetric: the liquid inlet is located on the right side of the computational domain, and the liquid and flue gas outlet is on the left. This reflects the design of real SCA, where a symmetric inlet and outlet configuration of the working media is impractical and even technically impossible. Consequently, the velocity field in the cell must also be asymmetric.

The mirror asymmetry of all velocity fields, except for that obtained using the k – ω, was noted above. This is explained by the asymmetric eddy structure of the flow (three eddies cannot be symmetrical relative to the vertical cross-sectional plane). This is due to the use of eddy-resolving models. Therefore, the feature in question is caused by the inlet and outlet flows, on the right and left, respectively.

The SST model is a combination of k – ω and k – ε, and the contribution of each part is determined by a weighting coefficient [26]. The qualitative similarity of the SST and k – ε modeling results is explained by the prevailing weight of the k – ε component.

The EARSM model uses a different Reynolds stress averaging procedure and can therefore serve as a benchmark for comparing RANS-class models. According to Figure 4h, the flow structure reproduced by the EARSM model is close to that of the k – ε and SST models, confirming the validity of the results obtained using them. It can be concluded that the k – ω model yields an unphysical result, since strict symmetry of the velocity field cannot exist in an asymmetric problem.

It was previously noted that the homogeneous DES model reproduces the maximum asymmetry of the free surface of the liquid. The obtained inclination angle (about 30º) does not correspond to the real hydrodynamic situation in the SCA; therefore, the authors consider this unphysical result to be an extreme case and deem it inappropriate to use a homogeneous DES model.

When discussing the reliability of various models, it is worth highlighting the similarity between most of them in both qualitative and quantitative results. Let us consider three extreme cases:

homogeneous k – ω (flow symmetry, minimum average liquid temperature);
heterogeneous RNG k – ε (maximum average liquid temperature);
homogeneous DES (maximum free surface deformation).

The remaining models are generally equivalent. However, for steady-state conditions, the SST model is recommended. It has demonstrated the greatest computational stability and allows for resolving the laminar-turbulent transition. For cases where process dynamics must be considered, the LES and heterogeneous SST models are recommended.

Conclusion. Numerical simulations of the SCA were performed using various methods for resolving turbulent fluctuations. The results of the numerical experiments suggest a similar flow pattern for all turbulence models. Similarities were found across four criteria:

formed gas-liquid jet;
the same range and diameter of the gas-liquid jet;
clear boundary of the free surface;
velocity field asymmetry (except k – ω).

All models, except the heterogeneous RNG k – ε, k – ω and DES, demonstrate insignificant scatter in integral parameters such as turbulent pulsation energy, average velocity, temperature, and liquid volume. These same models reproduce the asymmetry of the forced liquid flow in the apparatus. The k – ω model proved insensitive to asymmetry, and DES reproduced a free surface with an inclination angle of 30º, which does not correspond to the hydrodynamics of the SCA.

Numerical experiments demonstrated low computational stability in a problem with both free surface and dense bubble flow. Half of the calculations failed due to a solver error. However, it is found that the steady-state mode is more stable. This allows us to recommend the steady-state mode for cases where parameter changes over time are not of interest. This happens, for example, when studying a steady-state mode or when there is no need to take into account the accumulation of Lagrangian particles. Furthermore, SST can be recommended for SCA calculations. This model has demonstrated the greatest computational stability and also allows for resolving the laminar-turbulent transition.

Further research involves a detailed study of the reasons why turbulence models reproduce such extreme results as elevated mean temperature or significant curvature of the free surface.

1. Ansys CFX Reference Guide. Release 2021 R1. Canonsburg: ANSYS Inc; 2021. 450 p.

2. Ansys CFX Tutorials. Release 2024 R2. Canonsburg: Ansys Inc., 2024. 1002 p.

References

1. Haçat G, Yılmazoğlu M, Çıbık A, El Messoaudi N, Miyah Y, Knani S, et al. Advanced Turbulence Modeling Approaches for Water Treatment Processes: A Comprehensive Review. Journal of Water Process Engineering. 2025;80:109102. https://doi.org/10.1016/j.jwpe.2025.109102

2. Reynolds O. On the Dynamical Theory of Incompressible Viscous Fluids and the Determination of the Criterion. Philosophical Transactions of the Royal Society. 1895;186:46–50. URL: https://archive.org/details/philtrans05206033/mode/2up (дата обращения: 02.07.2026).

3. Sharfarets BP, Dmitriev SP. Modeling of Turbulent Fluid Motion Based on the Boussinesq Hypothesis. Review. Nauchnoe priborostroenie. 2018;28(3):101–107. https://doi.org/10.18358/np-28-3-i101108

4. Prandtl L. Uber die Ausgebildete Turbulenz. ZAMM. 1925;5:46–50. URL: https://link.springer.com/chapter/10.1007/978-3-662-11836-8_60 (дата обращения: 02.07.2026).

5. Kármán T. Mechanische Ähnlichkeit und Turbulenz. In: Sonderdrucke aus den Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen. Mathematisch-physische Klasse. Frankfurt am Main: Weidmannsche Buchhandlung wurde; 1930. 19 p. URL: https://www.researchgate.net/publication/266936962_Mechanische_Ahnlichkeit_und_Turbulenz/stats (дата обращения: 02.07.2026).

6. Spalart PR, Allmaras SR. A One-Equation Turbulence Model for Aerodynamic Flows. AIAA Paper. 1992;1:5–21. https://doi.org/10.2514/6.1992-439

7. Launder BE, Spalding DB. Mathematical Models of Turbulence. London: Academic Press; 1972. 169 p.

8. Wilcox D. Formulation of the k – ω Turbulence Model Revisited. AIAA Journal. 2008;46(11):2823–2838. https://doi.org/10.2514/1.36541

9. Durbin PA. Separated Flow Computations with the k – ε – ν² Model. AIAA Journal. 1995;33(4):659–664. https://doi.org/10.2514/3.12628

10. Ben Hassan Saïdi I, Schmelzer M, Cinnella P, Grasso F. CFD-Driven Symbolic Identification of Algebraic Reynolds-Stress Models. Journal of Computational Physics. 2022;457:111037. https://doi.org/10.1016/j.jcp.2022.111037

11. Wenzhi Huang, Yanhui Li, Rundong Zhang, Zhuohang Wu, Min Gao, Yu-Xin Ren. The Intrinsic Filters of the Finite Volume Methods and a New Interpretation of the Implicit Large Eddy Simulation. Journal of Computational Physics. 2026;549:114603. https://doi.org/10.1016/j.jcp.2025.114603

12. Annapeh HF, Kurushina VA. Analysis of the Drag-Reduction Ability of the Layout and Cross-Sectional Shapes of Subsea Structures in the Critical Flow Mode. Advanced Engineering Research (Rostov-on-Don). 2024;24(2):135–147. https://doi.org/10.23947/2687-1653-2024-24-2-135-147

13. Menter FR. Turbulence Modeling for Engineering Flows (a technical paper from Ansys Inc). Canonsburg, PA: Ansys; 2011. 25 p. URL: https://cfd.spb.ru/agarbaruk/doc/2011_Menter_Turbulence%20Modeling%20for%20Engineering%20Flows.pdf (дата обращения: 02.07.2026).

14. Kostyrya AV. Numerical Modeling of a Submerged Gas-Liquid Jet of Various Dispersity. Bulletin of Perm University. Physics. 2024;(3):5–12. https://doi.org/10.17072/1994-3598-2024-3-05-12

15. Nana Qi, Kai Zhang, Gang Xu, Yongping Yang, Hu Zhang. CFD-PBE Simulation of Gas-Phase Hydrodynamics in a Gas-Liquid-Solid Combined Loop Reactor. Petroleum Science. 2013;10:251–261. https://doi.org/10.1007/s12182-013-0274-5

16. Saad N Saleh, Ahmed A Mohammed, Farah K Al-Jubory, Shahzad Barghi. CFD Assessment of Uniform Bubbly Flow in a Bubble Column. Journal of Petroleum Science and Engineering. 2018;161:96–107. https://doi.org/10.1016/j.petrol.2017.11.002

17. Ostrovsky GM, Abiev RSh, Barabash VM, Bilenko LF, Bogdanov VV, Britov VP. A New Handbook for Chemists and Technologists. Processes and Equipment for Chemical Engineering, part I. Saint-Petersburg: Professional; 2004. 848 p. (In Russ.)

18. Peng Li, Xuhui Zhang, Xiaobing Lu. Three-Dimensional Eulerian Modeling of Gas – Liquid – Solid Flow with Gas Hydrate Dissociation in a Vertical Pipe. Chemical Engineering Science. 2019;196:145–165. https://doi.org/10.1016/j.ces.2018.10.053

19. Yimeng Du, Zhendong Jin, Chengjun Zhang, Yan Cui, Yulong Wang, Rongxuan Hu, et al. A Solver-Agnostic Lagrangian Approach for Efficient Particle Tracking in Unsteady Multiphase Flows. Computers & Fluids. 2026;305:106914. https://doi.org/10.1016/j.compfluid.2025.106914

20. Mengdie Wang, Chao Ni, Xiangning Bu, Yaoli Peng, Guangyuan Xie, Zhongchao Tan, et al. CFD‐PBM Simulation of the Column Flotation Unit of FCMC: Importance of Gas‐Liquid Interphase Forces Models. The Canadian Journal of Chemical Engineering. 2023;101(10):5925–5940. https://doi.org/10.1002/cjce.24828

21. Weiling Li, Wenqi Zhong. CFD Simulation of Hydrodynamics of Gas — Liquid — Solid Three-Phase Bubble Column. Powder Technology. 2015;286:766–788. https://doi.org/10.1016/j.powtec.2015.09.028

22. Gallier S. Aluminum Combustion in Strong Convective Flows. Combustion and Flame. 2023;249:112598. https://doi.org/10.1016/j.combustflame.2022.112598

23. Demin V, Kostyrya A. Dynamics of Gas — Liquid — Solid Three Phase Flow in Submerged Combustion Plant. Bulletin of the Saint Petersburg State Institute of Technology (Technical University). 2022;63(89):78–83. https://doi.org/10.36807/1998-9849-2022-63-89-78-83

24. Sadrehaghighi I. Flow Uncertainty. Quantification within CFD. Verification & Validation. Annapolis: CFD Open Series; 2025. 103 p. https://doi.org/10.13140/RG.2.2.26397.03042/7

25. Alabovskiy AN, Udyma PG. Submerged Combustion Devices. Moscow: Izdatelstvo MEI; 1994. 256 p. URL: https://www.c-o-k.ru/library/document/12790 (дата обращения: 02.07.2026).

26. Yang XL, Liu Y, Yang L. A Shear Stress Transport Incorporated Elliptic Blending Turbulence Model Applied to Near-Wall, Separated and Impinging Jet Flows and Heat Transfer. Computers & Mathematics with Applications. 2020;79(12):3257–3271. https://doi.org/10.1016/j.camwa.2020.01.024

About the Authors

V. A. Demin

Perm State National Research University; Perm National Research Polytechnic University
Russian Federation

Vitaly A. Demin, Dr.Sci. (Phys.-Math.), Head of the Theoretical Physics Department

ResearcherID: A-2841-2017

Scopus Author ID: 7101803084

SPIN-code: 9695-7830

15, Bukireva Str., Perm, 614990

A. V. Kostyrya

Perm National Research Polytechnic University
Russian Federation

Aleksey V. Kostyrya, Postgraduate student of the General Physics Department

ResearcherID: PDX-9947-2025

Scopus Author ID: 58976342400

SPIN-code: 9954-4762

15, Bukireva Str., Perm, 614990

The hydrodynamics of a submerged combustion vaporizer was studied. For the first time, eight turbulence models were compared within a single computational framework. Models that accurately reproduce velocity, temperature, and surface area were identified. Recommendations for selecting models for steady-state and transient conditions were provided. This approach improves the design accuracy of similar industrial installations. The results are applicable to brine heating and contaminated media purification.

Review

For citations:

JATS XML

This work is licensed under a Creative Commons Attribution 4.0 License.

ISSN 2687-1653 (Online)

Username
Password
	Remember me
Not a user? Register with this site Forgot your password?

User

Advanced Engineering Research (Rostov-on-Don)

Comparative Analysis of Turbulence Models for the Description of the Operation of a Submerged Combustion Apparatus

Full Text:

Abstract

Keywords

For citations:

References

About the Authors

Review

For citations:

Cookies policy