Energy Spectrum of Turbulent Velocity Pulsations at Arbitrary Values of Fluid Viscosity

Edward V. Teodorovich

Ishlinsky Institute for Problems in Mechanics, Russian Academy of Sciences, Moscow, Russia.

DOI: 10.4236/jmp.2020.1110092                

Abstract

The problem of calculating the energy spectrum of turbulent velocity pulsations in the case of homogeneous isotropic and stationary turbulence is considered. The domain of turbulent energy production is treated as “a black box” on which boundary the spectral energy flux is given. It is assumed that the spectrum is formatted due to intermodal interactions being local in the wave-number space that leads to a cascade mechanism of energy transfer along the wave-number spectrum and the possibility of using the renormalization-group method related to the Markovian features of the process under consideration. The obtained formula for energy spectrum is valid in a wide wave-number range and at arbitrary values of fluid viscosity. It is shown that in functional formulation of the statistical theory of turbulence, the procedure of separating local intermodal interactions, which govern energy transfer (straining effect), and filtering out nonlocal interactions, which have no influence on energy transfer (sweeping effect), is directly described without providing additional arguments or conjectures commonly used in the renormalization-group analysis of turbulent spectra.

Keywords

Spectral Flux, Spectral Energy Density, Local and Nonlocal Interactions, Inertial Term, Renormalization-Group Invariance, Epsilon-Expansion Procedure

1. Introduction

The problem of calculating the energy distribution over wave-numbers of turbulent fluid (spectral energy density) is a subject of numerous investigations. Even within the framework of the simplest model of homogeneous isotropic and stationary turbulence, solving this problem is complicated by the fact that the Navier-Stokes equations, which describe a turbulent fluid, are nonlinear and for obtaining a closed set of equations for energy one needs to express the quantity of third-order statistical moment in terms of the second-order statistical moment with the help of one or other hypotheses of closing (as a survey see [1]) or to construct a solution in the form of perturbation theory series with subsequent term-by-term averaging the series obtained [2]. In the last case, the technique of Feynman diagrams is commonly applied. The specific feature of the system under consideration consists in the fact that in a fully developed turbulence, a large number of various mode scales are excited and the effect of modes of all scales appears to be essential and has to be accounted for since a single act of intermodal interactions is only a link in a long cascade chain via the mechanism of energy transfer from the range of large-scale modes, where the turbulent energy of stochastic fluid velocity pulsations due to a development of instability of large-scale flows is produced, into the range of small-scale modes where the energy dissipates due to fluid viscosity.

Somewhat different approach to finding the spectral energy density beyond the scope of explicit applying the Navier-Stokes equations was proposed by A. N. Kolmogorov [3] who postulated that the energy transfer along the wave-number spectrum goes due to nonlinear intermodal interactions between the modes of close scales, whereas the interaction between the modes of essentially distinguished scales is realized through the cascade sequence of acts relevant to intermodal interactions between modes of intermediate scales (the Richardson-Kolmogorov cascade); in other words, it has a place “a locality in the wave-number space” of intermodal interactions that form a cascade process of energy transfer over wave-number spectrum. The question of a locality nature was discussed in detail in the surveys devoted to application of the renormalization-group method in turbulence theory [4] [5]. The locality relates to the fact that intermodal interaction between the modes with essentially different scales reduces to a primitive transfer of small-scale modes by large-scale ones without energy redistribution between modes (sweeping effect) [6]. In connection with this, when studying the spectrum form it arises the problem of selecting weak local interactions (dynamic interactions) forming the energy spectrum and acting against the ground of strong nonlocal (kinematic) interactions [7]. In the author’s paper [8] it was claimed and argued the statement that in the turbulence theory with applying the RG-method [4] [5] use of the ε-expansion procedure, well-known in the theory of critical phenomena, is a way to select the local intermodal interactions and filtering out the effects of non-local (distant) interactions.

Kolmogorov proposed to divide the wave-number spectrum into three parts.

1) The range of turbulent energy production in the domain of small wave numbers $k<k_g$ where the energy of turbulent pulsations is generated due to development of instability of large-scale flows. The energy production is simulated by the action of an external random force that is similar to the Langevin force used in the theory of random processes. The energy production range lies beyond the scope of our analyses and in our model it will be treated as “a black box” when one knows nothing concerning the processes in the box and only knows conditions on its boundary (spherical surface of radius $k_g$ in the wave-number space). As a boundary condition it is taken the spectral flux through the boundary surface $W\left(k_g\right)$, which is equal to the rate of energy pumping $\mathcal{E}$ received from the range $k<k_g$.

