Next Article in Journal
Study on An Effective Roadway Watering Scheme for Mitigating Pedestrian Thermal Comfort According to the Street Configuration
Next Article in Special Issue
The Decay of Energy and Scalar Variance in Axisymmetric Turbulence
Previous Article in Journal
Volatile Organic Compound Emissions from Indoor Fragrance Diffusers
Previous Article in Special Issue
Turbulence with Magnetic Helicity That Is Absent on Average
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Equilibrium Turbulent Transport in Convective Plumes Obtained from Closure Theory

by
Nobumitsu Yokoi
Institute of Industrial Science, University of Tokyo, Tokyo 153-8505, Japan
Visiting Researcher at Nordic Institute for Theoretical Physics (NORDITA).
Atmosphere 2023, 14(6), 1013; https://doi.org/10.3390/atmos14061013
Submission received: 26 April 2023 / Revised: 25 May 2023 / Accepted: 5 June 2023 / Published: 12 June 2023

Abstract

:
The non-equilibrium property of turbulence modifies the characteristics of turbulent transport. With the aid of response function formalism, such non-equilibrium effects in turbulent transport can be represented by the temporal variation of the turbulent energy (K) and its dissipation rate ( ε ) along the mean stream through the advective derivatives of K and ε . Applications of this effect to the turbulent convection with plumes are considered for the first time in this work. The non-equilibrium transport effects associated with plumes are addressed in two aspects. Firstly, the effect associated with a single plume is evaluated using data measured in the recent plume/jet experiments. The second argument is developed for the collective turbulent transport associated with multiple plumes mimicking the stellar convection zone. In this second case, for the purpose of capturing the plume motions into the advective derivatives, use has to be made of the time–space double-averaging procedure, where the turbulent fluctuations are divided into the coherent or dispersion component (which represents plume motions) and the incoherent or random component. With the aid of the transport equations of the coherent velocity stress and the incoherent counterpart, the interaction between the dispersion and random fluctuations are also discussed in the context of convective turbulent flows with plumes. It is shown from these analyses that the non-equilibrium effect associated with plume motions is of a great deal of relevance in the convective turbulence modeling.

1. Introduction

Turbulent flows encountered in astro/geophysical phenomena as well as laboratory plasma ubiquitously have huge Reynolds numbers. Since the direct numerical simulations (DNSs) of such flows are impossible in the foreseeable future, the turbulence or subgrid-scale (SGS) modeling approach provides a powerful tool for analyzing turbulent flows in these phenomena.
In conventional turbulence modeling, the simple expressions of turbulent fluxes such as the eddy viscosity and eddy diffusivity representations are very useful and have been extensively utilized [1]. For instance, the turbulent momentum flux, defined by the Reynolds stress, is modeled by the turbulent transport coefficients, eddy viscosity, coupled with the mean velocity strain. The eddy viscosity is expressed in terms of the quantities that characterize the dynamics and statistical properties of the turbulence. In the mixing-length formulation, the eddy viscosity ν T is expressed by the characteristic turbulent velocity v and length scale as ν T = ν T ( v , ) . In the standard K ε model, ν T is modeled in terms of the turbulent kinetic energy K and its dissipation rate ε as ν T = ν T ( K , ε ) . As compared with the mixing-length formulation, the K ε is more elaborated since the dynamics of the turbulent fields and their interaction with the mean fields are taken into account through the transport equations of K and ε [2,3].
In spite of the notable achievements of simple gradient-diffusion-type models for turbulent fluxes with the eddy viscosity and eddy diffusivity representations, several situations exist wherein such models would completely fail. One such deficiency appears when the eddy viscosity model is applied to a situation with cross-flow configuration. This failure is caused by the lack of vorticity or rotation effects in the eddy viscosity model. An alleviation of this deficiency of the eddy viscosity model was proposed with the aid of an analytical statistical theory of inhomogeneous turbulence applied to a non-reflectional symmetric system with an emphasis on the helicity effects [4,5]. Another kind of deficiency shows up in the case of turbulence with non-equilibrium properties. For instance, in homogeneous-shear turbulence, the evaluation of the eddy viscosity systematically fails because of the overestimates of the turbulent energy K and its dissipation rate ε . There, the energy injection increases with time, and the energy dissipation process cannot catch up with the injection at the energy-containing scale. This is caused by the non-equilibrium property of turbulence induced by the energy injection due to the flow shear. In order to alleviate this flaw in the framework of the K ε model with the eddy viscosity representation, the non-equilibrium effect has to be incorporated into the eddy viscosity representation [6].
Employing the response function formalism in turbulence closure studies provides a powerful way to construct a self-consistent turbulence theory without resorting to any parameter. There are several classical modern theories of turbulence closure [7,8,9,10]. Starting with the Liouville equation derived by Edwards [8], Herring [9,10] constructed a series of self-consistent field (SCF) theories for turbulence closure. In the formulation, the coupled equations of the two-time velocity correlation function and the response function are derived under the requirement that the probability distribution functions satisfy the self-consistency condition. The obtained equations for the correlation and response functions are similar to the those obtained with the direct-interaction approximation (DIA) by Kraichnan [7]. In the combination of the DIA with multiple-scale analysis, we address the convective turbulence with variable density and velocity shear. One important consequence of this multiple-scale DIA analysis lies in the point that it enables us to treat the non-equilibrium effect of turbulence through the Lagrangian derivative of the turbulence fields. If turbulence properties such as the turbulent energy K and its dissipation rate ε are spatially developing in the streamwise direction, the non-equilibrium effects are represented by the advective derivatives of K and ε , ( U · ) K and ( U · ) ε in the Lagrangian derivatives. In the model for the homogeneous-shear turbulent flow [6], the non-equilibrium properties of turbulence including the variations in the time and length scales of turbulence are taken into account in the expression of the eddy viscosity through the Lagrangian derivative along the mean velocity U of the turbulent energy and its dissipation rate as ν T = ν T ( K , ε , D K / D t , D ε / D t ) where D / D t = / t + U · .
There are several types of buoyant streams in nature. Some buoyant flows are purely driven by thermobuoyancy (plumes) and others are driven by a combination of thermobuoyancy and dynamical forcing (thermobuoyant jets, non-isothermal shear layers). In the buoyancy-dominated turbulent flows, the turbulent convection is well represented by large-scale plumes [11,12,13,14]. In order for a turbulence model to treat turbulent convective flows, the turbulence interaction with coherent structures such as a large-scale jet and plume has to be properly represented in the model. Because of the presence of jets and plumes, the turbulent fluxes deviate from the usual gradient diffusivity, viscosity, etc. [14]. The turbulent transport coefficients should properly reflect such coherent structure effects. Another example in astrophysics can be seen in supernova explosion studies. In the core-collapse supernovae (CCSNe) explosion problem, turbulent effects aid the relaunch of blast shock waves due to the neutrino heating [15,16,17,18,19,20,21]. In order to identify the turbulence effects most relevant to the explosion condition, Murphy and Meakin [21] extensively examined several types of closure models. By comparing the model simulation results with the 2D simulation counterparts, they found that the standard turbulent flux models based on the gradient-transport approximation failed to reproduce the turbulent transport of the mass, momentum and heat in the stellar interiors. This suggests that the buoyancy-dominated flows of core-collapse convection are characterized by large-scale plumes. The entrainment between the rising and sinking plumes strongly affects the transport properties of neutrino-driven turbulence. In particular, the turbulent dissipation, which is balanced by the buoyancy drive, is dominantly controlled by the entrainment of negatively buoyant plumes.
It has been theoretically and experimentally shown that the dissipation rates in jets and plumes follow the Kolmogorov equilibrium law only if their spreading rates do not vary with the streamwise distance [22,23]. Adopting the Lie symmetry group theory on the statistical turbulent model equations, Sunita and Layek [24] recently argued the case of non-equilibrium turbulent dissipation in buoyant axisymmetric plumes in unstratified stagnant ambient at the infinite-Reynolds-number limit. They considered the relationship between the dissipation rates of the turbulent kinetic energy and turbulent thermal fluctuation, and the entrainment coefficients in the Kolmogorov and non-Kolmogorov regions of turbulent dissipations. They showed that the well-known link between entrainment and dissipation in the Kolmogorov region of turbulence is also established in the non-Kolmogorov region.
As mentioned above, with the aid of the response function formalism, the non-equilibrium properties of plume motions can be expressed by the variations in the turbulent kinetic energy and its dissipation rate along the plume motions. Recent elaborated laboratory experiments with jets and negatively or positively buoyant plumes provide detailed data on the spatial distributions of the turbulent energy and its dissipation rate as well as the jet and plume velocity distributions [25,26,27,28]. Utilizing these data, we can evaluate the variation of the turbulent energy and its dissipation rate along the jet/plume motion. Such an evaluation enables us to examine how much non-equilibrium effect may be present in turbulent flows associated with the jet and plume. This is one of the novel points discussed in this work.
Negatively or positively buoyant plumes play a dominant role in mass, momentum and heat mixing in the stellar convection zone. In the stellar convection zone case, we have to simultaneously treat the effects of numerous plumes on the turbulent fluxes, such as the turbulent mass flux, momentum flux and heat or internal-energy flux, in the mean-field equations. The non-equilibrium effect associated with these numerous plumes is investigated with the aid of the numerical simulations of a simplified domain mimicking the stellar convection zone. A turbulence model incorporating the non-equilibrium effect into the expressions of the turbulence fluxes should be validated in comparison with the direct numerical simulations (DNSs) of the turbulent mass, momentum and heat or internal-energy fluxes. In order to properly capture the effects of plume motions, which are part of fluctuating motions in the framework of the mean-field modeling, we introduce some elaborated averaging procedures in this evaluation of the non-equilibrium turbulence model. This is a time–space double-averaging procedure, which has been recently applied to a stellar convection zone [29]. This point will also be discussed in this paper.
The organization of this paper is as follows. In Section 2, we present the mean-field equations with the turbulence correlations, which should be expressed in terms of the mean-field and transport coefficients representing the statistical properties of turbulence. In Section 3, closure models for the turbulent fluxes are promoted with the aid of their theoretical derivations, with special reference to the corrections to the conventional gradient-transport approximations. There, modifications due to the non-equilibrium effects are presented. In Section 4, the roles of the non-equilibrium effects are suggested in the context of recent plume/jet experiments. In Section 5, the role of the plumes in the turbulent transport is discussed in the context of the stellar convection. In addition, the interaction between the coherent and incoherent fluctuations is argued for using the transport equations of each component of fluctuation energy. Concluding remarks are given in Section 6.

