An ILLIAC program for the numerical simulation of homogeneous incompressible turbulence

An algorithm and ILLIAC computer program, developed for the simulation of homogeneous incompressible turbulence in the presence of an applied mean strain, are described. The turbulence field is represented spatially by a truncated triple Fourier series (spectral method) and followed in time using a fourth-order Runge-Kutta algorithm. These include: (1) transformation of variables suggested by Taylor's sudden-distortion theory; (2) implicit viscous diffusion by use of an integrating factor; (3) implicit pressure calculation suggested by Taylor's sudden-distortion theory, and (4) inexpensive control of aliasing by random and phased coordinate shifts

Nonlinear Evolution of Hydrodynamical Shear Flows in Two Dimensions

We examine how perturbed shear flows evolve in two-dimensional, incompressible, inviscid hydrodynamical fluids, with the ultimate goal of understanding the dynamics of accretion disks. To linear order, vorticity waves are swung around by the background shear, and their velocities are amplified transiently before decaying. It has been speculated that sufficiently amplified modes might couple nonlinearly, leading to turbulence. Here we show how nonlinear coupling occurs in two dimensions. This coupling is remarkably simple because it only lasts for a short time interval, when one of the coupled modes is in mid-swing. We focus on the interaction between a swinging and an axisymmetric mode. There is instability provided that k_{y,swing}/k_{x,axi} < omega/q, i.e., that the ratio of wavenumbers is less than the ratio of the axisymmetric mode's vorticity to the background vorticity. If this is the case, then when the swinging mode is in mid-swing it couples with the axisymmetric mode to produce a new leading swinging mode that has larger vorticity than itself; this new mode in turn produces an even larger leading mode, etc. Therefore all axisymmetric modes, regardless of how small in amplitude, are unstable to perturbations with sufficiently large azimuthal wavelength. We show that this shear instability occurs whenever the momentum transported by a perturbation has the sign required for it to diminish the background shear; only when this occurs can energy be extracted from the mean flow and hence added to the perturbation. For an accretion disk, this means that the instability transports angular momentum outwards while it operates.Comment: published versio

Energy transfer in isotropic turbulence at low Reynolds numbers

Detailed measurements were made of energy transfer among the scales of motion in incompressible turbulent fields at low Reynolds numbers generated by direct numerical simulation. It was observed that although the transfer resulted from triad interactions that were non-local in k space, the energy always transferred locally. The results are consistent with the notion of non-uniform advection of small weak eddies by larger and stronger ones, similar to transfer processes in the far dissipation range at high Reynolds numbers

Implementation of Lees-Edwards periodic boundary conditions for direct numerical simulations of particle dispersions under shear flow

A general methodology is presented to perform direct numerical simulations of particle dispersions in a shear flow with Lees-Edwards periodic boundary conditions. The Navier-Stokes equation is solved in oblique coordinates to resolve the incompatibility of the fluid motions with the sheared geometry, and the force coupling between colloidal particles and the host fluid is imposed by using a smoothed profile method. The validity of the method is carefully examined by comparing the present numerical results with experimental viscosity data for particle dispersions in a wide range of volume fractions and shear rates including nonlinear shear-thinning regimes

Numerical experiments in homogeneous turbulence

The direct simulation methods developed by Orszag and Patternson (1972) for isotropic turbulence were extended to homogeneous turbulence in an incompressible fluid subjected to uniform deformation or rotation. The results of simulations for irrotational strain (plane and axisymmetric), shear, rotation, and relaxation toward isotropy following axisymmetric strain are compared with linear theory and experimental data. Emphasis is placed on the shear flow because of its importance and because of the availability of accurate and detailed experimental data. The computed results are used to assess the accuracy of two popular models used in the closure of the Reynolds-stress equations. Data from a variety of the computed fields and the details of the numerical methods used in the simulation are also presented

Navier-Stokes Simulation of Homogeneous Turbulence on the CYBER 205

A computer code which solves the Navier-Stokes equations for three dimensional, time-dependent, homogenous turbulence has been written for the CYBER 205. The code has options for both 64-bit and 32-bit arithmetic. With 32-bit computation, mesh sizes up to 64 (3) are contained within core of a 2 million 64-bit word memory. Computer speed timing runs were made for various vector lengths up to 6144. With this code, speeds a little over 100 Mflops have been achieved on a 2-pipe CYBER 205. Several problems encountered in the coding are discussed

Pressure and higher-order spectra for homogeneous isotropic turbulence

The spectra of the pressure, and other higher-order quantities including the dissipation, the enstrophy, and the square of the longitudinal velocity derivative are computed using data obtained from direct numerical simulation of homogeneous isotropic turbulence at Taylor-Reynolds numbers R(sub lambda) in the range 38 - 170. For the pressure spectra we find reasonable collapse in the dissipation range (of the velocity spectrum) when scaled in Kolmogorov variables and some evidence, which is not conclusive, for the existence of a k(exp -7/3) inertial range where k = absolute value of K, is the modulus of the wavenumber. The power spectra of the dissipation, the enstrophy, and the square of the longitudinal velocity derivative separate in the dissipation range, but appear to converge together in the short inertial range of the simulations. A least-squares curve-fit in the dissipation range for one value of R(sub lambda) = 96 gives a form for the spectrum of the dissipation as k(exp 0)exp(-Ck eta), for k(eta) greater than 0.2, where eta is the Kolmogorov length and C is approximately equal to 2.5

Method permits mechanical and electrical checkout of piezoelectric transducers while installed in a system

Known dc voltage is applied and then removed suddenly in a method to permit checkout of the mechanical and electrical condition of piezoelectric transducers of the cantilever beam type, while installed in a system

On statistically stationary homogeneous shear turbulence

A statistically stationary turbulence with a mean shear gradient is realized in a flow driven by suitable body forces. The flow domain is periodic in downstream and spanwise directions and bounded by stress free surfaces in the normal direction. Except for small layers near the surfaces the flow is homogeneous. The fluctuations in turbulent energy are less violent than in the simulations using remeshing, but the anisotropy on small scales as measured by the skewness of derivatives is similar and decays weakly with increasing Reynolds number.Comment: 4 pages, 5 figures (Figs. 3 and 4 as external JPG-Files

The structure of intense vorticity in homogeneous isotropic turbulence

The structure of the intense vorticity regions is studied in numerically simulated homogeneous, isotropic, equilibrium turbulent flow fields at four different Reynolds numbers in the range Re(sub lambda) = 36-171. In accordance with previous investigators, this vorticity is found to be organized in coherent, cylindrical or ribbon-like, vortices ('worms'). A statistical study suggests that they are just especially intense features of the background, O(omega'), vorticity. Their radii scale with the Kolmogorov microscale and their lengths with the integral scale of the flow. An interesting observation is that the Reynolds number based on the circulation of the intense vortices, gamma/nu, increases monotonically with Re(sub lambda), raising the question of the stability of the structures in the limit of Re(sub lambda) approaching infinity. One and two-dimensional statistics of vorticity and strain are presented; they are non-gaussian, and the behavior of their tails depends strongly on the Reynolds number. There is no evidence of convergence to a limiting distribution in our range of Re(sub lambda), even though the energy spectra and the energy dissipation rate show good asymptotic properties in the higher Reynolds number cases. Evidence is presented to show that worms are natural features of the flow and that they do not depend on the particular forcing scheme