Next Article in Journal
Instability and Transition of a Boundary Layer over a Backward-Facing Step
Next Article in Special Issue
Numerical Investigation of Orifice Nearfield Flow Development in Oleo-Pneumatic Shock Absorbers
Previous Article in Journal
Characterization of Effective Diffusion within Viscoelastic Fluids with Elastic Instabilities
Previous Article in Special Issue
Combined Stereoscopic Particle Image Velocimetry Measurements in a Single Plane for an Impinging Jet around a Thin Control Rod
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Compressible Turbulence Model for Pressure—Strain

Laboratory of Aeronautical Structure and Control of System Flows, SACES (UR17DN01), School of Military Aviation, Borj el Amri 1142, Tunisia
*
Author to whom correspondence should be addressed.
Fluids 2022, 7(1), 34; https://doi.org/10.3390/fluids7010034
Submission received: 31 July 2021 / Revised: 18 August 2021 / Accepted: 21 August 2021 / Published: 14 January 2022
(This article belongs to the Special Issue Turbulent Flow)

Abstract

:
This work focuses on the performance and validation of compressible turbulence models for the pressure-strain correlation. Considering the Launder Reece and Rodi (LRR) incompressible model for the pressure-strain correlation, Adumitroaie et al., Huang et al., and Marzougui et al., used different modeling approaches to develop turbulence models, taking into account compressibility effects for this term. Two numerical coefficients are dependent on the turbulent Mach number, and all of the remaining coefficients conserve the same values as in the original LRR model. The models do not correctly predict the compressible turbulence at a high-speed shear flow. So, the revision of these models is the major aim of this study. In the present work, the compressible model for the pressure-strain correlation developed by Khlifi−Lili, involving the turbulent Mach number, the gradient, and the convective Mach numbers, is used to modify the linear mean shear strain and the slow terms of the previous models. The models are tested in two compressible turbulent flows: homogeneous shear flow and the newly developed plane mixing layers. The predicted results of the proposed modifications of the Adumitroaie et al., Huang et al., and Marzougui et al., models and of its universal versions are compared with direct numerical simulation (DNS) and experiment data. The results show that the important parameters of compressibility in homogeneous shear flow and in the mixing layers are well predicted by the proposal models.

1. Introduction

Compressible turbulence modelling is an essential element for many industrial problems. Abetter understanding of the compressibility effects is highly relevant in the design of aerospace, supersonic, and hypersonic flights; combustion field; and other engineering problems. Firstly, attention is paid to study the compressibility effects on homogeneous shear flow(the mean velocity is (Sx2, 0, 0), S = cte), which is a useful problem because this flow summarizes some of the important compressibility properties in a simplified setting. In addition, this flow has excessively been used in calibration and evaluating turbulence models. In this context, the DNS results of Blaisdell et al. [1] and Sarkar et al. [2] show that there are significant turbulence changes when the compressibility increases as the turbulent kinetic energy growth rate decrease with the increasing turbulent Mach number, M t = 2 K / a ¯ , a ¯ = γ R T ˜ . Both studies show that the dilatational terms, π d and ε c , on the R.H.S of the turbulent kinetic energy equation were found to be much smaller compared with the control compressibility effects. Sarkar [3], Simone et al. [4], and Hamba [5] developed DNS results and reached a similar conclusion concerning the role of the dilatational terms. It has been found in their DNS results that the structural compressibility effects affect the pressure field and then the pressure-strain, which is recognized as the main factor responsible for the strong changes in the magnitude of the Reynolds stress anisotropies, and thereafter the reduced trend of the growth rate of the turbulent kinetic energy when the compressibility increases. Similar conclusions are confirmed by the DNS results of Vreman et al. [6]; the experimental data of Goebel et al. [7] and Samimy et al. [8]; and, more recently, the DNS data of Pantano et al. [9] and Foysi et al. [10], in which it is reported that the compressibility effects on the pressure-strain are the main cause of the changes in the planar compressible mixing layers. Thus, we argue that the pressure-strain correlation is one of the mean terms contributing to the reduced growth rate and the changes of the Reynolds stress arising from the compressibility effects. Modeling the turbulent pressure-strain correlation occurs mainly at a high speed. Three independent compressible pressure-strain models, by Adumitroaie et al. [11], Huang et al. [12], and Marzougui et al. [13], are considered in this study. These models are derived by considering different variable density extensions of the Launder et al. LRR model [14], which account for the compressibility effects by using the turbulent Mach number.It has been shown that these models may be able to reproduce low and moderate compressibility effects. However, when the compressibility effects are more significant, the models do not correctly predict the decrease in the spreading rate of the mixing layers, as it is observed in [7,8,9,10], nor the reduction in the growth rate of turbulent kinetic energy [3,4,5]. The deficiencies of this closure are probably because LRR-compressibility correction models seem to be insufficient to induce important variation in calculations in accordance with the anisotropy turbulence strong changes when the compressibility is higher. Thus, one can see that in the models [11,12,13], the two coefficients that affect the linear term in relation to the Reynolds stress anisotropy and the mean stain are modified, which then become a function of the turbulent Mach number. All of the other remaining coefficients that affect the mean shear and the return to isotropy terms are conserved as in the LRR model [14], without any modification. However, modification of these models taking into account structural compressibility effects is needed for the pressure-strain correlation coefficient models. The present work focuses on this major issue. For this, more attention is paid to the results and the analysis of the DNS [3,4], in which some important compressibility discrepancies for homogeneous turbulent shear flows can be sound. It has been argued that the gradient Mach number, M g = S l / a ¯ ( S = U ˜ i , j U ˜ i , j , where S and l are the mean shear constant and integral length scale, respectively),is an appropriate parameter in addition to the turbulent Mach number for studying the structural compressibility effects and must be added to M t in the compressible modelling concept. Similar recommendations have been suggested in different experimental [7,8] and DNS [9] data, which identify the convective Mach number, M c = ( U 1 U 2 ) / ( a 1 + a 2 ) , where U 1 ,   a 1 and U 2 ,   a 2 denoting the velocity and the speed sound in the high speed stream and in the low speed stream, respectively, as an appropriate parameter in order to study the compressibility effects on the planar mixing layer as in [7,8,15].
In the present study, a revision ofthe models by Adumitroaie et al. [11], Huang et al. [12], and Marzougui et al. [13] for pressure-strain is considered, making the model coefficients as a function of M t ,   M g , and M c . The proposed models are tested in different compressible turbulent homogeneous shear flow and compressible mixing layers cases.

