Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Ray Traces Through Unsteady Jet Turbulence

Results of an ongoing effort to quantify the role turbulence in scattering sound in jets are reported. Using a direct numerical simulation database to provide the flow data, ray paths traced through the mean flow are compared with those traced through the actual time evolving turbulent flow. Significant scattering by the turbulence is observed. The most notable effect is that upstream traveling waves that are trapped in the potential core by the mean flow, which acts as a wave guide, easily escape in the turbulent flow. A crude statistical estimate based on ray number density suggests that directivity is modified by the turbulence, but no rigorous treatment of non-uniformities in the high-frequency approximation is attempted.

Freund, J. B.↗
DOE OSTI2020

Front propagation formulation for large eddy simulation of turbulent premixed flames

The front propagation formulation (FPF) for large eddy simulation (LES) of turbulent premixed flames is presented and discussed. The FPF method adopts a regularized Dirac delta function to minimize the spurious propagation of under-resolved filtered flame fronts in LES, while reproducing the reaction characteristics of premixed flame fronts. Here, the extension of the FPF method to LES and the relation to existing sub-filter combustion models are discussed. An algebraic model for the regularized Dirac delta function is presented and validated using direct numerical simulation (DNS) data for statistically planar premixed flames propagating in homogeneous isotropic turbulence. A posteriori validation for a filtered 1-D planar laminar flame is performed to demonstrate the characteristics of FPF in reproducing the propagation speed and the thickness of the under-resolved flame front. LES of a turbulent Bunsen flame is also performed to demonstrate the performance of the proposed method.

42 ENGINEERING↗
DOE OSTI2024

Reynolds-number scaling of wall-pressure–velocity correlations in wall-bounded turbulence

Wall-pressure fluctuations are a practically robust input for real-time control systems aimed at modifying wall-bounded turbulence. The scaling behaviour of the wall-pressure–velocity coupling requires investigation to properly design a controller with such input data so that it can actuate upon the desired turbulent structures. A comprehensive database from direct numerical simulations (DNS) of turbulent channel flow is used for this purpose, spanning a Reynolds-number range$Re_\tau \approx 550\unicode{x2013}5200$. Spectral analysis reveals that the streamwise velocity is most strongly coupled to the linear term of the wall pressure, at a Reynolds-number invariant distance-from-the-wall scaling of$\lambda _x/y \approx 14$(and$\lambda _x/y \approx 8$for the wall-normal velocity). When extending the analysis to both homogeneous directions in$x$and$y$, the peak coherence is centred at$\lambda _x/\lambda _z \approx 2$and$\lambda _x/\lambda _z \approx 1$for$p_w$and$u$, and$p_w$and$v$, respectively. A stronger coherence is retrieved when the quadratic term of the wall pressure is concerned, but there is only little evidence for a wall-attached-eddy type of scaling. An experimental dataset comprising simultaneous measurements of wall pressure and velocity complements the DNS-based findings at one value of$Re_\tau \approx 2$k, with ample evidence that the DNS-inferred correlations can be replicated with experimental pressure data subject to significant levels of (acoustic) facility noise. It is furthermore shown that velocity-state estimations can be achieved with good accuracy by including both the linear and quadratic terms of the wall pressure. An accuracy of up to 72 % in the binary state of the streamwise velocity fluctuations in the logarithmic region is achieved; this corresponds to a correlation coefficient of$\approx$0.6. This thus demonstrates that wall-pressure sensing for velocity-state estimation – e.g. for use in real-time control of wall-bounded turbulence – has merit in terms of its realization at a range of Reynolds numbers.

Mechanics↗
DOE OSTI2022

Aero-Optics of Hypersonic Turbulent Boundary Layers

