Next Article in Journal
Effect of Concrete Micro-Mechanical Properties under the Coupled Corrosion of Sulfate and High Water Head
Next Article in Special Issue
Application of Artificial Intelligence Technologies to Assess the Quality of Structures
Previous Article in Journal
Stochastic Modelling to Analyze the Impact of Electric Vehicle Penetration in Thailand
Previous Article in Special Issue
Multiparameter Optimization of Thrust Vector Control with Transverse Injection of a Supersonic Underexpanded Gas Jet into a Convergent Divergent Nozzle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Standard and Modified SST Models with the Consideration of the Streamline Curvature for Separated Flow Calculation in a Narrow Channel with a Conical Dimple on the Heated Wall

1
Research Department, Saint-Petersburg State Marine Technical University, Lotsmanskaya 3, 190121 Saint-Petersburg, Russia
2
Research Department, Saint-Petersburg State University of Civil Aviation, Pilots Str. 38, 196210 Saint-Petersburg, Russia
3
Research Department, Baltic State Technical University, 1-a Krasnoarmejskaya 1, 190005 Saint-Petersburg, Russia
4
Research Department, N.E. Zhukovskii Central Aerohydrodynamic Institute, Radio Str., 17, 107005 Moscow, Russia
*
Author to whom correspondence should be addressed.
Energies 2021, 14(16), 5038; https://doi.org/10.3390/en14165038
Submission received: 19 May 2021 / Revised: 2 August 2021 / Accepted: 3 August 2021 / Published: 17 August 2021
(This article belongs to the Special Issue Advanced Cyber-Physical System for Engines and Power System)

Abstract

:
The testing of the standard and modified SST models of the transfer of shear stresses was carried out on an example of calculating the heat transfer with an intense detached flow around a conical dimple with a slope angle of 45 ° on the heated wall of a narrow channel. It was shown that the standard turbulence model by Menter SST (MSST) of 2003, widely used in the packages Fluent, CFX, StarCCM+, etc., significantly underestimated the intensity of the return flow. A correction of this model was presented that took into account the influence of the curvature of streamlines within the framework of the Rodi-Leshziner-Isaev (RLI) approach for spatial separated flows. It was found that the predictions for the RLI MSST 2003 were close to the predictions for the original standard MSST 1993, in which the eddy viscosity was calculated using the vorticity modulus. At the same time, the predictions based on the modified one, following Smirnov-Menter (SM) MSST 2003, included in the ANSYS model catalog did not differ too much from the standard MSST 2003. The preference of the MSST modified within the RLI 2003 for calculating the heat transfer in intense separated flows was substantiated.

1. Introduction

Progress in computer systems related to the growth of their performance, development, and a wide application of software technologies in the world, as well as the improvement of differential-type turbulence models and vortex-resolving methods drastically changed the state and scope of engineering methods for the calculation of turbulent flows in comparison to the known approaches—for example, monographs [1,2,3]. The present article focuses on one of the most popular closure models of the Reynolds-averaged Navier–Stokes (RANS) equations—the Menter shear stress transport model (SST model) [4]. This model with two differential equations is one of numerous semi-empirical models presented in reference [5].
Constructing such models is associated with introducing semi-empirical constants into the equations. These constants are scaled with typical simple turbulent flows, including boundary and shear layers, as well as jets. Work [6] contained a set of experimental test benches for the determination of the semi-empirical constants of one of the most successful turbulence models—the dissipative two-parameter high-Reynolds number Launder–Spalding k ε model [7]. Improvements the standard SST model were carried out in works [8,9,10,11].
The present article pays attention to the actual fundamental and practical problems of how to model separated turbulent flows. As known, conventional turbulence models are constructed for fully developed turbulent flows. However separated flows combine different-scale mixing shear layers and zones of low-velocity backflows that cannot always be considered as fully developed turbulent flows. In these flows, stagnation (deceleration) zones can also appear at the places where near-wall flows are separated: with the regions of laminar flow and laminar-turbulent transition. It is obvious that the standard low-Reynolds number turbulence models constructed to show accelerating and decelerating channel flows with pressure gradients are not quite suitable to calculate complex flows. Moreover, 3D separated flows are characteristic of flow zones with intense swirling flows, as well as with localized vortex regions where the flow velocity is very moderate. It is also known that near-wall flows of the boundary layer type and free jets are characterized by different-value turbulent constants [1]. Thus, it is quite justified that the standard turbulence models constructed for typical flows are not always suitable for displaying complex vortex structures characteristic of 3D separated flows. They require corrections—in particular, those taking into account the streamline curvature influence. Therefore, for the Launder–Spalding k ε model [7] applied to two-dimensional twin parallel and concentric outflowing jets, Rodi and Leschziner [12,13] proposed the correction function:
f μ = 1 1 + C × R i t
that was multiplied by the semi-empirical constant C μ = 0.09 when determining the vortex viscosity coefficient. The additional constant C = 0.57 was determined analytically, and R i t was the turbulent Richardson number. It should be emphasized that the product ( C μ f μ ) was within the range of 0.02–0.15. In works [14,15,16,17,18,19,20,21,22], some results of the influence of curvature and rotation on the modeling of the vortex and separated flows were presented.
The Rodi–Leschziner approach to correcting the Launder–Spalding k ε model was used to solve a number of tasks [23,24,25,26,27,28] using physical analogs, but the additional constant C c had to be decreased to 0.1. External sub- and transonic separated flows around both a classical-shape circular cylinder and a sphere and tandems: two co-axial different-sized discs, a disc and a cylinder, and a stepped cylinder, as well as a Volkswagen car profile near a movable shield, were considered. It was shown that the numerical predictions of the integral and local force loads agreed well with the available experimental data.
F. Menter in his works published in the years 1993 and 2003 [29,30] proposed two SST models differing mainly in the vortex viscosity coefficient. The Menter 1993 SST model uses the vorticity modulus Ω and the advanced model—the strain rate tensor modulus S . It should be noted that the modified Menter 2003 SST model approaches all similar differential turbulence models containing the modulus S [7].
The Menter 1993 SST model was tested in detail in references [31,32,33,34,35,36,37,38,39,40,41,42,43,44,45], dealing with the calculations of 2D and 3D separated incompressible liquid and compressed gas flows. It was successfully tested for modeling turbulent recirculating flow in a square cavity with a movable cover (in comparison to other semi-empirical models) [31], a separated flow in a circular cavity on the walls of plane-parallel [32] and divergent [33] channels (S.V. Guvernyuk’s experiments, Research Institute of Mechanics, Moscow State University), a flow around a trench on a plane [46,47] (the experiments of A.I. Leontiev et al., N.E. Bauman Moscow State Technical University), discrete roughness elements in pipes [36] (G.A. Dreitser’s experiments, Moscow Aviation Institute), separated flow and heat transfer in a tube bundle [48] (A. Zhukauskas’ experiments), in the vicinity of spherical [38,49] and oval [50] dimples on a plate (the experiments of S.Z. Sapozhnikov et al., St. Petersburg Polytechnic University, as well as of G.I. Kiknadze et al.), convective heat transfer in the plane-parallel and divergent channels with a set of spherical [51] and conical [42] dimples placed on a wall (the experiments of Chyu et al., as well as of Banker et al.), and a supersonic axisymmetric flow around a tear-shaped body with a conical tapered flare [43] (A.N. Mikhalev’s experiments, A.F. Ioffe Physical-Technical Institute). The numerical study of unsteady turbulent flows around a circular cylinder under sub- and supercritical conditions in references [44,45] contained comparable available experimental data and numerical predictions obtained with the use of different software and turbulent models. The methodical tasks of how to develop the multiblock computational technique [52] for a solution of the Reynolds equations and their correct uses for the numerical simulations of the flow and heat transfer of bodies were discussed.
It should be mentioned that, when separated flows are calculated, the SST model with the vorticity modulus does not require vortex viscosity correction, like the k ε model. As mentioned in reference [53], constant vorticity cores appear in separated flow zones at high Reynolds numbers. These cores hinder the process of pumping vortex viscosity into them. At the same time, when the standard two-parameter dissipative model [23] is used, regions of abnormal high vortex viscosity are formed in the large-scale cores. To eliminate these regions, streamline curvature corrections must be introduced. The realization and testing of the modified Menter 2003 model [30] again raised the problem of the correct resolution of separated flows, since, in the modern semi-empirical turbulence models with the strain rate tensor modulus, false turbulence generation occurs in large-scale vortex cores. The additional constant C c associated with correcting the vortex viscosity in the shear stress transport model was redefined and taken as equal to 0.02. The scope of application of the Menter 2003 SST model was significantly expanded in the calculations of the trans-, super-, and hypersonic separated flows, as well as of the 3D jet-vortex structures. The approach to correcting the Menter 2003 SST model with consideration of the streamline curvature influence on the vortex viscosity coefficient was generalized by S.A. Isaev and was called the Rodi–Leschziner–Isaev approach (RLI) [54].
In general, the modified standard semi-empirical turbulence models allowing for the flow swirling and streamline curvature influence on turbulence characteristics are widely encountered in engineering practice. For example, Shur and Spalart [55,56] introduced the correction of the Spalart–Allmaras vortex viscosity with the consideration of the flow swirling in a one-parameter model. A similar approach was used by Smirnov and Menter (SM) [57] to modify the Menter 2003 SST model [4]. It is important to emphasize that, in this case, such an integral characteristic as the vortex viscosity does not undergo correction but only one of the constants in the equation for the specific dissipation rate.
Systematic comparisons of the predictions by means of two standard SST models and the modified Menter 2003 model [30] started with the use of the RLI and SM approaches that are in progress in this study.
Approving and testing the RLI approach when comparing the numerical predictions obtained by different semi-empirical models and software with the experimental data were made to solve the traditional tasks of recirculating the flows, as well as of the flow and heat transfer control due to forming intense large-scale vortex structures [58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77]. The first group of tasks included recirculating flows in a square cavity with a movable cover [53] in a round cavity on the wall of the plane-parallel [53,58] and U-shaped [59,60,78] channels for the periodic flow around the airfoil NACA0012 when the angle of attack was varied from 0° to 180° [61]. The predictions of the local and integral flow characteristics are in a good agreement with the experimental data of S.V. Guvernyuk [61] (Research Institute of Mechanics, Moscow State University) and Ian Castro [77] (University of Southampton).
The second group of tasks is characterized by interesting aerodynamic and thermophysical effects.
Recirculating flows in round vortex cells located on the plane-parallel channel wall [79] or on the backside of the 20%-thick airfoil MQ [62] are enhanced due to distributed suction (over the central body surface) and concentrated (slot) suction. In these cells, a maximum backflow velocity achieves a characteristic incoming flow velocity. It has been shown that vortex cells allow the carrying surface characteristics to be significantly increased at the Mach numbers from 0 to 0.55. A good agreement of the predicted and experimental profiles was obtained in a vortex cell on a channel wall over a wide range of the sucked air consumption coefficients [79] (S.V. Guvernyuk’s experiment, Research Institute of Mechanics, Moscow State University).
For the thin disc–cylinder tandem at the sub- and supersonic air incoming flow velocities, when the rational sizes of the disc and the gap between it and the cylinder end have been chosen, it is possible to reach the close-to-zero drag of the tandem layout of the body at low subsonic velocities and to significantly reduce the drag at supersonic flight speeds. The numerical predictions of the integral and local loads on the bodies in axisymmetric flow agree quite fairly with the results of the aeroballistic experiments [64,65] (A.N. Mikhalev, A.F., Ioffe Physical-Technical Institute) and the tube tests [80,81,82] (V.N. Konovalov, Aerospace Academy named after A.F. Mozhaiky) and A. Roshko [83] (CalTech).
The comparison of the predictions and experiments of V.I. Zapryagaev (S.A. Khristianovich Institute of Theoretical and Applied Mathematics SB RAS) revealed a fair agreement of the characteristics of a supersonic jet flowing into a submerged space [84].
The cycle of numerical and physical studies is concerned with the task of a subsonic periodic flow around a semi-circular cylinder at R e = 5 × 10 4 [69,70]. It appeared that the ultra-low coefficients of the lift force C y of the order of −1 [85] (it was established by Z.P. Sluchanovskaya, Research Institute of Mechanics, Moscow State University) were seen for such a body at the angles of attack close to zero. This effect was explained by the numerical predictions. A detailed numerical study was made using different software, semi-empirical models, and numerous grid structures. A physical experiment was also done to assess the influence of side washers on the integral characteristics, since these washers allow us to reach the quasi-two-dimensional regime of flow around a body [32,33] (S.V. Guvernyuk’s experiment, Research Institute of Mechanics, Moscow State University). It appeared that the standard near-wall SST models in the all considered packages (VP2/3 and Sigma-flow) predicted the significantly underestimated (up to 30%) power loads on a semi-circular cylinder, although the Menter 1993 SST model, unlike the Menter 2003 SST model, yielded somewhat better results, since it reproduced the separated and vortex flows with less errors. The streamline curvature corrections led to a significant decrease of the vortex viscosity level in the cores of the vortices and vortex flow intensification and contributed to a better agreement between the numerical predictions and experimental data. The RLI approach turns out to be preferred over the SM approach. A better agreement of the predictions obtained by the Menter 2003 SST model corrected according to the RLI approach with the use of different software and grids (multiblock and composite) was demonstrated.
A good agreement of the results of the numerical and physical modeling of convective heat transfer in a narrow channel with single spherical dimples of moderate [86] and average [87] depths (V.I. Terekhov’s experiment, S.S. Kutateladze Institute of Thermophysics SB RAS) and in a channel with an in-line spherical dimple bundle [42] (M.A. Gotovsky’s experiment, I. Polzunov Central Boiler and Turbine Institute) should also be mentioned. A recent publication [74] was concerned with the abnormal enhancement of a separated flow and heat transfer in an inclined oval-trench dimple in a stabilized area of the narrow channel with the use of different methods. A general quantitative agreement was found, and the disagreement of the zonal models was explained by the different turbulence models ( k ε Realizable, SA, and SST with the Durbin and Rodi–Leschziner–Isaev corrections) and by the applied software (StarCCM+ and VP2/3). Attention was paid to the low-pressure zone on the entrance spherical portion of a dimple with a dimensioned by double-hydraulic head pressure of −0.35. The negative relative friction in this zone attained a value of f / f p l −4.5, and the relative heat transfer N u / N u p l exceeded 5.
A good enough agreement of the predictions and experiments on the static pressure distribution was obtained by analyzing the 3D flow around models of high-rise buildings [75] (S.V. Guvernyuk’s experiment, Research Institute of Mechanics, Moscow State University).
The present study continues the cycle of calculations [53] of separated flows and heat transfers in the vicinity of a conical dimple with an inclination angle of 45 ° and a depth of 0.233 (in the spot diameter fractions) on the narrow channel wall—the computer analog of V.I. Terekhov’s setup [86] with the use of different SST turbulence models (S.S. Kutateladze Institute of Thermophysics SB RAS).