2) The inertial range is a domain of wave numbers where there is no energy production and the dissipation effects are negligible. In this range the spectral flux remains constant, and it goes the process of energy transfer from the range of small wave numbers into the range of large wave numbers via the cascade sequence of local intermodal interaction acts. In the case of high Reynolds numbers (very small fluid viscosity) the inertial range has a sufficiently long length. In the inertial range it takes a place the Kolmogorov formula for spectral energy density

$E\left(k\right)=C_K{\mathcal{E}}^{2/3}{k}^{-5/3},$ (1.1)

where $C_K$ is the Kolmogorov constant. This formula has been obtained only on the basis of dimensionality arguments without reference to the Navier-Stokes equations.

3) The dissipation range relates to the case when the Reynolds number is not high and dissipation effects are not neglectable. As the Reynolds number decreases (the fluid viscosity grows) the domain of inertial range existence tends to zero and the Kolmogorov formula appears to be inapplicable. Below we consider the problem of building the model that is true beyond the inertial range and obtain the formula for spectral energy distribution $E\left(k\right)$ with account for viscosity that is valid for all wave-numbers with the exception of ones from the energy production range.

2. The Problem of Calculating Energy Spectrum

If the fluid viscosity is accounted for, the spectral flux will depend on wave number and due to the locality of intermodal interactions the spectral energy density $E\left(k\right)$ will be determined by the value of spectral flux at given wave-number $W\left(k\right)$. This fact has to be accounted for when applying the dimensionality arguments. As the result, the formula for $E\left(k\right)$ may be written in the form of “generalized Kolmogorov formula”

$E\left(k\right)=C\left(k\right)W{\left(k\right)}^{2/3}{k}^{-5/3}$ (2.1)

Here $C\left(k\right)$ is a dimensionless function of dimensionless variables $W\left(k\right)/W\left(k_g\right)$, $C\left(k_g\right)$, $k/k_g$ and $H\left(k\right)=W\left(k\right)/{\nu }^3{k}^4$ ; $\nu$ is the fluid viscosity and $k_g$ is the upper boundary of energy production range.

Beyond the energy production range $k>k_d$ the equation of energy balance has the form

$\frac{\text{d}W\left(k\right)}{\text{d}k}\mathrm{<mo>\; =\; </mo>}-2\nu k^{2}E\left(k\right)$ (2.2)

that relates two unknown quantities $W\left(k\right)$ and $E\left(k\right)$. For solving the problem one may use the energy balance Equation (2.2) in combination with the “generalized Kolmogorov formula” (2.1), but in this case it remains one more unknown function $C\left(k\right)$, which is a functional analog of the Kolmogorov constant. Thus it arises the need to have one more equation followed from some additional considerations or hypotheses. One of them was proposed by Kovazhny [9], whose conjecture appears to be equivalent to the statement $C\left(k\right)=\mathrm{const}$. In this approximation Equation (2.2) gives

$W\left(k\right)=\mathcal{E}{\left[1+\frac{\nu C\left(k_g\right)}{2{\mathcal{E}}^{1/3}}\left(k^{4/3}-k_g^{4/3}\right)\right]}^3,\mathcal{E}=W(\; k\; g\; )$

$E\left(k\right)=C\left(k_g\right){\mathcal{E}}^{2/3}{\left[1+\frac{\nu C\left(k_g\right)}{2{\mathcal{E}}^{1/3}}\left(k^{4/3}-k_g^{4/3}\right)\right]}^2{k}^{-5/3}$ (2.3)

According to Equation (2.3) the spectral flux $W\left(k\right)$ decreases as the wave-number grows, however, at a certain value of wave-number

$k_d=k_g{\left[1+\frac{2{\mathcal{E}}^{1/3}}{\nu C\left(k_g\right)k_g^{4/3}}\right]}^{3/4}$

