Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homogeneous Isotropic Turbulence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Evaluation of the Use of Second Generation Wavelets in the Coherent Vortex Simulation Approach

The objective of this study is to investigate the use of the second generation bi-orthogonal wavelet transform for the field decomposition in the Coherent Vortex Simulation of turbulent flows. The performances of the bi-orthogonal second generation wavelet transform and the orthogonal wavelet transform using Daubechies wavelets with the same number of vanishing moments are compared in a priori tests using a spectral direct numerical simulation (DNS) database of isotropic turbulence fields: 256(exp 3) and 512(exp 3) DNS of forced homogeneous turbulence (Re(sub lambda) = 168) and 256(exp 3) and 512(exp 3) DNS of decaying homogeneous turbulence (Re(sub lambda) = 55). It is found that bi-orthogonal second generation wavelets can be used for coherent vortex extraction. The results of a priori tests indicate that second generation wavelets have better compression and the residual field is closer to Gaussian. However, it was found that the use of second generation wavelets results in an integral length scale for the incoherent part that is larger than that derived from orthogonal wavelets. A way of dealing with this difficulty is suggested.

Goldstein, D. E.↗

On the Lighthill relationship and sound generation from isotropic turbulence

In 1952, Lighthill developed a theory for determining the sound generated by a turbulent motion of a fluid. With some statistical assumptions, Proudman applied this theory to estimate the acoustic power of isotropic turbulence. Recently, Lighthill established a simple relationship that relates the fourth-order retarded time and space covariance of his stress tensor to the corresponding second-order covariance and the turbulent flatness factor, without making statistical assumptions for a homogeneous turbulence. Lilley revisited Proudman's work and applied the Lighthill relationship to evaluate directly the radiated acoustic power from isotropic turbulence. After choosing the time separation dependence in the two-point velocity time and space covariance based on the insights gained from direct numerical simulations, Lilley concluded that the Proudman constant is determined by the turbulent flatness factor and the second-order spatial velocity covariance. In order to estimate the Proudman constant at high Reynolds numbers, we analyzed a unique data set of measurements in a large wind tunnel and atmospheric surface layer that covers a range of the Taylor microscale based on Reynolds numbers 2.0 x 10(exp 3) less than or equal to R(sub lambda) less than or equal to 12.7 x 10(exp 3). Our measurements demonstrate that the Lighthill relationship is a good approximation, providing additional support to Lilley's approach. The flatness factor is found between 2.7 - 3.3 and the second order spatial velocity covariance is obtained. Based on these experimental data, the Proudman constant is estimated to be 0.68 - 3.68.

Zhou, YE↗

The nonlinear Galerkin method: A multi-scale method applied to the simulation of homogeneous turbulent flows

Using results of Direct Numerical Simulation (DNS) in the case of two-dimensional homogeneous isotropic flows, the behavior of the small and large scales of Kolmogorov like flows at moderate Reynolds numbers are first analyzed in detail. Several estimates on the time variations of the small eddies and the nonlinear interaction terms were derived; those terms play the role of the Reynolds stress tensor in the case of LES. Since the time step of a numerical scheme is determined as a function of the energy-containing eddies of the flow, the variations of the small scales and of the nonlinear interaction terms over one iteration can become negligible by comparison with the accuracy of the computation. Based on this remark, a multilevel scheme which treats differently the small and the large eddies was proposed. Using mathematical developments, estimates of all the parameters involved in the algorithm, which then becomes a completely self-adaptive procedure were derived. Finally, realistic simulations of (Kolmorov like) flows over several eddy-turnover times were performed. The results are analyzed in detail and a parametric study of the nonlinear Galerkin method is performed.

Debussche, A.↗

Determining the alpha dynamo parameter in incompressible homogeneous magnetohydrodynamic turbulence

Alpha, an important parameter in dynamo theory, is proportional to either the kinetic, current, magnetic, or velocity helicity of the fluctuating magnetic field and fluctuating velocity field. The particular helicity to which alpha is proportional depends on the assumptions used in deriving the first order smoothed equations that describe the alpha effect. In two cases, when alpha is proportional to either the magnetic helicity or velocity helicity, alpha is determined experimentally from two point measurements of the fluctuating fields in incompressible, homogeneous turbulence having arbitrary symmetry. For the other two possibilities, alpha is determined if the turbulence is isotropic.