2. Problem Statement

Turbulent heat transfer in the air flow in a narrow channel with a single conical dimple at an inclination angle of 45 ° on a heated wall (Figure 1) was calculated. The stabilized area of flow in the channel, corresponding to the experiments [86], was considered. The profile of the fully developed turbulent flow was assigned at the computational domain input. The dimple spot diameter D and the bulk flow velocity U were chosen as characteristic parameters. The Reynolds number was 104. The channel height was equal to 0.33, and the channel width was 2.5. The conical dimple depth h / D = 0.233 . The edge rounding radius was equal to 0.025. The turbulence level was assigned close to the experimental one (1.5%). The turbulence scale was of the order of the characteristic sizes. The Prandtl number was taken as equal to 0.7.
The heat and hydrodynamic tasks (the natural convection effects were not taken into account) were solved separately. The flow at the channel inlet was isothermal and had a temperature T r e f = 293 K . The lower streamlined channel wall with a dimple was heated. The supplied constant heat flux q was recalculated in a dimensionless form by the formula:
q ¯ = q λ P r R e T r e f D
where λ is the thermal conductivity. The side channel walls are adiabatic. The upper wall is isothermal and has a temperature T r e f taken as the dimensionalization scale. The soft boundary conditions are assigned at the channel outlet for T . The order of q ¯ is 3.408 · 10−5.
The task of convective heat transfer in 3D steady fully developed turbulent incompressible fluid flow in the channel with a single conical dimple on the heated wall was solved by means of the mathematical model based on the system of the Reynolds-averaged steady Navier–Stokes equations (RANS) and the energy equation. To close these equations, the modified Menter shear stress transport (SST) models, which proved themselves well in applications for typical near-wall flows, including separated flows [29,30], were used. The k ω model proposed by Menter is a generalization of two turbulence models: the Launder–Spalding k ε model for shear flow zones far from the wall [7] and the Saffman–Wilcox k ω model for the near-wall region [5]. The ideas realized in the Johnson–King turbulence model [88] were used to construct the zonal shear stress transport model. As already mentioned, to determine the vortex viscosity, the Menter 1993 model [29] used the vorticity modulus Ω , and the modified Menter 2003 model [30], as most of the differential-type semi-empirical models, included the strain rate tensor modulus S into the expression for vortex viscosity. It is important to emphasize that the semi-empirical models are mainly scaled with near-wall flows. Therefore, the necessity appears to adjust them to high-intensity separated flows. The Rodi–Leschziner–Isaev approach to correcting the Menter 2003 SST model means that the vortex viscosity coefficient μ t is subjected to the correction function:
f c = 1 1 + C c × R i t
where the additional constant C c = 0.02 (index c stands for corrected) and the product f c × C μ is to be within the range 0.02 < f c × C μ < 0.15 (in the standard k ω model [89], the semi-empirical constant C μ = 0.09 is used in the expression for vortex viscosity). It is fundamentally important to determine the turbulent Richardson number ( R i ) for 3D flows (for 2D flows, it depends on a local flow velocity and curvature radius [12]). Such a formula is proposed:
R i t = ( 1 C μ ω ) 2 ( × V ) ( V | V | × ( V r V | V | ) )
where r is the radius vector, V is the local flow velocity vector, and C μ = 0.09 ; ω is the specific dissipation rate.
The near-wall conditions for SST models were formulated so that the normal derivative to the wall for turbulence energy was equal to zero. The specific dissipation rate of the turbulence energy in the near-wall cell was determined according to reference [30].
The methodology of solving the Reynolds-averaged steady Navier–Stokes equations is based on the concept of process splitting and is realized in the pressure correction procedure SIMPLEC [90,91,92,93], able to simulate convective heat transfer on multiblock overlapping structured grids. The developed methodology has been applied to cluster (multicore and multiprocessor) systems and generalized to cases of unstructured grids. Like reference [90], the generalized transport equation is written in increments of dependent variables. The explicit hand-side of the equation is discretized according to the second-order approximation schemes (QUICK Leonard’s quadratic upwind interpolation [94] for the convective terms of the momentum equations, the TVD scheme [95] for the convective terms of the equations for turbulence characteristics, and the central-difference scheme for the diffusion terms). The implicit hand-side is presented with the use of the upwind scheme with one-sided differences. The use of the centered grids with the dependent variables located in the centers of computational cells leads to the need to monotonize the pressure fields by introducing Rhie–Chow correction [96,97]. The solution methods of algebraic equations are the preconditioned BiCGSHAB [98] with the AMG preconditioner from Demidov’s library (amgl) [99] for pressure correction and the ILU0 preconditioner for other variables.
The developed VP2/3 code (Velocity–Pressure, 2D/3D) uses the multiblock computational technique (MCT) described in reference [90]. This technique is based on a set of different-scale, tier, and overlapping structured grids around the structural elements of the physical tasks of the relevant scales. In the two rows of the near-boundary cells of each of the overlapping or overset grids, the parameters are determined through linear interpolation [90].
The computational domain is the space bounded by the walls of the narrow channel with a length of 6.7. The computational domain input is located at a distance of 1.4 from the projection of the dimple center onto the lower wall plane taken as the center of the Cartesian coordinate ( x , y , z ) system.
To solve the task, the basic multiblock computational grid containing 2.187 million cells is used. The grid consists of four different-scale fragments: (1) a rectangular channel grid CG (with a minimum step of 0.08 along the longitudinal coordinate and of 0.07 across the transverse coordinate) with densely packed nodes to the wall. It is meant for displaying flow in the channel; (2) a detailed grid meant for a detailed description of the flow in the near wake of the dimple, a Cartesian grid WG covering the 3.3 × 0.1 × 1.7 vicinity of the dimple (the front boundary is at a distance of 0.44 from the center); (3) a cylindrical grid near the dimpled curvilinear wall adjacent to the spherical dimple surface, an elliptical grid Ring with densely packed nodes to the wall; and (4) an oblique grid covering the axis of the dimpled cylindrical grid—“0.2 × 0.2 patches”. The near-wall step is chosen equal to 10−4. The channel grid contains 99 × 83 × 87 cells, and the detailed rectangular grid around the spherical dimple—132 × 35 × 59 cells. For detailed grid 2, the minimum value of the longitudinal and transverse steps is equal to 0.03. The cylindrical elliptical grid around the dimple contains 104 × 83 × 143, cells and the oblique grid covering the axis—19 × 83 × 19 cells.
The present work compares the local and integral characteristics of flow and heat transfer in the dimpled channel for different SST models.
The tables contain the integral values of the heat transfer N u m m from the control area of the dimpled wall and from the surface inside the dimple, as well as the hydraulic loss coefficients ζ on the sections of the dimpled channel. The thermal performance N u m m / N u m m p l (index pl stands for plate) and the relative hydraulic loss coefficients ζ / ζ p l on the 2.5 × 1.5 control area around the dimple with a 0.5 center shift downstream are determined by the methods of references [51,55]. The thermal performance ( N u m m d / N u m m d p l ) (index d stands for dimple) of the surface bounded by the dimple and the relative hydraulic loss coefficient ( ζ d / ζ d p l ) on the area of the channel with cross-sections at the leading and trailing edges of the dimple are also determined by the above methods. The thermal-hydraulic performance is determined as the ratio of the thermal performance to the relative hydraulic loss:
T H P = N u m m / N u m m p l ζ / ζ p l
The characteristics of N u m m / N u m m p l and ζ / ζ p l and N u m m d / N u m m d p l and ζ d / ζ d p l are presented in Table 1 and Table 2 as a fraction in small numbers.
The compared parameters include the distributions of the pressure drop p p p l on the streamlined wall with the dimple and on the plane-parallel channel wall, the relative friction f / f p l , and the wall temperatures in the middle longitudinal and transverse sections of the dimple. The distributions of the relative Nusselt numbers, both local N u / N u p l in the middle section of the dimple and integrated N u m / N u m p l over the transverse strips along the longitudinal and transverse coordinates of the control section, are compared. The dependence of the Cartesian flow velocity components u , v   and   w ; the turbulence coefficient k ; the vortex viscosity μ t ; and the wall temperature T w on coordinate y are compared.