the flux becomes zero and at $k>k_d$ it becomes negative whereas the spectral energy grows that corresponds to the transport of the energy produced by a certain fictitious source from small-scale flows to large-scale ones. This “nonphysical result” points out to the fact that the approximation $C\left(k\right)=\mathrm{const}=C\left(k_g\right)$ is unsatisfactory and needs for a refinement. In particular, in author’s paper [10] it was proposed to treat the result (2.3) as two first terms of the series expansion in fluid viscosity powers of exponent that reproduces the results obtained by Pao [11] (see also [12]). Below it will be shown that this proposal appears to be true.

3. Cascade Mechanism of Turbulent Energy Transfer and Renormalization-Group Method

The required additional equation, which allows one to find $C\left(k\right)$, may be obtained by taking into account the cascade mechanism of energy transfer by local in the wave-number space intermodal interactions.

The locality of intermodal interaction acts, which are treated as links of cascade chain, manifest itself in the fact that there is no a certain scale which stands out of another scales (the equal role in cascade chain of all links with given scale). In this case the characteristics of all links in cascade chain $W\left(k\right)$ and $C\left(k\right)$ are defined only by the characteristics $W_0=W\left(k_0\right)$, $C_0=C\left(k_0\right)$ of the link with $k=k_0$ selected as initial one and are independent of how this link was formatted (the independence of a previous history). This means that the process of energy transfer over wave-number spectrum is a Markovian process. If the energy transfer be a Markovian process, we will have

$W\left(k\right)=W\left(k;W_0,C_0,k_0\right),C\left(k\right)=C\left(k;W_0,C_0,k_0\right)$

And the dimensionality arguments enable one to represent the functions desired in the form of dimensionless functions of dimensionless variables

$W\left(k\right)=W_0{\phi }_1\left(k/k_0,H\left(k_0\right),C_0\right),C\left(k\right)=C_0{\phi }_2\left(k/k_0,H\left(k_0\right),C_0\right),$

$H\left(k\right)=W\left(k\right)/{\nu }^3{k}^4$ (3.1)

From Equation (3.1) it follows that the dimensionless functions ${\phi }_{\mathrm{1,2}}\left(x\mathrm{,}h\mathrm{,}c\right)$ of dimensionless variables $x=k/k_0$, $h=H\left(k_0\right)$ and $c=C_0$ obey the normalization conditions

${\phi }_{1,2}\left(1,h,c\right)=1$ (3.2)

The functions ${\phi }_{\mathrm{1,2}}\left(x\mathrm{,}h\mathrm{,}c\right)$ possess a certain additional symmetry related to an ambiguity in a way of setting the boundary conditions that are reduced to a choice of a certain link with wave-number $k_0$ as an initial link in cascade chain and specifying the characteristics of this link, namely, $W_0$ and $C_0$ ; in what follows the value $k_0$ will be referred to as the normalization point.

If another link with the wave-number $k=k_1$ be taken as initial one and the values of parameters of this link $W\left(k_1\right)=W_1$ and $C\left(k_1\right)=c_1$ be given, the form of functions ${\phi }_{\mathrm{1,2}}$ remains unchanged; this means that the following relationships have to be satisfied

$W\left(k\right)=W_0{\phi }_1\left(\frac{k}{k_0},h,c\right)=W_1{\phi }_1\left(\frac{k}{k_1},h_1,c_1\right)$

$C\left(k\right)=c{\phi }_2\left(\frac{k}{k_0},h,c\right)=c_1{\phi }_2\left(\frac{k}{k_1},h_1,c_1\right)$ (3.3)

Due to the presence of ambiguity in a choice of the value $k_0$ (unit of scale) it follows that the functions ${\phi }_{\mathrm{1,2}}\left(x\mathrm{,}h\mathrm{,}c\right)$ are invariant with respect to the operation of scale transformation $k_0\to k_1$ and relevant change (renormalization) of the parameters $W_0\to W_1=W\left(k_1\right)$ and $c\to c_1=C\left(k_1\right)$.

The totality of above pointed operators of scale transformations in combination with renormalization of governing parameters obeys the group composition law, contains the operators of identical and inverse transformations, i.e. it made up a group called the renormalization group (RG-group), and the invariance of the function forms under RG-transformations referred to as RG-invariance. In a special case when RG-transformations and RG-invariance are related to change in putting initial or boundary conditions (our case) the term “functional self-similarity” is used [13].

Putting $k=k_1$ in Equation (3.3) and using the normalization condition (3.2)

$W_1=W_0{\phi }_1\left(\alpha ,h,c\right),c_1=c{\phi }_2\left(\alpha ,h,c\right),\alpha =k_1/k_0$ (3.4)

