Scale disparity and spectral transfer in anisotropic numerical turbulence

Estimating the effective Reynolds number in implicit large-eddy simulation

In implicit large-eddy simulation (ILES), energy-containing large scales are resolved, and physics capturing numerics are used to spatially filter out unresolved scales and to implicitly model subgrid scale effects. From an applied perspective, it is highly desirable to estimate a characteristic Reynolds number (Re)-and therefore a relevant effective viscosity-so that the impact of resolution on predicted flow quantities and their macroscopic convergence can usefully be characterized. We argue in favor of obtaining robust Re estimates away from the smallest scales of the simulated flow-where numerically controlled dissipation takes place and propose a theoretical basis and framework to determine such measures. ILES examples include forced turbulence as a steady flow case, the Taylor-Green vortex to address transition and decaying turbulence, and simulations of a laser-driven reshock experiment illustrating a fairly complex turbulence problem of current practical interest.

Nonlocal interactions in hydrodynamic turbulence at high Reynolds numbers: the slow emergence of scaling laws

We analyze the data stemming from a forced incompressible hydrodynamic simulation on a grid of 2048(3) regularly spaced points, with a Taylor Reynolds number of R(lambda) ~ 1300. The forcing is given by the Taylor-Green vortex, which shares similarities with the von Kàrmàn flow used in several laboratory experiments; the computation is run for ten turnover times in the turbulent steady state. At this Reynolds number the anisotropic large scale flow pattern, the inertial range, the bottleneck, and the dissipative range are clearly visible, thus providing a good test case for the study of turbulence as it appears in nature. Triadic interactions, the locality of energy fluxes, and longitudinal structure functions of the velocity increments are computed. A comparison with runs at lower Reynolds numbers is performed and shows the emergence of scaling laws for the relative amplitude of local and nonlocal interactions in spectral space. Furthermore, the scaling of the Kolmogorov constant, and of skewness and flatness of velocity increments is consistent with previous experimental results. The accumulation of energy in the small scales associated with the bottleneck seems to occur on a span of wave numbers that is independent of the Reynolds number, possibly ruling out an inertial range explanation for it. Finally, intermittency exponents seem to depart from standard models at high R(lambda), leaving the interpretation of intermittency an open problem.

Large-scale flow effects, energy transfer, and self-similarity on turbulence

The effect of large scales on the statistics and dynamics of turbulent fluctuations is studied using data from high resolution direct numerical simulations. Three different kinds of forcing, and spatial resolutions ranging from 256(3) to 1024(3), are being used. The study is carried out by investigating the nonlinear triadic interactions in Fourier space, transfer functions, structure functions, and probability density functions. Our results show that the large scale flow plays an important role in the development and the statistical properties of the small scale turbulence. The role of helicity is also investigated. We discuss the link between these findings and intermittency, deviations from universality, and possible origins of the bottleneck effect. Finally, we briefly describe the consequences of our results for the subgrid modeling of turbulent flows.

Nonlocal modulation of the energy cascade in broadband-forced turbulence

Classically, large-scale forced turbulence is characterized by a transfer of energy from large to small scales via nonlinear interactions. We have investigated the changes in this energy transfer process in broadband forced turbulence where an additional perturbation of flow at smaller scales is introduced. The modulation of the energy dynamics via the introduction of forcing at smaller scales occurs not only in the forced region but also in a broad range of length scales outside the forced bands due to nonlocal triad interactions. Broadband forcing changes the energy distribution and energy transfer function in a characteristic manner leading to a significant modulation of the turbulence. We studied the changes in this transfer of energy when changing the strength and location of the small-scale forcing support. The energy content in the larger scales was observed to decrease, while the energy transport power for scales in between the large and small scale forcing regions was enhanced. This was investigated further in terms of the detailed transfer function between the triad contributions and observing the long-time statistics of the flow. The energy is transferred toward smaller scales not only by wave numbers of similar size as in the case of large-scale forced turbulence, but by a much wider extent of scales that can be externally controlled.

Shell-to-shell energy transfer in magnetohydrodynamics. I. Steady state turbulence