3. Grid Convergence Justification

The influence of the computational grids on the uncertainty of the numerical predictions of characteristics of the separated flow and heat transfer in the narrow channel with a conical dimple is analyzed for the SST model corrected within the RLI approach.
Three multiblock different-scale grids are considered: simple (S), simplified (A), basic (B), and refined (C), containing approx. 984,000 (S), 1,561,000 (A), 2,187,000 (B), and 2,667,000 (C) cells.
The simple multiblock cell (S) is rough. Unlike the basic cell, it has a rather rough rectangular grid WG with minimum steps in the span-wise and streamwise directions Δ x = Δ z = 0.05 . The space grid is very rarefied inside the dimple and has a number of computational cells that is almost an order of magnitude smaller than that in the basic cell. The simplified grid (A) differs from the basic one by a smaller number of cells in the curvilinear O-type cell in the area of the conical dimple—745,000 cells instead of 1,211,000 cells.
The grid (C), refined in comparison to the basic one, has a more detailed grid WG in the channel with the minimum steps Δ x = Δ z = 0.05 in the streamwise and spanwise directions.
Some of the analysis results of the grid convergence and the numerical predictions are listed in Table 2 and Table 3.
A noticeable refinement of the predictions of the extreme characteristics is seen at a switch from the calculations on the rough multiblock grid (S) to the calculations on the grid (A) with an increasing number of cells in the area of the conical grid. The differences of p m a x , p m i n , u m a x , u m i n , v m a x , v m i n , w m a x , w m i n , k m a x , and μ t m a x on the considered grids are of the order of 5–10%. The switch from the basic grid (B), at which the turbulence models are compared, to the improved grid (C) slightly influences the characteristics.
The integral thermal and hydraulic characteristics of the control section of the channel with the conical dimple undergo less influence of the multiblock grids in comparison to the local extreme separated flow parameters. A small difference is seen in these characteristics when the rough grid (S) is replaced with the refined grid (A). A small change in these characteristics is seen at the switch from the rough grid (S) to the refined grid (A) within the dimple.
The data of Table 2 and Table 3 demonstrate the multiblock grid convergence results with an increase in the number of cells and the packing of cells in the dimple region.

