Numerical simulations of fully developed turbulent channel flow at three Reynolds numbers up to Reτ=590 are reported. It is noted that the higher Reynolds number simulations exhibit fewer low Reynolds number effects than previous simulations at Reτ=180. A comprehensive set of statistics gathered from the simulations is available on the web at http://www.tam.uiuc.edu/Faculty/Moser/channel.

Direct numerical simulation of turbulent channel flow up to Reτ=590

The structure of energy-containing turbulence in the outer region of a zero-pressure-
 gradient boundary layer has been studied using particle image velocimetry (PIV) to measure the instantaneous velocity fields in a streamwise-wall-normal plane. Experiments performed at three Reynolds numbers in the range 930 0) that occur on a locus inclined at 30–60° to the wall.In the outer layer, hairpin vortices occur in streamwise-aligned packets that propagate with small velocity dispersion. Packets that begin in or slightly above the buffer layer are very similar to the packets created by the autogeneration mechanism (Zhou, Adrian & Balachandar 1996). Individual packets grow upwards in the streamwise direction at a mean angle of approximately 12°, and the hairpins in packets are typically spaced several hundred viscous lengthscales apart in the streamwise direction. Within the interior of the envelope the spatial coherence between the velocity fields induced by the individual vortices leads to strongly retarded streamwise momentum, explaining the zones of uniform momentum observed by Meinhart & Adrian (1995). The packets are an important type of organized structure in the wall layer in which relatively small structural units in the form of three-dimensional vortical structures are arranged coherently, i.e. with correlated spatial relationships, to form much longer structures. The formation of packets explains the occurrence of multiple VITA events in turbulent ‘bursts’, and the creation of Townsend's (1958) large-scale inactive motions. These packets share many features of the hairpin models proposed by Smith (1984) and co-workers for the near-wall layer, and by Bandyopadhyay (1980), but they are shown to occur in a hierarchy of scales across most of the boundary layer.In the logarithmic layer, the coherent vortex packets that originate close to the wall frequently occur within larger, faster moving zones of uniform momentum, which may extend up to the middle of the boundary layer. These larger zones are the induced interior flow of older packets of coherent hairpin vortices that originate upstream and over-run the younger, more recently generated packets. The occurence of small hairpin packets in the environment of larger hairpin packets is a prominent feature of the logarithmic layer. With increasing Reynolds number, the number of hairpins in a packet increases.

Vortex organization in the outer region of the turbulent boundary layer

▪ Abstract We review the direct numerical simulation (DNS) of turbulent flows. We stress that DNS is a research tool, and not a brute-force solution to the Navier-Stokes equations for engineering problems. The wide range of scales in turbulent flows requires that care be taken in their numerical solution. We discuss related numerical issues such as boundary conditions and spatial and temporal discretization. Significant insight into turbulence physics has been gained from DNS of certain idealized flows that cannot be easily attained in the laboratory. We discuss some examples. Further, we illustrate the complementary nature of experiments and computations in turbulence research. Examples are provided where DNS data has been used to evaluate measurement accuracy. Finally, we consider how DNS has impacted turbulence modeling and provided further insight into the structure of turbulent boundary layers.

/pdf/direct-numerical-simulation-a-tool-in-turbulence-research-2a650vh1lf.pdf

DIRECT NUMERICAL SIMULATION: A Tool in Turbulence Research

▪ Abstract The single-point statistics of turbulence in the ‘roughness sub-layer’ occupied by the plant canopy and the air layer just above it differ significantly from those in the surface layer. The mean velocity profile is inflected, second moments are strongly inhomogeneous with height, skewnesses are large, and second-moment budgets are far from local equilibrium. Velocity moments scale with single length and time scales throughout the layer rather than depending on height. Large coherent structures control turbulence dynamics. Sweeps rather than ejections dominate eddy fluxes and a typical large eddy consists of a pair of counter-rotating streamwise vortices, the downdraft between the vortex pair generating the sweep. Comparison with the statistics and instability modes of the plane mixing layer shows that the latter rather than the boundary layer is the appropriate model for canopy flow and that the dominant large eddies are the result of an inviscid instability of the inflected mean velocity profi...

Turbulence in plant canopies

The accuracy of numerical simulation algorithms is one of main concerns in modern Computational Fluid Dynamics. Development of new and more accurate mathematical models requires an insight into the problem of numerical errors. In order to construct an estimate of the solution error in Finite Volume calculations, it is first necessary to examine its sources. Discretisation errors can be divided into two groups: errors caused by the discretisation of the solution domain and equation discretisation errors. The first group includes insufficient mesh resolution, mesh skewness and non-orthogonality. In the case of the second order Finite Volume method, equation discretisation errors are represented through numerical diffusion. Numerical diffusion coefficients from the discretisation of the convection term and the temporal derivative are derived. In an attempt to reduce numerical diffusion from the convection term, a new stabilised and bounded second-order differencing scheme is proposed. Three new methods of error estimation are presented. The Direct Taylor Series Error estimate is based on the Taylor series truncation error analysis. It is set up to enable single-mesh single-run error estimation. The Moment Error estimate derives the solution error from the cell imbalance in higher moments of the solution. A suitable normalisation is used to estimate the error magnitude. The Residual Error estimate is based on the local inconsistency between face interpolation and volume integration. Extensions of the method to transient flows and the Local Residual Problem error estimate are also given. Finally, an automatic error-controlled adaptive mesh refinement algorithm is set up in order to automatically produce a solution of pre-determined accuracy. It uses mesh refinement and unrefinement to control the local error magnitude. The method is tested on several characteristic flow situations, ranging from incompressible to supersonic flows, for both steady-state and transient problems.

