Aspects of supercritical turbulence: direct numerical simulation of an Osub2/Hsub2 temporal mixing layer
Explore the source record for details and available documents.
supercritical turbulence↗
SEARCH · Search NASA
Results for “Direct numerical simulation”
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.
Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing
A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.
tubulent multicompenent-species↗
Investigation of Numerical Schemes for Direct Numerical Simulations of Supersonic Turbulent Boundary Layers
No abstract available
Morris, Christopher I.↗
Direct Numerical Simulation of Turbulent Couette-Poiseuille Flow With Zero Skin Friction
The near-wall scaling of mean velocity U(yw) is addressed for the case of zero skin friction on one wall of a fully turbulent channel flow. The present DNS results can be added to the evidence in support of the conjecture that U is proportional to the square root of yw in the region just above the wall at which the mean shear dU=dy = 0.
Coleman, Gary N.↗
Direct Numerical Simulation of Hypersonic Boundary Layer Transition in a Digital Wind Tunnel at Mach 8
Explore the source record for details and available documents.
boundary layer transition, acoustics↗
Multi-dimensional high order essentially non-oscillatory finite difference methods in generalized coordinates
The nonlinear stability of compact schemes for shock calculations is investigated. In recent years compact schemes were used in various numerical simulations including direct numerical simulation of turbulence. However to apply them to problems containing shocks, one has to resolve the problem of spurious numerical oscillation and nonlinear instability. A framework to apply nonlinear limiting to a local mean is introduced. The resulting scheme can be proven total variation (1D) or maximum norm (multi D) stable and produces nice numerical results in the test cases. The result is summarized in the preprint entitled 'Nonlinearly Stable Compact Schemes for Shock Calculations', which was submitted to SIAM Journal on Numerical Analysis. Research was continued on issues related to two and three dimensional essentially non-oscillatory (ENO) schemes. The main research topics include: parallel implementation of ENO schemes on Connection Machines; boundary conditions; shock interaction with hydrogen bubbles, a preparation for the full combustion simulation; and direct numerical simulation of compressible sheared turbulence.
Shu, Chi-Wang↗
Computation of the sound generated by isotropic turbulence
The acoustic radiation from isotropic turbulence is computed numerically. A hybrid direct numerical simulation approach which combines direct numerical simulation (DNS) of the turbulent flow with the Lighthill acoustic analogy is utilized. It is demonstrated that the hybrid DNS method is a feasible approach to the computation of sound generated by turbulent flows. The acoustic efficiency in the simulation of isotropic turbulence appears to be substantially less than that in subsonic jet experiments. The dominant frequency of the computed acoustic pressure is found to be somewhat larger than the dominant frequency of the energy-containing scales of motion. The acoustic power in the simulations is proportional to epsilon (M(sub t))(exp 5) where epsilon is the turbulent dissipation rate and M(sub t) is the turbulent Mach number. This is in agreement with the analytical result of Proudman (1952), but the constant of proportionality is smaller than the analytical result. Two different methods of computing the acoustic power from the DNS data bases yielded consistent results.
Sarkar, S.↗
A numerical study of the turbulent Ekman layer
The three-dimensional time-dependent turbulent flow in a neutrally stratified Ekman layer over a smooth flat surface was numerically simulated by directly solving the Navier-Stokes equations. Issues addressed using the direct numerical simulation (DNS) fields include the presence or absence of large-scale coherent structures ('longitudinal' or 'roll' vortices) in neutrally stratified Ekman-layer turbulence, the effects of the horizontal component of the angular velocity vector (i.e., latitude), and implications for models of the PBL. Experimental and DNS profiles are compared.
Coleman, G. N.↗
Numerical simulations of compressible mixing layers
Direct numerical simulations of two-dimensional temporally growing compressible mixing layers are presented. The Kelvin-Helmholtz instability is initially excited by a white-noise perturbation superimposed onto a hyperbolic tangent meanflow profile. The linear regime is studied at low resolution in the case of two flows of equal temperatures, for convective Mach numbers from 0.1 to 1 and for different values of the Reynolds number. At higher resolution, the complete evolution of a two-eddy mixing layer between two flows of different temperatures is simulated at moderate Reynolds number. Similarities and differences between flows of equal convective Mach numbers are discussed.
Normand, Xavier↗
Coherent structures in a simulated turbulent mixing layer
A direct numerical simulation of a plane turbulent mixing layer has been performed. The simulation was initialized using two turbulent velocity fields obtained from direct numerical simulation of a turbulent boundary layer at momentum thickness Reynolds number 300 (Spalart, 1988). The mixing layer is allowed to evolve long enough for self-similar linear growth to occur, with the visual thickness Reynolds number reaching 14,000. The simulated flow is examined for evidence of the coherent structures expected in a mixing layer (rollers and rib vortices). Before the onset of self-similar growth, such structures are present with properties similar to the corresponding laminar or transitional structures. In the self-similar growth regime, however, only the rollers are present with no indication of rib vortices and no indication of conventional pairing. This results in a reduction of mixing and layer growth.
Moser, Robert D.↗
The Effect of Two-Dimensional Forcing on Self-Similar Turbulent Plane Wakes
Three direct numerical simulations of time-evolving turbulent plane wakes with velocity deficit Reynolds numbers of about 2,000 have been simulated using a spectral numerical method with up to 600 x 260 x 160 modes. The initial conditions for the simulations are generated from direct numerical simulations of a turbulent boundary layer (momentum thickness Reynolds number of 670), and varying amounts of additional two- dimensional, forcing. In order to preserve the self-similar flow evolution, the forcing is implemented by multiplying all the two-dimensional modes in the initial condition by a constant factor. In the "natural" case no additional forcing is used; in the "forced" and "heavily forced" cases this factor is 5 and 20, respectively. The wake spreading rate Is increased by factors of 1.7 and 7.1 for the two forced cases. The Reynolds stresses are also increased by a similar or even larger factor. These results indicate that the plane wake is much more sensitive to initial forcing than the plane mixing layer. As in the plane mixing layer, two-dimensional forcing promotes more organized large-scale vortical flow structures and these structures axe sometimes separated by "braid regions" containing streamwise "rib" vortices, unlike in the unforced wake.
Rogers, M. M.↗
LES, DNS and RANS for the analysis of high-speed turbulent reacting flows
The objective of this research is to continue our efforts in advancing the state of knowledge in Large Eddy Simulation (LES), Direct Numerical Simulation (DNS), and Reynolds Averaged Navier Stokes (RANS) methods for the analysis of high-speed reacting turbulent flows. In the first phase of this research, conducted within the past six months, focus was in three directions: RANS of turbulent reacting flows by Probability Density Function (PDF) methods, RANS of non-reacting turbulent flows by advanced turbulence closures, and LES of mixing dominated reacting flows by a dynamics subgrid closure. A summary of our efforts within the past six months of this research is provided in this semi-annual progress report.
Givi, Peyman↗
Some Results Relevant to Statistical Closures for Compressible Turbulence
For weakly compressible turbulent fluctuations there exists a small parameter, the square of the fluctuating Mach number, that allows an investigation using a perturbative treatment. The consequences of such a perturbative analysis in three different subject areas are described: 1) initial conditions in direct numerical simulations, 2) an explanation for the oscillations seen in the compressible pressure in the direct numerical simulations of homogeneous shear, and 3) for turbulence closures accounting for the compressibility of velocity fluctuations. Initial conditions consistent with small turbulent Mach number asymptotics are constructed. The importance of consistent initial conditions in the direct numerical simulation of compressible turbulence is dramatically illustrated: spurious oscillations associated with inconsistent initial conditions are avoided, and the fluctuating dilatational field is some two orders of magnitude smaller for a compressible isotropic turbulence. For the isotropic decay it is shown that the choice of initial conditions can change the scaling law for the compressible dissipation. A two-time expansion of the Navier-Stokes equations is used to distinguish compressible acoustic and compressible advective modes. A simple conceptual model for weakly compressible turbulence - a forced linear oscillator is described. It is shown that the evolution equations for the compressible portions of turbulence can be understood as a forced wave equation with refraction. Acoustic modes of the flow can be amplified by refraction and are able to manifest themselves in large fluctuations of the compressible pressure.
Ristorcelli, J. R.↗
Dynamics of Aerosol Particles in Stationary, Isotropic Turbulence
A detailed study of the dynamics of sub-Kolmogorov-size aerosol particles in stationary isotropic turbulence has been performed. The study combined direct numerical simulations (DNS; directed by Prof. Collins) and high-resolution experimental measurements (directed by Prof. Meng) under conditions of nearly perfect geometric and parametric overlap. The goal was to measure the accumulation of particles in low-vorticity regions of the flow that arises from the effect commonly referred to as preferential concentration. The grant technically was initiated on June 13, 2000; however, funding was not available until July 11, 2000. The grant was originally awarded to Penn State University (numerical simulations) and SUNY-Buffalo (experiments); however, Prof. Collins effort was moved to Cornell University on January 2002 when he joined that university. He completed the study there. A list of the specific tasks that were completed under this study is presented.
Collins, Lance R.↗
An experimental study of the effects of rapid rotation on turbulence
Experiments, large eddy simulations and direct numerical simulations all show that rapid rotation has a striking effect on homogeneous nearly isotropic turbulence. The cascade process is effectively inhibited by rotation, and thus dissipation is greatly reduced. Some attempts have been made to incorporate this effect in turbulence models. Numerical simulation further showed the somewhat surprising result that anisotropic homogeneous turbulence subjected to rotation tended towards an isotropic state; however, the residual anisotropy was not zero. Reynolds (1989) performed a detailed analysis using Rapid Distortion Theory (RDT) and showed that a reduction in the anisotropy is indeed expected and if the anisotropy is produced by irrotational strain then the anisotropy tensor b(sub ij) is asymptotically driven to half its initial value. Our objective is to extend the work of the experiments mentioned above to lower turbulent Rossby numbers and to confirm some of the results obtained by Reynolds (1989) for anisotropic turbulence.
Veeravalli, S. V.↗
LES, DNS and RANS for the analysis of high-speed turbulent reacting flows
The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES), direct numerical simulation (DNS), and Reynolds averaged Navier Stokes (RANS) methods for the computational analysis of high-speed reacting turbulent flows. In the second phase of this work, covering the period 1 Sep. 1993 - 1 Sep. 1994, we have focused our efforts on two research problems: (1) developments of 'algebraic' moment closures for statistical descriptions of nonpremixed reacting systems, and (2) assessments of the Dirichlet frequency in presumed scalar probability density function (PDF) methods in stochastic description of turbulent reacting flows. This report provides a complete description of our efforts during this past year as supported by the NASA Langley Research Center under Grant NAG1-1122.
Givi, Peyman↗
LES, DNS and RANS for the analysis of high-speed turbulent reacting flows
The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES), direct numerical simulation (DNS), and Reynolds averaged Navier Stokes (RANS) methods for the computational analysis of high-speed reacting turbulent flows. In the second phase of this work, covering the period 1 Aug. 1994 - 31 Jul. 1995, we have focused our efforts on two programs: (1) developments of explicit algebraic moment closures for statistical descriptions of compressible reacting flows and (2) development of Monte Carlo numerical methods for LES of chemically reacting flows.
Adumitroaie, V.↗