4. Results and Discussions

Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 and Table 2 and Table 3 present some of the results of the methodical study.
The local and integral characteristics of the separated flow and heat transfer in the channel with a conical dimple of depth 0.233 and an inclination angle of 45 ° are calculated with the use of four SST turbulence models. The Menter 1993 [29] and 2003 SST models [4], as well as the two Menter 2003 models with the consideration of the streamline curvature influence corrected within the framework of the RLI [54] and SM [57] approaches, are considered. All the calculations used the basic multiblock grid.
Figure 2 illustrates the separated flow streamlines in the conical dimple and the computer visualization of the jet-vortex structure in it. They were obtained by means of the Menter 2003 SST model corrected according to the RLI approach. The calculation results obtained with the use of the other models are in qualitative agreement with the presented flow streamlines.
The streamlines in Figure 2a are the trajectories of the liquid particles in the layer of computational cells adjacent to the curvilinear wall with indication of their direction. The calculated streamlines are the experimental visualization of the air flow along the wall obtained with the use of soot–oil coatings. They allow analyzing the vortex structures of 3D separated near-wall flows.
Figure 2b shows the 3D trajectories of liquid particles at the singularity points of source and sink type in the source streamline. The coordinates of the points are: the focus on the left slope (−0.144, −0.183, and −0.282) is the red strip; the sink on the right slope (−0.144, −0.183, and −0.282) is the blue strip; the focus in the center of the bottom (−0.088, −0.23, and 0.005) is the magenta strip; and the focus at the bottom (−0.128, −0.23, and −0.168) is the green strip.
As follows from Figure 2, the separated steady flow zone covers the entire inner region of the dimple. The flow separation and attachment lines coincide with the leeward and windward rounded edges of the dimple, respectively. For deep dimples, the asymmetrical flow with an inclined jet-vortex structure inside the dimple is similar to the flow in the vicinity of the spherical dimple on the plate and on the narrow channel wall [100,101]. In the considered case, the near-wall air flow enters the conical dimple through the gap between the separation and attachment lines on the left side of the dimple. On the right side, air flows into the gap between the separation line, deviating from the dimple edge and the sink line of the interacting flows on the windward and side slopes. On the inner surface of the dimple, the rather complex separated flow is divided into a flow around the inclined slopes and a flow at the flat bottom similar to the case when the flow is formed inside the cylindrical dimple [102]. On the left slope near the bottom, the source streamline shows the singularity point of focus (source) type, indicating the place where the inclined tornado-like vortex (the red strip in Figure 2b) is generated. Air enters this place mainly when it flows along the dimple slopes, bypassing the bottom. The external flow entering the dimple from above interacts with the windward slope of the dimple.
Air is sucked into the spiral vortex not only from the surface of the left side wall but also when it flows on the right and leeward inclined slopes. The source streamline on the right side of the slope shows the occurrence of the sink-type singularity point, where the external flow turns in the vicinity of the windward edge of the dimple (the blue strip in Figure 2b). The separated air flow at the conical dimple bottom is divided into two zones adjacent to the windward and leeward slopes. The separation line of the zones illustrates the occurrence of two focus-type singularity points, in the vicinity of which, column-shaped spiral vortices are formed. They are shown by green and magenta colors (Figure 2b).
Figure 3, Figure 4 and Figure 5 compare the pressure fields, the relative friction, and the Nusslet number values on the streamlined surface of the control area of the heated wall with a conical dimple. The air source streamlines on the dimpled wall are plotted in the fields of the characteristics.
Figure 3 shows a sharp change in the pressure fields in the zone on the left side of the windward rounded edge of the conical dimple, where the flow stagnation region is also located. The low-pressure region is located behind the dimple in the vicinity of the surface curvature discontinuity line. By the maximum absolute value, the low pressure exceeds the maximum stagnation pressure at the edge. Generally, the static pressure changes inside and near the dimple are small, based on the color palette. All the pressure fields and the source streamlines obtained for the considered SST turbulence models are slightly different from each other.
In Figure 4, the fields of the relative x-component friction f x in the projection onto the control plane, coinciding with the channel heated wall, refer to the same characteristics at the appropriate points on the non-dimpled wall of the plane-parallel channel. The air source streamlines (top view) are superimposed on the calculated fields of f x / f x p l . The general features of the presented streamlines are the sharp flow acceleration in the low-pressure narrow regions just behind the dimples, where f x / f x p l reaches a value of 8.
The relative negative friction regions are formed on the windward slopes and at the dimple bottom. The minimum values of f x / f x p l decrease very strongly up to −1.5. The difference in the fields obtained with the use of the Menter 1993 and 2003 SST models corrected according to RLI and in the fields of f x / f x p l calculated with the use of the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM are significant. At the same time, the predictions with the use of the mentioned pairs of the turbulence models appear to be pretty close.
In Figure 5, the fields of the relative Nusselt numbers demonstrate the narrow zones of the maximum heat transfer at the windward edge of the conical dimple. N u / N u p l reaches a value of 2.1. These high-heat zones correlate with the high-gradient regions, i.e., with the maximum values of f x / f x p l .
It should be noted that the low heat load regions with N u / N u p l of the order of 0.3–0.5 appear in the vicinity of the places where tornado-like vortex structures are generated on the left slopes and on the dimple bottom parts adjacent to them, as well as on the leeward slopes.
The above-mentioned separated flow splitting at the dimple bottom into the zones adjacent to the windward and leeward slopes, where the transverse separating streamline is formed, correlates with different heat transfer regions.
The intense air flow on the windward slope and on the bottom part adjacent to it are characteristic of the higher relative Nusselt numbers N u / N u p l than on the leeward slope. It should be noted that the fields of N u / N u p l calculated with the use of the Menter 1993 model and the Menter 2003 model corrected according to RLI are close.
Figure 6 and Figure 7 compare the distributions of the characteristics of the flow and heat transfer over the surface ( p p p l , f x / f x p l , T w , relative Nusselt number N u / N u p l , and transverse strip-integrated averaged Nusselt numbers N u m / N u m p l ) in the middle longitudinal and transverse sections of the control area with the conical dimple calculated with the use of the different SST models. Their analysis allowed us to compare the standard SST models and the models corrected with consideration of the streamline curvature, as well as to determine a quantitative difference between the predictions obtained with their use.
In Figure 6a,b, the pressure drop p p p l is determined in the projections of the dimpled wall points onto the control plane, coinciding with the dimpled heated flat wall of the channel, and at the appropriate points of the plane-parallel channel. As already mentioned in the analysis of the pressure fields on the dimpled wall of the channel, the maximum pressure is seen in the vicinity of the windward edge of the dimple. The static pressure drop in the middle longitudinal section is also observed on the rounded surface of the dimple edge, achieving a value of 0.25. A sharp pressure drop to a value of approx. −0.25 is seen behind the dimple edge at the place of curvature discontinuity. On the windward side, p p p l decreases rather quickly when it moves away from the edge and somewhat increases near the bottom, up to a value of 0.1. While approaching the dimple center on the windward side, the pressure drop takes a value of −0.04 at the flat bottom, i.e., a slight rarefaction is observed in the dimple center. A small negative pressure drop close to a constant one is maintained on the bottom part adjacent to the leeward slope of the dimple, and on the slope itself, it is close to zero. A small static pressure drop up to a value of −0.02 is seen in the vicinity of the leading edge of the dimple. Generally, the pressure drop predictions obtained with the use of different SST models in the middle longitudinal section are close enough. Slight differences can be seen in the center of the dimple and at its trailing edge. The rarefaction and maximum stagnation pressure levels, predicted with the use of the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM, appear to be underestimated. The changes in the predictions obtained by the different SST models for the pressure drop in the middle transverse section of the dimple appear to be quite visible in the background of the not-too-significant extreme values from −0.038 to 0.025. The differences are especially big in the dimple center. They reach a value of 0.02, which is almost 100% of the minimum pressure drop value determined by the Menter 2003 SST model. Thus, we should emphasize a good quantitative agreement of the p p p l predictions, obtained by the Menter 1993 SST model and the Menter 2003 SST model corrected according to RLI, and a significant deviation from the predictions, obtained by the Menter 2003 model incorporated into the widely used software CFX, Fluent, StarCCM+, OpenFOAM, etc., and from those obtained by this model corrected according to SM.
In Figure 6c,d, the relative friction distributions f x / f x p l in the middle longitudinal and transverse sections of the control area of the dimpled wall of the channel demonstrate that an inhomogeneous separated flow is formed inside the dimple. As already noted, the high-gradient narrow zone with multiplying the increasing f x / f x p l arises at the windward rounded edge. The narrow zone with a decreasing relative friction up to a value of 0.7 is also seen behind the above-mentioned high-gradient zone. The separated flow zone in the longitudinal section of the dimple can be divided into three areas: an intense backflow area on the windward side of the dimple with the minimum value of the ratio f x / f x p l equal to −1.4, a backflow area on the bottom flat part adjacent to the leeward side of the dimple with the minimum value of the ratio f x / f x p l equal to −1.6, and a secondary flow area on the bottom flat part adjacent to the leeward side of the dimple with the maximum value of the ratio f x / f x p l equal to 0.3. Attention should be paid to the large maximum absolute values of the relative friction of backflow in the dimple that are caused by developing an intense inclined swirling flow inside the separated flow zone. As in the analysis of the static pressure drop distribution predictions obtained with the use of different SST models, the longitudinal distributions of the relative friction appear to be close, although there is some difference in the minimum values of f x / f x p l for the SST models: −1.1 for the Menter 2003 model and −1.6—for the Menter 1993 model. In the transverse section, the minimum value of f x / f x p l equal to −1.1 is in the vicinity of the junction of the left slope of the dimple and its bottom, i.e., it is strongly shifted from the center. The maximum value of f x / f x p l inside the separated flow zone is seen at the dimple bottom near the right slope. The extreme values of f x / f x p l obtained with the use of the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM are underestimated. A maximum deviation is very significant and is almost 50% of the extreme value. Therefore, for the Menter 2003 SST model, f x / f x p l = 0.7 .
The longitudinal and transverse temperature distributions in the middle sections of the control area of the dimpled wall of the channel (Figure 6e,f) are very illustrative. In front of the dimple, the wall temperature is equal to 1.0055. In the longitudinal section on the leeward slope of the dimple, the temperature rises sharply, reaching a value of 1.011. The wall temperature decreases significantly to 1.008 at the junction of the leeward slope and the flat wall of the bottom. Further moving to the dimple center, the temperature rises rapidly to its maximum value ( T w = 1.0135 ), corresponding to the zero friction point and located in the region of the practically stagnation secondary separated flow of air. On the remaining part of the inner surface of the dimple, the temperature decreases rather quickly and has a small local maximum value at the junction of the bottom and the windward slope, achieving a minimum value of 1.003. It should be mentioned that the minimum value of T w at the place of dimple curvature discontinuity at the switch to the smooth wall of the channel is significantly smaller than the wall temperature in front of the dimple. The differences in the predictions of the extreme wall temperature values inside the dimple obtained with the use of different turbulence models are very noticeable. The value of T w m i n on the leeward slope estimated by the Menter 2003 model is 1.01 and appears to be re-estimated with an error of 0.004. The values of T w m a x on the dimple bottom obtained with the use of all SST models are close, but the maximum values, predicted by the Menter 2003 SST model and the Meter 2003 SST model corrected according to SM, appear to be shifted to the leeward slope. The most distinct differences in the temperature predictions obtained with the use of different SST wall models are seen in the middle transverse section of the dimple. A sharp wall temperature increase on the left slope of the dimple to a value of 1.0115, obtained by the Menter 1993 SST model and the Menter 2003 SST model corrected according to RLI, appears to be somewhat smaller in comparison to the maximum values of T w , obtained by the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM. The greatest differences in the temperature predictions are seen on the bottom flat part adjacent to the left slope of the dimple. The minimum wall temperature value equal to 1.007 is obtained using two SST models: the Menter 1993 SST model and the Menter 2003 SST model corrected according to RLI. According to the Menter 2003 SST model, the temperature in this zone is 1.0085, and according to the SM correction, it is 1.008. The wall temperature behavior on the right slope and on the bottom flat part adjacent to it is almost the same for all the SST models, with the exception of the disagreement of the maximum temperature values on the slope near the rounded edge, where the predictions by the Menter 2003 SST model and its correction according to SM are somewhat overestimated in comparison to T w m a x = 1.009 . It should be noted that there is one local maximum temperature value of approx. 1.009 on the dimple bottom at the junction of the bottom and the right slope.
Generally, the predictions of the considered characteristics obtained with the use of the Menter 1993 model and the Menter 2003 model corrected according to RLI correlate well with each other. At the same time, there are differences in the predictions obtained with the use of the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM that are mainly seen in the field of attaining extreme values. The greatest deviations are characteristic of the Menter 2003 SST model.
Figure 7 illustrates the calculation relations of the relative local Nusselt numbers N u / N u p l and the transverse strip-integrated Nusselt numbers N u m / N u m p l along the x- and transverse z-coordinates in the middle section of the control area of the heated wall with a conical dimple for the different SST turbulence models.
As already discussed, for the pressure drops and the relative x-component friction, the longitudinal distributions of the local and integral heat loads obtained by the standard and modified SST models are rather close, with the exception of the rather narrow zones in the fields attaining extreme values. The maximum values of ( N u / N u p l ) m a x are seen on the windward slope, on the bottom flat part adjacent to the windward slope, and on the windward rounded edge of the dimple. The values of ( N u / N u p l ) m a x are equal to 0.7, 1.1, and 2.05, respectively. The minimum values of ( N u / N u p l ) m i n are observed in the vicinity of the leeward rounded edge, on the dimple bottom in the region of the interacting backward and secondary flows, and at the junction of the bottom and the windward slope. The values of ( N u / N u p l ) m i n . are equal to 0.5, 0.4, and 0.9, respectively. The maximum difference in the predictions of N u / N u p l is of the order of 0.15, which is approx. 20% of the extreme value.
The relative Nusselt numbers N u m / N u m p l | x integrated over the transverse strips of the control area of the heated wall with a dimple achieved a maximum value equal to 1.36 in the strip, corresponding to the vicinity of the windward rounded edge, and a minimum value equal to 0.68 in the strip, passing through the junction of the dimple bottom and the leeward slope. The maximum prediction of N u m / N u m p l obtained by the means of the different SST models is of the order of 0.05.
The dependence of N u / N u p l | z for the different SST models have slight differences in the vicinity of the extreme values and substantial differences in the central region of the dimple bottom. The local minimum values of N u / N u p l , equal to 0.5 and 0.65, are observed on the side left and right slopes of the conical dimple in the vicinity of the rounded edges. The local minimum value equal to 0.66 is at the junction of the bottom and the right slope. The local maximum value of N u / N u p l , equal to 0.9, corresponds to the right slope of the dimple. The greatest differences between the predictions of N u / N u p l are observed in the central region of the dimple at the place, where column-shaped tornado-like structures are generated. The deviation of the prediction obtained by the Menter 2003 SST model from that found by the Menter 1993 SST model is 0.17 when the value of N u / N u p l is equal to 0.86.
The relative Nusselt numbers N u m / N u m p l | z averaged over the longitudinal strips of the control area of the heated wall with a dimple have a local maximum value equal to 1.12 in the strip corresponding to the rounded edge of the left slope of the dimple. The maximum value of N u m / N u m p l is seen in the center of the control area corresponding to the center of the dimple and the gradient flow zone in the wake behind it. This value predicted by means of the Menter 1993 SST model is 1.15 and appears to be much higher than that obtained by means of the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM ( N u m / N u m p l = 1.1 ). Generally, the predictions with the use of the Menter 2003 SST model corrected according RLI agree well with those obtained by the Menter 1993 SST model. At the same time, the predictions of the Menter 2003 SST model are significantly different from them. Unlike the correction according to RLI, the correction according SM has no strong influence on the results.
Figure 8 compares the calculated profiles of the Cartesian velocity components of the separated flow, temperature, and turbulence characteristics in the conical dimple center obtained with the use of different SST models. The predicted values of u ( y ) are typical for separated flows in deep dimples [72,76]. The maximum backflow velocity reaches a value of 0.4. The prediction obtained with the use of the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM somewhat underestimates this value (up to 0.3). In the channel flow core, all models yield the same profiles. The velocity profiles of the secondary swirling flow w ( y ) are different for the different SST models. The maximum value of the positive velocity w near the wall is 0.27, and the minimum value of w in the shear layer is −0.15. In the channel flow core, approximately the same velocity w = 0.05 is maintained, showing that the right–left shear flow is formed under the dimple. The use of the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM underestimates the secondary flow intensity. w m a x decreases up to 0.2 and w m i n increases up to −0.11. Thus, the maximum deviation of w is 30%. The profiles of the vertical velocity v ( y ) show that a rather weak ascending flow with a maximum velocity of 0.055 is formed near the dimple bottom. This flow meets the descending flow at the maximum intensity in the shear layer of the order of v m i n = 0.04 . Above the dimple in the flow core, the vertical velocity component is negative and decreases linearly while approaching the upper wall. The use of the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM underestimates the intensity of the ascending and descending flows. The maximum deviation of the predictions of v does not exceed 0.01, which is 20–25% of the extreme value of v .
The temperature profiles T ( y ) in the separated flow in the dimple are divided into two prediction groups that agree well with each other. The first group uses the predictions obtained with the use of the Menter 1993 SST model and the Menter 2003 SST model corrected according RLI, and the second group uses the predictions obtained with the use of the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM. In the first group, the S-shaped profile with a minimum value of 1.04 and a maximum value of 1.05 is formed inside the dimple. In the second group, the maximum value is practically absent, and in the swirling flow core in the dimple, the temperature is approximately the same at a level of 1.04.
The turbulence energy profiles k ( y ) show that the maximum value of k is realized in the shear layer under the dimple (in the vicinity of y = 0 ). Inside the separated flow zone, the predictions of k are strongly different. In the above-mentioned first group, the swirling flow core is lowly turbulent, and the value of k is equal to 0.017–0.02. It should be emphasized that the maximum value of k , in this case, is 0.04. The second group models predict higher turbulence levels in the separated flow. Therefore, in the swirling flow inside the dimple, the turbulence energy of the fluctuations increases, and the maximum value of k in the shear layer is within 0.044–0.047. It is interesting to mention that, in the channel flow core, the profiles k ( y ) coincide.
As described in references [53,54,63], when the Menter 2003 SST model is used to calculate the separated flows, the effect of artificial pumping of the vortex viscosity occurs in the cores of large-scale vortices. As seen from Figure 8f, the predictions by means of the Menter 2003 SST model are characterized by a bell-shaped distribution when the maximum vortex viscosity value in the swirling flow core is equal to 0.0028. It should be noted that this maximum value hardly exceeds the maximum value of μ t in the shear layer (0.0022). The SM correction somewhat decreases the level of the vortex viscosity pumping (up to 0.0018), but it does not eliminate it. In the shear layer, μ t reaches a value of 0.002. The first group models actually eliminate the false vortex viscosity pumping. The profile μ t ( y ) inside the dimple becomes a double peak with a maximum value of μ t in the swirling flow core within 0.001–0.0015. The second peak corresponds to the shear layer with a maximum value of µt = 0.0018. Outside the dimple in the channel flow core, the predictions of μ t ( y ) by the means of all the SST models coincide.
The comparison of the extreme characteristics of the flow and turbulence in Table 4 complements the analysis of the local characteristic profiles and the pressure distributions over the dimple surface. It is easy to see that the predictions of the characteristics confirm that the considered SST models are divided into two groups. The first group deals with the Menter 1993 SST model and the Menter 2003 SST model corrected according to RLI, and the second deals with the Menter 2003 SST model and the Menter 2003 SST model corrected according to SM. Only the maximum velocities in the dimpled channel appear to be close to each other. The predictions with the use of the second group models significantly decrease the intensity of the backflow and secondary flow and swirling flows, as well as ascending and descending flows. This is largely due to the overestimation of turbulent extreme characteristics.
The same conclusion can be made by comparing the thermal and thermal-hydraulic characteristics of the dimpled channel (Table 2). The predictions of the heat transfer, thermal performance, and relative hydraulic losses of the control area of the channel and the surface bounded by the conical dimple contour appear to be rather close when the first group of the SST turbulence models is used. The Menter 2003 SST model yields the underestimated predictions of the thermal and thermal-hydraulic performances of the control area of the dimpled channel. It should be noted that the relative heat transfer of the channel area appears to be higher than that of an individual dimple. This significantly contributes to the thermal performance of the heat transfer in the wake behind the dimple.