we arrive at the functional RG equations

${\phi }_{1,2}\left(x,h,c\right)={\phi }_{1,2}\left(\alpha ,h,c\right){\phi }_{1,2}\left(x,h,c\right)\left(\frac{x}{\alpha },\frac{h{\phi }_1\left(\alpha ,h,c\right)}{{\alpha }^4},c{\phi }_2\left(\alpha ,h,c\right)\right)$ (3.5)

here $\alpha$ is an arbitrary dimensionless parameter.

These equations are similar to the Kolmogorov-Chapman semi-group equations in the theory of Markovian random processes (the Einstein-Smolukhovsky equations in physics of Brownian motion).

By differentiating the functional RG-equations with respect to $\alpha$ and next putting $\alpha =1$ we obtain the RG-differential equations

$r_{1,2}\left(h,c\right){\phi }_{1,2}\left(x,h,c\right)+\left\{-x\frac{\partial }{\partial x}+\left[r_1\left(h,c\right)-4\right]h\frac{\partial }{\partial h}+r_2\left(h,c\right)c\frac{\partial }{\partial c}\right\}{\phi }_{1,2}\left(x,h,c\right)=0$ (3.6)

here

$r_{1,2}\left(h,c\right)={\frac{\partial {\phi }_{1,2}\left(x,h,c\right)}{\partial x}|}_{x=1}$ (3.7)

The functions $r_{\mathrm{1,2}}$ are similar to the operators of infinitesimal transformations in the Lie theory of continuous groups.

In terms of the functions ${\phi }_{\mathrm{1,2}}$ the balance energy Equation (2.2) takes the form

$\frac{\partial {\phi }_1\left(x,h,c\right)}{\partial x}=-2c{h}^{-1/3}{x}^{1/3}{\phi }_1^{2/3}\left(x,h,c\right){\phi }_2\left(x,h,c\right)$ (3.8)

from which it follows

$r_1\left(h\mathrm{,}c\right)=-2c{h}^{-1/3}$ (3.9)

However, the differential RG Equations (3.6) contain a new unknown quantity $r_2\left(h\mathrm{,}c\right)$ determined by the function $C\left(k\right)$ in the generalized Kolmogorov Equation (2.10). Knowledge of this function is necessary for solving the RG differential equation. However, in our analyses it is enough to know this quantity in the lowest perturbation theory approximation, namely, the first-order term of a series expansion in the fluid viscosity. This is in agreement with the procedure of improving the perturbation theory by applying the RG method proposed by Bogoyubov and Shirkov in the quantum field theory [14].

To the zero-order approximation of perturbation theory it is assigned the case when $C\left(k\right)=\mathrm{const}$, i.e. ${\phi }_2=1$ and $r_2=0$, that reproduce the Kovazhny theory containing “nonphysical singularity”. The first-order approximation contains a term being proportional to fluid viscosity of first degree. An account for dimensionality arguments gives the following representation for the function ${\phi }_2\left(x\mathrm{,}h\mathrm{,}c\right)$

${\phi }_2\left(x\mathrm{,}h\mathrm{,}c\right)=1+\mu c{h}^{-1/3}\left(x^{4/3}-1\right)$

where $\mu$ is a certain dimensionless parameter calculated by using statistical solving the Navier-Stokes equations within the framework of perturbation theory [2] at given external random force. This lies beyond the scope of the model under consideration. From that it follows $r_2\left(h\mathrm{,}c\right)=\mu r_1\left(h\mathrm{,}c\right)$ and we will analyze the relevant solution

${\phi }_1\left(x\mathrm{,}h\mathrm{,}c\right)={\left\{1-\frac{1-3\mu }{2}c{h}^{-1/3}\left(x^{4/3}-1\right)\right\}}^{3/\left(1-3\mu \right)}\mathrm{,}$

${\phi }_2\left(x\mathrm{,}h\mathrm{,}c\right)={\phi }_1^{\mu }\left(x\mathrm{,}h\mathrm{,}c\right)$ (3.10)

From Equation (3.10) one can see that in the range $0<\mu <1/3$ the solution contains “nonphysical” singularity and the position of a singular point goes from $k_d$ to infinity as $\mu$ goes from zero to 1/3. Since the transition to $\mu \to 1/3$ is weakly defined, we consider this case separately by direct solving the differential RG Equation (3.6).