Matthaeus, W. H.↗

Validation of a Pseudo-Sound Theory for the Pressure-Dilatation in DNS of Compressible Turbulence

The results of an asymptotic theory for statistical closures for compressible turbulence are explored and validated with the direct numerical simulation of the isotropic decay and the homogeneous shear. An excellent collapse of the data is seen. The slow portion is found to scale, as predicted by the theory, with the quantity M(sub t)(sup 2) and epsilon(sub s). The rapid portion has an unambiguous scaling with alpha(sup 2)M(sub t)(sup s)epsilon(sub s)[P(sub k)/epsilon - l](Sk/epsilon)(sup 2). Implicit in the scaling is a dependence, as has been noted by others, on the gradient Mach number. A new feature of the effects of compressibility, that of the Kolmogorov scaling coefficient, alpha, is discussed. It is suggested that alpha may contain flow specific physics associated with large scales that might provide further insight into the structural effects of compressibility.

Ristorcelli, J. R.↗

On the Two-Dimensionalization of Homogeneous Rotating Turbulence

Large-eddy simulation of the incompressible Navier-Stokes equations has been used to examine the long-time development of initially isotropic turbulence subjected to solid-body rotation. The simulations were carried out using a pseudo-spectral method with 128 x 128 x 512 collocation points in a computational domain that is four times larger along the rotation axis than in the other directions; subgrid-scale motions were parameterized using a spectral eddy viscosity model modified for system rotation. Simulation results show that the correlation length along the rotation am's of velocities orthogonal to the rotation vector exhibits rapid growth while the integral length-scale of velocities aligned with the rotation axis is relatively unaffected by rotation. Examination of the energy spectrum of two-dimensional, two-component motions indicates the presence of an inverse cascade of energy. System rotation also causes an alignment of vorticity along the rotation axis with relatively stronger cyclonic vorticity than anticyclonic. The onset of anisotropic effects are well characterized by Rossby numbers defined in terms of both macroscopic and microscopic quantities.

Squires, K. D.↗

The alpha dynamo parameter and measurability of helicities in magnetohydrodynamic turbulence

Alpha, an important parameter in dynamo theory, is shown to be proportional to either the kinetic, current, magnetic, or velocity helicities of the fluctuating magnetic field and fluctuating velocity field. The particular helicity to which alpha is proportional depends on the assumptions used in deriving the first-order smoothed equations that describe the alpha effect. In two cases, viz., when alpha is proportional to either the magnetic helicity or velocity helicity, alpha can be determined experimentally from two-point measurements of the fluctuating fields in incompressible, homogeneous turbulence with arbitrary rotational symmetry. For the other two possibilities, alpha can be determined if the turbulence is isotropic.

Matthaeus, W. H.↗

Infrared properties of an anisotropically stirred fluid

A renormalization group is developed for the Navier-Stokes equations driven by an anisotropically correlated random stirring force. The stirring force generates homogeneous turbulence with a preferred direction. The force correlation is the sum of a small anisotropic perturbation and an isotropic correlation chosen, so that the fixed point of renormalization group has a k exp -5/3 energy spectrum. Fixed points for the anisotropic correlation are found near this isotropic fixed point. Two types of anisotropy are analyzed. when the additional stirring is in the plane perpendicular to the preferred direction, the renormalized viscosity is increased. When it is aligned with the preferred direction, the viscosity is decreased. A possible connection with the inverse energy cascade of two-dimensional turbulence is discussed.

Rubinstein, Robert↗

The structure of correlation tensors in homogeneous anisotropic turbulence