5. Conclusions

A methodical study of a rather intense separated flow and heat transfer in a narrow channel with a conical dimple located on a heated wall with a slope angle of 45°, associated with testing the standard and modified models of shear stress transfer, widely used in the popular packages Fluent-CFX, StarCCM+, etc. of SST turbulence models was carried out. At present, the Menters 2003 SST model (MSST2003) is actively used in engineering practices, in which, as in many other semi-empirical turbulence models, the strain rate tensor module is used. It should be noted that the original Menters 1993 SST model (MSST1993) uses the vorticity modulus to determine the eddy viscosity. When simulating intense separated and swirling flows, it is necessary to correct the Menters 2003 SST model for the effect of streamline curvature, as the model stimulates the false pumping of the vortex viscosity into the cores of large-scale vortices. The MSST2003 correction was presented and described in detail within the framework of the Rodi-Leshziner-Isaev approach, associated with the direct correction of the eddy viscosity (MSST2003RLI) in the calculations of spatial separated flows. The well-known MSST2003 with Smirnov-Menter correction, included in the catalog of models of the CFX package, was also tested, with the introduction of a correction function instead of a constant in the equation for the specific rate of dissipation of the turbulent energy (MSST2003SM). For a sufficiently intense separated flow in a rather deep conical dimple on the wall of a narrow channel with a return flow rate of about 40% of the average mass flow rate in the channel, a significant spread in the uncertainty of the numerical predictions was obtained. It turned out that the standard MSST2003 significantly (by 25%) underestimates the intensity of a separated flow in a conical dimple on a channel wall. At the same time, it is known that the standard MSST1993 quite satisfactorily displays the vortex and separated flows. The predictions for the flow and heat transfer in the conical dimple obtained from it and from the modified MSST2003 RLI turned out to be very close. The predictions for the modified MSST2003SM did not differ much from the predictions using the MSST2003. Thus, the correction of the widely used standard MSST of 2003 within the framework of the Rodi-Leshziner-Isaev (RLI) approach was substantiated, which made it possible to significantly refine the modeling of heat transfer in intense separated and swirling flows.

Author Contributions