The equation for $ln{\phi }_1\left(x\mathrm{,}h\mathrm{,}c\right)$ can be written in the form

$\left\{-x\frac{\partial }{\partial x}-4h\frac{\partial }{\partial h}+\left[r_1\left(h\mathrm{,}c\right)h\frac{\partial }{\partial h}+r_2\left(h\mathrm{,}c\right)c\frac{\partial }{\partial c}\right]\right\}ln{\phi }_1\left(x\mathrm{,}h\mathrm{,}c\right)=-r_1\left(h\mathrm{,}c\right)$ (3.11)

From the easily verified identity

$\left[r_1\left(h\mathrm{,}c\right)h\frac{\partial }{\partial h}+r_2\left(h\mathrm{,}c\right)c\frac{\partial }{\partial c}\right]r_1\left(h\mathrm{,}c\right)=\left[-r_1\left(h\mathrm{,}c\right)/3+r_2\left(h\mathrm{,}c\right)\right]r_1\left(h\mathrm{,}c\right)$

it follows that the expression in square brackets is zero when $r_2=r_1/3$. Thus one can seek for the solution to Equation (3.10) in the form

$ln{\phi }_1\left(x\mathrm{,}h\mathrm{,}c\right)=r_1\left(h\mathrm{,}c\right)F(\; x\; )$

where the function $F\left(x\right)$ obeys the equation

$1+\left[-x\frac{\text{d}}{\text{d}x}+\frac{4}{3}\right]F\left(x\right)=0,F\left(1\right)=0,{\frac{\text{d}F\left(x\right)}{\text{d}x}|}_{x=1}=1$

that gives

$F\left(x\right)=\frac{3}{4}\left(x^{4/3}-1\right)$

As the result the formulas for spectral flux and spectral energy density take the forms

$W\left(k\right)=\mathcal{E}exp\left\{-\frac{3}{2}\frac{c\nu }{{\mathcal{E}}^{1/3}}\left(k^{4/3}-k_g^{4/3}\right)\right\}$ (3.12)

$E\left(k\right)=c{\mathcal{E}}^{2/3}{k}^{-5/3}exp\left\{-\frac{3}{2}\frac{c\nu }{{\mathcal{E}}^{1/3}}\left(k^{4/3}-k_g^{4/3}\right)\right\}$ (3.13)

Here we choose the wave-number $k_g$, corresponding to the upper boundary of turbulent energy production range, as a normalization point, and use the notation $\mathcal{E}=W\left(k_g\right)$ and $c=C\left(k_g\right)$. In particular case when $k_g=0$ the results (3.12)-(3.13) appear to be identical to formulas proposed by Pao [11] [12], who put forward the hypothesis on proportionality of energy density $E\left(k\right)$ to spectral flux $W\left(k\right)$ without any argumentation.

If the parameter $\mu$ exceeds one third the spectral characteristics monotonically decrease as k tends to infinity and no singularities arise. In this case the formulas for energy flux and spectral energy density take the forms

$W\left(k\right)=\mathcal{E}{\left\{-\frac{3}{2}\frac{c\nu }{{\mathcal{E}}^{1/3}}\left(k^{4/3}-k_g^{4/3}\right)\right\}}^{-3/\left(3\mu -1\right)}$ (3.14)

$E\left(k\right)=c{\mathcal{E}}^{2/3}{k}^{-5/3}{\left\{-\frac{3}{2}\frac{c\nu }{{\mathcal{E}}^{1/3}}\left(k^{4/3}-k_g^{4/3}\right)\right\}}^{-\left(3\mu +2\right)/\left(3\mu -1\right)}$ (3.15)

4. Local and Nonlocal Interactions

Our analysis is based on the statement that the energy spectrum is formed due to intermodal interactions being local in the wave-number space, and this enables one to tell about a cascade mechanism of energy transfer over wave-number spectrum and a Markovian character of the process. Precisely this property of the process gives a possibility to assume that the solution describing the cascade chain depends uniquely on numerical parameters of the link in the chain treated as initial one and is independent of the fact how this link was formatted (independence on previous history). The property of independence on previous history of formatting the initial and boundary conditions (functional self-similarity) is an inherent one for differential equations that do not contain integral terms. Solutions to such equations possess the property of invariance with respect to the way of setting additional (initial and boundary) conditions. Namely, under the shift of the hyper-surface on which additional conditions are given and relevant changing (renormalization) numerical parameters specifying additional conditions the solution remains to be unchanged. The balance energy Equation (2.2) relates to the class of such equations.