Error analysis and estimation for the finite volume method with applications to fluid flows

Spectral methods for the Navier-Stokes equations with one infinite and two periodic directions

Three direct numerical simulations of incompressible turbulent plane mixing layers have been performed. All the simulations were initialized with the same two velocity fields obtained from a direct numerical simulation of a turbulent boundary layer with a momentum thickness Reynolds number of 300 computed by Spalart (J. Fluid Mech. 187, 61, 1988). In addition to a baseline case with no additional disturbances, two simulations were begun with two-dimensional disturbances of varying strength in addition to the boundary layer turbulence. After a development stage, the baseline case and the case with weaker additional two-dimensional disturbances evolve self-similarly, reaching visual thickness Reynolds numbers of up to 20 000. This self-similar period is characterized by a lack of large-scale organized pairings, a lack of streamwise vortices in the 'braid' regions, and scalar mixing that is characterized by 'marching' Probability Density Functions (PDFs). The case begun with strong additional two-dimensional disturbances only becomes approximately self-similar, but exhibits sustained organized large-scale pairings, clearly defined braid regions with streamwise vortices that span them, and scalar PDFs that are 'nonmarching.' It is also characterized by much more intense vertical velocity fluctuations than the other two cases. The statistics and structures in several experiments involving turbulent mixing layers are in better agreement with those of the simulations that do not exhibit organized pairings.

Direct Simulation of a Self-Similar Turbulent Mixing Layer

The structures of the vorticity fields in several homogeneous irrotational straining flows and a homogeneous turbulent shear flow were examined using a database generated by direct numerical simulation of the unsteady Navier-Stokes equations. In all cases, strong evidence was found for the presence of coherent vortical structures. The initially isotropic vorticity fields were rapidly affected by imposed mean strain and the rotational component of mean shear and developed accordingly. In the homogeneous turbulent shear-flow cases, the roll-up of mean vorticity into characteristic hairpin vortices was clearly observed, supporting the view that hairpin vortices are an important vortical structure in all turbulent shear flows; the absence of mean shear in the homogeneous irrotational straining flows precludes the presence of hairpin-like vortices.

The structure of the vorticity field in homogeneous turbulent flows

The Kelvin Helmholtz roll up of three dimensional, temporally evolving, plane mixing layers were simulated numerically. All simulations were begun from a few low wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity profile. The spanwise disturbance wavelength was taken to be less than or equal to the streamwise wavelength associated with the Kelvin Helmholtz roll up. A standard set of clean structures develop in most of the simulations. The spanwise vorticity rolls up into a corrugated spanwise roller, with vortex stretching creating strong spanwise vorticity in a cup shaped region at the vends of the roller. Predominantly streamwise rib vortices develop in the braid region between the rollers. For sufficiently strong initial three dimensional disturbances, these ribs collapse into compact axisymmetric vortices. The rib vortex lines connect to neighboring ribs and are kinked in the opposite direction of the roller vortex lines. Because of this, these two sets of vortex lines remain distinct. For certain initial conditions, persistent ribs do not develop. In such cases the development of significant three dimensionality is delayed. When the initial three dimensional disturbance energy is about equal to, or less than, the two dimensional fundamental disturbance energy, the evolution of the three dimensional disturbance is nearly linear (with respect to the mean and the two dimensional disturbances), at least until the first Kelvin Helmholtz roll up is completed.

The three-dimensional evolution of a plane mixing layer - The Kelvin-Helmholtz rollup

The evolution of three-dimensional temporally evolving plane mixing layers through as many as three pairings has been simulated numerically. All simulations were begun from a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. Three-dimensional perturbations were used with amplitudes ranging from infinitesimal to large enough to trigger a rapid transition to turbulence. Pairing is found to inhibit the growth of infinitesimal three-dimensional disturbances, and to trigger the transition to turbulence in highly three-dimensional flows. The mechanisms responsible for the growth of three-dimensionality and onset of transition to turbulence are described. The transition to turbulence is accompanied by the formation of thin sheets of spanwise vorticity, which undergo secondary rollups. The post-transitional simulated flow fields exhibit many properties characteristic of turbulent flows.

https://ntrs.nasa.gov/api/citations/19930014228/downloads/19930014228.pdf

Michael M. Rogers

Papers

Spectral methods for the Navier-Stokes equations with one infinite and two periodic directions

Direct Simulation of a Self-Similar Turbulent Mixing Layer

The structure of the vorticity field in homogeneous turbulent flows

The three-dimensional evolution of a plane mixing layer - The Kelvin-Helmholtz rollup

The three-dimensional evolution of a plane mixing layer: pairing and transition to turbulence