Search NASA⌕ Search

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.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Asynchronicity in opposed-piston RCMs: Does it matter?

Rapid Compression Machines (RCMs) are widely utilized to study combustion phenomena at engine-relevant conditions, and significant efforts are typically made to create a quiescent environment, particularly for investigations of autoignition chemistry. Opposed-piston configurations can be advantageous due to shorter compression times and reduced surface area to volume ratios. Each side must be actuated simultaneously, but this can be challenging in practice. These devices, like most RCMs, utilize hydraulics for actuation, speed control and arrestation of the piston at the end of the stroke; there is no mechanical control or linkage of the two piston trajectories. To quantify the magnitudes and effects of piston asynchronous behavior, this work employs both detailed experimental measurements and, for the first time, high-fidelity, Direct Numerical Simulation (DNS). The boundary conditions are carefully considered applying insight from high-resolution linear variable differential transformer (LVDT) measurements of the piston trajectory and a zero-dimensional kinematics model of the piston-shaft assembly. Sufficient resolution in the piston crevice region is used. The complicated fluid dynamical behavior that can evolve during piston compression and the ensuing delay processes due to offset timings from t offset = 0-10 ms is elucidated. It is found that near t offset = 6 ms and beyond, the boundary layer on the face of the first-seating piston can be sufficiently perturbed, due initially to reemergence of gas from the crevice of the firstseating piston, so that the adiabatic core can become degraded at long ignition delay times. Substantial mixing of colder gas into the interior of the reaction chamber can alter the measurements, similar to effects previously observed for improper piston crevice configuration. In conclusion, experimental techniques to mitigate asynchronous behavior are discussed and demonstrated.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS2019

Nonlinear Görtler Vortices and Their Secondary Instability in a Hypersonic Boundary Layer

Nonlinear development of the Görtler instability over a concave surface gives rise to a highly distorted inflectional flow field in the boundary layer that leads to both wall-normal and spanwise gradients in the flow. Such nonlinear structures are susceptible to strong, high-frequency secondary instabilities that may lead to the onset of laminar-turbulent transition. The present numerical study uses direct numerical simulations and linear secondary instability theory to investigate finite amplitude Görtler vortices and their secondary instability characteristics, respectively, in the hypersonic flow over an axisymmetric cone with a concave aft body. To complement previous studies in the literature wherein the Görtler instability was usually studied for a flat plate and initiated at some upstream location by imposing an eigenfunction as the inflow condition or by blowing and suction at the wall, the present investigation is focused on fully realizable Görtler instability that is excited by an azimuthally periodic array of surface protuberances. Furthermore, while the previous work had mostly focused on the secondary instability of Görtler vortices with cross-plane velocity contours that resembled bell-shaped structures, the present results confirm that fully developed mushroom structures also exist in the hypersonic regime when the Görtler vortex amplitude is sufficiently large. Computations further reveal that the dominant modes of secondary instability in these mushroom-shaped structures correspond to an antisymmetic (i.e., sinuous) “stem” mode that concentrates within the strong, nearly wall-normal internal shear layers surrounding the stem regions underneath the caps of the mushroom structures. Additionally, there exist a multitude of other significantly unstable secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes originating from the Mack mode instability play an important role during the nonlinear breakdown process.

Laminar turbulent transition↗
NASA NTRS2017

Shed Vortex Structure and Phase-Averaged Velocity Statistics in Symmetric/Asymmetric Turbulent Flat Plate Wakes

The near wake of a flat plate is investigated via direct numerical simulations (DNS). Many earlier experimental investigations have used thin plates with sharp trailing edges and turbulent boundary layers to create the wake. This results in large theta divided by D (sub TE) values (theta is the boundary layer momentum thickness towards the end of the plate and D (sub TE) is the trailing edge thickness). In the present study the emphasis is on relatively thick plates with circular trailing edges (CTE) resulting in theta divided by D values less than one (D is the plate thickness and the diameter of the CTE), and vigorous vortex shedding. The Reynolds numbers based on the plate length and D are 1.255 x 10 (sup 6) and 10,000, respectively. Two cases are computed; one with turbulent boundary layers on both the upper and lower surfaces of the plate (statistically the same, symmetric wake, Case TT) and, a second with turbulent and laminar boundary layers on the upper and lower surfaces, respectively (asymmetric case, Case TL). The data and understanding obtained is of considerable engineering interest, particularly in turbomachinery where the pressure side of an airfoil can remain laminar or transitional because of a favorable pressure gradient and the suction side is turbulent. Shed-vortex structure and phase-averaged velocity statistics obtained in the two cases are compared here. The upper negative shed vortices in Case TL (turbulent separating boundary layer) are weaker than the lower positive ones (laminar separating boundary layer) at inception (a factor 1.27 weaker in terms of peak phase-averaged spanwise vorticity at first appearance of a peak). The upper vortices weaken rapidly as they travel downstream. A second feature of interest in Case TL is a considerable increase in the peak phase-averaged, streamwise normal intensity (random component) with increasing streamwise distance (x divided by D) that occurs nears the positive vortex cores. This behavior is observed for a few diameters in the near wake. This is counter to Case TT where the peak value essentially decreases with increasing x divided by D. Both these effects are examined in detail and the important contributors are identified.