The study of turbulence with spatially homogeneous but anisotropic statistical properties has applications in space physics and laboratory plasma physics. The first step in the systematic study of such fluctuations is the elucidation of the kinematic properties of the relevant statistical objects, which are the correlation tensors. The theory of isotropic tensors, developed by Robertson, Chandrasekhar and others, is reviewed and extended to cover the general case of turbulence with a pseudo-vector preferred direction, without assuming mirror reflection invariance. Attention is focused on two point correlation functions and it is shown that the form of the decomposition into proper and pseudo-tensor contributions is restricted by the homogeneity requirement. It is also shown that the vector and pseudo-vector preferred direction cases yield different results. An explicit form of the two point correlation tensor is presented which is appropriate for analyzing interplanetary magnetic fluctuations. A procedure for determining the magnetic helicity from experimental data is presented.

Matthaeus, W. H.↗

Statistical analysis of turbulence data from the NASA Marshall Space Flight Center atmospheric boundary layer tower array facility

Statistical properties of atmospheric turbulence near the Earth's surface are presented. Emphasis is placed on the probability density distribution two point spatial correlation, length scale, and two point and single point spectrum. Comparison of the data with isotropic homogeneous models is made. In general, agreement with the models is poor. For the design of aircraft during operations in the lower levels associated with approach, takeoff, and terrain following, it appears necessary to improve existing models or develop new nonisotropic turbulence models.

Frost, W.↗

Study of the velocity gradient tensor in turbulent flow

The behavior of the velocity gradient tensor, A(ij)=delta u(i)/delta x(j), was studied using three turbulent flows obtained from direct numerical simulation The flows studies were: an inviscid calculation of the interaction between two vortex tubes, a homogeneous isotropic flow, and a temporally evolving planar wake. Self-similar behavior for each flow was obtained when A(ij) was normalized with the mean strain rate. The case of the interaction between two vortex tubes revealed a finite sized coherent structure with topological characteristics predictable by a restricted Euler model. This structure was found to evolve with the peak vorticity as the flow approached singularity. Invariants of A(ij) within this structure followed a straight line relationship of the form: gamma(sup 3)+gammaQ+R=0, where Q and R are the second and third invariants of A(ij), and the eigenvalue gamma is nearly constant over the volume of this structure. Data within this structure have local strain topology of unstable-node/saddle/saddle. The characteristics of the velocity gradient tensor and the anisotropic part of a related acceleration gradient tensor H(ij) were also studied for a homogeneous isotropic flow and a temporally evolving planar wake. It was found that the intermediate principal eigenvalue of the rate-of-strain tensor of H(ij) tended to be negative, with local strain topology of the type stable-node/saddle/saddle. There was also a preferential eigenvalue direction. The magnitude of H(ij) in the wake flow was found to be very small when data were conditioned at high local dissipation regions. This result was not observed in the relatively low Reynolds number simulation of homogeneous isotropic flow. A restricted Euler model of the evolution of A(ij) was found to reproduce many of the topological features identified in the simulations.

Cheng, Wei-Ping↗

Modeling of pressure spectra in a turbulent shear flow

The development of the pressure spectrum has been followed from its separation into turbulent-turbulent and turbulent-mean square contributions to a unified spectral model. The model is obtained by Fourier transforming the integral solution to Poisson's equation for an isotropic, homogeneous, constant mean shear flow. The variations in the asymptotic form for the mean-shear spectrum is discussed. It is noted that the spectra are associated with the particular velocity field selected, and that by an appropriate choice of spectral forms, the mean-square pressure of the flow may be calculated.

Beuther, P.↗

Simulation of rotation effect in turbulence dissipation rate equation

It was observed that the Reynolds stress transport equations using the usual isotropic dissipation rate transport equations do not accurately predict the correct evolution of the kinetic energy in both homogeneous shear and isotropic decay. This same effect carries over to more practical flows such as those with strong irrotational strains and/or adverse pressure gradients. The deficiencies in these predictions were attributed to the need for modifications in the isotropic dissipation rate equation. Pope and Hanjalic and Launder have proposed modifications to this equation. These modifications were tested in the cases of homogeneous shear and iso-decay and their performance was analyzed. The results are compared to the rotating flow models of Bardina and Raj.

Raj, Rishi S.↗

Modeling the turbulent kinetic energy equation for compressible, homogeneous turbulence