However, in the statistical theory of turbulence, the chain of equations that relates statistical moments of various orders arises due to nonlinearity of the Navier-Stokes equation, and this relation necessary contains integral terms. In particular, the equation for the Fourier-transform of the second-order statistic moment of velocity field $B\left(k\mathrm{,}\omega \right)$ has the form

$L^{\left(0\right)}\left(k\mathrm{,}\omega \right)B\left(k\mathrm{,}\omega \right)+T\left(k\mathrm{,}\omega \right)=D^{\left(0\right)}\left(k\mathrm{,}\omega \right)G\left(k\mathrm{,}\omega \right)$ (4.1)

where the second-rank tensor $T\left(k\mathrm{,}\omega \right)$, also referred to as the inertial term, is expressed via the integral of the third-order statistical moment of turbulent velocity pulsations (Equation (A.4) in Appendix). The inertial term describes the processes of momentum and energy redistribution due to mixing induced by velocity pulsations. The relevant intermodal interactions are obviously nonlocal, and the question of validity the form of balance energy Equation (2.2) and an ability to use various symmetry properties like the RG-invariance arises. Thus we arrive at the problem of filtering out these interactions. In the Yakhot-Orszag renormalization-group theory of turbulence [15], this problem is solved by using the ε-expansion procedure when one first calculates the quantity desired in the low-order approximation in ε near the point $\epsilon =0$ and next extends the result obtained into the point $\epsilon =4$. This procedure is similar to the t’Hooft-Veltman dimensional regularization method in quantum field theory.

An account for the effect of turbulent velocity pulsations leads to a correction to fluid viscosity (a turbulent viscosity), and as a result the representation for the Fourier-transform of reverse Green’s function takes the form

$G^{-1}\left(k\mathrm{,}\omega \right)={\left[G^{\left(0\right)}\left(k\mathrm{,}\omega \right)\right]}^{-1}-\Sigma \left(k\mathrm{,}\omega \right)$ (4.2)

${\left[G^{\left(0\right)}\left(k,\omega \right)\right]}^{-1}=L^{\left(0\right)}\left(k,\omega \right)=-i\omega +\nu k^2$ is the reverse Green function of the linear problem. (This equation can be rewritten in the form of Dyson’s equation well-known in quantum-field-theory

$G\left(k,\omega \right)=G^{\left(0\right)}\left(k,\omega \right)+G^{\left(0\right)}\left(k,\omega \right)\Sigma \left(k,\omega \right)G\left(k,\omega \right)$

where $\Sigma \left(k,\omega \right)$ is referred to as the self-energy operator.)

Another result of account for a mixing by turbulent velocity pulsations consists in appearance of effective random force which variance is written as

$D\left(k,\omega \right)=D^{\left(0\right)}\left(k,\omega \right)+D^{\left(1\right)}\left(k,\omega \right)$ (4.3)

Here $D^{\left(0\right)}\left(k,\omega \right)$ is the variance of external random forces that simulate the emergence of stochasticity due to development of instability of large-scale flows and is similar to Langevin forces in the theory of random processes; the second summand $D^{\left(1\right)}\left(k,\omega \right)$ arises due to account for transport phenomena produced by mixing processes.

The quantity $D\left(k,\omega \right)$ enters into the equation for the second-order statistic moment of velocity field $B\left(k,\omega \right)$

$B\left(k,\omega \right)=G\left(k,\omega \right)D\left(k,\omega \right)G\left(-k\mathrm{,}-\omega \right)$ (4.4)

first obtained by Schwinger [16] when building the theory of quantized fields beyond the scope of perturbation methods. The Schwinger approach is based on statistical description of quantized fields in terms of characteristic (generating) functional. In the statistical theory of turbulence this equation was obtained by Wyld [2] with the help of summing up the perturbation theory series and applying the technique of Feynman diagrams.

In the space-time variables the inertial term $T\left(\mathrm{1,2}\right)$ contains the third-order statistical moment of velocity field (see Appendix, Equation (A4)) and the question arises whether the Fourier-transform of inertial term will contain integral terms and how to close the set of equations by excluding the third-order statistical moment.