Conceptualization, S.I. and A.S.; methodology, S.I.; software, A.S.; validation, N.T.; investigation, L.I.; resources, V.K.; data curation, A.U.; writing—original draft preparation, S.I.; writing—review and editing, N.T.; visualization, A.E. and N.T.; project administration, D.N. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the Ministry of Education and Science of the Russian Federation within the framework of the Program of the World-Class Science Center “Advanced Digital Technologies” (agreement dated 16 November 2020, No. 075-15-2020-903).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hunt, J.C.R. Reynolds Turbulent flows in engineering. J. Fluid Mech. 1975, 70, 414–416. [Google Scholar] [CrossRef]
  2. Spalart, P.R. Strategies for turbulence modelling and simulations. Int. J. Heat Fluid Flow 2000, 21, 252–263. [Google Scholar] [CrossRef]
  3. Leschziner, M.A.; Drikakis, D. Turbulence modelling and turbulent-flow computation in aeronautics. Aeronaut. J. 1968 2002, 106, 349–384. [Google Scholar]
  4. Menter, F.; Kuntz, M.; Langtry, R. Ten Years of Industrial Experience with the SST Turbulence Model. Turbul. Heat Mass Transf. 2003, 4, 625–632. [Google Scholar]
  5. Wilcox, D.C. Solutions Manual: Turbulence Modeling for CFD, 3rd ed.; DCW Industries: La Canada, CA, USA, 2006; ISBN 9781928729099 1928729096. [Google Scholar]
  6. Chen, C.-J.; Jaw, S.-Y. Fundamentals of Turbulence Modeling; Taylor & Francis: London, UK, 1999. [Google Scholar]
  7. Launder, B.E.; Spalding, D.B. The numerical computation of turbulent flows. Comput. Methods Appl. Mech. Eng. 1974, 3, 269–289. [Google Scholar] [CrossRef]
  8. Hellsten, A. Some improvements in Menter’s k-omega SST turbulence model. In 29th AIAA, Fluid Dynamics Conference; Fluid Dynamics and Co-located Conferences; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 1998. [Google Scholar]
  9. Hellsten, A. New Advanced k-w Turbulence Model for High-Lift Aerodynamics. AIAA J. 2005, 43, 1857–1869. [Google Scholar] [CrossRef]
  10. Spalart, P.R.; Garbaruk, A.V. The Predictions of Common Turbulence Models in a Mature Vortex. Flow Turbul. Combust. 2019, 102, 667–677. [Google Scholar] [CrossRef]
  11. Matyushenko, A.A.; Garbaruk, A.V. Non-linear correction for the k-ω SST turbulence model. J. Phys. Conf. Ser. 2017, 929, 012102. [Google Scholar] [CrossRef]
  12. Leschziner, M.A.; Rodi, W. Calculation of Annular and Twin Parallel Jets Using Various Discretization Schemes and Turbulence-Model Variations. J. Fluids Eng. 1981, 103, 352–360. [Google Scholar] [CrossRef]
  13. Leschziner, M.A.; Rodi, W. Computation of strongly swirling axisymmetric free jets. AIAA J. 1984, 22, 1742–1747. [Google Scholar] [CrossRef]
  14. Dong, S.K.; Chung, M.K. Curvature corrections to Reynolds stress model for computation of turbulent recirculating flows. AIAA J. 1992, 30, 2968–2970. [Google Scholar] [CrossRef]
  15. Bradshaw, P. The analogy between streamline curvature and buoyancy in turbulent shear flow. J. Fluid Mech. 1969, 36, 177–191. [Google Scholar] [CrossRef]
  16. Rodi, W.; Scheuerer, G. Calculation of curved shear layers with two-equation turbulence models. Phys. Fluids 1983, 26, 1422–1436. [Google Scholar] [CrossRef]
  17. Launder, B.E.; Priddin, C.H.; Sharma, B.I. The Calculation of Turbulent Boundary Layers on Spinning and Curved Surfaces. J. Fluids Eng. 1977, 99, 231–239. [Google Scholar] [CrossRef]
  18. Knight, D.; Saffman, P. Turbulence model predictions for flows with significant mean streamline curvature. In Proceedings of the 16th Aerospace Sciences Meeting, Huntsville, AL, USA, 16–18 January 1978; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 1978. [Google Scholar]
  19. Wallin, S.; Johansson, A.V. Modelling streamline curvature effects in explicit algebraic Reynolds stress turbulence models. Int. J. Heat Fluid Flow 2002, 23, 721–730. [Google Scholar] [CrossRef]
  20. Shur, M.; Strelets, M.; Travin, A.; Spalart, P. Two numerical studies of trailing vortices. In Proceedings of the 36th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 12–15 January 1998; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 1998. [Google Scholar]
  21. Burr, R.C.; Menter, F.R.; Grotjans, H.; Frühauf, H.H. Development of a Curvature Correction Method for the k-w Turbulence Model. In Computational Fluid Dynamics ’98; Wiley: Hoboken, NJ, USA, 1998; pp. 1105–1110. [Google Scholar]
  22. Zhang, Q.; Yang, Y. A new simpler rotation/curvature correction method for Spalart–Allmaras turbulence model. Chin. J. Aeronaut. 2013, 26, 326–333. [Google Scholar] [CrossRef] [Green Version]
  23. Bobyshev, V.K.; Isaev, S.A. Numerical study of the effects of the current turbulence on the flow along cylinder with a front disk situated. J. Eng. Phys. 1990, 58, 556–572. [Google Scholar]
  24. Isaev, S.A.; Suprun, V.M.; Shul’zhenko, O.A. Numerical and physical simulation of axisymmetric streamlining of a staged cylinder with a low-velocity flow of air. J. Eng. Phys. 1991, 60, 342–347. [Google Scholar] [CrossRef]
  25. Isaev, S.A. Numerical simulation of the axisymmetric low-velocity flow around a cylinder with coaxial disks. J. Eng. Phys. Thermophys. 1995, 68, 16–21. [Google Scholar] [CrossRef]
  26. Isaev, S.A. Numerical investigation of the mechanism of drag reduction of a body with a leading separation zone. J. Eng. Phys. Thermophys. 1995, 68, 786–792. [Google Scholar] [CrossRef]
  27. Isaev, S.A.; Kudryavtsev, N.A.; Sudakov, A.G. Numerical Modeling of a Turbulent Incompressible Viscous Flow Along Bodies of a Curvilinear Shape in the Presence of a Mobile Shield. J. Eng. Phys. Thermophys. 1998, 71, 613–626. [Google Scholar] [CrossRef]
  28. Bobyshev, V.K.; Isaev, S.A. Numerical Investigation of the Effect of Compressibility on the Mechanism of Decreasing the Motion Drag of a Cylinder with Organized Stall Regions in a Turbulent Flow of Viscous Gas. J. Eng. Phys. Thermophys. 1998, 71, 600–606. [Google Scholar] [CrossRef]
  29. Menter, F. Zonal Two Equation k-w Turbulence Models For Aerodynamic Flows. In Proceedings of the 23rd Fluid Dynamics, Plasmadynamics, and Lasers Conference, Orlando, FL, USA, 6–9 July 1993; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 1993. [Google Scholar]
  30. Menter, F.; Ferreira, J.C.; Esch, T.; Konno, B. The SST Turbulence Model with Improved Wall Treatment for Heat Transfer Predictions in Gas Turbines. In Proceedings of the International Gas Turbine Congress, Tokyo, Japan, 2–7 November 2003. [Google Scholar]
  31. Isaev, S.A.; Baranov, P.A.; Kudryavtsev, N.A.; Lysenko, D.A.; Usachov, A.E. Complex analysis of turbulence models, algorithms, and grid structures at the computation of recirculating flow in a cavity by means of VP2/3 and fluent packages. Part 2. Estimation of models adequacy. Thermophys. Aeromech. 2006, 13, 55–65. [Google Scholar] [CrossRef]
  32. Isaev, S.A.; Guvernyuk, S.V.; Zubin, M.A.; Prigorodov, Y.S. Numerical and physical modeling of a low-velocity air flow in a channel with a circular vortex cell. J. Eng. Phys. Thermophys. 2000, 73, 337–344. [Google Scholar] [CrossRef]
  33. Isaev, S.A.; Baranov, P.A.; Guvernyuk, S.V.; Zubin, M.A. Numerical and Physical Modeling of Turbulent Flow in a Divergent Channel with a Vortex Cell. J. Eng. Phys. Thermophys. 2002, 75, 269–276. [Google Scholar] [CrossRef]
  34. Isaev, S.A.; Baranov, P.A.; Kudryavtsev, N.A.; Usachov, A.E. Analysis of vortex heat transfer in a transverse flow past a trench on a plane using multiblock computation technologies and different semi-empirical models of turbulence. J. Eng. Phys. Thermophys. 2004, 77, 1236–1246. [Google Scholar] [CrossRef]
  35. Isaev, S.A.; Leontiev, A.I.; Kudryavtsev, N.A. Numerical simulation of hydrodynamics and heat transfer under conditions of turbulent transverse flow past a “trench” on a plane surface. High Temp. 2005, 43, 89–102. [Google Scholar] [CrossRef]
  36. Dreitser, G.A.; Isaev, S.A.; Lobanov, I.E. Calculation of convective heat transfer in a pipe with periodically arranged surface vortex generators. High Temp. 2005, 43, 214–221. [Google Scholar] [CrossRef]
  37. Isaev, S.A.; Baranov, P.A.; Kudryavtsev, N.A. Numerical Simulation of Heat Transfer under Conditions of Turbulent Separated Flow in Tube Banks. High Temp. 2004, 42, 290–301. [Google Scholar] [CrossRef]
  38. Baranov, P.A.; Isaev, S.A.; Leontiev, A.I.; Mityakov, A.V.; Mityakov, V.Y.; Sapozhnikov, S.Z. Experimental and numerical modelling of vortex heat transfer in turbulent flow past spherical dimple on a plane. Thermophys. Aeromech. 2002, 9, 497–508. [Google Scholar]
  39. Isaev, S.A.; Leont’ev, A.I.; Kiknadze, G.I.; Kudryavtsev, N.A.; Gachechiladze, I.A. Comparative Analysis of the Vortex Heat Exchange in Turbulent Flows around a Spherical Hole and a Two-Dimensional Trench on a Plane Wall. J. Eng. Phys. Thermophys. 2005, 78, 749–761. [Google Scholar] [CrossRef]
  40. Isaev, S.A.; Leont’ev, A.I.; Mityakov, A.V.; Pyshnyi, I.A.; Usachov, A.E. Intensification of Tornado Turbulent Heat Exchange in Asymmetric Holes on a Plane Wall. J. Eng. Phys. Thermophys. 2003, 76, 266–270. [Google Scholar] [CrossRef]
  41. Isaev, S.A.; Leontiev, A.I.; Baranov, P.A.; Pyshnyi, I.A.; Usachov, A.E. Numerical analysis of the vortex intensification of heat transfer in a channel with a set of deep spherical dimples on one of the walls. Dokl. Phys. 2002, 47, 755–757. [Google Scholar] [CrossRef]
  42. Banker, R.; Belenkiy, M.Y.; Gotovskii, M.; Fokin, B.S. Experimental and Computational Investigation of the Hydrodynamics and Heat Transfer in a Flat Channel of Variable Width for Smooth and Intensified Surfaces. Heat Transf. Res. 2004, 35, 10. [Google Scholar] [CrossRef]
  43. Isaev, S.A.; Mikhalev, A.N.; Sudakov, A.G.; Usachev, A.E. Simulation of turbulent flow around a tear-shaped body with a tapered flare. Tech. Phys. 2007, 52, 991–997. [Google Scholar] [CrossRef]
  44. Isaev, S.A.; Baranov, P.A.; Kudryavtsev, N.A.; Zhukova, Y.V. Numerical simulation of unsteady heat exchange at a turbulent flow around a circular cylinder. Part 1. Methodic study. Thermophys. Aeromech. 2005, 12, 27–38. [Google Scholar]
  45. Isaev, S.A.; Baranov, P.A.; Kudryavtsev, N.A.; Lysenko, D.A.; Usachev, A.E. Comparative analysis of the calculation data on an unsteady flow around a circular cylinder obtained using the VP2/3 and fluent packages and the Spalart-Allmaras and Menter turbulence models. J. Eng. Phys. Thermophys. 2005, 78, 1199–1213. [Google Scholar] [CrossRef]
  46. Leontiev, A.I.; Kiselev, N.A.; Vinogradov, Y.A.; Strongin, M.M.; Zditovets, A.G.; Burtsev, S.A. Experimental investigation of heat transfer and drag on surfaces coated with dimples of different shape. Int. J. Therm. Sci. 2017, 118, 152–167. [Google Scholar] [CrossRef]
  47. Leontiev, A.I.; Zditovets, A.G.; Kiselev, N.A.; Vinogradov, Y.A.; Strongin, M.M. Experimental investigation of energy (temperature) separation of a high-velocity air flow in a cylindrical channel with a permeable wall. Exp. Therm. Fluid Sci. 2019, 105, 206–215. [Google Scholar] [CrossRef]
  48. Zhukauskas, A.A.; Ulinskas, R. Heat Transfer in Tube Banks in Crossflow; Hemisphere Pub. Corp.: New York, NY, USA, 1988; ISBN 0-89116-685-8. [Google Scholar]
  49. Mityakov, V.Y.; Mityakov, A.V.; Sapozhnikov, S.Z.; Isaev, S.A. Local heat fluxes on the surfaces of dimples, ditches, and cavities. Therm. Eng. 2007, 54, 200–203. [Google Scholar] [CrossRef]
  50. Kiknadze, G.I.; Gachechiladze, I.A.; Gorodkov, A.Y. Self-Organization of Tornado-Like Jets in Flows of Gases and Liquids and the Technologies Utilizing This Phenomenon. In Proceedings of the HT2009, San Francisco, CA, USA, 19–23 July 2009; Volume 3, pp. 547–560. [Google Scholar]
  51. Chyu, M.K.; Yu, Y.; Ding, H. Heat Transfer Enhancement in Rectangular Channels with Concavities. J. Enhanc. Heat Transf. 1999, 6, 429–439. [Google Scholar] [CrossRef]
  52. Isaev, S.A.; Zhdanov, V.L.; Niemann, H.-J. Numerical study of the bleeding effect on the aerodynamic characteristics of a circular cylinder. J. Wind Eng. Ind. Aerodyn. 2002, 90, 1217–1226. [Google Scholar] [CrossRef]
  53. Isaev, S.A.; Baranov, P.A.; Zhukova, Y.V.; Usachov, A.E.; Kharchenko, V.B. Correction of the Shear-Stress-Transfer Model with Account of the Curvature of Streamlines in Calculating Separated Flows of an Incompressible Viscous Fluid. J. Eng. Phys. Thermophys. 2014, 87, 1002–1015. [Google Scholar] [CrossRef]
  54. Isaev, S.A. Experience of application of SST-model-2003 with correction on streamline curvature according to Rodi-Leshziner-Isaev approach for (U)RANS calculations of separated and vortex sub-and supersonic flows. AIP Conf. Proc. 2018, 2027, 020015. [Google Scholar]
  55. Spalart, P.R.; Shur, M. On the sensitization of turbulence models to rotation and curvature. Aerosp. Sci. Technol. 1997, 1, 297–302. [Google Scholar] [CrossRef]
  56. Shur, M.L.; Strelets, M.K.; Travin, A.K.; Spalart, P.R. Turbulence Modeling in Rotating and Curved Channels: Assessing the Spalart-Shur Correction. AIAA J. 2000, 38, 784–792. [Google Scholar] [CrossRef]
  57. Smirnov, P.E.; Menter, F.R. Sensitization of the SST Turbulence Model to Rotation and Curvature by Applying the Spalart–Shur Correction Term. J. Turbomach. 2009, 131, 041010. [Google Scholar] [CrossRef]
  58. Baranov, P.A.; Guvernyuk, S.V.; Zubin, M.A.; Isaev, S.A.; Usachov, A.E. Application of various models of turbulence for calculation of incompressible internal flows. TsAGI Sci. J. 2017, 48, 31–42. [Google Scholar] [CrossRef]
  59. Isaev, S.A.; Baranov, P.A.; Usachov, A.E.; Zhukova, Y.V.; Vysotskaya, A.A.; Malyshkin, D.A. Simulation of the Turbulent air Flow Over a Circular Cavity with a Variable Opening Angle in an U-Shaped Channel. J. Eng. Phys. Thermophys. 2015, 88, 902–917. [Google Scholar] [CrossRef]
  60. Isaev, S.A.; Kalinin, E.I.; Tereshkin, A.A.; Usachov, A.E. Modeling a decrease in hydraulic losses during turbulent flow in a U-bend channel with a circular cavern with a large opening angle. Tech. Phys. Lett. 2015, 41, 298–300. [Google Scholar] [CrossRef]
  61. Baranov, P.A.; Guvernyuk, S.V.; Isaev, S.A.; Soudakov, A.G.; Usachov, A.E. Simulation of periodic vortical structures in the airfoil wake. TsAGI Sci. J. 2014, 45, 273–292. [Google Scholar] [CrossRef]
  62. Isaev, S.; Baranov, P.; Popov, I.; Sudakov, A.; Usachov, A. Improvement of aerodynamic characteristics of a thick airfoil with a vortex cell in sub- and transonic flow. Acta Astronaut. 2017, 132, 204–220. [Google Scholar] [CrossRef]
  63. Isaev, S.A.; Baranov, P.A.; Sudakov, A.G.; Popov, I.A. Verification of the standard model of shear stress transport and its modified version that takes into account the streamline curvature and estimation of the applicability of the Menter combined boundary conditions in calculating the ultralow profile drag for an optimally configured cylinder–coaxial disk arrangement. Tech. Phys. 2016, 61, 1152–1161. [Google Scholar]
  64. Isaev, S.A.; Lipnitskii, Y.M.; Mikhalev, A.N.; Panasenko, A.V.; Usachov, A.E. Simulation of the supersonic turbulent flow around a cylinder with coaxial disks. J. Eng. Phys. Thermophys. 2011, 84, 827. [Google Scholar] [CrossRef]
  65. Isaev, S.A.; Baranov, P.A.; Mikhalev, A.N.; Sudakov, A.G. Modeling the effect of head drag reduction for a cylinder with a protruding disk at high mach numbers. Tech. Phys. Lett. 2014, 40, 996–999. [Google Scholar] [CrossRef]
  66. Isaev, S.A.; Lipnitskii, Y.M.; Baranov, P.A.; Panasenko, A.V.; Usachov, A.E. Simulation of a turbulent supersonic underexpanded jet flowing into a submerged space with the help of a shear stress transfer model. J. Eng. Phys. Thermophys. 2012, 85, 1357–1371. [Google Scholar] [CrossRef]
  67. Isaev, S.A.; Miau, J.-J.; Sudakov, A.G.; Usachov, A.E. Analysis of extremal lift behavior of a semicircular airfoil in a turbulent airflow at a near-zero angle of attack. Tech. Phys. Lett. 2015, 41, 737–739. [Google Scholar] [CrossRef]
  68. Isaev, S.A.; Baranov, P.A.; Zhukova, Y.V.; Kalinin, E.I.; Miau, J.J. Verification of the Shear-Stress Transfer Model and its Modifications in the Calculation of a Turbulent Flow Around a Semicircular Airfoil with a Zero Angle of Attack. J. Eng. Phys. Thermophys. 2016, 89, 73–89. [Google Scholar] [CrossRef]
  69. Isaev, S.; Baranov, P.; Popov, I.; Sudakov, A.; Usachov, A.; Guvernyuk, S.; Sinyavin, A.; Chylunin, A.; Mazo, A.; Kalinin, E. Ensuring safe descend of reusable rocket stages – Numerical simulation and experiments on subsonic turbulent air flow around a semi-circular cylinder at zero angle of attack and moderate Reynolds number. Acta Astronaut. 2018, 150, 117–136. [Google Scholar] [CrossRef]
  70. Isaev, S.; Baranov, P.; Popov, I.; Sudakov, A.; Usachov, A.; Guvernyuk, S.; Sinyavin, A.; Chulyunin, A.; Mazo, A.; Demidov, D.; et al. Numerical simulation and experiments on turbulent air flow around the semi-circular profile at zero angle of attack and moderate Reynolds number. Comput. Fluids 2019, 188, 1–17. [Google Scholar] [CrossRef] [Green Version]
  71. Isaev, S.; Leontiev, A.; Chudnovsky, Y.; Nikushchenko, D.; Popov, I.; Sudakov, A. Simulation of Vortex Heat Transfer Enhancement in the Turbulent Water Flow in the Narrow Plane-Parallel Channel with an Inclined Oval-trench Dimple of Fixed Depth and Spot Area. Energies 2019, 12, 1296. [Google Scholar] [CrossRef] [Green Version]
  72. Isaev, S.A.; Kornev, N.V.; Leontiev, A.I.; Hassel, E. Influence of the Reynolds number and the spherical dimple depth on turbulent heat transfer and hydraulic loss in a narrow channel. Int. J. Heat Mass Transf. 2010, 53, 178–197. [Google Scholar] [CrossRef]
  73. Isaev, S.A.; Leontiev, A.I. Problems of simulating tornado-like heat transfer in turbulent flow past a dimpled relief on a narrow channel wall. J. Eng. Phys. Thermophys. 2010, 83, 783–793. [Google Scholar] [CrossRef]
  74. Isaev, S.A.; Chulyunin, A.Y.; Nikushchenko, D.V.; Sudakov, A.G.; Usachov, A.E. Analysis of the anomalous intensification of a separate flow and heat transfer in a stabilized section of a narrow channel with single-row, inclined, oval-trench dimples with the use of various grids and turbulence models. High Temp. 2021, 59, 106–114. [Google Scholar]
  75. Isaev, S.A.; Baranov, P.A.; Zhukova, Y.V.; Tereshkin, A.A.; Usachov, A.E. Simulation of the Wind Effect on an Ensemble of High-Rise Buildings by means of Multiblock Computational Technologies. J. Eng. Phys. Thermophys. 2014, 87, 112–123. [Google Scholar] [CrossRef]
  76. Isaev, S.A.; Schelchkov, A.V.; Leontiev, A.I.; Gortyshov, Y.F.; Baranov, P.A.; Popov, I.A. Vortex heat transfer enhancement in the narrow plane-parallel channel with the oval-trench dimple of fixed depth and spot area. Int. J. Heat Mass Transf. 2017, 109, 40–62. [Google Scholar] [CrossRef]
  77. Savelsberg, R.; Castro, I. Vortex flows in open cylindrical-section cavities. Exp. Fluids 2009, 46, 485–497. [Google Scholar] [CrossRef]
  78. Saqr, K.M.; Aly, H.S.; Kassem, H.I.; Sies, M.M.; Wahid, M.A. Computations of shear driven vortex flow in a cylindrical cavity using a modified k-ε turbulence model. Int. Commun. Heat Mass Transf. 2010, 37, 1072–1077. [Google Scholar] [CrossRef]
  79. Isaev, S.A.; Guvernyuk, S.V.; Zubin, M.A.; Baranov, P.A.; Ermakov, A.M. Numerical and Physical Simulation of the Low-Velocity Air Flow in a Diffuser with a Circular Cavity in the Case of Suction of the Air from the Central Cylindrical Body Positioned in the Cavity. J. Eng. Phys. Thermophys. 2015, 88, 186–199. [Google Scholar] [CrossRef]
  80. Belov, I.A.; Isaev, S.A.; Konovalov, V.N.; Mitin, A.Y. Application of the concept of an ideal fluid for calculating the separation flow past blunt bodies with account for the turbulent shear layer at the boundary of the separation region. Pisma Zh Tekh Fiz 1984, 10, 1217–1220. [Google Scholar]
  81. Belov, I.A.; Isaev, S.A.; Konovalov, V.N.; Mitin, A.Y. Estimation of the wave resistance of bodies of revolution with a front separation zone in a supersonic flow. Izv Akad Nauk SSSR 1985, 4, 47–51. [Google Scholar]
  82. Belov, I.A.; Isaev, S.A.; Konovalov, V.N.; Mitin, A.Y. Modeling of large-scale vortex structures in a turbulent supersonic flow past a blunt body. Izv Akad Nauk SSSR 1987, 15, 101–107. [Google Scholar]
  83. Klamo, J.T.; Leonard, A.; Roshko, A. The effects of damping on the amplitude and frequency response of a freely vibrating cylinder in cross-flow. Bluff Body Wakes Vortex Induc. Vib. BBVIV-4 2006, 22, 845–856. [Google Scholar] [CrossRef]
  84. Zapryagaev, V.I.; Rybak, S.P.; Kavun, I.N.; Kiselev, N.P.; Gubanov, D.A.; Sobolev, A.V. Research of supersonic jet flows applicable to problems of aerospace technology. AIP Conf. Proc. 2021, 2351, 020007. [Google Scholar]
  85. Sluchanovskaya, Z.P. Pressure distributions on the surfaces of a rectangular, a trihedral, and a semicircular cylinders and their aerodynamic coefficients. Tr. Inst. Mekhaniki 1973, 52–60. [Google Scholar]
  86. Terekhov, V.I.; Kalinina, S.V.; Mshvidobadze, Y.M. Heat Transfer Coefficient and Aerodynamic Resistance on a Surface with a Single Dimple. J. Enhanc. Heat Transf. 1997, 4, 131–145. [Google Scholar] [CrossRef]
  87. Terekhov, V.I.; Kalinina, S.V.; Mshvidobadze, Y.M. Pressure field and resistance of a single cavity with sharp and rounded edges. J. Appl. Mech. Tech. Phys. 1993, 34, 337–338. [Google Scholar] [CrossRef]
  88. Johnson, D.A.; King, L.S. A mathematically simple turbulence closure model for attached and separated turbulent boundary layers. AIAA J. 1985, 23, 1684–1692. [Google Scholar] [CrossRef]
  89. Wilcox, D.C. Reassessment of the scale-determining equation for advanced turbulence models. AIAA J. 1988, 26, 1299–1310. [Google Scholar] [CrossRef]
  90. Isaev, S.A.; Baranov, P.A.; Usachov, A.E. Multiblock Computational Technologies in the VP2/3 Package on Aerothermo-Dynamics; LAP LAMBERT Academic Publishing: Saarbrucken, Germany, 2013; 316p. (In Russian) [Google Scholar]
  91. Ferziger, J.H.; Perić, M. Computational Methods for Fluid Dynamics, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 1999. [Google Scholar]
  92. Patankar, S.V. Numerical Heat Transfer and Fluid Flow, 1st ed.; Hemisphere Pub. Corp.: New York, NY, USA, 1980; ISBN 0-07-048740-5. [Google Scholar]
  93. Zheng, Y.; Liou, M.-S. A novel approach of three-dimensional hybrid grid methodology: Part 1. Grid generation. Comput. Methods Appl. Mech. Eng. 2003, 192, 4147–4171. [Google Scholar] [CrossRef] [Green Version]
  94. Leonard, B.P. A stable and accurate convective modelling procedure based on quadratic upstream interpolation. Comput. Methods Appl. Mech. Eng. 1979, 19, 59–98. [Google Scholar] [CrossRef]
  95. van Leer, B. Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method. J. Comput. Phys. 1979, 32, 101–136. [Google Scholar] [CrossRef]
  96. Rhie, C.M.; Chow, W.L. Numerical study of the turbulent flow past an airfoil with trailing edge separation. AIAA J. 1983, 21, 1525–1532. [Google Scholar] [CrossRef]
  97. Pascau, A.; Garcia-Polanco, N. Consistency of simplec scheme in collocated grids. In Proceedings of the V European Conference on Computational Fluid Dynamics, Lisbon, Portugal, 14–17 June 2010. [Google Scholar]
  98. Saad, Y. Iterative Methods for Sparse Linear Systems, 2nd ed.; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2003; ISBN 0-89871-534-2. [Google Scholar]
  99. Demidov, D. AMGCL—A C++ library for efficient solution of large sparse linear systems. Softw. Impacts 2020, 6, 100037. [Google Scholar] [CrossRef]
  100. Isaev, S.A.; Leont’ev, A.I.; Baranov, P.A.; Pyshnyi, I.A. Numerical Analysis of the Influence of the Depth of a Spherical Hole on a Plane Wall on Turbulent Heat Exchange. J. Eng. Phys. Thermophys. 2003, 76, 61–69. [Google Scholar] [CrossRef]
  101. Isaev, S.A.; Schelchkov, A.V.; Leontiev, A.I.; Baranov, P.A.; Gulcova, M.E. Numerical simulation of the turbulent air flow in the narrow channel with a heated wall and a spherical dimple placed on it for vortex heat transfer enhancement depending on the dimple depth. Int. J. Heat Mass Transf. 2016, 94, 426–448. [Google Scholar] [CrossRef]
  102. Isaev, S.A.; Leontiev, A.I.; Baranov, P.A.; Popov, I.A.; Shchelchkov, A.V.; Gabdrakhmanov, I.P. Numerical Simulation of the Intensification of the Heat Exchange in a Plane-Parallel Channel with a Cylindrical Shallow Dimple on the Heated Wall. J. Eng. Phys. Thermophys. 2016, 89, 1186–1201. [Google Scholar] [CrossRef]