Aero-optics refers to optical distortions due to index-of-refraction gradients that are induced by aerodynamic density gradients. At hypersonic flow conditions, the bulk velocity is many times the speed of sound and density gradients may originate from shock waves, compressible turbulent structures, acoustic waves, thermal variations, etc. Due to the combination of these factors, aero-optic distortions are expected to differ from those common to sub-sonic and lower super-sonic speeds. This report summarizes the results from a 2019-2022 Laboratory Directed Research and Development (LDRD) project led by Sandia National Laboratories in collaboration with the University of Notre Dame, New Mexico State University, and the Georgia Institute of Technology. Efforts extended experimental and simulation methodologies for the study of turbulent hypersonic boundary layers. Notable experimental advancements include development of spectral de-aliasing techniques for highspeed wavefront measurements, a Spatially Selective Wavefront Sensor (SSWFS) technique, new experimental data at Mach 8 and 14, a Quadrature Fringe Imaging Interferometer (QFII) technique for time-resolved index-of-refraction measures, and application of QFII to shock-heated air. At the same time, model advancements include aero-optic analysis of several Direct Numerical Simulation (DNS) datasets from Mach 0.5 to 14 and development of wall-modeled Large Eddy Simulation (LES) techniques for aero-optic predictions. At Mach 8 measured and predicted root mean square Optical Path Differences agree within confidence bounds but are higher than semi-empirical trends extrapolated from lower Mach conditions. Overall, results show that aero-optic effects in the hypersonic flow regime are not simple extensions from prior knowledge at lower speeds and instead reflect the added complexity of compressible hypersonic flow physics.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Some characteristics of small-scale turbulence in a turbulent duct flow

The fine-scale structure of turbulence in a fully developed turbulent duct flow is examined by considering the 3D velocity derivative field obtained from direct numerical simulations at two relatively small Reynolds numbers. The magnitudes of all mean-square derivatives (normalized by wall variables) increase with the Reynolds number, the increase being largest at the wall. These magnitudes are not consistent with the assumption of local isotropy except perhaps near the duct center-line. When the assumption of local isotropy is relaxed to one of local axisymmetry, or invariance with respect to rotation about a coordinate axis (here chosen in the streamwise direction), satisfactory agreement is indicated by the data outside the wall region. Support for axisymmetry is demonstrated by anisotropy invariant maps of the dissipation and vorticity tensors.

Antonia, R. A.↗
NASA NTRS2002

Computation of High Speed Jet Noise

The objective of this work was to use direct numerical simulation (DNS) techniques to study the physics of noise generation by a high-speed turbulent jet. A Mach 0.9, Reynolds number 3,600 jet was selected because of available experimental data. New numerical methods for generating disturbances at the nozzle and computing far-field sound were developed and reported in the course of this work. Over 25 million mesh points were used in the simulations which ran for over 50,000 timesteps and required over 50,000 processor hours on state-of-the-art parallel computer systems to complete. Figures show a visualization of the jet and sound field, a comparison of the mean flow development with the experiment, a directivity comparison with the experiment, and time spectrum comparison with the experiment. Agreement is seen to be excellent. These are fully document in the attached references. Full details of the work, detailed achievements and conclusions are discussed in appendices, which are copies of publications that resulted from this work. We have studied noise mechanisms in supersonic jets, the refraction of sound by turbulence in subsonic jets, and noise sources in conjunction with a DNS of a Mach 0.9 jet.

Freund, Jonathan B.↗
NASA NTRS2000

A Study of Wake Development and Structure in Constant Pressure Gradients

Motivated by the application to high-lift aerodynamics for commercial transport aircraft, a systematic investigation into the response of symmetric/asymmetric planar turbulent wake development to constant adverse, zero, and favorable pressure gradients has been conducted. The experiments are performed at a Reynolds number of 2.4 million based on the chord of the wake generator. A unique feature of this wake study is that the pressure gradients imposed on the wake flow field are held constant. The experimental measurements involve both conventional LDV and hot wire flow field surveys of mean and turbulent quantities including the turbulent kinetic energy budget. In addition, similarity analysis and numerical simulation have also been conducted for this wake study. A focus of the research has been to isolate the effects of both pressure gradient and initial wake asymmetry on the wake development. Experimental results reveal that the pressure gradient has a tremendous influence on the wake development, despite the relatively modest pressure gradients imposed. For a given pressure gradient, the development of an initially asymmetric wake is different from the initially symmetric wake. An explicit similarity solution for the shape parameters of the symmetric wake is obtained and agrees with the experimental results. The turbulent kinetic energy budget measurements of the symmetric wake demonstrate that except for the convection term, the imposed pressure gradient does not change the fundamental flow physics of turbulent kinetic energy transport. Based on the turbulent kinetic energy budget measurements, an approach to correct the bias error associated with the notoriously difficult dissipation estimate is proposed and validated through the comparison of the experimental estimate with a direct numerical simulation result.