2. Equations of Compressible Hydrodynamic Turbulence

The material of this work is plumes in compressible convective turbulence. The turbulent fluxes, such as the turbulent mass flux, the Reynolds stress, the turbulent heat, or internal-energy flux, in the mean-field equations are analyzed with the aid of the two-scale direct-interaction approximation (TSDIA) formulation. In this method, the response or Green’s function is used to incorporate the non-equilibrium properties of the turbulence along the mean or grid-scale (GS) flow into the turbulent transport coefficients. By applying the present formulation to the statistical analysis of the turbulent transports in the presence of numerous plumes, a time–space double-averaging procedure is adopted. Since there is no general way of determining the time-domain filtering, we have to make use of the physical considerations on the plume statistics and dynamics including its lifetime for determining the filter in the time domain. This point should be further explored beyond the present formulation described below.

2.1. Mean-Field Equations and Turbulent Correlations

We adopt the conventional Reynolds decomposition procedure
f = F + f , F = f
with
f = ( ρ , u , p , e )
F = ( ρ , U , P , E ) ,
f = ( ρ , u , p , e ) ,
where ρ is density, u is the velocity, p is the pressure and e is the internal energy of fluid. With the Reynolds averaging (1), the mean-field equations are obtained in the following form, and the density equation:
ρ t + · ρ U = · ρ u ,
the mean velocity equation:
t ρ U α + x a ρ U a U α = γ 1 x α ρ E + x a μ S a α x a ρ u a u α + U a ρ u α + U α ρ u a + ρ g α ,
and the mean internal-energy equation:
t ρ E + · ρ U E = · κ C V E · ρ e u + E ρ u + U ρ e γ 1 ρ E · U + ρ e · u + E ρ · u .
In (6), S α β is the deviatric part of the mean velocity strain defined by
S α β = U β x α + U α x β 2 3 · U δ α β ,
and g = { g α } is the gravitational acceleration.
In the mean-field Equations (5)–(7), we have several second-order turbulence correlations. They are the turbulent mass flux: ρ u , the Reynolds stress: u u , the turbulent internal-energy flux: e u , the internal-energy dilatation correlation: e · u , the density dilatation correlation: ρ · u , and the density–internal-energy correlation: ρ e . In the conventional turbulence modeling approach, the expressions of these turbulent correlations are modeled in a heuristic manner. In marked contrast, in a self-consistent turbulence modeling approach based on an analytical theory, the expressions of these correlations are obtained from a closure theory of inhomogeneous turbulence. In the present work, we use the analytical results from a multiple-scale renormalization perturbation expansion method; two-scale direct-interaction approximation (TSDIA) analysis of the fluctuation equations [3,30,31]. The analysis of incompressible and strongly compressible turbulence with a magnetic field using the TSDIA formalism can be seen, for example, in [32,33,34].

2.2. Simplest Modeling Turbulent Fluxes

In the turbulence modeling approaches, in order to reduce the huge burden of treating complexity of turbulent flows, it is often required that the turbulent fluxes be modeled in the simplest possible manner. A heuristic approach of modeling is to adopt the gradient-diffusion approximations for the turbulent fluxes appearing in (5)–(7). There, the turbulent momentum, mass and energy fluxes are assumed to be written in the form:
u α u β D = ν T S α β ,
ρ u = ν T σ ρ ρ ,
e u = ν T σ e E ,
where the suffix D denotes the deviatoric or traceless part of a tensor as A D α β = A α β ( 1 / 3 ) A a a δ α β . Here, the transport coefficient ν T is the turbulent or eddy viscosity, and σ ρ and σ e are the turbulent Prandtl numbers for the density and internal energy, respectively. The turbulent transport coefficients represent the statistical properties of turbulence. For example, the turbulent viscosity is expressed in terms of the intensity of turbulence or turbulent energy (per mass) K ( = u 2 / 2 ) .
In the mixing-length theory (MLT) model, the turbulent viscosity is written as
ν T = v C
apart from the numerical factor, where v is the characteristic intensity of turbulent velocity and C is the characteristic length in the energy-containing eddies. In practical applications of the model to the geophysical and astrophysical turbulent flows, depending on the problem we consider, the pressure or density-scale height is often adopted as the mixing length. In terms of the turbulent energy K, Equation (12) is expressed as
ν T = K 1 / 2 C
(the Kolmogorov–Prandtl expression).
This model expression for the turbulent viscosity ν T can be obtained from the consideration of the equilibrium energy spectrum of turbulence as follows. The Kolmogorov scaling of turbulence is based on the local equilibrium between the production rate P K and the dissipation rate ε K of the turbulent energy as
P K ε K .
Hereafter, the dissipation rate of turbulent energy, ε K , is simply denoted as ε . Assuming the equilibrium scaling of Kolmogorov:
E ( k ) = K o ε 2 / 3 k 5 / 3
( K o ( 1.5 ) : Kolmogorov constant), the turbulent energy is expressed as
K = d k E ( k ) = k = k C d k K o ε 2 / 3 k 3 / 5 = C K 0 ε 2 / 3 C 2 / 3 ,
where C ( = 2 π / k C ) is the largest eddy size. Here, we simply related it to the infrared cut-off wavenumber k C , and C K 0 is the model constant linked to the Kolmogorov constant K o . With (13), the turbulent or eddy viscosity ν T in (9)–(11) is evaluated as
ν T = K 1 / 2 C = ε 1 / 3 C 4 / 3 = K 2 / ε
apart from the numerical factors. The final equality is the well-known eddy viscosity expression in the standard K ε model. With this turbulent viscosity expression, Equation (9) constitutes the eddy viscosity representation of the Reynolds stress or turbulent momentum transport.

3. Non-Equilibrium Effect

Non-Equilibrium Effect on Eddy Viscosity