We investigate the transfer of energy from large scales to small scales in fully developed forced three-dimensional magnetohydrodynamics (MHD) turbulence by analyzing the results of direct numerical simulations in the absence of an externally imposed uniform magnetic field. Our results show that the transfer of kinetic energy from large scales to kinetic energy at smaller scales and the transfer of magnetic energy from large scales to magnetic energy at smaller scales are local, as is also found in the case of neutral fluids and in a way that is compatible with the Kolmogorov theory of turbulence. However, the transfer of energy from the velocity field to the magnetic field is a highly nonlocal process in Fourier space. Energy from the velocity field at large scales can be transferred directly into small-scale magnetic fields without the participation of intermediate scales. Some implications of our results to MHD turbulence modeling are also discussed.

Onset of turbulence in accelerated high-Reynolds-number flow

A new criterion, flow drive time, is identified here as a necessary condition for transition to turbulence in accelerated, unsteady flows. Compressible, high-Reynolds-number flows initiated, for example, in shock tubes, supersonic wind tunnels with practical limitations on dimensions or reservoir capacity, and high energy density pulsed laser target vaporization experimental facilities may not provide flow duration adequate for turbulence development. In addition, for critical periods of the overall flow development, the driving background flow is often unsteady in the experiments as well as in the physical flow situations they are designed to mimic. In these situations transition to fully developed turbulence may not be realized despite achievement of flow Reynolds numbers associated with or exceeding stationary flow transitional criteria. Basically our transitional criterion and prediction procedure extends to accelerated, unsteady background flow situations the remarkably universal mixing transition criterion proposed by Dimotakis [P. E. Dimotakis, J. Fluid Mech. 409, 69 (2000)] for stationary flows. This provides a basis for the requisite space and time scaling. The emphasis here is placed on variable density flow instabilities initiated by constant acceleration Rayleigh-Taylor instability (RTI) or impulsive (shock) acceleration Richtmyer-Meshkov instability (RMI) or combinations of both. The significant influences of compressibility on these developing transitional flows are discussed with their implications on the procedural model development. A fresh perspective for predictive modeling and design of experiments for the instability growth and turbulent mixing transitional interval is provided using an analogy between the well-established buoyancy-drag model with applications of a hierarchy of single point turbulent transport closure models. Experimental comparisons with the procedural results are presented where use is made of three distinctly different types of acceleration driven instability experiments: (1) classical, relatively low speed, constant acceleration RTI experiments; (2) shock tube, shockwave driven RMI flow mixing experiments; (3) laser target vaporization RTI and RMI mixing experiments driven at very high energy density. These last named experiments are of special interest as they provide scaleable flow conditions simulating those of astrophysical magnitude such as shock-driven hydrodynamic mixing in supernova evolution research.

Energy transfer in Rayleigh-Taylor instability

The spatial structure and energy budget for Rayleigh-Taylor instability are examined using results from a 512 x 512 x 2040 point direct numerical simulation. The outer-scale Reynolds number of the flow follows a rough t(3) power law and reaches a final value of about 5500. Taylor microscales and Reynolds numbers are plotted to characterize anisotropy in the flow and document progress towards the mixing transition. A mixing parameter is defined which characterizes the relative rates of entrainment and mixing in the flow. The spectrum of each term in the kinetic energy equation is plotted, at regular time intervals, as a function of the inhomogeneous direction and the two-dimensional wave number for the homogeneous directions. The energy spectrum manifests the beginning of an inertial range by the latter stages of the simulation. The production and dissipation spectra become increasingly opposite and separate in wave space as the flow evolves. The transfer spectrum depends strongly on the inhomogeneous direction, with the net transfer being from large to small scales. Energy transfer at the bubble/spike fronts is strictly positive. Extensive cancellation occurs between the pressure and advection terms. The dilatation term produces negligible energy transfer, but its overall effect is to move energy from high to low density regions.

Nature of the energy transfer process in compressible turbulence

Using a two-point closure theory, the eddy-damped-quasinormal-Markovian approximation, we have investigated the energy transfer process and triadic interactions of compressible turbulence. In order to analyze the compressible mode directly, the Helmholtz decomposition is used. The following issues were addressed: (1) What is the mechanism of energy exchange between the solenoidal and compressible modes, and (2) is there an energy cascade in the compressible energy transfer process? It is concluded that the compressible energy is transferred locally from the solenoidal part to the compressible part. It is also found that there is an energy cascade of the compressible mode for high turbulent Mach number. Since we assume that the compressibility is weak, the magnitude of the compressible (radiative or cascade) transfer is much smaller than that of the solenoidal cascade. These results are further confirmed by studying the triadic energy transfer function, the most fundamental building block of the energy transfer.