direct numerical simulations↗
NASA NTRS2022

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary-layer transition↗
NASA NTRS1994

Navier-Stokes simulations of the effects of suction holes on a flat plate boundary layer

Direct numerical simulations (DNS) of the Navier-Stokes equations are employed to explore the effects of suction holes on transition in a laminar flat plate boundary layer. The Navier-Stokes equations are cast in vorticity-velocity form. Periodicity is imposed in spanwise direction; all other spatial derivatives are discretized with fourth order compact differences. An explicit fourth order Runge-Kutta scheme is employed for the time-integration of the vorticity transport equations. Suction is applied through a row of holes aligned in spanwise direction. For low suction strengths, each hole generates a pair of stable streamwise vortices. When the suction strength exceeds a critical value, the vortices become unstable. For high suction strengths, vortex shedding occurs right at the suction holes. Our numerical findings agree well with experimental observations.

Meitz, Hubert L.↗
DOE OSTI2024

a priori uncertainty quantification of reacting turbulence closure models using Bayesian neural networks

While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulations (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. Here, in this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN, which can be used for out-of-distribution query detection. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Direct simulations of chemically reacting turbulent mixing layers

The results of direct numerical simulations of chemically reacting mixing layers are presented. The reaction considered is a binary, irreversible reaction with no heat release, so that only the effect of the turbulence on the chemical reaction is investigated. The simulation results are shown to be consistent with similarity theory, and are found to be in approximate agreement with laboratory data, even though there are no adjustable constants in the method.

Riley, J. J.↗
NASA NTRS1987

Angular distribution of turbulence in wave space

An alternative to the one-point closure model for turbulence, the large eddy simulation (LES), together with its more exact relative, direct numerical simulation (DNS) are discussed. These methods are beginning to serve as partial substitutes for turbulence experiments. The eddy damped quasi-normal Markovian (EDQNM) theory is reviewed. Angular distribution of the converted data was examined in relationship to EDQNM.

Coleman, G.↗
NASA NTRS1992

Towards understanding turbulent scalar mixing

In an effort towards understanding turbulent scalar mixing, we study the effect of molecular mixing, first in isolation and then by accounting for the effects of the velocity field. The chief motivation for this approach stems from the strong resemblance of the scalar probability density function (PDF) obtained from the scalar field evolving from the heat conduction equation that arises in a turbulent velocity field. However, the evolution of the scalar dissipation is different for the two cases. We attempt to account for these differences, which are due to the velocity field, using a Lagrangian frame analysis. After establishing the usefulness of this approach, we use the heat-conduction simulations (HCS), in lieu of the more expensive direct numerical simulations (DNS), to study many of the less understood aspects of turbulent mixing. Comparison between the HCS data and available models are made whenever possible. It is established that the beta PDF characterizes the evolution of the scalar PDF during mixing from all types of non-premixed initial conditions.

Girimaji, Sharath S.↗
NASA NTRS1991

Towards large eddy and direct simulation of complex turbulent flows

Recent advances in the methodology for direct numerical simulation of turbulent flows and some of the current applications are reviewed. It is argued that high-order finite difference schemes yield solutions with comparable accuracy to the spectral methods with the same number of degrees of freedom. The effects of random inflow conditions on the downstream evolution of turbulence are discussed.

Moin, Parviz↗
NASA NTRS1995

Simulation of crossflow instability on a supersonic highly swept wing

A direct numerical simulation (DNS) algorithm has been developed and validated for use in the investigation of crossflow instability on supersonic swept wings, an application of potential relevance to the design of the High-Speed Civil Transport (HSCT). The algorithm is applied to the investigation of stationary crossflow instability on an infinitely long 77-degree swept wing in Mach 3.5 flow. The results of the DNS are compared with the predictions of linear parabolized stability equation (PSE) methodology. In-general, the DNS and PSE results agree closely in terms of modal growth rate, structure, and orientation angle. Although further validation is needed for large-amplitude (nonlinear) disturbances, the close agreement between independently derived methods offers preliminary validation of both DNS and PSE approaches.

Pruett, C. David↗
NASA NTRS2002

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗
NASA NTRS2005

A Priori Analysis of Subgrid-Scale Models for Large Eddy Simulations of Supercritical Binary-Species Mixing Layers

Models for large eddy simulation (LES) are assessed on a database obtained from direct numerical simulations (DNS) of supercritical binary-species temporal mixing layers. The analysis is performed at the DNS transitional states for heptane/nitrogen, oxygen/hydrogen and oxygen/helium mixing layers. The incorporation of simplifying assumptions that are validated on the DNS database leads to a set of LES equations that requires only models for the subgrid scale (SGS) fluxes, which arise from filtering the convective terms in the DNS equations. Constant-coefficient versions of three different models for the SGS fluxes are assessed and calibrated. The Smagorinsky SGS-flux model shows poor correlations with the SGS fluxes, while the Gradient and Similarity models have high correlations, as well as good quantitative agreement with the SGS fluxes when the calibrated coefficients are used.

turbulence↗
NASA NTRS1988

Direct simulations of chemically reacting turbulent mixing layers, part 2

The results of direct numerical simulations of chemically reacting turbulent mixing layers are presented. This is an extension of earlier work to a more detailed study of previous three dimensional simulations of cold reacting flows plus the development, validation, and use of codes to simulate chemically reacting shear layers with heat release. Additional analysis of earlier simulations showed good agreement with self similarity theory and laboratory data. Simulations with a two dimensional code including the effects of heat release showed that the rate of chemical product formation, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release. Subsequent three dimensional simulations showed similar behavior, in agreement with laboratory observations. Baroclinic torques and thermal expansion in the mixing layer were found to produce changes in the flame vortex structure that act to diffuse the pairing vortices, resulting in a net reduction in vorticity. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers were shown to result from vorticity generation by baroclinic torques.

Metcalfe, Ralph W.↗
NASA NTRS1996

Dynamic Models for LES of Turbulent Front Propagation With a Spectral Method

Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.

Im, H. G.↗
NASA NTRS1993

Numerical simulation of stability and stability control of high speed compressible rotating couette flow

The nonlinear temporal evolution of disturbances in compressible flow between infinitely long, concentric cylinders is investigated through direct numerical simulations of the full, three-dimensional Navier-Stokes and energy equations. Counter-rotating cylinders separated by wide gaps are considered with supersonic velocities of the inner cylinder. Initially, the primary disturbance grows exponentially in accordance with linear stability theory. As the disturbances evolve, higher harmonics and subharmonics are generated in a cascading order eventually reaching a saturation state. Subsequent highly nonlinear stages of the evolution are governed by the interaction of the disturbance modes, particularly the axial subharmonics. Nonlinear evolution of the disturbance field is characterized by the formation of high-shear layers extending from the inner cylinder towards the center of the gap in the form of jets similar to the ejection events in transitional and turbulent wall-bounded shear flows.

Biringen, Sedat↗
NASA NTRS1984

Direct simulations of chemically reacting turbulent mixing layers

The report presents the results of direct numerical simulations of chemically reacting turbulent mixing layers. The work consists of two parts: (1) the development and testing of a spectral numerical computer code that treats the diffusion reaction equations; and (2) the simulation of a series of cases of chemical reactions occurring on mixing layers. The reaction considered is a binary, irreversible reaction with no heat release. The reacting species are nonpremixed. The results of the numerical tests indicate that the high accuracy of the spectral methods observed for rigid body rotation are also obtained when diffusion, reaction, and more complex flows are considered. In the simulations, the effects of vortex rollup and smaller scale turbulence on the overall reaction rates are investigated. The simulation results are found to be in approximate agreement with similarity theory. Comparisons of simulation results with certain modeling hypotheses indicate limitations in these hypotheses. The nondimensional product thickness computed from the simulations is compared with laboratory values and is found to be in reasonable agreement, especially since there are no adjustable constants in the method.

Riley, J. J.↗
NASA NTRS1993

Time accurate simulations of compressible shear flows

The objectives of this research are to employ direct numerical simulation (DNS) to study the phenomenon of mixing (or lack thereof) in compressible free shear flows and to suggest new means of enhancing mixing in such flows. The shear flow configurations under investigation are those of parallel mixing layers and planar jets under both non-reacting and reacting nonpremixed conditions. During the three-years of this research program, several important issues regarding mixing and chemical reactions in compressible shear flows were investigated.

Givi, Peyman↗