2. Basic Equations

In general, compressible turbulent flow is described by continuity, Navier–Stokes, energy, and state equations. It is well known that the basic equations of the mean quantities used in describing turbulence closure schemes are essentially those using the Favre average. These equations are formally similar to those governed by incompressible turbulent flows. Obviously, this technique gives reason to the extension of the incompressible turbulent models to study compressible turbulent flows. This is one of the essential advantages provided by the density weighting technique for modeling compressible turbulence. For this study, the Favre averaged continuity, momentum, and specific internal energy equations are respectively written as follows [11,16]:
t ρ ¯ + x i ( ρ ¯ U ˜ i )   =   0 ,
t ( ρ ¯ U ˜ i ) + x i ( ρ ¯ U ˜ i   U ˜ j ) = x j ( τ ˜ i j + τ i j ¯ x j ρ u i u j ¯ p ¯ δ i j ) ,
t ρ ¯ c ¯ v T ˜ + x j ρ ¯ c ¯ v T ˜ U ˜ j = ϕ e + π d x j ( c v ρ u j T ¯ )   .
where ϕ e = p ¯ x i ( U ˜ i + u i ¯   ) + x i ( κ x i T ¯ ) + τ i j u i , j ¯   ,   τ ˜ i j = 2 μ ¯ S ˜ i j 2 3 μ ¯ U ˜ k , k δ i j

3. Turbulence Models

The turbulence models used in this work are closely related to the standard Reynolds stress model from which the Favre averaged Reynolds stress, ρ u i u j ¯ / ρ ¯ , is described by the following equations [11,16]:
t ( ρ ¯ R i j ) + x m ( ρ ¯ U ˜ m R i j ) = Pr i j + D i j + P i j + ε i j + V i j .
where the symbols Pr i j , D i j , P i j , ε i j , and V i j represent the turbulent production, turbulent diffusion, pressure-strain correlation, turbulent dissipation, and the mass flux variation, respectively.
Pr i j = ρ ¯ R j m U ˜ i , m ρ ¯ R i m U ˜ j , m , D i j = ( ρ u i u j u m ¯ + p u j ¯ δ i m + p u i ¯ δ j m τ i m u j ¯ τ j m u i ¯ ) , m , P i j = p ( u i , j + u j , i ) ¯ = P i j + 2 / 3 p u k , k ¯ δ i j , ε i j = τ i m u j , m ¯ τ j m u i , m ¯   , V i j = p ¯ , j u i ¯ + p ¯ , i u j ¯ + τ ˜ i m , m u j ¯ + τ ˜ j m , m u i ¯   .
The turbulent dissipation in the compressible turbulence was proposed in [2,17], as follows:
ε = ε s + ε c .
where, for homogeneous shear flow turbulence, ρ ¯ ε s = μ ¯ ω i ω i ¯ ( ω i is the fluctuating vorticity) and ε c = 4 / 3 μ ¯ d 2 ¯ are the incompressible and dilatational (or compressible) parts of the turbulent dissipation rate, respectively. The authors argued that the incompressible part of the dissipation can be modeled by using the incompressible equation model [16], namely:
t ( ρ ¯ ε s ) + x k ( ρ ¯ ε s U ˜ k ) = ρ ¯ ε s K ( C ε 1 R k m x m U ˜ k C ε 2 ε s ) x k ( C ε 3 ρ ¯ K ε s R k m x m ε s )

4. Compressible Turbulence Model for the Pressure—Strain