In the standard K ε model, the turbulent viscosity is expressed in terms of the turbulent energy K and its dissipation rate ε as
ν T = C ν K 2 ε ,
where C ν is model constant usually taken as C ν = 0.09 . However, it has long been known that the standard K ε model with the turbulent viscosity (18) applied to practical flows has some obvious drawbacks. The representative situations where the standard K ε model fails include the cross-flow configuration with the circumferential or rotating velocity component as well as the axial one, which impinges flow configuration. Another example is the overestimate of the turbulent energy K and its dissipation rate ε in the homogeneous-shear turbulence. The results of the standard K ε show a tremendously overestimated turbulent energy K as compared with the counterpart in the direct numerical simulations (DNSs) as schematically depicted in Figure 1. Even in one of the simplest flow geometries such as the homogeneous shear, a simple turbulent viscosity model fundamentally fails to reproduce turbulent flow characteristics.
In order to rectify such drawbacks of the standard K ε model, several improvements to the model have been proposed. A possible improvement of the model is the implementation of the non-equilibrium effect on the basis of a statistical theory of inhomogeneous turbulence. The non-equilibrium effect on turbulent viscosity was investigated with the two-scale direct-interaction approximation (TSDIA) [30,35]. In the framework of the TSDIA, the formal solution of the lowest-order turbulent velocity arising from the large-scale inhomogeneities in space and time is given in terms of the homogeneous isotropic basic field u 0 ( k ; τ 1 ) and the response function G ( k ; τ , τ 1 ) as
u 1 α ( k ; τ ) = U b X a τ d τ 1 G α b ( k ; τ , τ 1 ) u 0 a ( k ; τ 1 ) τ d τ 1 G α a ( k ; τ , τ 1 ) D u 0 a ( k ; τ 1 ) D T + 2 M d a b ( k ) δ ( k p q ) d p d q τ d τ 1 G α d ( k ; τ , τ 1 ) q b q 2 u 0 a ( p ; τ 1 ) u 0 c ( q ; τ 1 ) X c M a b c d ( k ) δ ( k p q ) d p d q τ d τ 1 G α d ( k ; τ , τ 1 ) X c ( u 0 a ( p ; τ 1 ) u 0 b ( q ; τ 1 ) ) .
Here, the second term on the right-hand side related to the Lagrange derivative represents the non-equilibrium properties of the turbulent field along the flow. The turbulent kinetic energy K with the non-equilibrium effect is written as
K = u 2 / 2 = 1 2 d k u ( k ; τ ) u ( k ; τ ) / δ ( k + k ) = ( u 0 u 0 + u 0 u 1 + u 1 u 0 + ) / 2 = C K 1 ε 2 / 3 C 2 / 3 C K 2 ε 3 / 2 C 4 / 3 D ε D t C K 3 ε 1 / 3 C 1 / 3 D C D t ,
where C is the length scale of the largest eddy and D / D t ( = / t + U · ) is the Lagrange or material derivative along the mean-flow velocity. Here, C K n ( n = 1–3) are the model constants whose magnitudes are related to those of the energy spectrum and the time scale of turbulence ( σ 0 in (A5) and ω 0 in (A6)). A brief derivation of (20) is given in Appendix A.
We solve (20) with respect to C by an iteration method with the first term of (20) as the starting term:
K = C K 1 ε 2 / 3 C 2 / 3 .
Then, we have C as
C = C 1 K 3 / 2 ε 1 ,
where C 1 is the model constant. With (13), this lowest-order expression of the length scale corresponds to the simplest eddy viscosity model (18).
Substituting (18) into (20), we have the first iteration result as
C = C 1 K 3 / 2 ε 1 + C 2 K 3 / 2 ε 2 D K D t C 3 K 5 / 2 ε 3 D ε D t
( C 2 and C 3 are model constants). The second and third terms of (23) represent the non-equilibrium effect. Equation (23) shows that the length scale of turbulence, as well as the time scale and energy of turbulence, is corrected by the D / D t -related non-equilibrium effect terms. Equation (23) can be approximated as
C = E 1 C N 1 K D D t K 2 ε ,
where E is the equilibrium characteristic length defined by
E = K 3 / 2 / ε ,
and C N is the non-equilibrium-effect-related coefficient. Depending on the streamwise development of the turbulent energy K and its dissipation rate ε , the length scale of turbulence changes as compared to the equilibrium cases. At the same time, the time scale of the non-equilibrium turbulence may be expressed as
τ NS = NS K 1 / 2 C S 3 / 2 K ε 1 C N 1 K D D t K 2 ε .
This indicates that the effective time scale may change depending on the D / D t effect.
The non-equilibrium turbulent viscosity ν T can be expressed as
ν T = C ν K 2 ε C ν K K 2 ε 2 D K D t + C ν ε K 3 ε 3 D ε D t ,
where C ν , C ν K and C ν ε are the model constants whose values can be evaluated from theoretical analyses. Equation (27) is rewritten as
ν T = C ν K 2 ε 1 C N K 1 ε D K D t C N ε K ε 2 D ε D t ,
or approximately
ν T = C ν K 2 ε 1 C N 1 K D D t K 2 ε
( C N K , C N ε and C N : model constants). Here, we should note that this non-equilibrium effect does not necessarily require compressibility. This effect shows up even in incompressible cases.
In practical situations wherein turbulent flows are simulated, a large positive D / D t term may occur locally; in such a case, (29) may lead to a physically undesirable negative turbulent viscosity. In order to alleviate such an unphysical model result to enhance the applicability of the model expression, we use a simple Padé approximation to (29) and have
ν T = ν TE 1 + C N 1 K D D t K 2 ε 1 for D D t K 2 ε > 0 , ν TE 1 C N 1 K D D t K 2 ε for D D t K 2 ε < 0 ,
where the equilibrium eddy viscosity ν TE is defined in the same form as (18) by
ν TE = C ν K 2 ε .
Equation (30) indicates that, depending on the sign of the Lagrangian or material derivative ( D / D t ) ( K 2 / ε ) , the turbulent viscosity is effectively decreased or increased by the non-equilibrium effect. In the case of ( D / D t ) ( K 2 / ε ) > 0 , ν T is decreased, while in the case of ( D / D t ) ( K 2 / ε ) < 0 , ν T is increased. This implies that, in this case, the non-equilibrium effect works in the direction to stabilize the deviation from the equilibrium state (return to equilibrium).
The non-equilibrium model expression for the eddy viscosity (30) has been applied to the homogeneous-shear turbulence, where the standard K ε model with the conventional eddy viscosity (31) leads to an overestimated turbulent energy K and its dissipation rate ε as compared to the DNS results (Figure 1). It was reported that the non-equilibrium model (30) with C N = 0.8 leads to the results consistent with the DNS ones [6]. In this sense, the validity of the non-equilibrium effects in turbulent mixing has been already confirmed in a simple geometrical flow.

4. Non-Equilibrium Effects in Jet and Plume Experiments

As was referred to in Section 1, in the buoyancy-dominated turbulent flows, the turbulent convection is characterized by large-scale plumes. As this consequence, the turbulent fluxes are deviated from the usual gradient-transport approximation. It has been recognized that entrainment dominates the evolution of buoyant plumes [11]. As was shown in a previous study of a supernova convection [21], global turbulence models introducing the entrainment parameters provide the best results for describing the evolution of buoyant plumes. However, such models highly depend on the geometrical configurations of the simulation, and are difficult to apply to a general 3D situation. In this sense, it is preferable to develop a local model without resorting to global parameters, provided that the local model works well enough for reproducing the basic behavior of the convection turbulence with plumes. In order to validate the model, here we consider turbulent flows with much simpler geometrical configurations, namely jets and plumes, utilizing recent experimental results.
The non-equilibrium effect may help develop such a local model. In the presence of a plume, we can expect that the turbulence properties vary in the plume-flow direction. The non-equilibrium effect associated with the plume motion can be incorporated through the variation in the turbulent energy K and its dissipation rate ε along the plume motion in terms of the Lagrangian or advective derivatives, D K / D t = [ ( / t ) + U · ] K and D ε / D t = [ ( / t ) + U · ] ε (Figure 2). The effect of non-equilibrium turbulence can be implemented into the model of turbulent viscosity as (30) with the equilibrium turbulent viscosity given by (31).
Recently, a series of detailed experimental studies of the budgets of turbulent energy, Reynolds stresses and dissipation rate has been performed in single- and multi-phase jets/plumes, including a turbulent round jet discharged into an initially stationary ambient [27], a buoyant multi-phase bubble plume [28], a negatively buoyant multi-phase particle plume [25] and a variable density round jet with co-flow [26].

4.1. Round Jets