Figure 1. Area of the plane-parallel channel with a conical dimple on the heated wall (a) with the the Cartesian x,y,z coordinate system and multiblock grids in the vicinity of the dimple (b). 1—channel grid (CG); 2—refined Cartesian grid (wake grid—WG) covering the dimple, 3—O-shaped grid Ring inside the dimple, and 4—grid Rec—axial “patch”. The upper wall is removed.
Figure 1. Area of the plane-parallel channel with a conical dimple on the heated wall (a) with the the Cartesian x,y,z coordinate system and multiblock grids in the vicinity of the dimple (b). 1—channel grid (CG); 2—refined Cartesian grid (wake grid—WG) covering the dimple, 3—O-shaped grid Ring inside the dimple, and 4—grid Rec—axial “patch”. The upper wall is removed.
Energies 14 05038 g001
Figure 2. Air source streamlines on the wall surfaces of the channel with a conical dimple (a) and the vortex structure (b).
Figure 2. Air source streamlines on the wall surfaces of the channel with a conical dimple (a) and the vortex structure (b).
Energies 14 05038 g002
Figure 3. Calculated pressure fields and the source streamlines in the control area of the channel wall with the dimples for the considered standard SST models and the modified Menter SST models: (a) Menter 1993 SST model, (b) Menter 2003 SST model, (c) Menter SST model corrected according to RLI, and (d) Menter SST model corrected according to SM.
Figure 3. Calculated pressure fields and the source streamlines in the control area of the channel wall with the dimples for the considered standard SST models and the modified Menter SST models: (a) Menter 1993 SST model, (b) Menter 2003 SST model, (c) Menter SST model corrected according to RLI, and (d) Menter SST model corrected according to SM.
Energies 14 05038 g003
Figure 4. Calculated fields of the relative friction f x / f x p l and the source streamlines in the control area of the channel walls with the dimples of different standard SST models and the modified Menter SST models: (a) Menter 1993 SST model, (b) Menter 2003 SST model, (c) Menter SST model corrected according to RLI, and (d) Menter SST model corrected according to SM.
Figure 4. Calculated fields of the relative friction f x / f x p l and the source streamlines in the control area of the channel walls with the dimples of different standard SST models and the modified Menter SST models: (a) Menter 1993 SST model, (b) Menter 2003 SST model, (c) Menter SST model corrected according to RLI, and (d) Menter SST model corrected according to SM.
Energies 14 05038 g004
Figure 5. Calculated fields of the relative Nusselt numbers N u / N u p l and the source streamlines in the control area of the channel wall with the dimples for the different standard SST models and the modified Menter SST models: (a) Menter 1993 SST model, (b) Menter 2003 SST model, (c) Menter SST model corrected according to RLI, and (d) Menter SST model corrected according to SM.
Figure 5. Calculated fields of the relative Nusselt numbers N u / N u p l and the source streamlines in the control area of the channel wall with the dimples for the different standard SST models and the modified Menter SST models: (a) Menter 1993 SST model, (b) Menter 2003 SST model, (c) Menter SST model corrected according to RLI, and (d) Menter SST model corrected according to SM.
Energies 14 05038 g005
Figure 6. Distributions of the static pressure drop p p p l (a,b), the relative friction f x / f x p l (c,d), and the temperature T w (e,f) in the middle longitudinal (a,c,e) and transverse (b,d,f) sections of the dimpled wall of the channel for the different standard SST models and the modified Menter SST models. Red line—Menter 1993 SST model, Blue line—Menter 2003 SST model, Green line—Menter SST model corrected according to RLI, and Purple line—Menter SST model corrected according to SM.
Figure 6. Distributions of the static pressure drop p p p l (a,b), the relative friction f x / f x p l (c,d), and the temperature T w (e,f) in the middle longitudinal (a,c,e) and transverse (b,d,f) sections of the dimpled wall of the channel for the different standard SST models and the modified Menter SST models. Red line—Menter 1993 SST model, Blue line—Menter 2003 SST model, Green line—Menter SST model corrected according to RLI, and Purple line—Menter SST model corrected according to SM.
Energies 14 05038 g006aEnergies 14 05038 g006b
Figure 7. Distributions of the relative Nusselt numbers N u / N u p l (a,b) and the transverse strip-averaged relative Nusselt numbers N u m / N u m p l (c,d) in the middle longitudinal (a,c) and transverse (b,d) sections of the dimpled walls for the different standard SST models and the modified Menter SST models. Red line—Menter 1993 SST model, Blue line—Menter 2003 SST model, Green line—Menter SST model corrected according to RLI, and Purple line—Menter SST model corrected according to SM.
Figure 7. Distributions of the relative Nusselt numbers N u / N u p l (a,b) and the transverse strip-averaged relative Nusselt numbers N u m / N u m p l (c,d) in the middle longitudinal (a,c) and transverse (b,d) sections of the dimpled walls for the different standard SST models and the modified Menter SST models. Red line—Menter 1993 SST model, Blue line—Menter 2003 SST model, Green line—Menter SST model corrected according to RLI, and Purple line—Menter SST model corrected according to SM.
Energies 14 05038 g007aEnergies 14 05038 g007b
Figure 8. Profiles of the longitudinal (a), transverse (b), and vertical (c) Cartesian velocity components; temperature (d); turbulence energy (e); and vortex viscosity (f) at the central point of the conical dimples for the different standard SST models and the modified Menter SST models. Red line—Menter 1993 SST model, Blue line—Menter 2003 SST model, Green line—Menter SST model corrected according to RLI, and Purple line—Menter SST model corrected according to SM.
Figure 8. Profiles of the longitudinal (a), transverse (b), and vertical (c) Cartesian velocity components; temperature (d); turbulence energy (e); and vortex viscosity (f) at the central point of the conical dimples for the different standard SST models and the modified Menter SST models. Red line—Menter 1993 SST model, Blue line—Menter 2003 SST model, Green line—Menter SST model corrected according to RLI, and Purple line—Menter SST model corrected according to SM.
Energies 14 05038 g008aEnergies 14 05038 g008b
Table 1. Thermal and thermal-hydraulic performances of the control area and the hydraulic losses of different density grids.
Table 1. Thermal and thermal-hydraulic performances of the control area and the hydraulic losses of different density grids.
Grid N u m m d 10 3 ζ d T H P d N u m m 10 3 ζ T H P
(S)36.719.04090.820/1.326
0.618
44.867.83571.026/1.150
0.892
(A)37.478.77460.837/1.287
0.650
45.067.83281.031/1.150
0.897
(B)37.489.04050.837/1.326
0.631
45.067.81801.031/1.148
0.898
(C)37.529.02170.838/1.323
0.633
45.057.79901.031/1.145
0.900
Table 2. Thermal and thermal-hydraulic performances of the control area and the hydraulic losses for the standard and modified SST models.
Table 2. Thermal and thermal-hydraulic performances of the control area and the hydraulic losses for the standard and modified SST models.
MSST N u m m d 10 3 ζ d T H P d N u m m 10 3 ζ T H P
199344.766.8198143.716.81081
200338.419.39630.858/1.378
0.622
45.937.97501.051/1.171
0.900
RLI34.498.90450.771/1.306
0.657
44.198.09001.011/1.188
0.851
SM37.489.04050.837/1.326
0.631
45.067.81801.031/1.148
0.898
Table 3. Extreme flow parameter values and the turbulence characteristics for different density grids.
Table 3. Extreme flow parameter values and the turbulence characteristics for different density grids.
Grid p m a x p m i n u m a x u m i n v m a x v m i n w m a x w m i n 10 k m a x 10 μ t m a x
(S)0.146−0.4031.288−0.4090.437−0.2790.252−0.3470.7230.300
(A)0.159−0.4251.290−0.4350.467−0.2890.261−0.3650.7180.272
(B)0.159−0.4281.287−0.4340.467−0.2890.264−0.3650.7200.272
(C)0.161−0.4321.288−0.4350.470−0.2900.267−0.3670.7190.271
Table 4. Extreme flow parameter values and the turbulence characteristics for the standard and modified SST models.
Table 4. Extreme flow parameter values and the turbulence characteristics for the standard and modified SST models.
MSST p m a x p m i n u m a x u m i n v m a x v m i n w m a x w m i n 10 k m a x 10 μ t m a x
19930.169−0.4411.288−0.4330.465−0.2820.363−0.2560.9040.313
20030.125−0.3751.283−0.3590.417−0.2640.290−0.2290.7700.377
RLI0.159−0.4281.287−0.4340.467−0.2890.365−0.2640.7200.272
SM0.137−0.3921.285−0.3780.437−0.2670.315−0.2350.8240.354
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Isaev, S.; Nikushchenko, D.; Sudakov, A.; Tryaskin, N.; Egorova, A.; Iunakov, L.; Usachov, A.; Kharchenko, V. Standard and Modified SST Models with the Consideration of the Streamline Curvature for Separated Flow Calculation in a Narrow Channel with a Conical Dimple on the Heated Wall. Energies 2021, 14, 5038. https://doi.org/10.3390/en14165038