Within the framework of statistical description of turbulence in terms of characteristic (generating) functional it can be obtained the formula for inertial term

$T\left(k,\omega \right)=-\Sigma \left(k,\omega \right)B\left(k,\omega \right)-G\left(k,\omega \right)D^{\left(1\right)}\left(k,\omega \right)$ (4.5)

firstly pointed out in author’s paper [17]. It should be noted that this formula is an exact one since no approximations or additional conjectures were used in its derivation.

The quantities $\Sigma$ and $D^{\left(1\right)}$ are defined by solving the Navier-Stokes equations, but these quantities can be excluded by using Equation (4.2) and Equation (4.4) that gives

$T\left(k,\omega \right)=-L^{\left(0\right)}\left(k,\omega \right)B\left(k,\omega \right)+G\left(k,\omega \right)D^{\left(0\right)}\left(k,\omega \right)$ (4.6)

For the quantity $T\left(k\right)={\displaystyle \int T\left(k\mathrm{,}\omega \right)\text{d}\omega }/2\pi$ we obtain

$T\left(k\right)=-\nu k^{2}B\left(k\right)+G\left(k\right)D^{\left(0\right)}(\; k\; )$

If to put $T\left(k\right)=\text{d}W\left(k\right)/\text{d}k$, we arrive at the balance energy equation that coincides with Equation (2.2) beyond the energy production range. The above presented procedure is another way of excluding nonlocal intermodal interactions when building the theory of turbulent spectra.

5. Conclusion

The theory of spectral energy distribution is based on the Kolmogorov conjecture that the energy spectrum is formatted by intermodal interactions being local in wave-number space. From this, it follows a cascade mechanism of energy transfer along the wave-number space, the Markovian properties of the process, as well as an ability to apply various similarity arguments such as the property of renormalization-group invariance (functional self-similarity). In this case, the problem of separating local intermodal interactions (straining effects) and filtering out nonlocal ones (sweeping effects) arises. An account for mixing processes by turbulent velocity pulsations (swimming effects) reveals in the form of an addition $D^{\left(1\right)}$ to a variance of external random forces (turbulent random force) and addition Σ to viscous term in the Navier-Stokes equations (turbulent viscosity). In Section 4 it was shown that in the functional formulation of statistical description of turbulence, these quantities prove to be excluded from the balance energy equation. As the result, the problem of selecting local intermodal interactions and filtering out nonlocal ones appears to be solved exactly without applying other methods such as the ε-RG (see more recent survey [18]).

Acknowledgements

The work was carried out within the framework of the State Task (project No AAAA-A17-117021310375-7).

Appendix. Basic Equations in Space-Time Variables

To simplify the writing in formulas we will use digital notation for space-time variables and the index of vector components $\left\{r_1\mathrm{,}{t}_1\mathrm{,}{\alpha }_1\right\}\equiv 1$ according to which $u_{{\alpha }_1}\left(r_1\mathrm{,}{t}_1\right)\equiv u\left(1\right)$. Also it will be implied the integration over space-time coordinates and the summing over component indexes for coinciding digital numbers (the Einstein rule), i.e.

$u\left(1\right)v\left(1\right)\equiv {\displaystyle \int u\left(r_1\mathrm{,}{t}_1\right)\cdot v\left(r_1\mathrm{,}{t}_1\right)\text{d}{r}_1\text{d}{t}_1}$

The Navier-Stokes equations (NSE) with the presence of the external random force $X\left(1\right)$ and the external regular force $f(\; 1\; )$

$L^{\left(0\right)}\left(\mathrm{1,2}\right)u\left(2\right)+\frac{1}{2}V\left(\mathrm{1\; <mo>\; |\; </mo>2,3}\right)u\left(2\right)u\left(3\right)=X\left(1\right)+f\left(1\right)$ (A.1)

(for more details related to notation see [16]).

Variance of external random force

$D^{\left(0\right)}\left(\mathrm{1,2}\right)=\langle X\left(1\right)X\left(2\right)\rangle$ (A.2)

The equation for $B\left(\mathrm{1,2}\right)=\langle u\left(1\right)u\left(2\right)\rangle$