Compressible turbulence modeling of the pressure-strain correlation is well known as an important problem in many interesting engineering applications related to environment, combustion flows, and hypersonic flights. During the last decade, different experiments and numerical simulations have been developed to understand the compressibility effects on pressure-strain. Despite the complexity of these problems, the last few years have been marked by an abundance of numerical results and a lot of advancement in second order closure. At first, researchers concentrated on the roles of the dilatational terms: pressure dilatation and the compressible part of the turbulent dissipation. In this context, Sarkar [3], Simone et al. [4], Hamba [5], and Blaisdell [18], using their DNS results, showed that the dilatational terms represent nearly 12% of turbulent kinetic energy production. They concluded that the notable decrease inthe growth rate of the turbulent kinetic energy arising from compressibility was due to the inhibited turbulent production and not to the dilatational terms. In addition, these conclusions were confirmed by Adumitroaie et al. [11], Huang et al. [12], Goebel et al. [7], Foysi et al. [10], and Pantano et al. [9], who reported that the dilatational terms did not affect the compressible mixing layers. On the other hand, based on their research, they insisted that compressibility has a significant effect on the pressure-strain correlation via the pressure field. Consequently, the pressure-strain correlation requires a topic motivation in the second order model, which based on the averaged Navier−Stokes (RANS) equations. Sarkar [3] performed two series of DNS, cases A1 to A4, for varying the initial gradient Mach number M g 0 at a constant initial turbulent Mach number, M t 0 = 0.4 , and cases B1 to B3 for constant M g 0 , ( M g 0 = 0.22 ) with varying M t 0 . The reason behind the variety of DNS cases is to intercept the source changes in the turbulence structure at different compressibility levels, which are necessary for accurate modeling. The primary remarks coming from these DNS concern the behavior of the turbulent Mach number and the gradient Mach number in relation to the strong changes of the turbulent anisotropy. One can see that when the initial gradient Mach number is low, after an initial slight increase with the dimensionless time, S t ( S t 5 ), and M g shows a trend to become asymptotically constant, contrary to Mt, which grows constantly with S t . In the case of A4, where it increases very highly, similar results are found for DNS [4]. Obviously, this gives a reason for several researchers, such as Huang et al. [12], Hamba [5], and Park et al. [19], to distinguish the DNS results [3] into three levels of compressibility: low, moderate, and high compressibility. In addition, the DNS results show that the behavior of M g seems to be similar to the turbulence anisotropy, as it can be clearly seen for the shear stress anisotropy component b 12 . Consequently, the authors concluded that M g is an appropriate parameter to measure the strength of the structural compressibility, and it should be added to M t in order to study the homogeneous highly sheared effects on the structural compressibility. So, a major challenge is to develop a model for the pressure-strain that is able to correctly capture the different compressibility levels. Thus, some compressible models have been derived for the pressure-strain correlation, most of which are derived from a simple extension of the incompressible LRR model [14], which reads as below:
P i j * = C 1 ρ ¯ ε s b i j +   C 2 ρ ¯ K   ( S ˜ i j 1 3 S ˜ l l δ i j ) + C 3 ρ ¯ K [ b i k S ˜ j k + b j k S ˜ i k 2 3 b m l S ˜ m l δ i j ] + C 4 ρ ¯ K [ b i k Ω ˜ j k + b j k Ω ˜ i k ]
S ˜ i j = 0.5 ( U ˜ i , j + U ˜ j , i ) , Ω ˜ i j = 0.5 ( U ˜ i , j U ˜ j , i ) and b i j = R i j / 2 K δ i j / 3 .
The coefficients models are: C 1 = 3 , C 2 = 0.8 , C 3 = 1.75 and C 4 = 1.31 .
In this study, three compressible models for the pressure-strain correlation are considered. The models were developed by Adumitroaie et al. [11], Huang et al. [12], and Marzougui et al. [13]. The authors used different modeling approaches to modify the LRR incompressible model of the pressure-strain, making the coefficients Ci (i = 1, 2, 3, and 4) in Equation (7) a function of the turbulent Mach number, as seen in Table 1.
Adumitroaie et al. [11] derived a compressible model taking into account the compressibility effects. Their model for the pressure-strain is given by the following:
p i j * = C 1 ρ ¯ ε s b i j + ( 4 5 + 2 5 d 1 ) ρ ¯ K   ( S ˜ i j 1 3 S ˜ l l δ i j ) + 2 ρ ¯ K ( 1 C 3 + 2 d 2 ) [ b i k S ˜ j k + b j k S ˜ i k 2 3 b m l S ˜ m l δ i j ] ρ ¯ K ( 1 C 4 2 d 2 ) [ b i k Ω ˜ j k + b j k Ω ˜ i k 4 3 d 2 S ˜ k k b i j ]
The compressible coefficients d 1 and d 2 are determined from some compressible closures for the pressure-dilatation correlation. Here, we use the model of Sarkar et al. [2] to determine the coefficients d 2 = 0.15 M t and d 1 = 0 .
Huang et al. [12] assumed that the incompressible modeling approach of the pressure-strain can be used to develop turbulent models taking into account compressibility effects. The authors used a damping function to modify the LRR model for the pressure-strain as follows:
P i j * = C 1 ρ ¯ ε s b i j +   C 2 ρ ¯ K   ( S ˜ i j 1 3 S ˜ l l δ i j ) + C 3 ρ ¯ K [ b i k S ˜ j k + b j k S ˜ i k 2 3 b m l S ˜ m l δ i j ] + C 4 ρ ¯ K [ b i k Ω ˜ j k + b j k Ω ˜ i k ]
C 1 = 3.6 , C 2 = 0.8 , C 3 =   1 . 2 + 0.25 exp ( 0.05 / M t 3 ) and C 4 =   1 . 2 - 0.25 exp ( 0.05 / M t 3 ) .
Marzougui et al. [13] used the concept of the turbulent kinetic energy growth rate to introduce a compressibility correction on the LRR model coefficients [14], which became a polynomial functions of the turbulent Mach number.
P i j * = C 1 ρ ¯ ε s b i j +   C 2 ρ ¯ K   ( S ˜ i j 1 3 S ˜ l l δ i j ) + C 3 ρ ¯ K [ b i k S ˜ j k + b j k S ˜ i k 2 3 b m l S ˜ m l δ i j ] + C 4 ρ ¯ K [ b i k Ω ˜ j k + b j k Ω ˜ i k ]
C 1 = 3   ( 1.0 0.8 M t 2 ) , C 2 = 0.8 , C 3 =   1 . 75   ( 1 1.4 M t 2 ) and C 4 =   1 . 31   ( 1 0.5 M t ) .
The application of these models on compressible homogeneous shear flow has shown predictions that are in acceptable agreement with the DNS of Sarkar [3] for cases A1, A2 and A3, who did not report the case A4. Particularly, in this case, the model [13] predictions deviated excessively from the DNS results for all of the Reynolds stress anisotropy components, and the shear stress anisotropy b 12 predictions were opposite to the DNS results. This can be attributed to the excessive reduction observed in the predicted behavior of the pressure-strain shear component. On the other hand, the model gave a high pressure-strain contribution in the normal stress anisotropy.

Proposal Model

The starting point of the proposed compressibility model is from some analyses and developments concerning the models developed by Adumitroaie et al. [11], Huang et al. [12], and Marzougui et al. [13]. They adopted different modeling approaches to modify the coefficients of the LRR model, which became a function of the turbulent Mach number (see Table 1). However, as can be seen, this modification is solely concerned with the coefficients C3 and C4, which affect the polynomial linear term of the Reynolds stress anisotropy and the mean strain rate; the other coefficients, C2, which affects the mean strain rate, and C1 of the return to isotropy model, are conserved as in the LRR model, without any compressibility correction. On the other hand, different analyses have been carried to show the influence of the pressure-strain on the Reynolds stress behavior. Hamba [5] presented a fine analysis for the compressible homogeneous shear flow case, and confirmed that the reduction of the transverse component P 22 of the pressure-strain correlation principally caused the reduction of the transverse Reynolds stress R 22 , which in turn induced a systematic reduction of the shear Reynolds stress, the streamwise component P 11 of the pressure-strain, and then the growth rate of the turbulent kinetic energy. Thus, the compressibility correction of the coefficients C3 and C4 seem to be sufficient to capture compressibility effects. In this context, according to Park et al. [18] and Huang et al. [12], in addition to the compressibility correction of the coefficients C3 and C4, the coefficient C2 should be corrected with compressible parameters, such as M t and M g , or others. One can see that C2 directly affects the shear component P 12 of the pressure-strain, which has an evident contribution in the transport equation for the Reynolds shear stress, R 12 . On the other hand, the reduction of P 12 , which works as a sink term in the transport equation for R 12 , leads to an increase in the growth rate of the turbulent kinetic energy via the growth of R 12 . This is not suitable using model. So, more attention should be paid to the modeling for P 12 . Khlifi et al. [20] considered an equation of the dilatation fluctuation to modify the incompressible C2 and the retour to isotropy C1-coefficients [14], as follows:
C 2 = 0.8 ( 1 + 0.45 M t 4 ) exp ( 0.015 M g ) ,
C 1 = 3 ( 10.7 M t 2 )
This is the first point on which the present study is based in order to revise the previous models. For this revision, the Khlifi et al. [20] model was chosen to modify the C1 and C2 coefficients in [11,12,13]. The reason behind this choice is that model [20] involves Mt and Mg, as suggested by different studies cited previously, and it is linked to M t 4 to distinguish between low- M t and high- M t regimes. Thus, all of the coefficients of models [11,12,13] are expressed as a function of the turbulent Mach number and the gradient Mach number. Considering Equations (11) and (12), the proposal coefficients models are summarized in Table 2.