Using a set of stereoscopic particle image velocimetry (SPIV) measurements of a turbulent round water jet, the budget of the turbulent energy equation was investigated in Lai and Socolofsky [27]. In this experiment, it was reported that the magnitudes of the normal Reynolds stresses R i i decreased by 30 % over 27 x / D 37 in the downstream direction (D: exit diameter of the jet nozzle). In particular, the streamwise variation of the axial velocity fluctuation ( u x ) 2 / ε at the jet centerline shows a decreasing tendency at the major region of the axial locations 31 x / D 37 (Table 1). Although in Lai and Socolofsky [27], these experimental results are argued in the context of the approximate achievement of the self-preserving state in the jet [36]; this decrease tendency may be interpreted as a deviation from the self-preserving state due to the non-equilibrium effect. Since the magnitude of axial velocity fluctuation ( u x ) 2 is similar to and just slightly larger than the horizontal counterparts of the velocity fluctuations, ( u y ) 2 and ( u z ) 2 , and since the energy dissipation rate ε is well represented by that of the isotropic turbulence, ε iso , the evolution of squared axial velocity fluctuation, ( u x ) 2 / ε iso , provides a reasonable surrogate for the basic behavior of K 2 / ε :
K 2 ε ( u x ) 2 ε iso 2 .
The decrease in ( u x ) 2 / ε iso in the downstream direction indicates that the sign of the non-equilibrium effect term is negative as
D D t K 2 ε < 0 ( in turbulent round jets ) .
It follows from (30) that the turbulent viscosity is enhanced by the non-equilibrium effect as
ν NE = ν E 1 C N 1 K D D t K 2 ε > ν E .
In order to evaluate the non-equilibrium effect, we introduced the non-equilibrium correction factor Λ by
Λ = C N 1 K D D t K 2 ε ,
where C N is the model constant estimated as the order of unity ( C N = 0.8 ) through applications to the homogeneous-shear flow turbulence [6]. This factor is approximated in terms of the Lagrangian derivative as
Λ C N 1 K ( U · ) K 2 ε .
Using the data of Lai and Socolofsky [27], Λ is estimated as
Λ C N U c ( u x ) 2 2 1 U c Δ { [ ( u x ) 2 / ε iso ] 2 } D Δ ( x / D ) 1.7
with C N = 1 (Table 1). Here, Δ { [ ( u x ) 2 / ε iso ] 2 } in the nominator and Δ ( x / D ) in the denominator are the increments in [ ( u x ) 2 / ε iso ] 2 and x / D along the jet streamwise direction.
In Lai and Socolofsky [27], the budgets of the turbulent energy and dissipation rate equations were also investigated. Using the balance in the budgets, the coefficients of each term of the K ε model were examined. They suggested the befitting directions for the change in model coefficients. The eddy viscosity constant C ν was increased from the standard value of C ν = 0.09 to 0.09–0.27, whereas a reduction of C ν to 0.07 is suggested from the matching for the jet spreading rate. This implies that it is inadequate to use a constant model coefficient C ν for the turbulent viscosity. The enhancement of eddy viscosity due to the non-equilibrium effect is basically in the correct direction to the evaluation of C ν , and the estimated value of Λ (37) is consistent with the suggested change of C ν . The degree of enhancement/suppression should be determined by the non-equilibrium property of the turbulence, which depends on the spatial location.

4.2. Buoyant Bubble Plumes

In Lai and Socolofsky [28], the turbulent energy budgets in bubble plume were experimentally investigated. There, the turbulent energy evolution was measured both in the adjustment phase of the plume dependent on the source conditions (source size and geometry) and in the asymptotic phase independent of the source condition.
These experimental results show that, in the asymptotic phase (Case A), the turbulent energy K, its dissipation rate ε and K 2 / ε in the plume core region ( r b g ) decrease in the axial or downstream direction along the plume velocity, while in the adjustment phase (Case B) K, ε and K 2 / ε increase (Table 2). However, the relative variation in the turbulent kinetic energy at the downstream to the upstream is much higher in the adjustment phase (Case B) than in the asymptotic phase (Case A). This means that the turbulent energy inhomogeneity along the plume flow, which may lead to the non-equilibrium effect is more prominent in the adjustment phase than in the asymptotic phase. This is natural since the non-equilibrium property is expected to be much higher in the adjustment phase than in the asymptotic phase. In the adjustment phase, the Lagrangian derivative of K 2 / ε is positive as
D D t K 2 ε > 0 ( in the adjustment phase of buoyant bubble plume ) .
In this case, from (30)
ν NE = ν E 1 + Λ 1 < ν E ( in the adjustment phase of buoyant bubble plume ) .
In the experiment of Lai and Socolofsky [28], Λ is estimated as
Λ C N 1 K W c Δ ( K 2 / ε ) D Δ ( z / D ) 0.20
with C N = 1 (Table 2). Here, Δ ( K 2 / ε ) and Δ ( z / D ) are the increments of K 2 / ε and z / D along the vertical streamwise height, respectively. This suggests that due to the non-equilibrium effect, the turbulent mixing in the buoyant bubble plumes is expected to be suppressed by a few 10% as compared to the standard equilibrium turbulent transport.

4.3. Variable Density Jets

There are several studies on the buoyancy-driven variable density turbulence with jets and plumes [26,37]. In recent experiments by Charonko and Prestridge [26], variable density effects in turbulent round jets with co-flow at high- and low-density ratios are experimentally investigated. There, the axial and radial mixing mechanisms of jet are compared between the cases with low and high variable density contrast. It was confirmed that the Reynolds stress is suppressed by the turbulent mass flux in the increased variable density case. It was also reported that the radial transport of momentum and the energy cascade down to small scales are significantly suppressed in the variable density case. These suppressions mostly take place in locations where the density fluctuation is large.
At the upstream region ( x 1 / d 0 < 16.3 , x 1 : axial location; d 0 : inner diameter), the turbulent energy—regardless of both whether averaging is mass-weighted or the Reynolds and whether turbulent energy is normalized by excess momentum or the initial kinetic energy of the mean flow—develops with a downstream distance along the centerline as
D K D t > 0 ( at upstream region of variable density jets ) .
This increase in turbulent energy along the mean flow takes place in both the small- and large-density variation cases. Since the evolution of the dissipation rate ε with the axial distance is not explicitly presented in Charonko and Prestridge [26], that of K 2 / ε cannot be accurately evaluated. However, the analysis of several flows later in Section 5.1 implies that the trend of the K 2 / ε evolution along the mean flow or the sign of K 2 / ε can be surmised from the evolution of K itself. From the non-equilibrium effect D ( K 2 / ε ) / D t [ ( K / ε ) D K / D t ] in (30), we see that the turbulent transport is suppressed at the upstream region.
At the same time, Charonko and Prestridge [26] showed that the variable density effects decrease both the peak and asymptotic values of turbulent energy. By considering the budget in the evolution equation of turbulent energy, they suggested that in the high variable density case, the density fluctuation contributes to the suppression of energy transfer into small-scale fluctuations and preserving the mean-flow structures further downstream. These features can be well understood from the viewpoint of a density-variance effect; the presence of density variance K ρ ( = ρ 2 ) leads to the suppression of turbulent transport as in (30). Charonko and Prestridge [26] concluded that the variable density effect must be modeled to accurately capture dynamical mixing in jets and other variable density flow phenomena. The combination of the non-equilibrium and density-variance effects proposed in the present work is a candidate for modeling the variable density effect in these flows.

5. Plumes in Stellar Convection Zone

In the previous section, we argued the non-equilibrium effect associated with a single plume and jet. In modeling such a situation, the velocity of the single plume/jet is treated as the local mean velocity, and the non-equilibrium effect associated with the mean velocity can be captured by the Lagrangian or advective derivative D / D t = / t + U · [ U ( u ) : mean velocity]. However, in the practical application of the turbulence model to geo/astrophysical flow phenomena, we often encounter situations where the effect of numerous plumes and jets on the global or averaged flow dynamics must be properly evaluated.
One example is the turbulent transport in the stellar convection zone. In some flow configurations, such as in the convective motion in a closed domain, the mean velocity formulated with a simple ensemble average or averaging in the statistically homogeneous plane is very small because of the statistical cancellations due to symmetry. Furthermore, such an averaged quantity cannot capture localized flow structures such as plumes. Because of the small mean flow ( u U 0 ), the non-equilibrium effect represented by the mean flow advection, ( U · ) K 2 / ε in D K / D t , is not suitable for capturing the plume effects. With a simple ensemble or space-averaging procedure, the direct impact due to the non-equilibrium effect might be negligible. However, the presence of grid-scale (GS) or large-eddy velocity component suggests that the counterpart of the non-equilibrium effect at GS may arise as a result of the alternation of the subgrid-scale (SGS) turbulent transport.

5.1. Non-Equilibrium Effect in the Stellar Convection