$L^{\left(0\right)}\left(\mathrm{1,}{1}^{\prime }\right)B\left(1^{\prime }\mathrm{,2}\right)+T\left(\mathrm{1,2}\right)=G\left(\mathrm{1,}{1}^{\prime }\right)D^{\left(9\right)}\left(1^{\prime }\mathrm{,2}\right)$ (A.3)

Inertial term

$T\left(\mathrm{1,2}\right)=\frac{1}{2}V\left(\mathrm{1\; <mo>\; |\; </mo>3,4}\right)B\left(\mathrm{3,4,2}\right)\mathrm{,}B\left(\mathrm{1,2,3}\right)=\langle u\left(1\right)u\left(2\right)u\left(3\right)\rangle$ (A.4)

Exact representation for inertial term

$T\left(\mathrm{1,2}\right)=-B\left(\mathrm{1,}{1}^{\prime }\right)\Sigma \left(1^{\prime }\mathrm{,2}\right)-G\left(\mathrm{1,}{1}^{\prime }\right)D^{\left(1\right)}\left(1^{\prime }\mathrm{,2}\right)$ (A.5)

The Schwinger-Wyld formula

$B\left(\mathrm{1,2}\right)=G\left(\mathrm{1,}{1}^{\prime }\right)G\left(\mathrm{2,}{2}^{\prime }\right)D\left(1^{\prime }\mathrm{,}{2}^{\prime }\right)$ (A.6)

$D\left(\mathrm{1,2}\right)$ is the variance of effective random forces.

Another form of Equation (A.6)

$B\left(\mathrm{1,}{1}^{\prime }\right)G^{-1}\left(1^{\prime }\mathrm{,2}\right)=G\left(\mathrm{1,}{1}^{\prime }\right)D\left(1^{\prime }\mathrm{,2}\right)$ (A.7)

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1]	Monin, A.S. and Yaglom, A.M. (1975) Statistical Fluid Mechanics, Vol. 2. MIT Press, Cambridge.
[2]	Wyld, H.W. (1961) Annals of Physics, 14, 143-165. https://doi.org/10.1016/0003-4916(61)90056-2
[3]	Kolmogorov, A.N. (1941) Doklady Akademii Nauk SSSR, 30, 299-303.
[4]	Teodorovich, E.V. (1993) Izvestiya, Russian Academy of Sciences: Atmospheric and Oceanic Physics, 29, 149-163.
[5]	Zhou, Y. and Vahala, G. (1992) Physical Review A, 46, 1136-1139. https://doi.org/10.1103/PhysRevA.46.1136
[6]	Kadomtsev, B.B. (1964) Plasma Turbulence. In: Problems in Theory of Plasmas, Vol. 4, Atomizdat, Moscow, 188. (In Russian)
[7]	Teodorovich, E.V. (1990) Fluid Dynamics, 25, 522-527. https://doi.org/10.1007/BF01049856
[8]	Zakharov, V.E. and L’vov, V.S. (1975) Izvestiya Vuzov. Radiofizika, 28, 1470-1487. https://doi.org/10.1007/BF01040337
[9]	Kovasznay, L.S.G. (1948) Journal of the Aeronautical Sciences, 15, 745-753. https://doi.org/10.2514/8.11707
[10]	Teodorovich, E.V. (2017) Applied Mathematics and Mechanics, 81, 450-454. https://doi.org/10.1016/j.jappmathmech.2018.03.013
[11]	Pao, Y.-H. (1965) Physics of Fluids, 8, 1063-1075. https://doi.org/10.1063/1.1761356
[12]	Leslie, D.C. (1973) Developments in the Theory of Turbulence. Clarendon Press, Oxford.
[13]	Shirkov, D.V. (1982) Soviet Physics—Doklady, 27, 197-199.
[14]	Bogolyubov, N.N. and Shirkov, D.V. (1984) Introduction to the Theory of Quantized Fields. Nauka, Moscow. (In Russian)
[15]	Yakhot, V. and Orszag, S.A. (1986) Journal of Scientific Computing, 1, 3-51. https://doi.org/10.1007/BF01061452
[16]	Schwinger, J. (1949) Physical Review, 75, 651-679. https://doi.org/10.1103/PhysRev.75.651
[17]	Teodorovich, E.V. (2013) Journal of Modern Physics, 4, 56-63. https://doi.org/10.4236/jmp.2013.41010
[18]	Zhou, Y. (2010) Physics Reports, 488, 1-49. https://doi.org/10.1016/j.physrep.2009.04.004