5. Applications

5.1. Simulation of Compressible Homogeneous Shear Flow

For compressible homogeneous shear flow, the mean velocity gradient is given by the following:
U ˜ i , j = S δ i 1 δ j 2 .
For homogeneous shear flow, ρ ¯ = c t e and T ˜ = T ˜ ( t ) is related to the Reynolds-average of the pressure using the state equation for ideal gas:
P ¯ = ρ ¯ R T ˜ .
The Favre averaged basic second order model equations are as follows:
ρ ¯ d d t R i j = Pr i j + P i j * 2 3 ρ ¯ ( ε s + ε c ) δ i j + 2 3 p d ¯ δ i j ,
ρ ¯ d d t ε s = C ε 1 ρ ¯ ε s K R k m x m U ˜ k C ε 2 ρ ¯ ε s 2 K .
Assuming that the mean specific heat is constant, the equivalently temperature equation for the Reynolds averaged energy may be written in a simplified form [21], namely:
ρ ¯ C ¯ v d d t T ˜ = ρ ¯ ( ε s + ε c ) - π d .
Contraction i = j in Equation (12) leads to an equation for the Favre-averaged turbulent kinetic energy, K = 0.5 ρ u i u i ¯ / ρ ¯ ,as follows:
ρ ¯ d d t K = P ρ ¯ ( ε s + ε c ) + p d ¯ .
where P = ρ ¯ R i j U ˜ i , j is the turbulent production.
The transport equation for the turbulent Mach number, M t = 2 K / γ R T ˜ , can be obtained from combining Equations (17) and (18) as follows [21]:
d d t M t = M t 2 ρ ¯ K ( 1 + 0.5 γ ( γ 1 ) M t 2 ) ( p d ¯ ρ ¯ ε ) + M t 2 K P .
where γ = c p / c v .

5.2. Simulation of Compressible Mixing Layers

The flow is governed by the averaged Navier-Stokes equations associated with those that describe the energy, Reynolds stress, and turbulent dissipation. The simplest resulting continuity, momentum, and energy equation for stationary mixing layers can be written as follows [16]:
x i ρ ¯ U ˜ i = 0 ,
x i ( ρ ¯ U ˜ i   U ˜ j ) = x j ( ρ u i u j ¯ ) ,
x j ( ρ ¯ C v T ˜ U ˜ j ) = x j ( C v ρ u j T ¯ ) + ρ ¯ ( ε s + ε c ) - p u i , i ¯ .
The Reynolds stress is the solution of the follow equation
x m ( ρ ¯ U ˜ m R i j ) = ( R i m U ˜ j , m + R j m U ˜ i , m ) + x m (   ρ u i u j u m ¯ ) + φ i j * + 2 3 p u i , i ¯ δ i j 2 3 ε δ i j .
The turbulent solenoidal dissipation rate shall is calculated from the classical model equation, namely:
x k ( ρ ¯ ε s U ˜ k ) = ρ ¯ ε s K ( C ε 1 R k m x m U ˜ k C ε 2 ε s ) x k ( C ε 3 ρ ¯ K ε s R k m x m ε s ) .
In the above mentioned transport equations, different terms should be modeled, and the gradient diffusion hypothesis is used to represent the following:
  • The turbulent heat flux [16]:
ρ u i T ¯ = C T K ε ρ u i u m ¯ x m T ˜ .
  • The diffusion term [16]
ρ u i u j u m ¯ = C s K ρ ¯ ε ρ u i u m ¯ x m ρ u j u m ¯ .

6. Results and Discussion

The ability of the models of Adumitroaie et al. [11], Huang et al. [12], and Marzougui et al. [13] models, as well as their proposal modified forms, called the Adumitroaie modified, Huang modified, and Marzougui modified models, respectively, for the pressure-strain correlation (see Table 1 and Table 2), in order to predict compressible homogeneous turbulent shear and planar mixing layer turbulent flows will now be considered.

6.1. Homogeneous Shear Flow