Recently, in the context of the stellar convection, the non-equilibrium effects associated with the convective plume motion were investigated in the framework of double averaging [29]. For the purpose of exploring the non-equilibrium effect associated with the plumes, fluid dynamics in the stellar convection zone was investigated with the aid of direct numerical simulations (DNSs). In the setup of the simulations, the diving plumes are driven by the cooling layer at the upper surface of the stellar convection zone. Because of the cooling, the surface layer is convectively unstable.
In the non-locally driven convection case, there exist a lot of downward plumes in the shallow region of the convection zone as in the lower plot of Figure 3b.
In order to capture the dynamics of plume motion, the time–space double-averaging procedure was adopted there. A field quantity f is divided into
f = f ¯ + f ˜ f ¯ f ¯ coherent fluctuation f ¯ + f f f ¯ incoherent fluctuation ,
where f ¯ denotes the time average, f denotes the space or ensemble average and f ˜ ( = f ¯ f ¯ ) denotes the residual of the temporal average subtracted by the ensemble average part. This part may be called the dispersion or coherent fluctuation of f, f ( = f f ¯ ) , the fluctuation around the time average f ¯ . As for the time averaging, the resolved part of a quantity f, f ¯ , is defined by the relation
f ¯ ( x ; t ) = f ( x ; s ) G ( t s ) d s
with the simplest top-hat time filter function G ( t s ) defined by
G ( t s ) = 1 / T ( | t s | T / 2 ) , 0 ( otherwise ) .
Here, T is the averaging time window, which should be put in the range:
τ T Ξ ,
where τ is the eddy turnover time of turbulence and Ξ is the time scale of the mean-field evolution.
In this formalism, the fluctuation around the space or ensemble average, f , is divided into the coherent and incoherent parts as
f = f ˜ + f .
With this double-averaging procedure, the time-averaged velocity u ¯ is divided into the space-averaged velocity u ¯ and the dispersion velocity u ˜ as
u ¯ = u ¯ + u ˜ .
The non-equilibrium effect along the plume motion is represented by the factor
Λ ˜ D = ( u ˜ · ) u 2 ¯ .
This means that, if the turbulent energy is inhomogeneous along the plume motion u ˜ , the turbulent fluxes are enhanced or suppressed depending on the sign of Λ ˜ D .
Realizations of the non-equilibrium factor Λ ˜ D = ( u ˜ · ) u 2 ¯ in the stellar convection simulations are plotted in Figure 4. In this case, statistically, the plume velocity is in the downward direction ( u ˜ z < 0 ) while the turbulent energy u 2 ¯ decreases in the downward direction ( u 2 ¯ / z > 0 ). As a result, the non-equilibrium factor Λ ˜ D is statistically negative as
Λ ˜ D = ( u ˜ · ) u 2 ¯ = u ˜ z u 2 ¯ z < 0 .
The spatial distribution of Λ ˜ D is concentrated in the region near the surface where the diving plume structures are prominently observed. The magnitude of Λ ˜ D shows a much larger value in the case with the smaller averaging time window T defined in (44). This clearly shows that the non-equilibrium effect associated with the plume motion can be captured by the Lagrangian derivative of the turbulent energy along the plume motion.
In the framework of the time–space-averaging procedure, the non-equilibrium effect is incorporated into the turbulent internal-energy flux. On the basis of the eddy viscosity expression (30) in the usual ensemble averaging procedure, the turbulent internal-energy diffusivity κ T with the non-equilibrium effect in the time–space-averaging procedure may be formulated as
κ T = κ TE 1 C ˜ ρ ¯ 1 / 3 Λ ˜ D ,
where κ TE is the equilibrium turbulent internal-energy diffusivity, ρ ¯ is the horizontally averaged density and C ˜ is the model constant. The mean density dependence ρ ¯ 1 / 3 in (50) is obtained from the argument that the coherent dissipation rate ε ˜ is proportional to the cubic root of the buoyancy flux associated with a plume. See Section 6.3 and Appendix C of Yokoi, Masada and Takiwaki [29] for the detailed arguments on this dependence.
The turbulent internal-energy flux e u z obtained from the DNSs in the locally and non-locally driven convection cases, which, respectively, correspond to (a) and (b) of Figure 3, are plotted in Figure 5. We see from the figure that the spatial profiles of these two cases are fairly different. First, the magnitude of e u z is much larger in the non-locally driven convection case. At the same time, turbulent transport is significantly localized near the surface region in the non-locally driven case. These prominent properties of turbulent transport in the non-locally driven convection cannot be properly reproduced by the usual gradient-diffusion model with the mixing-length theory.
The spatial distribution of e u z calculated by the model with the non-equilibrium effect (50):
e u z = κ T E
which is also plotted in Figure 5 (E: mean internal energy). The comparison with the DNS results shows that the spatial distribution of e u z in the non-locally driven convection is well reproduced by the present model with the non-equilibrium effect.

5.2. Interaction between Coherent and Incoherent Fluctuations

In this double-averaging formulation, the fluctuation around the space or ensemble average, f , is divided into the coherent and incoherent components as f = f ˜ + f . These two components are not independent but interact with each other. For instance, let us consider the Reynolds stress in the time–space double-averaging procedure. The evolutions of the Reynolds stress of the coherent velocity fluctuation, u ˜ u ˜ ¯ , are subject to
t + u ¯ · u ˜ i u ˜ j ¯ = P ˜ i j + Π ˜ i j ε ˜ i j + T ˜ i j x + P ˜ i j ,
where P ˜ ( = { P ˜ i j } ) , Π ˜ ( = { Π ˜ i j } ) , ε ˜ ( = { ε ˜ i j } ) and · T ˜ ( = { T ˜ i j / x } ) are the production, re-distribution, dissipation and the transport rates of the coherent components, respectively. They are defined as
P ˜ i j = u ˜ j u ˜ u i x u ˜ i u ˜ u j x ,
Π ˜ i j = 1 ρ 0 p ˜ u ˜ i x j + u ˜ j x i ,
ε ˜ i j 2 ν u ˜ i x u ˜ j x ,
T ˜ i j x = + x u ˜ u ˜ i u ˜ j + p ˜ u ˜ i δ j + p ˜ u ˜ j δ i + ν x u ˜ i u ˜ j .
On the other hand, the evolution of the Reynolds stress of the incoherent velocity fluctuation, u u ¯ , is subject to
t + u ¯ · u i u j ¯ = P i j + Π i j ε i j + T i j x + P i j ,
where P ( = { P i j } ) , Π ( = { Π i j } ) , ε ( = { ε i j } ) and · T ( = { T i j / x } ) are defined as
P i j = u j u u i x u i u u j x ,
Π i j = 1 ρ 0 p u i x j + u j x i 2 ν u i x u j x ,
ε i j = 2 ν u i x u j x ,
T i j x = + x u u i u j + p u i δ j + p u j δ i + ν x u i u j .
The terms (53)–(56) and (58)–(61) are similar to the counterparts of the Reynolds stress equation without using the double-averaging procedure. Actually, for example, the addition of P ˜ i j and P i j just gives the usual production term of u i u j . On the other hand, the final terms in (52) and (57) originate from the double-averaging procedure. These terms, P ˜ i j and P i j , are defined by
P i j = u u j ˜ u ˜ i x u u i ˜ u ˜ j x = P ˜ i j ,
and represent the transfer rates between the coherent and incoherent components. Here, the dispersion part of the incoherent Reynolds stress is defined by
u i u j ˜ = u i u j ¯ u i u j ¯ .
If we add (52) and (57), these two terms cancel each other out, and will not contribute to the budget of the total Reynolds stress u u but to the transfer between the coherent and incoherent components. The dispersion part of the incoherent Reynolds stress u u j ˜ (63) coupled with the spatial structure of the coherent motion, u ˜ j / x , the turbulent stress are transferred between the coherent and incoherent components.
Taking the contraction of i and j in (52) and (57), we obtain the evolution equations of the coherent and incoherent energies as
t + u ¯ · 1 2 u ˜ 2 = P ˜ ε ˜ + · T ˜ + u u m ˜ u ˜ m x ,
t + u ¯ · 1 2 u 2 = P ε + · T u u m ˜ u ˜ m x ,
where P ˜ and P are the production rates, ε ˜ and ε are the dissipation rates, · T ˜ and · T is the transport rate of the coherent and incoherent turbulent energies, respectively. The final terms in (64) and (65) represent the energy transfer rates between the coherent and incoherent turbulent energies. They are denoted by
P K = u u m ˜ u ˜ m x = P K ˜ .
If we approximate the dispersion part of the Reynolds stress by the eddy viscosity-type model by
u u m ˜ ν ˜ u ˜ m x ,
the energy transfer rate from the coherent to incoherent fluctuations, P K , is expressed as
P K = u u m ˜ u ˜ m x + ν ˜ u ˜ m x 2 > 0 .
As (66) shows, the sink of the coherent fluctuation energy, P K ˜ < 0 , corresponds to the production of the incoherent fluctuation energy, P K > 0 . In this case, the kinetic energy of the plume or coherent fluctuation motion driven by the surface cooling is transferred towards the kinetic energy of the random or incoherent fluctuation motions as schematically depicted in Figure 6.
In this picture, the diving plumes driven by the surface cooling enhances the turbulent transport through the non-equilibrium effect along the plume motions. During this process, the energy of plumes (coherent fluctuations) is transferred to the energy of random noises (incoherent fluctuations) through the inhomogeneous motions of plumes u ˜ coupled with the incoherent fluctuations u u ˜ . Of course, there are several other mechanisms such as the production P , dissipation ε , etc. in (65), that contribute to the evolution of the incoherent fluctuation energy u 2 . In this sense, the present picture is a fairly simplified one.
With the non-equilibrium effect, the energy transfer rate to the incoherent component of the fluctuation is written as
P K = u u m ˜ u ˜ m x + ν ˜ 1 Λ ˜ D u ˜ m x 2 .
As the DNSs of the surface cooling case show (Figure 4), the non-equilibrium-effect coefficient Λ ˜ D is negative:
Λ ˜ D < 0 .
In this case, the energy transfer to the incoherent or random fluctuation is enhanced by the non-equilibrium turbulence effect along the plume motion. This matches the tendency of the surface cooling-driven convection case. The presence of the diving plumes strongly enhances the turbulence mixing through the non-equilibrium effect. This picture is based on several simplifications, but still suggests the importance of the non-equilibrium effect associated with the plume motions in the turbulent transport in convection.