Thomas, Flint O.↗
DOE OSTI2024

Multiscale simulation of spatially correlated microstructure via a latent space representation

When deformation gradients act on the scale of the microstructure of a part due to geometry and loading, spatial correlations and finite-size effects in simulation cells cannot be neglected. We propose a multiscale method that accounts for these effects using a variational autoencoder to encode the structure–property map of the stochastic volume elements making up the statistical description of the part. In this paradigm the autoencoder can be used to directly encode the microstructure or, alternatively, its latent space can be sampled to provide likely realizations. Furthermore, we demonstrate the method on three examples using the common additively manufactured material AlSi10Mg in: (a) a comparison with direct numerical simulation of the part microstructure, (b) a push forward of microstructural uncertainty to performance quantities of interest, and (c) a simulation of functional gradation of a part with stochastic microstructure.

Elastoplasticity↗
NASA NTRS1990

A perspective on coherent structures and conceptual models for turbulent boundary layer physics

Direct numerical simulations of turbulent boundary layers have been analyzed to develop a unified conceptual model for the kinematics of coherent motions in low Reynolds number canonical turbulent boundary layers. All classes of coherent motions are considered in the model, including low-speed streaks, ejections and sweeps, vortical structures, near-wall and outer-region shear layers, sublayer pockets, and large-scale outer-region eddies. The model reflects the conclusions from the study of the simulated boundary layer that vortical structures are directly associated with the production of turbulent shear stresses, entrainment, dissipation of turbulence kinetic energy, and the fluctuating pressure field. These results, when viewed from the perspective of the large body of published work on the subject of coherent motions, confirm that vortical structures may be considered the central dynamic element in the maintenance of turbulence in the canonical boundary layer. Vortical structures serve as a framework on which to construct a unified picture of boundary layer structure, providing a means to relate the many known structural elements in a consistent way.

Robinson, Stephen K.↗
NASA NTRS1991

Evolution and dynamics of shear-layer structures in near-wall turbulence

Near-wall flow structures in turbulent shear flows are analyzed, with particular emphasis on the study of their space-time evolution and connection to turbulence production. The results are obtained from investigation of a database generated from direct numerical simulation of turbulent channel flow at a Reynolds number of 180 based on half-channel width and friction velocity. New light is shed on problems associated with conditional sampling techniques, together with methods to improve these techniques, for use both in physical and numerical experiments. The results clearly indicate that earlier conceptual models of the processes associated with near-wall turbulence production, based on flow visualization and probe measurements need to be modified. For instance, the development of asymmetry in the spanwise direction seems to be an important element in the evolution of near-wall structures in general, and for shear layers in particular. The inhibition of spanwise motion of the near-wall streaky pattern may be the primary reason for the ability of small longitudinal riblets to reduce turbulent skin friction below the value for a flat surface.

Johansson, Arne V.↗
NASA NTRS2015

Direct Simulation of Flow over A Wall-Mounted Hump

Simulating separated flows at high Reynolds numbers using Reynolds Averaged Navier-Stokes (RANS) modeled equations remains a challenge in aeronautics. The main hindrance to progress stems from the lack of extended de- tailed data pertinent to the root cause of failure of RANS models, as well as the little progress in RANS modeling innovations in the past several decades. The goal of the current effort is to generate data for separated flow at a Reynolds numbers where conventional models are challenged. We use Direct Numerical Simulations to model turbulent flow over the wall-mounted hump configuration to investigate the physics of flow separation and boundary layer recovery, as well as provide data relevant to the modeling community. A chord-based Reynolds number of Re(sub c) = 47,500 is considered with a turbulent inflow profile of Re(sub θ) = 1,400 (θ /c = 3%). We use FDL3DI, a code that solves the compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard re- cycling/rescaling method of generating turbulent boundary layers as inflow to the computational domain. Two differ- ent configurations of the upper-wall are analyzed for two sets of boundary conditions (slip and no-slip). The results are compared with the available higher Re(sub c) (= 936,000, Re(sub θ) = 7,200, θ /c = 0.77%) experiment for major flow features. The simulated lower Rec allows for DNS-like mesh resolutions, and adequately wide spans. The results from these simulations show earlier separation and delayed reattachment compared to Re(sub c) = 936,000, and significantly higher skin friction in the forebody of the hump. We also find that the upper-wall shape and boundary condition influence pressure distribution over the hump, whereas skin friction is only influenced by the boundary condition.