The averaged transport for compressible homogeneous turbulence, given by Equations (15)–(19), is solved numerically using a fourth-order Runge−Kutta numerical scheme. Figures presented in this paper show the comparisons between the predictions obtained by the proposed models, called the Adumitroaie modified, Huang modified, and Marzougui modified models, and those from their corresponding original models, developed by Adumitroaie et al. [11], Huang et al. [12], and Marzougui et al. [13], respectively, and with the DNS results [3]. The reason for considering various models tests is to verify the proposal modifications validation and its added value in the prediction of the compressibility effects; especially when the proposed corrections are applied to three models build on different approaches, as shown in the indicated Section 4. In the same context, the DNS results considered here correspond to the different initial conditions listed in Table 3, from a low compressibility in case A1, to case A4, for which the compressibility is higher. From all of the figures, it is clear that both models provide a qualitative performance to reproduce the DNS results for all cases. However, one can see the accuracy of the proposal models for the prediction of the major characteristic parameters for compressible homogeneous shear flow, which can be expected from the model’s revision.
Figure 1 shows the normalized dissipation esk (esk = ε s / S K ) for cases A1 to A4 from DNS [3]. The results show that the proposal models better predict the trend of the decrease of ε s / S K when M g 0 increases, as compressibility effects cause a significant reduction in the turbulent production from numerical simulation cases A1 to A4. In fact, ( ε s / S K = 2 b 12 ( ε s / P ) , according to DNS [3], ( ε s / P ) appeared to be insensitive to compressibility effects and showed little differences between cases A1 to A4. This implies that the strong reduction in the shear stress anisotropy, b 12 , is responsible for the lowered gradient Mach number increases.
Figure 2 presents the behavior of the turbulent Mach number for cases A1 to A4. It is clearly seen that for both models, the original models and their modified versions are similar in their prediction of the correct trend of, an increase with creasing initial M g 0 for little dimensionless time ( S t 5 ). For ( S t 10 ), the difference between both models predictions are smaller in cases A1 and A2. On the other hand, in case A3 and particularly in case A4, the models become different; this can be seen in Figure 2d, which corresponds to case A4. One can see that the modified models predict the asymptotic trend of M t better than the original models.
Figure 3, Figure 4 and Figure 5 show the non-dimensional time S t variation of the Reynolds stress anisotropies of b 11 , b 12 , and b 22 , respectively, for cases A1 to A4. For case A1, which corresponds to low compressibility, the results are shown in Figure 3a, Figure 4a and Figure 5a, with both models being nearly similar, and the difference between their predictions being smaller. The models provide an acceptable performance for reproducing the DNS [3] results of this case. The results for case A2 and A3 are shown in Figure 3b,c, Figure 4b,c and Figure 5b,c. The Adumitroaie et al. [11] and Huang et al. [12] models appear to be insensitive to the increase of the streamwise b 11 and the transverse b 22 of the Reynolds stress anisotropies when the compressibility increases; the models’ results are in disagreement with the DNS data. However, these models have an acceptable performance for the prediction of the correct behavior of the reduced shear stress anisotropy b 12 . For the Marzougui et al. [13] model, the results are in good agreement with the DNS data, the model predicts the asymptotic behaviors of the shear stress anisotropy and the normal components b 11 and b 22 well.
Figure 3d, Figure 4d and Figure 5d show the result predictions for case A4, which corresponds to a high compressibility. The DNS data show that there is a strong amplification of energy arising from the compressibility and then the normal stress components become stronger. The Adumitroaie et al. [11] and Huang et al. [12] models predicted the same anisotropy levels as in previous cases; these models were still unable to predict the changes in the magnitude of the normal Reynolds stress anisotropy when the compressibility was higher. As can be seen in the previous figures, the models [11,12] appear to be insensitive to the increase of the streamwise b 11 and the transverse b 22 Reynolds stress anisotropies when the compressibility increases. However, these models are in qualitative agreement with the DNS data for the shear stress anisotropy b 12 .
The Marzougui et al. [13] model predictions deviated from the DNS data for the all of the Reynolds stress anisotropies components. The authors attributed the deficiency of their model to the fact that A4 falls within the rapid distortion limited for which the non-linear effects of the pressure-strain are unimportant, and they recommend an eventual revision of the model; this is what has been investigated by Khlifi et al. [20]. On the other hand, from Figure 3, Figure 4 and Figure 5, one can see the improvements brought about by the proposed modifications of the existent models [11,12,13] for the pressure-strain correlation in the prediction of the changes in the magnitude of the compressible Reynolds stress anisotropies. The modified models give results for the normal components b 11 and b 22 , which are in reasonable agreement with the DNS data and predict the asymptotic shear stress anisotropy well for all cases.
Figure 6, Figure 7 and Figure 8 present the behavior of the pressure-strain correlation P i j = P * i j / 2 S K , ( i ,   j = 1 ,   2 ) for cases A1 to A4 from DNS [3].The results for case A1 are shown in Figure 6a, Figure 7a and Figure 8a; all of the models are in qualitative acceptable agreement with the DNS results. The results for cases A2 to A4 are shown in Figure 6b–d, Figure 7b–d and Figure 8b–d; the modified models yielded reasonable agreement with the DNS data, and clearly predicted the reduced pressure-strain components better than their original versions. This was concluded from the DNS results [3,4], in which it is argued that the pressure fluctuations and all components of the pressure-strain correlation showed a monotone decrease with the increasing gradient Mach number.

6.2. Mixing Layers

In this study, thecomputation of two free streams of a fully developed compressible mixing layer (see Figure 9) is examined. The flows are characterized typically by the parameters s = ρ 2 / ρ 1 and r = U 2 / U 1 , which are the density and velocity ratios, respectively. The experiment conditions of Goebel et al. [7] are listed in Table 4.
The basic Equations (20)–(26) on which the second order model for the stationary compressible mixing layers is based are solved using a finite difference scheme. The grid of the computational physical domain, which is a rectangular box defined by the set of point (x, y), has 6666 × 41 points. The initial profiles for ε s , ρ ¯ and, T ˜ which are not available in the experiment of Goebel et al. [7], are generated as follows:
  • The initial profile of the turbulent dissipation is determined from the turbulent viscosity model.
ε s = C μ ρ ¯ K 2 ρ u v ¯ y U ˜ C μ = 0.09 ,
  • The initial profile of the temperature is obtained from the following similarity
U ˜ U 2 ( U 1 U 2 ) = T ˜ T 2 ( T 1 T 2 ) ,
  • The state equation of perfect gas is used to determine the initial profile of the density.
The values of the constant models used in the present simulation are as follows:
C ε 1 = 1.4 , C ε 2 = 1.8 , C μ = 0.09 , C ε = 0.25 , C T = 0.26 .
According to Sarkar [3], homogeneous shear flow is closely related to the mixing layers, this allows M g to be connected to M c . Thus, according to Khlifi et al. [20], the coefficients Ci of the Adumitroaie et al. [11] model are expressed as a function of the turbulent Mach number and the convective Mach number, as shown in Table 5.
In this work, the simulation was limited to the Adumitroaie et al. [11] and Adumitroaie modified (referred to presently) models, as listed in Table 1 and Table 5. Figure 10, Figure 11, Figure 12 and Figure 13 compare the computed results from the Adumitroaie et al. model and the present models with the experiments of Goebel et al. [7] for three cases: M c = 0.2 , M c = 0.69 , and M c = 0.86 , for which the important characteristic parameter of the developed planar mixing layers is represented in relation to the similarity variable y * = ( y y c ) / δ , where y is the local cross stream coordinate and y c is the cross-stream coordinate corresponding to U * = 0.5 . The normalized stream mean velocity U * = ( U ˜ U 2 ) / ( U 1 U 2 ) , is plotted in Figure 10, the calculated profiles for both of the two models are in good agreement with the experiment results for a low convective Mach number M c = 0.2 , and the largest convective Mach number M c = 0.86 , as shown in Figure 10a,b, respectively. The expected effects of the two Mach numbers M t and M c on the changes in pressure-strain correlation are clearly seen in Figure 11, Figure 12 and Figure 13, which compare the Reynolds similarity intensity profiles as follows: the streamwise intensity R 11 = ρ u 2 ¯ / ρ ¯ ( U 1 U 2 ) 2 , the transverse intensity R 22 = ρ v 2 ¯ / ρ ¯ ( U 1 U 2 ) 2 , and the shear stress R 12 = ρ u 2 v 2 ¯ / ρ ¯ ( U 1 U 2 ) 2 obtained by the present model and by Adumitroaie et al. [11] model with the experiment results of Goebel et al. [7]. The two models have nearly similar behaviors for small values of convective Mach numbers ( M c = 0.2 ), as observed in Figure 11a,b, Figure 12a,b and Figure 13a,b. When the compressibility is more significant, M c = 0.86 , Figure 11c, Figure 12c and Figure 13c show that the present model results are in reasonable agreement with the experiment data compared with those obtained by the Adumitroaie et al. [11]. The compressibility effects via the Mc-influence on the mixing layers anisotropy can be more clearly seen in Figure 14, which compares the variation of the maximum values of the Reynolds intensities: ( R 11 ) max , ( R 22 ) max , and ( R 12 ) max versus the convective Mach number with different experiment and DNS data [7,9,10]. We found that the present model better predicts the decrease in the maximum values of all Reynolds stress intensities with increasing the initial value of M c . However, the influence of M c on the variation of ( R 11 ) max is not significant (Figure 14a) as can be seen for the other components in Figure 14b,c. This behavior can be found in the analysis of [19], who concluded that the reduction of the pressure-strain directly caused a reduction of R 22 . However, the reduction of R 11 involved turbulent production and the pressure-strain terms.
Figure 15 compares the centerline values of the Reynolds stress anisotropy components: ( b 11 ) max , ( b 22 ) max and ( b 12 ) max versus convective Mach number computed by the present model and those given by Adumitroaie et al. model with the DNS datas of Freund et al. [22]. It is found that the present model predicts the decreasing trend of shear stress anisotropy b 12 and the increasing trend of all the normal stress anisotropies b 11 and b 22 with increasing convective Mach number well. The Adumitroaie et al. model shows that such trends occurred, but incorrectely followed the DNS data; thus, the peaks of b i j were not relatively affected by M c . This gives reason to imply M c in addition to M t for the eventual modeling of the pressure-strain correlation, which controls the anisotropy changes arising from the compressibility effects. Regarding Figure 14a, one can remark that both models give nearly similar peaks for R 11 , which decrease less than the other components R 22 and R 12 when increasing M c . This leads to a monotone increase of all the normal stress anisotropies with M c , as can be seen in Figure 15a,c; this behavior was found by the DNS of Vreman et al. [6] and Freund et al. [22]. Obviously, this act comes from compressibility via a reduction in the magnitude of all components of the pressure-strain correlation with M c .
From the previous results, one can conclude that the pressure-strain closure involving the parameter M g appears to be suited to study compressibility effects on highly sheared homogeneous turbulence. In addition, the convective Mach number M c seems to be an appropriate parameter to study such effects on mixing layers. There are differences in the accuracy between the models in the prediction compressibility effects on homogeneous shear flow and mixing layers. These differences can be attributed to the modeling approaches on which the previous models are built.

7. Conclusions

In this study, the Favre Reynolds stress model is used for the prediction compressibility effects in two important turbulent flow cases—homogeneous shear flow and the spatially planar mixing layers. Evaluation of the density extension turbulence models of Adumitroaie et al. [11], Huang et al. [12], and Marzougui et al. [13] for the pressure-strain correlation were examined. A revision of these models by considering the Khlifi et al. [20] model, lead to expressing all of the coefficients models [11,12,13] in function of the gradient and convective Mach numbers, in addition to the turbulent Mach number. The application of the proposed models referred to as Adumitroaie modified, Huang modified, and Marzougui modified models, to predict the previous flows shows satisfactory agreement with the available DNS data and the experiment results. The proposed models appear to be able to accurately predict the structural compressibility effects on homogeneous shear flow as a significant decrease in the magnitude of the Reynolds shear stress and a reduction in the pressure-strain components with increasing initial values of the gradient Mach number. In addition, the proposed model, referred to as the Adumitroaie modified model, successfully predicted the changes in the similarity of the Reynolds shear stress arising from the compressibility effects on the mixing layers. From the previous results, it is clear that the proposed models’ predictions are better than those obtained by their original models. Therefore, as a priority, the blending between different compressible models is seen as an important issue in the modeling of the pressure-strain correlation.

Author Contributions

Writing-original draft preparation, H.K.; writing-Review editing, H.K and A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

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.

Nomenclature

pPressure
TTemperature
tTime
aSpeed of sound
bijReynolds stress anisotropy
cpSpecific heat at constant pressure
cvSpecific heat at constant volume
RIdeal gas constant
RijReynolds stress
MtTurbulent mach number
MgGradient mach number
McConvective mach number
KTurbulent kinetic energy
u i Velocity in the direction xi
d’Fluctuation of the dilatation
(.),iXi-derivative
Greek symbols
γ Specific heat ratio
ε Turbulent dissipation
ε s Solenoidal dissipation
ε c Compressible dissipation
ρ Density
μ Viscosity coefficient
κ Thermal conductivity coefficient
π d = p d ¯ Pressure−dilatation correlation
P i j * Deviator of the pressure-strain tensor
δ i j Kronecker delta
τ i j Viscous stress tensor
Statistic symbols
(.)″Favre fluctuation
(.)′Reynolds fluctuation
(   )   ˜ Favre averaged
( . ) ¯ Reynolds averaged

References

  1. Blaisdell, G.A.; Mansour, N.N.; Reynolds, W.C. Compressibility effects on the growth and structure of homogeneous turbulent shear flow. J. Fluid Mech. 1983, 256, 443–485. [Google Scholar] [CrossRef]
  2. Sarkar, S.; Erlebacher, G.; Hussaini, M.Y.; Kress, H.O. The analysis and modeling of dilatational terms in compressible turbulence. J. Fluid Mech. 1991, 227, 473–493. [Google Scholar] [CrossRef]
  3. Sarkar, S. The stabilizing effects of compressibility in turbulent shear flows. J. Fluid Mech. 1995, 282, 163–186. [Google Scholar] [CrossRef] [Green Version]
  4. Simone, S.; Coleman, G.N.; Cambon, C. The effect of compressibility on turbulent shear flow: A rapid distorsion-theory and direct numerical simulation study. J. Fluid Mech. 1997, 330, 307–338. [Google Scholar] [CrossRef] [Green Version]
  5. Fujihiro, H. Effets of pressure fluctuations on turbulence growth compressible homogeneous shear flow. Phys. Fluid. 1999, 11, 1623–1635. [Google Scholar]
  6. Vreman, A.W.; Sandham, N.D.; Luo, K.H. Compressible mixing layer growth rate and turbulence characteristics. J. Fluid Mech. 1996, 330, 235–258. [Google Scholar] [CrossRef] [Green Version]
  7. Goebel, S.G.; Dutton, J.C. Experimental study of compressible mixing layers. AIAA J. 1991, 29, 538–546. [Google Scholar] [CrossRef]
  8. Samimy, M.; Elliot, G.S. Effects of compressibility on the characteristics of the free shear layers. AIAA J. 1990, 89, 439–445. [Google Scholar] [CrossRef] [Green Version]
  9. Pantano, C.; Sarkar, S. A study of compressibility effects in the high speed turbulent shear layer direct simulation. J. Fluid Mech. 2002, 451, 329–371. [Google Scholar] [CrossRef] [Green Version]
  10. Foysi, H.; Sarkar, S. The compressible mixing layers: An LES study. Theor. Comput. Fluid Dyn. 2010, 24, 565–588. [Google Scholar] [CrossRef] [Green Version]
  11. Adumitroaie, V.; Ristorcelli, J.R.; Taulbee, D.B. Progress in Favre Reynolds stress closures for compressible flows. Phys. Fluids 1999, A11, 2696–2719. [Google Scholar] [CrossRef] [Green Version]
  12. Hung, S.; Song, F. Modelling of pressure strain correlation compressible turbulent flow. Acta Mech. Sin. 2008, 24, 37–43. [Google Scholar] [CrossRef]
  13. Hamed, M.; Hechmi, K.; Taieb, L. Extension of the Launder Reece and Rodi on Homogeneous shear flow. Eur. Phys. J. B 2005, 45, 147–154. [Google Scholar]
  14. Launder, B.E.; Reece, G.J.; Rodi, W. Progress in the development of a Reynolds-stress turbulence closure. J. Fluid Mech. 1975, 68, 537–566. [Google Scholar] [CrossRef] [Green Version]
  15. Hechmi, K.; Taieb, L. On the compressibility effects in mixing layers. Therm. Sci. 2016, 20, 1473–1484. [Google Scholar]
  16. Speziale, C.G.; Sarkar, S. Second Order Closure Models for Supersonic Turbulent Flows; ICASE Report; NASA Langley Research Center: Hampton, VA, USA, 1991.
  17. Zeman, O. Dilatational dissipation, the concept and application in modeling compressible mixing layers. Phys. Fluids 1990, A2, 178–188. [Google Scholar] [CrossRef]
  18. Blaisdell, G.A.; Sarkar, S. Investigation of the pressure-strain correlation in compressible homogeneous turbulent shear flow. ASME-Publ.-FED 1993, 151, 133–138. [Google Scholar]
  19. Park, C.H.; Park, S.O. Compressible turbulence model for the pressure strain correlation. J. Turbul. 2005, 6, 1–24. [Google Scholar] [CrossRef]
  20. Hechmi, K.; Taieb, L. A compressibility correction of the pressure strain correlation model in turbulent flow. Comptes Rendus Mécanique 2013, 341, 567–581. [Google Scholar]
  21. Speziale, C.G.; Abid, R.; Mansour, N. Evaluation of Reynolds stress turbulence closures in compressible homogeneous shear flow. J. Appl. Math. Phys. 1995, 46, 5717–5736. [Google Scholar]
  22. Freund, J.B.; Lele, S.K.; Monin, P. Compressibility effects in a turbulent annular mixing layer. J. Fluid Mech. 2000, 421, 229–267. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Time evolution of the normalized dissipation, (esk = ε s / S K in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Figure 1. Time evolution of the normalized dissipation, (esk = ε s / S K in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Fluids 07 00034 g001
Figure 2. Time evolution of the turbulent Mach number, in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Figure 2. Time evolution of the turbulent Mach number, in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Fluids 07 00034 g002
Figure 3. Time evolution of the streamwise Reynolds stress anisotropy, in the following cases: (a) 1; (b) A2; (c) A3; (d) A4.
Figure 3. Time evolution of the streamwise Reynolds stress anisotropy, in the following cases: (a) 1; (b) A2; (c) A3; (d) A4.
Fluids 07 00034 g003aFluids 07 00034 g003b
Figure 4. Time evolution of the shear Reynolds stress anisotropy, b 12 , in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Figure 4. Time evolution of the shear Reynolds stress anisotropy, b 12 , in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Fluids 07 00034 g004
Figure 5. Time evolution of the transverse Reynolds stress anisotropy, in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Figure 5. Time evolution of the transverse Reynolds stress anisotropy, in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Fluids 07 00034 g005
Figure 6. Time evolution of the pressure-strain component, in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Figure 6. Time evolution of the pressure-strain component, in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Fluids 07 00034 g006aFluids 07 00034 g006b
Figure 7. Time evolution of the pressure-strain component, P 22 in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Figure 7. Time evolution of the pressure-strain component, P 22 in the following cases: (a) A1; (b) A2; (c) A3; (d) A4.
Fluids 07 00034 g007
Figure 8. Time evolution of the pressure-strain component in the following cases: (a) A1, (b) A2; (c) A3; (d) A4.
Figure 8. Time evolution of the pressure-strain component in the following cases: (a) A1, (b) A2; (c) A3; (d) A4.
Fluids 07 00034 g008
Figure 9. Turbulent mixing layer.
Figure 9. Turbulent mixing layer.
Fluids 07 00034 g009
Figure 10. Similarity profiles of the mean velocity, in the following cases: (a,b).
Figure 10. Similarity profiles of the mean velocity, in the following cases: (a,b).
Fluids 07 00034 g010
Figure 11. Similarity profiles of R 11 in the cases: (a) M c = 0.2 , (b) M c = 0.69 , and (c) M c = 0.86 .
Figure 11. Similarity profiles of R 11 in the cases: (a) M c = 0.2 , (b) M c = 0.69 , and (c) M c = 0.86 .
Fluids 07 00034 g011
Figure 12. Similarity profiles of in the cases: (ac).
Figure 12. Similarity profiles of in the cases: (ac).
Fluids 07 00034 g012
Figure 13. Similarity profiles of R 12 in the cases: (a) M c = 0.2 , (b) M c = 0.69 , and (c) M c = 0.86 .
Figure 13. Similarity profiles of R 12 in the cases: (a) M c = 0.2 , (b) M c = 0.69 , and (c) M c = 0.86 .
Fluids 07 00034 g013
Figure 14. The influence of the convective Mach number on the maximum Reynolds stress: (a) ( R 11 ) max , (b) ( R 22 ) max and, (c) ( R 12 ) max .
Figure 14. The influence of the convective Mach number on the maximum Reynolds stress: (a) ( R 11 ) max , (b) ( R 22 ) max and, (c) ( R 12 ) max .
Fluids 07 00034 g014
Figure 15. The influence of the convective Mach number on the maximum Reynolds stress anisotropies: (a) ( b 11 ) max , (b) ( b 22 ) max , and (c) ( b 12 ) max .
Figure 15. The influence of the convective Mach number on the maximum Reynolds stress anisotropies: (a) ( b 11 ) max , (b) ( b 22 ) max , and (c) ( b 12 ) max .
Fluids 07 00034 g015
Table 1. Numerical coefficients of the pressure-strain model.
Table 1. Numerical coefficients of the pressure-strain model.
ModelC1C2C3C4
Adumitroaie [11]30.8 1.75 + 0.15 M t 1.3 0.15 M t
Huang [12]3.60.8 1.2 + 0.25 exp ( 0.05 / M t 3 ) 1.2 0.25 exp ( 0.05 / M t 3 )
Marzougui [13] 3   ( 1.0 0.8 M t 2 ) 0.8 1 . 75   ( 1 1.4 M t 2 ) 1 . 31   ( 1 0.5 M t )
Table 2. Numerical coefficients of the pressure-strain model.
Table 2. Numerical coefficients of the pressure-strain model.
Model C1C2C3C4
Adumitroaie modified 3 . ( 1.0 0.7 M t 2 ) 0.8 ( 1 + 0.45 M t 4 ) e 0.015 M g 1.75 + 0.15 M t 1.3 0.15 M t
Huang modified 3 . ( 1.0 0.7 M t 2 ) 0.8 ( 1 + 0.45 M t 4 ) e 0.015 M g 1.2 + 0.25 e 0.05 M t 1.2 0.25 e 0.05 M t
Marzougui modified 3 . ( 1.0 0.7 M t 2 ) 0.8 ( 1 + 0.45 M t 4 ) e 0.015 M g C 3 =   1 . 75 ( 1 1.4 M t 2 ) C 4 =   1 . 31 ( 1 0.5 M t )
Table 3. Initial conditions for homogeneous shear flow: DNS [3].
Table 3. Initial conditions for homogeneous shear flow: DNS [3].
Case Mt0 ( ε s / S K ) 0 Mg0b11b22b12
A10.41.80.22000
A20.43.60.44000
A30.45.40.66000
A40.410.81.32000
B10.135.40.22000
B20.23.60.22000
B30.41.80.22000
Table 4. Experiment of Goebel and Dutton.
Table 4. Experiment of Goebel and Dutton.
M c r = U 2 U 1 s = ρ 2 ρ 1
0.20.760.78
0.4600
0.6900
0.860.160.6
Table 5. Numerical coefficients of the pressure-strain model.
Table 5. Numerical coefficients of the pressure-strain model.
ModelC1C2C3C4
Adumitroaie modified 3 ( 1.0 2.5 M t 2 ) 0.8 ( 1 + 4.5 M t 4 ) exp ( 0.00022 M c ) 1.75 + 0.15 M t 1.3 0.15 M t
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khlifi, H.; Bourehla, A. A Compressible Turbulence Model for Pressure—Strain. Fluids 2022, 7, 34. https://doi.org/10.3390/fluids7010034

AMA Style

Khlifi H, Bourehla A. A Compressible Turbulence Model for Pressure—Strain. Fluids. 2022; 7(1):34. https://doi.org/10.3390/fluids7010034

Chicago/Turabian Style

Khlifi, Hechmi, and Adnen Bourehla. 2022. "A Compressible Turbulence Model for Pressure—Strain" Fluids 7, no. 1: 34. https://doi.org/10.3390/fluids7010034

APA Style

Khlifi, H., & Bourehla, A. (2022). A Compressible Turbulence Model for Pressure—Strain. Fluids, 7(1), 34. https://doi.org/10.3390/fluids7010034

Article Metrics

Back to TopTop