AMA Style

Isaev S, Nikushchenko D, Sudakov A, Tryaskin N, Egorova A, Iunakov L, Usachov A, Kharchenko V. Standard and Modified SST Models with the Consideration of the Streamline Curvature for Separated Flow Calculation in a Narrow Channel with a Conical Dimple on the Heated Wall. Energies. 2021; 14(16):5038. https://doi.org/10.3390/en14165038

Chicago/Turabian Style

Isaev, Sergey, Dmitry Nikushchenko, Alexandr Sudakov, Nikita Tryaskin, Ann Egorova, Leonid Iunakov, Alexandr Usachov, and Valery Kharchenko. 2021. "Standard and Modified SST Models with the Consideration of the Streamline Curvature for Separated Flow Calculation in a Narrow Channel with a Conical Dimple on the Heated Wall" Energies 14, no. 16: 5038. https://doi.org/10.3390/en14165038

APA Style

Isaev, S., Nikushchenko, D., Sudakov, A., Tryaskin, N., Egorova, A., Iunakov, L., Usachov, A., & Kharchenko, V. (2021). Standard and Modified SST Models with the Consideration of the Streamline Curvature for Separated Flow Calculation in a Narrow Channel with a Conical Dimple on the Heated Wall. Energies, 14(16), 5038. https://doi.org/10.3390/en14165038

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