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
is expressed by the characteristic turbulent velocity
v and length scale
ℓ as
. In the standard
K–
model,
is modeled in terms of the turbulent kinetic energy
K and its dissipation rate
as
. 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
,
and
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
of the turbulent energy and its dissipation rate as
where
.
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.
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,
and
(
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
decreased by 30 % over
in the downstream direction (
D: exit diameter of the jet nozzle). In particular, the streamwise variation of the axial velocity fluctuation
at the jet centerline shows a decreasing tendency at the major region of the axial locations
(
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
is similar to and just slightly larger than the horizontal counterparts of the velocity fluctuations,
and
, and since the energy dissipation rate
is well represented by that of the isotropic turbulence,
, the evolution of squared axial velocity fluctuation,
, provides a reasonable surrogate for the basic behavior of
:
The decrease in
in the downstream direction indicates that the sign of the non-equilibrium effect term is negative as
It follows from (
30) that the turbulent viscosity is enhanced by the non-equilibrium effect as
In order to evaluate the non-equilibrium effect, we introduced the non-equilibrium correction factor
by
where
is the model constant estimated as the order of unity (
) through applications to the homogeneous-shear flow turbulence [
6]. This factor is approximated in terms of the Lagrangian derivative as
Using the data of Lai and Socolofsky [
27],
is estimated as
with
(
Table 1). Here,
in the nominator and
in the denominator are the increments in
and
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
was increased from the standard value of
to 0.09–0.27, whereas a reduction of
to
is suggested from the matching for the jet spreading rate. This implies that it is inadequate to use a constant model coefficient
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
, and the estimated value of
(
37) is consistent with the suggested change of
. 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
in the plume core region (
) decrease in the axial or downstream direction along the plume velocity, while in the adjustment phase (Case B)
K,
and
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
is positive as
In the experiment of Lai and Socolofsky [
28],
is estimated as
with
(
Table 2). Here,
and
are the increments of
and
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 (
,
: axial location;
: 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
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
cannot be accurately evaluated. However, the analysis of several flows later in
Section 5.1 implies that the trend of the
evolution along the mean flow or the sign of
can be surmised from the evolution of
K itself. From the non-equilibrium effect
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
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 [: 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 (), the non-equilibrium effect represented by the mean flow advection, in , 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
where
denotes the time average,
denotes the space or ensemble average and
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,
, the fluctuation around the time average
. As for the time averaging, the resolved part of a quantity
f,
, is defined by the relation
with the simplest top-hat time filter function
defined by
Here,
T is the averaging time window, which should be put in the range:
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,
, is divided into the coherent and incoherent parts as
With this double-averaging procedure, the time-averaged velocity
is divided into the space-averaged velocity
and the dispersion velocity
as
The non-equilibrium effect along the plume motion is represented by the factor
This means that, if the turbulent energy is inhomogeneous along the plume motion , the turbulent fluxes are enhanced or suppressed depending on the sign of .
Realizations of the non-equilibrium factor
in the stellar convection simulations are plotted in
Figure 4. In this case, statistically, the plume velocity is in the downward direction (
) while the turbulent energy
decreases in the downward direction (
). As a result, the non-equilibrium factor
is statistically negative as
The spatial distribution of
is concentrated in the region near the surface where the diving plume structures are prominently observed. The magnitude of
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
with the non-equilibrium effect in the time–space-averaging procedure may be formulated as
where
is the equilibrium turbulent internal-energy diffusivity,
is the horizontally averaged density and
is the model constant. The mean density dependence
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
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
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
calculated by the model with the non-equilibrium effect (
50):
which is also plotted in
Figure 5 (
E: mean internal energy). The comparison with the DNS results shows that the spatial distribution of
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,
, is divided into the coherent and incoherent components as
. 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,
, are subject to
where
,
,
and
are the production, re-distribution, dissipation and the transport rates of the coherent components, respectively. They are defined as
On the other hand, the evolution of the Reynolds stress of the incoherent velocity fluctuation,
, is subject to
where
,
,
and
are defined as
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
and
just gives the usual production term of
. On the other hand, the final terms in (
52) and (
57) originate from the double-averaging procedure. These terms,
and
, are defined by
and represent the transfer rates between the coherent and incoherent components. Here, the dispersion part of the incoherent Reynolds stress is defined by
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
but to the transfer between the coherent and incoherent components. The dispersion part of the incoherent Reynolds stress
(
63) coupled with the spatial structure of the coherent motion,
, 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
where
and
are the production rates,
and
are the dissipation rates,
and
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
If we approximate the dispersion part of the Reynolds stress by the eddy viscosity-type model by
the energy transfer rate from the coherent to incoherent fluctuations,
, is expressed as
As (
66) shows, the sink of the coherent fluctuation energy,
, corresponds to the production of the incoherent fluctuation energy,
. 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
coupled with the incoherent fluctuations
. Of course, there are several other mechanisms such as the production
, dissipation
, etc. in (
65), that contribute to the evolution of the incoherent fluctuation energy
. 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
As the DNSs of the surface cooling case show (
Figure 4), the non-equilibrium-effect coefficient
is negative:
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.