banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Forcing for statistically stationary compressible isotropic turbulence
Rent this article for


Image of FIG. 1.
FIG. 1.

Evolution of , which is the dilatational to solenoidal dissipation ratio, for unsplit and split forcing. When a single forcing term is used [Eq. (7)] there is no control over , and it often continues to grow throughout the simulation. With a split solenoidal/dilatational forcing [Eq. (17)] quickly adjusts to the imposed value (solid horizontal lines). For clarity, single forcing term data are averaged over a 5 s window. The light gray is unaveraged data, and shows high variability in the single term forcing results.

Image of FIG. 2.
FIG. 2.

The ratio of dilatational to solenoidal kinetic energy reaches the equilibrium value about the same time as ; however, the specific value is different than and depends on the turbulent Mach number.

Image of FIG. 3.
FIG. 3.

Evolution of the total dissipation for unsplit and split forcing. The forcing coefficients can be cast in a form that allows the specification of the dissipation (solid horizontal lines) and, thus, the Kolmogorov microscale. The split forcing method adheres closely to the imposed dissipation (solid horizontal lines) unlike the single term forcing method.

Image of FIG. 4.
FIG. 4.

Time variation of the turbulent kinetic energy for unsplit and split forcing. The increase in for unsplit forcing causes the turbulent kinetic energy to increase as well.

Image of FIG. 5.
FIG. 5.

The energy flux is equal to the dissipation for a range of wavenumbers for low- forcing but not for full spectrum linear forcing, which requires larger resolutions to develop an inertial range. The results correspond to the parameters from run 1c except as noted.

Image of FIG. 6.
FIG. 6.

The energy content at low wavenumbers is lower when full spectrum linear forcing is used, compared to low- forcing, which results in a lower overall Reynolds number. The results correspond to the parameters from run 1c.

Image of FIG. 7.
FIG. 7.

Compensated kinetic energy spectra obtained from low- split forced simulations at different and values: (a) total kinetic energy, (b) solenoidal, and (c) dilatational parts of the kinetic energy.

Image of FIG. 8.
FIG. 8.

Third-order structure function from low-wavenumber split linear forcing (thick lines, black on-line) and full spectrum split linear forcing (thin lines, red on-line) for series 1c. At a particular resolution, low- forcing produces a higher peak than full spectrum forcing.

Image of FIG. 9.
FIG. 9.

Structure function curves for the high resolution, nearly incompressible simulation 2a, with grid cells and , using low- split forcing. All three curves peak near the theoretically expected values (horizontal lines). Bottom: 4/15 law, ; middle: 4/5 law, ; and top: 4/3 law, .

Image of FIG. 10.
FIG. 10.

Isotropy relation at third order for simulation 2a, with grid cells and using low- split forcing. Solid line (red on-line): , dashed line (blue on-line): .

Image of FIG. 11.
FIG. 11.

Error in average statistics when using separately each of three realizations (thin lines) or the average of the three (thick line). The horizontal axis is the total model time required to compute the average. This comparison shows that computing time is better spent measuring statistics over one simulation, rather than taking the average of an ensemble of simulations.

Image of FIG. 12.
FIG. 12.

Error in the modified wavenumber when the sixth-order compact finite difference scheme is on the same grid as a spectral method (a) and on a finer grid (b) with grid spacing contracted by a factor of .


Generic image for table
Table I.

Parameter values of simulations, where series 1 was used for long-time simulations on a mesh (Figs. 1–6) and series 2 was used for high-resolution simulations on a mesh (Figs. 5–10). All simulations use and . For the split forcing described below, the initial values of and target are close to their values in the stationary state.

Generic image for table
Table II.

Low- forcing achieves nearly double the Taylor Reynolds number as that of full spectrum forcing at the same resolution with the same criterion for all runs.


Article metrics loading...


Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Forcing for statistically stationary compressible isotropic turbulence