6. Conclusions

Linear gradient-transport-type turbulence models for the turbulent fluxes with the eddy viscosity and diffusivity representation work poorly in some turbulent flows. For example, the K ε turbulence model based on the eddy viscosity ν T = C ν K 2 / ε cannot be applied to the homogeneous-shear turbulence; it leads to the overestimation of K and ε . Convective turbulence with plumes/jets and supersonic turbulence accompanying shocks constitute other examples. Simulations with a simple linear model provide substantially over- or underestimated turbulent transports (turbulent viscosity, mass flux, heat flux, etc.). To rectify this deficiency of linear gradient-transport models, in this work, the non-equilibrium effect mediated by the streamwise variation in turbulence are incorporated into the turbulence model. A turbulence model with these effects implemented in combination was proposed. In the model, in addition to the usual eddy viscosity and diffusivity representations for the turbulent transports, deviations from the simple representations due to the non-equilibrium effect are taken into account.
Relevance of the non-equilibrium effects was argued with utilizing recent elaborated experiments in turbulent jets and plumes. It was pointed out that the expected results of the non-equilibrium effect on turbulent transports in the turbulent jets and plumes are transport enhancement or suppression due to the non-equilibrium effect (depending on the turbulence variation along the jets/plumes). These evaluations of turbulent transports seem to be in preferable directions to correct the discrepant results obtained from simple conventional gradient-transport-type models.

Funding

This research was funded by the Japan Society for the Promotion of Science (JSPS) Grants-in-Aid for Scientific Research: JP18H01212 and JP23H01199. The Isaac Newton Institute (INI) for Mathematical Sciences program: “Frontier in dynamo theory: From the Earth to the stars” (DYT2) was supported by EPSRC grant no EP/R014604/1. The Nordic Institute for Theoretical Physics (NORDITA) program: “Toward a comprehensive model for the galactic magnetic field” was partly supported by NordForsk and the Royal Astronomical Society.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data from the numerical simulations and analysis presented in this paper are available from the author upon reasonable request.

Acknowledgments

The author would like to thank Youhei Masada and Tomoya Takiwaki for providing data from numerical simulations and valuable discussions on density variance and non-equilibrium effects in stellar convection simulations. Part of this work was started during the author’s stay at the INI, Cambridge, during the DYT2 program from September to December 2022. The author also would like to thank the NORDITA for support and hospitality during the program: Towards a comprehensive model for the galactic magnetic field, in April 2023. The author would like to cordially dedicate this work to the memory of Jack Herring, who had kept encouraging the author in working turbulence closure theory and convection modeling.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Non-Equilibrium Turbulent Energy from the Two-Scale Direct-Interaction Approximation Analysis

From the two-scale direct-interaction approximation (TSDIA) analysis, the turbulent energy K ( u 2 / 2 ) is known to be expressed in terms of as [30]
K = d k Q ( k ; τ , τ ) d k τ d τ 1 G ( k , X ; τ , τ 1 ) D Q ( k , X ; τ , τ 1 ) D T ,
where Q ( k , X ; τ , τ , T ) and G ( k , X ; τ , τ , T ) are the energy spectral function and the response function defined by
Q ( k , X ; τ , τ , T ) = σ ( k , X ; T ) exp ω ( k , X ; T ) | τ τ | ,
G ( k , X ; τ , τ , T ) = H ( k , τ τ ) exp ω ( k , X ; T ) ( τ τ ) ,
respectively. Here the Heaviside step function defined by
H ( x ) = 1 for x 0 , 0 for x < 0 ,
and
σ ( k , X ; T ) = σ 0 ε ( X , T ) 2 / 3 k 11 / 3 ,
ω ( k , X ; T ) = ω 0 ε ( X , T ) 1 / 3 k 2 / 3 .
As we see below, Equation (A5) corresponds to the Kolmogorov scaling of homogeneous and isotropic equilibrium turbulence.
Substituting (A2)–(A6) into Equation (A1), we have
K = d k σ ( k , X ; T ) exp [ ω ( k , X ; T ) | τ τ | ] d k τ d τ 1 H ( τ τ 1 ) exp [ ω ( k , X ; T ) ( τ τ 1 ) ] × D D T σ ( k , X ; T ) exp [ ω ( k , X ; T ) | τ τ 1 | ]
The first term of (A7) denoted by K E is evaluated as
K E = d k σ ( k , X ; T ) exp [ ω ( k , X ; T ) | τ τ | ] = 4 π k C d k k 2 σ 0 ε 2 / 3 k 11 / 3 = 4 π σ 0 k C d k ε 2 / 3 k 5 / 3 ,
where
k C = 2 π / C
is the infrared cut-off wave number representing the largest eddy size of turbulence, C . We further introduce the scale transformation
s = k / k C
to express the relative scale based on the infrared cut-off wave number. With this length scale, the equilibrium energy part can be expressed as
K E = 4 π σ 0 ε 2 / 3 k C s 1 d s s 5 / 3 = 6 π σ 0 ε 2 / 3 k C 2 / 3 = 3 ( 2 π ) 1 / 3 σ 0 ε 2 / 3 C 2 / 3 .
The second term in (A7) represents the non-equilibrium effect on the turbulent energy. This part is evaluated as
K N = d k τ d τ 1 H ( τ τ 1 ) exp [ ω ( k , X ; T ) ( τ τ 1 ) ] × D D T σ ( k , X ; T ) exp [ ω ( k , X ; T ) | τ τ 1 | ] = d k τ d τ 1 D σ ( k , X ; T ) D T exp [ 2 ω ( k , X ; T ) ( τ τ 1 ) ] + d k τ d τ 1 ( τ τ 1 ) σ ( k , X ; T ) D ω D T exp [ 2 ω ( k , X ; T ) ( τ τ 1 ) ] = 4 π k 2 d k 1 2 ω D σ D T exp [ 2 ω ( k , X , T ) ( τ τ 1 ) ] τ 1 = τ + 4 π d k k 2 σ 2 ω 2 D ω D T exp [ 2 ω ( k , X ; T ) ( τ τ 1 ) ] τ 1 = τ = 2 π d k k 2 1 ω D σ D T σ 2 ω 2 D ω D T
Here, in integration with respect to τ 1 , we integrated by parts, and used the fact that neither σ nor ω explicitly depend on τ 1 . With Equations (A5) and (A6), (A12) is calculated as
K N = π σ 0 ω 0 k C 4 / 3 ε 2 / 3 D ε D T + π 6 σ 0 ω 0 ε 1 / 3 k C 7 / 3 D k C D T .
It follows from Equation (A9) that the non-equilibrium part of the energy is expressed in terms of ε and C as
K N = 1 2 ( 2 π ) 1 / 3 σ 0 ω 0 ε 2 / 3 C 4 / 3 D ε D T 1 12 ( 2 π ) 4 / 3 σ 0 ω 0 ε 1 / 3 C 1 / 3 D C D T .
With the addition of (A11), the turbulent energy is expressed as
K = 3 ( 2 π ) 1 / 3 σ 0 ε 2 / 3 C 2 / 3 1 2 ( 2 π ) 1 / 3 σ 0 ω 0 ε 2 / 3 C 4 / 3 D ε D T 1 12 ( 2 π ) 4 / 3 σ 0 ω 0 ε 1 / 3 C 1 / 3 D C D T ,
or with model constants C K 1 , C K 2 and C K 3 , the expression for the turbulent kinetic energy is written as
K = C K 1 ε 2 / 3 C 2 / 3 C K 2 ε 3 / 2 C 4 / 3 D ε D t C K 3 ε 1 / 3 C 1 / 3 D C D t ,
where the model constants are defined as
C K 1 = 3 ( 2 π ) 1 / 3 σ 0 , C K 2 = 1 2 ( 2 π ) 1 / 3 σ 0 ω 0 , C K 3 = 1 12 ( 2 π ) 4 / 3 σ 0 ω 0 .
Equation (A16) is Equation (20) in Section 3.