Wall-Mounted↗
NASA NTRS1999

Large Eddy Simulation in the Computation of Jet Noise

Noise can be predicted by solving Full (time-dependent) Compressible Navier-Stokes Equation (FCNSE) with computational domain. The fluctuating near field of the jet produces propagating pressure waves that produce far-field sound. The fluctuating flow field as a function of time is needed in order to calculate sound from first principles. Noise can be predicted by solving the full, time-dependent, compressible Navier-Stokes equations with the computational domain extended to far field - but this is not feasible as indicated above. At high Reynolds number of technological interest turbulence has large range of scales. Direct numerical simulations (DNS) can not capture the small scales of turbulence. The large scales are more efficient than the small scales in radiating sound. The emphasize is thus on calculating sound radiated by large scales.

Mankbadi, R. R.↗
DOE OSTI2025

Calibrating a finite-strain phase-field model of fracture for bonded granular materials with uncertainty quantification

To study the mechanical behavior of mock high explosives, an experimental and simulation program was developed to calibrate, with quantified uncertainty, a material model of the bonded granular material Idoxuridine and nitroplasticized Estane-5703. This paper reports on the efficacy of such a framework as a generalizable methodology for calibrating material models against experimental data with uncertainty quantification. Additionally, this paper studies the effect of two manufacturing temperatures and three initial granular configurations on the unconfined compressive behavior of the resulting bonded granular materials. In each of these cases, the same calibration framework was used; in that, hundreds of high-fidelity direct numerical simulations using a new, graphics processing unit-enabled, high-performance finite element method software, Ratel, were run to calibrate a finite-strain phase-field fracture model against experimental data. It was found that manufacturing temperature influenced the elastic response of the mock high explosives, with higher temperatures yielding a stiffer response. By contrast, it was found that the initial configuration of the grains had a negligible impact on the overall behavior of the mock high explosives though it remains possible that local damage accumulation within the specimens could be altered by the initial configurations. Overall, the calibration framework was successful at creating well-calibrated models, showing its usefulness as an engineering and scientific tool.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1993

Direct simulation of isothermal-wall supersonic channel flow

The motivation for this work is the fact that in turbulent flows where compressibility effects are important, they are often poorly understood. A few examples of such flows are those associated with astrophysical phenomena and those found in combustion chambers, supersonic diffusers and nozzles, and over high-speed airfoils. For this project, we are primarily interested in compressibility effects near solid surfaces. Our main objective is an improved understanding of the fundamentals of compressible wall-bounded turbulence, which can in turn be used to cast light upon modeling concepts such as the Morkovin hypothesis and the Van Driest transformation. To this end, we have performed a direct numerical simulation (DNS) study of supersonic turbulent flow in a plane channel with constant-temperature walls. All of the relevant spatial and temporal scales are resolved so that no sub grid scale or turbulence model is necessary. The channel geometry was chosen so that finite Mach number effects can be isolated by comparing the present results to well established incompressible channel data. Here the fluid is assumed to be an ideal gas with constant specific heats, constant Prandtl number, and power-law temperature-dependent viscosity. Isothermal-wall boundary conditions are imposed so that a statistically stationary state may be obtained. The flow is driven by a uniform (in space) body force (rather than a mean pressure gradient) to preserve stream wise homogeneity, with the body force defined so that the total mass flux is constant.

Coleman, Gary N.↗
NASA NTRS2000

Numerical Simulations of Homogeneous Turbulence Using Lagrangian-Averaged Navier-Stokes Equations

The Lagrangian-averaged Navier-Stokes (LANS) equations are numerically evaluated as a turbulence closure. They are derived from a novel Lagrangian averaging procedure on the space of all volume-preserving maps and can be viewed as a numerical algorithm which removes the energy content from the small scales (smaller than some a priori fixed spatial scale alpha) using a dispersive rather than dissipative mechanism, thus maintaining the crucial features of the large scale flow. We examine the modeling capabilities of the LANS equations for decaying homogeneous turbulence, ascertain their ability to track the energy spectrum of fully resolved direct numerical simulations (DNS), compare the relative energy decay rates, and compare LANS with well-accepted large eddy simulation (LES) models.

Mohseni, Kamran↗