The turbulent kinetic energy transport equation, which is the basis of turbulence models, is investigated for homogeneous, compressible turbulence using direct numerical simulations performed at CTR. It is shown that the partition between dilatational and solenoidal modes is very sensitive to initial conditions for isotropic decaying turbulence but not for sheared flows. The importance of the dilatational dissipation and of the pressure-dilatation term is evidenced from simulations and a transport equation is proposed to evaluate the pressure-dilatation term evolution. This transport equation seems to work well for sheared flows but does not account for initial condition sensitivity in isotropic decay. An improved model is proposed.

Aupoix, B.↗

Particle transport in Markov stochastic mixtures with spatial gradients

Markov geometries are a prototype class of stochastic media that are widely used to model complex disordered systems. In a series of recent works, we have shown that spatially homogeneous (but possibly non-isotropic) Poisson tessellations can be conveniently used in order to generate three-dimensional realizations of Markov media and thus characterize the features of particle transport e.g. in fragmented fuel elements following severe nuclear accidents or Rayleigh-Taylor turbulent mixing for inertial confinement fusion. In this paper we show that Poisson tessellations can be extended to take into account the presence of spatial gradients, which commonly occur in real-world applications due to material stratification. For this purpose, we will provide an explicit construction of spatially-non-homogeneous Poisson tessellations, fully characterize their statistical properties, and finally illustrate the behaviour of particle transport in such random media. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of turbulence in the orthonormal wavelet representation

The usefulness of the wavelet transform for the analysis of turbulent flow fields is explored by examining the wavelet transform properties of a decomposition of turbulent velocity fields into modes that exhibit the localization in a wavenumber and physical space. The calculations are performed on 3D fields from direct numerical simulations of isotropic flow and homogeneous shear flow, and from measurements in two laboratory wind-tunnel experimental velocity signals (boundary layer and wake behind a circular cylinder). The analysis confirmed that there is strong spatial intermittency in nonlinear quantities; their mean spectral behavior results from a delicate balance between large positive and negative excursions. The wavelet analysis is a way to quantify these observations in a standardized fashion by using 'flow-independent eddies' to decompose the velocity field.

Meneveau, Charles↗

Linearly Forced Isotropic Turbulence

Stationary isotropic turbulence is often studied numerically by adding a forcing term to the Navier-Stokes equation. This is usually done for the purpose of achieving higher Reynolds number and longer statistics than is possible for isotropic decaying turbulence. It is generally accepted that forcing the Navier-Stokes equation at low wave number does not influence the small scale statistics of the flow provided that there is wide separation between the largest and smallest scales. It will be shown, however, that the spectral width of the forcing has a noticeable effect on inertial range statistics. A case will be made here for using a broader form of forcing in order to compare computed isotropic stationary turbulence with (decaying) grid turbulence. It is shown that using a forcing function which is directly proportional to the velocity has physical meaning and gives results which are closer to both homogeneous and non-homogeneous turbulence. Section 1 presents a four part series of motivations for linear forcing. Section 2 puts linear forcing to a numerical test with a pseudospectral computation.

Lundgren, T. S.↗

Eulerian and Lagrangian second-order statistics of superfluid 4 He grid turbulence

Here, we use particle-tracking velocimetry to study Eulerian and Lagrangian second-order statistics of superfluid 4 He grid turbulence. The Eulerian energy spectra at scales larger than the mean distance between quantum vortex lines behave classically with close to Kolmogorov-1941 scaling and are almost isotropic. The Lagrangian second-order structure functions and frequency power spectra, measured at scales comparable with the intervortex distance, demonstrate a sharp transition from nearly classical behavior to a regime dominated by the motion of quantum vortex lines. Employing the homogeneity of the flow, we verify a set of relations that connect various second-order statistical objects that stress different aspects of turbulent behavior, allowing a multifaceted analysis. We use the two-way bridge relations between Eulerian energy spectra and second-order structure functions to reconstruct the energy spectrum from the known second-order velocity structure function and vice versa. The Lagrangian frequency spectrum reconstructed from the measured Eulerian spectrum using the Eulerian-Lagrangian bridge differs from the measured Lagrangian spectrum in the quasiclassical range, which calls for further investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