Time-ordered fluctuation-dissipation relation for incompressible isotropic turbulence

Statistical Dynamics of Mean Flows Interacting with Rossby Waves, Turbulence, and Topography

Abridged statistical dynamical closures, for the interaction of two-dimensional inhomogeneous turbulent flows with topography and Rossby waves on a beta–plane, are formulated from the Quasi-diagonal Direct Interaction Approximation (QDIA) theory, at various levels of simplification. An abridged QDIA is obtained by replacing the mean field trajectory, from initial-time to current-time, in the time history integrals of the non-Markovian closure by the current-time mean field. Three variants of Markovian Inhomogeneous Closures (MICs) are formulated from the abridged QDIA by using the current-time, prior-time, and correlation fluctuation dissipation theorems. The abridged MICs have auxiliary prognostic equations for relaxation functions that approximate the information in the time history integrals of the QDIA. The abridged MICs are more efficient than the QDIA for long integrations with just two relaxation functions required. The efficacy of the closures is studied in 10-day simulations with an easterly large-scale flow impinging on a conical mountain to generate rapidly growing Rossby waves in a turbulent environment. The abridged closures closely agree with the statistics of large ensembles of direct numerical simulations for the mean and transients. An Eddy Damped Markovian Inhomogeneous Closure (EDMIC), with analytical relaxation functions, which generalizes the Eddy Dampened Quasi Normal Markovian (EDQNM) to inhomogeneous flows, is formulated and shown to be realizable under the same circumstances as the homogeneous EDQNM.

Eulerian field-theoretic closure formalisms for fluid turbulence

The formalisms of Wyld [Ann. Phys. 14, 143 (1961)] and Martin, Siggia, and Rose (MSR) [Phys. Rev. A 8, 423 (1973)] address the closure problem of a statistical treatment of homogeneous isotropic turbulence (HIT) based on techniques primarily developed for quantum field theory. In the Wyld formalism, there is a well-known double-counting problem, for which an ad hoc solution was suggested by Lee [Ann. Phys. 32, 292 (1965)]. We show how to implement this correction in a more natural way from the basic equations of the formalism. This leads to what we call the Improved Wyld-Lee Renormalized Perturbation Theory. MSR had noted that their formalism had more vertex functions than Wyld's formalism and based on this felt Wyld's formalism was incorrect. However a careful comparison of both formalisms here shows that the Wyld formalism follows a different procedure to that of the MSR formalism and so the treatment of vertex corrections appears in different ways in the two formalisms. Taking that into account, along with clarifications made to both formalisms, we find that they are equivalent and we demonstrate this up to fourth order.

Direct-numerical-simulation-based measurement of the mean impulse response of homogeneous isotropic turbulence

A technique for measuring the mean impulse response function of stationary homogeneous isotropic turbulence is proposed. Such a measurement is carried out here on the basis of direct numerical simulation (DNS). A zero-mean white-noise volume forcing is used to probe the turbulent flow, and the response function is obtained by accumulating the space-time correlation between the white forcing and the velocity field. This technique to measure the turbulent response in a DNS numerical experiment is a research tool in that field of spectral closures where the linear-response concept is invoked either by resorting to renormalized perturbations theories or by introducing the well-known fluctuation-dissipation relation (FDR). Although the results obtained in the present work are limited to relatively low values of the Reynolds number, a preliminary analysis is possible. Both the characteristic form and the time scaling properties of the response function are investigated in the universal subrange of dissipative wave numbers; a comparison with the response approximation given by the FDR is proposed through the independent DNS measurement of the correlation function. Very good agreement is found between the measured response and Kraichnan's description of random energy-range advection effects.

Eulerian spectral closures for isotropic turbulence using a time-ordered fluctuation-dissipation relation

Procedures for time-ordering the covariance function, as given in a previous paper [K. Kiyani and W. D. McComb, Phys. Rev. E 70, 066303 (2004)], are extended and used to show that the response function associated at second order with the Kraichnan-Wyld perturbation series can be determined by a local (in wave number) energy balance. These time-ordering procedures also allow the two-time formulation to be reduced to time-independent form by means of exponential approximations and it is verified that the response equation does not have an infrared divergence at infinite Reynolds number. Last, single-time Markovianized closure equations (stated in our previous paper) are derived and shown to be compatible with the Kolmogorov distribution without the need to introduce an ad hoc constant.