References

  1. Launder, B.E.; Spalding, D.B. Lectures in Mathematical Models of Turbulence; Academic Press: London, UK, 1972. [Google Scholar]
  2. Rodi, W. Turbulence Models and Their Application in Hydraulics: A State-of-the-Art Review, 3rd ed.; Taylor and Francis: London, UK, 2000. [Google Scholar]
  3. Yokoi, N. Turbulence, transport and reconnection. In Topics in Magnetohydrodynamic Topology, Reconnection and Stability Theory: CISM International Centre for Mechanical Sciences, 591; MacTaggart, D., Hillier, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2020; pp. 177–265. [Google Scholar]
  4. Yokoi, N.; Yoshizawa, A. Statistical analysis of the effects of helicity in inhomogeneous turbulence. Phys. Fluids A 1993, 5, 464–477. [Google Scholar] [CrossRef]
  5. Yokoi, N.; Brandenburg, A. Large-scale flow generation by inhomogeneous helicity. Phys. Rev. E 2016, 93, 033125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Yoshizawa, A.; Nisizima, S. A non-equilibrium representation of the turbulent viscosity based on the two-scale turbulence theory. Phys. Fluids A 1993, 5, 3302–3304. [Google Scholar] [CrossRef]
  7. Kraichnan, R.H. The structure of isotropic turbulence at very high Reynolds numbers. J. Fluid Mech. 1959, 5, 497–543. [Google Scholar] [CrossRef] [Green Version]
  8. Edwards, S.F. The statistical dynamics of homogeneous turbulence. J. Fluid Mech. 1964, 18, 239–273. [Google Scholar] [CrossRef]
  9. Herring, J.R. Self-consistent-field approach to turbulence theory. Phys. Fluids 1965, 8, 2219–2225. [Google Scholar] [CrossRef]
  10. Herring, J.R. Self-consistent-field approach to nonstationary turbulence. Phys. Fluids 1966, 9, 2106–2110. [Google Scholar] [CrossRef]
  11. Spruit, H.C. Convection in stellar envelopes: A changing paradigm. Mem. Soc. Astron. Ital. 1997, 68, 397–413. [Google Scholar]
  12. Rast, M.P. Compressible plume dynamics and stability. J. Fluid Mech. 1998, 369, 125–149. [Google Scholar] [CrossRef]
  13. Cossette, J.-F.; Rast, M.P. Supergranulation as the largest buoyant driven convective scale of the Sun. Astrophys. J. Lett. 2016, 829, L17. [Google Scholar] [CrossRef] [Green Version]
  14. Brandenburg, A. Stellar mixing length theory. Astrophys. J. 2016, 832, 6. [Google Scholar] [CrossRef] [Green Version]
  15. Bethe, H.A.; Wilson, J.R. Revival of stalled supernova shock by neutrino heating. Astrophys. J. 1985, 295, 14–23. [Google Scholar] [CrossRef]
  16. Couch, S.M.; Warren, M.L.; O’Connor, E.P. Simulating turbulence-aided neutrino-driven core-collapse supernova explosions in one dimension. Astrophys. J. 2020, 890, 127. [Google Scholar] [CrossRef]
  17. Janka, H.-T.; Müller, E. Neutrino heating, convection, and mechanism of Type-II supernova explosions. Astron. Astrophys. 1996, 306, 167–198. [Google Scholar]
  18. Mabanta, J.W.; Murphy, Q.A. How turbulence enables core-collapse supernova explosion. Astrophys. J. 2018, 856, 22. [Google Scholar] [CrossRef] [Green Version]
  19. Marek, A.; Janka, H.-T. Delayed neutrino-driven supernova explosions aided by the standing accretion-shock instability. Astrophys. J. 2009, 694, 664–696. [Google Scholar] [CrossRef] [Green Version]
  20. Müller, B. A critical assessment of turbulence models for 1D core-collapse supernova simulations. Mon. Not. R. Astron. Soc. 2019, 487, 5304–5323. [Google Scholar] [CrossRef] [Green Version]
  21. Murphy, Q.A.; Meakin, C. A global turbulence model for neutrino-driven convection in core-collapse supernovae. Astrophys. J. 2011, 742, 74. [Google Scholar] [CrossRef]
  22. Layek, G.C.; Sunita, V. Non-Kolmogorov scaling and dissipation laws in planar turbulent plume. Phys. Fluids 2018, 30, 115105. [Google Scholar] [CrossRef]
  23. Cafiero, G.; Vassilicos, J.C. Nonequilibrium turbulence scalings and self-similarity in turbulent planar jets. Proc. R. Soc. Lond. A 2019, 475, 20190038. [Google Scholar]
  24. Sunita, V.; Layek, G.C. Nonequilibrium turbulent dissipation in buoyant axisymmetric plume. Phys. Rev. Fluids 2021, 6, 104602. [Google Scholar] [CrossRef]
  25. Bordoloi, A.D.; Lai, C.C.K.; Clark, L.; Veliz, G.; Variano, E. Turbulence statistics in a negatively buoyant multiphase plume. J. Fluid Mech. 2020, 896, A19. [Google Scholar] [CrossRef]
  26. Charonko, J.J.; Prestridge, K. Variable-density mixing in turbulent jets with coflow. J. Fluid Mech. 2017, 825, 887–921. [Google Scholar] [CrossRef]
  27. Lai, C.C.K.; Socolofsky, S.A. Budgets of turbulent kinetic energy, Reynolds stresses, and dissipation in a turbulent round jet discharged into a stagnant ambient. Environ. Fluid Mech. 2019, 19, 349–377. [Google Scholar] [CrossRef]
  28. Lai, C.C.K.; Socolofsky, S.A. The turbulent kinetic energy budget in a bubble plume. J. Fluid Mech. 2019, 865, 993–1041. [Google Scholar] [CrossRef]
  29. Yokoi, N.; Masada, Y.; Takiwaki, T. Modelling stellar convective transport with plumes—I. Non-equilibrium turbulence effect in double-averaging formulation. Mon. Not. R. Astron. Soc. 2022, 516, 2718–2735. [Google Scholar] [CrossRef]
  30. Yoshizawa, A. Statistical analysis of the deviation of the Reynolds stress from its eddy-viscosity representation. Phys. Fluids 1984, 27, 1377–1387. [Google Scholar] [CrossRef]
  31. Yoshizawa, A. Hydrodynamic and Magnetohydrodynamic Turbulent Flows: Modeling and Statistical Theory; Kluwer: Dordrecht, The Netherlands, 1998; pp. 265–303. [Google Scholar]
  32. Yokoi, N. Cross helicity and related dynamo. Geophys. Astrophys. Fluid Dyn. 2013, 107, 114–184. [Google Scholar] [CrossRef] [Green Version]
  33. Yokoi, N. Electromotive force in strongly compressible magnetohydrodynamic turbulence. J. Plasma Phys. 2018, 84, 735840501. [Google Scholar] [CrossRef]
  34. Yokoi, N. Mass and internal-energy transports in strongly compressible magnetohydrodynamic turbulence. J. Plasma Phys. 2018, 84, 735840603. [Google Scholar] [CrossRef]
  35. Okamoto, M. Theoretical investigation of an eddy-viscosity-type representation of the Reynolds stress. J. Phys. Soc. Jpn. 1994, 63, 2101–2122. [Google Scholar] [CrossRef]
  36. Tennekes, H.; Lumley, J.L. A First Course in Turbulence; The MIT Press: Cambridge, MA, USA, 1972. [Google Scholar]
  37. Livescu, D.; Ristorcelli, J.R. Variable-density mixing in buoyancy-driven turbulence. J. Fluid Mech. 2008, 605, 145–180. [Google Scholar] [CrossRef]
Figure 1. Schematically depicted evolution of the turbulent energy in the simulations of homogeneous-shear flow. •: DNS, : standard K ε model simulation; and ―: non-equilibrium model simulation. Onset of the model simulations is set at a scaled time S t 0 using the DNS data (S: shear rate).
Figure 1. Schematically depicted evolution of the turbulent energy in the simulations of homogeneous-shear flow. •: DNS, : standard K ε model simulation; and ―: non-equilibrium model simulation. Onset of the model simulations is set at a scaled time S t 0 using the DNS data (S: shear rate).
Atmosphere 14 01013 g001
Figure 2. The non-equilibrium effect associated with a plume motion. A plume is thermobuoyantly formed above a heat source. The non-equilibrium effect is taken into account through the Lagrangian or advective derivatives of the turbulent energy K and its dissipation rate ε , D K / D t = ( / t + U · ) K and D ε / D t = ( / t + U · ) ε , where U is the plume velocity.
Figure 2. The non-equilibrium effect associated with a plume motion. A plume is thermobuoyantly formed above a heat source. The non-equilibrium effect is taken into account through the Lagrangian or advective derivatives of the turbulent energy K and its dissipation rate ε , D K / D t = ( / t + U · ) K and D ε / D t = ( / t + U · ) ε , where U is the plume velocity.
Atmosphere 14 01013 g002
Figure 3. Entropy distributions in the direct numerical simulations (DNSs) for the locally driven case (a) and the non-locally driven case (b). The horizontal cross-sections of the entropy fluctuation s ( = s s ) at the top surface (Top). In the non-locally or cooling-driven case, the horizontal extension of the cell structures is much more limited than the counterpart in the locally driven case. The vertical cross-sections of the entropy fluctuation s from the horizontal mean (bottom). In the non-locally driven case, the low-entropy down flow or plume structures produced at the surface are prominent in the upper region.
Figure 3. Entropy distributions in the direct numerical simulations (DNSs) for the locally driven case (a) and the non-locally driven case (b). The horizontal cross-sections of the entropy fluctuation s ( = s s ) at the top surface (Top). In the non-locally or cooling-driven case, the horizontal extension of the cell structures is much more limited than the counterpart in the locally driven case. The vertical cross-sections of the entropy fluctuation s from the horizontal mean (bottom). In the non-locally driven case, the low-entropy down flow or plume structures produced at the surface are prominent in the upper region.
Atmosphere 14 01013 g003
Figure 4. Spatial distributions of the non-equilibrium effect factor Λ ˜ D = ( u ˜ · ) u 2 ¯ . The left axis 1 ( z / d ) is the depth from the surface ( z = d ). The magnitude of the non-equilibrium effect factor | Λ ˜ D | depends on the time averaging window T.
Figure 4. Spatial distributions of the non-equilibrium effect factor Λ ˜ D = ( u ˜ · ) u 2 ¯ . The left axis 1 ( z / d ) is the depth from the surface ( z = d ). The magnitude of the non-equilibrium effect factor | Λ ˜ D | depends on the time averaging window T.
Atmosphere 14 01013 g004
Figure 5. Spatial distributions of the turbulent internal-energy flux e u z obtained by the direct numerical simulations (DNSs) and the model with the non-equilibrium effect. The left axis 1 ( z / d ) is the depth from the surface ( z = d ). The DNS results for the locally driven convection case (⋯), for the non-locally driven convection case ( ) and the result from the model with the non-equilibrium effect for the non-locally driven convection case (―).
Figure 5. Spatial distributions of the turbulent internal-energy flux e u z obtained by the direct numerical simulations (DNSs) and the model with the non-equilibrium effect. The left axis 1 ( z / d ) is the depth from the surface ( z = d ). The DNS results for the locally driven convection case (⋯), for the non-locally driven convection case ( ) and the result from the model with the non-equilibrium effect for the non-locally driven convection case (―).
Atmosphere 14 01013 g005
Figure 6. Schematic picture of the interaction between the coherent and incoherent fluctuations. Plumes (coherent fluctuations, depicted by thick black curved lines) are driven by the surface cooling. The energy of the plume motions are transferred to the energy of the random noise (incoherent fluctuation, depicted by gray circle eddies) by the interaction P K if P K > 0 .
Figure 6. Schematic picture of the interaction between the coherent and incoherent fluctuations. Plumes (coherent fluctuations, depicted by thick black curved lines) are driven by the surface cooling. The energy of the plume motions are transferred to the energy of the random noise (incoherent fluctuation, depicted by gray circle eddies) by the interaction P K if P K > 0 .
Atmosphere 14 01013 g006
Table 1. Streamwise variation of normalized turbulent intensities, turbulent energy and turbulent axial velocity fluctuation normalized by the square root of the dissipation rate at the jet centerline. x / D : jet streamwise length, D: jet nozzle exit diameter ( D = 11 mm ), U c : centerline jet velocity ( U c = 50 mm s 1 ). Data are taken from Figures 7 and 8 of Lai and Socolofsky [27].
Table 1. Streamwise variation of normalized turbulent intensities, turbulent energy and turbulent axial velocity fluctuation normalized by the square root of the dissipation rate at the jet centerline. x / D : jet streamwise length, D: jet nozzle exit diameter ( D = 11 mm ), U c : centerline jet velocity ( U c = 50 mm s 1 ). Data are taken from Figures 7 and 8 of Lai and Socolofsky [27].
x / D ( u x ) 2 / U c ( u y ) 2 / U c ( u z ) 2 / U c K / U c 2 ( u x ) 2 / ε iso
310.260.200.190.1437 3.25 × 10 2
370.260.200.190.1437 2.65 × 10 2
(−20%)
Table 2. Evolutions of turbulent energy and its dissipation rate along the plume flow. Case A is the asymptotic phase and Case B is the adjustment phase. z / D : vertical streamwise height; D: dynamic length scale ( D = 6.8 cm for Case A and D = 20.4 cm for Case B); W c : plume centerline velocity; and b g : Gaussian plume radius. Data are taken from Figures 20 and 25 of Lai and Socolofsky [28].
Table 2. Evolutions of turbulent energy and its dissipation rate along the plume flow. Case A is the asymptotic phase and Case B is the adjustment phase. z / D : vertical streamwise height; D: dynamic length scale ( D = 6.8 cm for Case A and D = 20.4 cm for Case B); W c : plume centerline velocity; and b g : Gaussian plume radius. Data are taken from Figures 20 and 25 of Lai and Socolofsky [28].
Case A (Asymptotic Phase)
z / D W c ( cm s 1 ) b g K / W c 2 ( ε / W c 3 ) b g × 10 3 K ( cm 2 s 2 ) ε ( cm 2 s 3 ) K 2 / ε ( cm 2 s 1 )
6.620.645.460.1350.0457.5164.4251.34
11.017.908.000.1650.08652.8761.6545.34
(−8.1%)(−4.3%)(−11.7%)
Case B (Adjustment Phase)
z / D W c ( cm s 1 ) b g K / W c 2 ( ε / W c 3 ) b g × 10 3 K ( cm 2 s 2 ) ε ( cm 2 s 3 ) K 2 / ε ( cm 2 s 1 )
2.1928.505.270.150.035121.84153.7496.56
3.6625.148.530.2250.0896142.20165.78121.97
(+16.7%)(+7.8%)(+26.3%)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yokoi, N. Non-Equilibrium Turbulent Transport in Convective Plumes Obtained from Closure Theory. Atmosphere 2023, 14, 1013. https://doi.org/10.3390/atmos14061013

AMA Style

Yokoi N. Non-Equilibrium Turbulent Transport in Convective Plumes Obtained from Closure Theory. Atmosphere. 2023; 14(6):1013. https://doi.org/10.3390/atmos14061013

Chicago/Turabian Style

Yokoi, Nobumitsu. 2023. "Non-Equilibrium Turbulent Transport in Convective Plumes Obtained from Closure Theory" Atmosphere 14, no. 6: 1013. https://doi.org/10.3390/atmos14061013

APA Style

Yokoi, N. (2023). Non-Equilibrium Turbulent Transport in Convective Plumes Obtained from Closure Theory. Atmosphere, 14(6), 1013. https://doi.org/10.3390/atmos14061013

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