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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

On the Universality of the Kolmogorov Constant in Numerical Simulations of Turbulence

Motivated by a recent survey of experimental data, we examine data on the Kolmogorov spectrum constant in numerical simulations of isotropic turbulence, using results both from previous studies and from new direct numerical simulations over a range of Reynolds numbers (up to 240 on the Taylor scale) at grid resolutions up to 512(exp 3). It is noted that in addition to k(exp -5/3) scaling, identification of a true inertial range requires spectral isotropy in the same wavenumber range. We found that a plateau in the compensated three-dimensional energy spectrum at k(eta) approx. = 0.1 - -0.2, commonly used to infer the Kolmogorov constant from the compensated three-dimensional energy spectrum, actually does not represent proper inertial range behavior. Rather, a proper, if still approximate, inertial range emerges at k(eta) approx. = 0.02 - 0.05 when R(sub lambda) increases beyond 140. The new simulations indicate proportionality constants C(sub 1) and C in the one- and three-dimensional energy spectra respectively about 0.60 and 1.62. If the turbulence were perfectly isotropic then use of isotropy relations in wavenumber space (C(sub 1) = 18/55 C) would imply that C(sub 1) approx. = 0.53 for C = 1.62, in excellent agreement with experiments. However the one- and three-dimensional estimates are not fully consistent, because of departures (due to numerical and statistical limitations) from isotropy of the computed spectra at low wavenumbers. The inertial scaling of structure functions in physical space is briefly addressed. Since DNS is still restricted to moderate Reynolds numbers, an accurate evaluation of the Kolmogorov constant is very difficult. We focus on providing new insights on the interpretation of Kolmogorov 1941 similarity in the DNS literature and do not consider issues pertaining to the refined similarity hypotheses of Kolmogorov (K62).

Yeung, P. K.↗
NASA NTRS1988

Studying Turbulence Using Numerical Simulation Databases, 2. Proceedings of the 1988 Summer Program

The focus of the program was on the use of direct numerical simulations of turbulent flow for study of turbulence physics and modeling. A special interest was placed on turbulent mixing layers. The required data for these investigations were generated from four newly developed codes for simulation of time and spatially developing incompressible and compressible mixing layers. Also of interest were the structure of wall bounded turbulent and transitional flows, evaluation of diagnostic techniques for detection of organized motions, energy transfer in isotropic turbulence, optical propagation through turbulent media, and detailed analysis of the interaction of vortical structures.

Source record↗
DOE OSTI2024

High-Fidelity Simulation of the Light-to-Dense Stratification Transient in the HiRJET Facility

Density stratification in a large enclosure is a crucial phenomenon to heat transfer and sustainable passive heat removal of a sodium fast reactor during reactivity transients. However, engineering turbulence models were identified to have unsatisfactory performance in predicting propagation of a stratified front. Yet, the scarcity of high-resolution data for stratification hampers the development of models. To explor e applications of leveraging direct numerical simulation (DNS) data to support turbulence model development, this work conducted DNS using NekRS to study a long stratification transient in the High-Resolution Jet (HiRJET) experimental facility. Here, this work considers an experiment run where light fluid is injected into a tank containing a denser fluid with a relative density difference of 1.5%. Formation of the stratified layer is identified as impingement of the buoyant jet promoting mixing of the two fluids. Based on the transient statistics, transport of the concentration can be characterized by regions with dominating effects of turbulent mixing, buoyant dissipation, and molecular diffusion, respectively, as moving away from the elevation of jet impingement. Concentration near the stratified front also exhibits oscillation at Brunt-Väisälä frequency. Preliminary validation of the simulation showed encouraging agreement of the concentration distribution with the reference experiment.

buoyant mixing↗
DOE OSTI2020

Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays

The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Heat transfer and transport property contrast effects on the compressible Rayleigh-Taylor instability

In extreme environments, the Rayleigh-Taylor instability (RTI) may occur under large variations in density and temperature and with fluid transport properties strongly dependent on temperature. Direct numerical simulations of the 3D fully compressible RTI are conducted, examining the idealized configuration of a hotter, less dense fluid pushing against a colder, denser fluid. Various temperature ratios and transport property configurations are explored to examine how heat conduction, large variations in transport properties, and sudden changes in transport properties can affect the evolution of the mixing layer. Nonuniform fluid expansion and contraction induced by heat transfer can significantly affect local density differences and overall instability growth, causing profile asymmetries about the initial interface for flow and mixing statistics. The departures from classical self-similar development of the instability along with misalignment between regions of mixing and regions of most intense turbulent activity caused by both heat transfer and transport property contrasts are examined. After sudden changes in fluid transport properties, which may occur as a result of rapid heating (e.g., in inertial confinement fusion), the flow quickly responds and begins to relax towards quasi-self-similar late-time evolution. For many dynamical quantities such as vorticity and dissipation, this late-time evolution resembles that of the configuration that already started with the final transport property magnitudes, suggesting that these quantities depend only on the transport properties and not on past flow history, provided that the density field distributions for the flows remain similar. On the other hand, the mixing evolution after the transport property change is unique, implying that both property magnitudes and previous history are impactful on the mixing. In conclusion, these simulations demonstrate how various temperature-related effects are extremely important to consider in compressible RTI flows with large temperature variations.

42 ENGINEERING↗
NASA NTRSDate not supplied

Assessment of UMUSCL Scheme for DNS of Turbulent Flows

Direct numerical simulation (DNS) is performed using an unstructured node-based finite-volume methodology for canonical turbulent flow configurations to assess the effect of the upwinding parameter (κ) in the unstructured monotonic upstream scheme for conservation laws (UMUSCL) reconstruction scheme. The NASA FUN3D solver is used, and four turbulent validation cases are considered: periodic Taylor-Green vortex, (2) flow through a plane channel, (3) flow through a channel with a constriction, and (4) flow over a flat plate. The results are presented for κ = 0.5, 0.9 and 1 on hexahedral grids, with the last value corresponding to a zero-numerical dissipation scheme. Detailed comparisons with available reference data indicate that for the smallest length-scales simulated, κ = 0.5 is too dissipative while κ = 1 could lead to an energy build-up. Overall, κ =0.9 yields satisfactory results in terms of accuracy and robustness for all the cases considered in this study.

Direct Numerical Simulation↗
NASA NTRS1995

Linear Wave Motion from Concentrated Harmonic Sources in Blasius Flow

The motion of individual linear instability waves in shear flows is well described by existing theoretical and numerical methods. However, naturally occuring sources produce coherent wave motions with broadband spanwise wavenumber and frequency spectra, and the different spectral components interact both linearly and nonlinearily. This paper describes a series of calculations for the parameters of three different experiments using locally parallel linear stability theory (LST), the Parabolized Stability Equations (PSE), and Direct Numerical Simulation (DNS). The calculations illustrate the strengths and weaknesses of the different methods, the extent to which the methods agree or disagree, and, finally the extent to which agreement with the measurements can be attained, given that the experiments also have their own difficulties.

linear↗
NASA NTRS2019

Laminar-Turbulent Transition Upstream of the Entropy-Layer Swallowing Location in Hypersonic Boundary Layers

Numerical and experimental studies have demonstrated that modal growth of planar Mack modes is responsible for laminar-turbulent transition on sharp cones at hypersonic speeds. However, the physical mechanisms that lead to transition onset upstream of the entropy-layer swallowing location over sufficiently blunt geometries are not well understood as yet. Modal amplification is too weak or nonexistent to initiate transition at moderate-to-large bluntness values. Nonmodal analysis shows that, with increasing nose bluntness, both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification. However, because of the relatively weak signature of the nonmodal traveling disturbances within the boundary-layer region, the route to transition onset subsequent to the nonmodal growth remains unclear. Thus, nonlinear parabolized stability equations (NPSE) and direct numerical simulations (DNS) have been used to investigate the potential transition mechanisms over a 7-degree blunt cone that was tested in the AFRL Mach-6 high-Reynoldsnumber facility. Computations are performed to separately follow the nonlinear development of two classes of inflow disturbances, namely, a pair of oblique traveling waves with equal but opposite angles with respect to the mean flow direction and a planar traveling wave. Results in both cases show an excellent agreement between the NPSE and DNS predictions, establishing that the NPSE is an accurate and efficient technique for predicting the nonlinear development for these particular nonmodal traveling disturbances. Computations reveal that the oblique mode interactions lead to the generation of stationary streaks inside the boundary layer that, in turn, facilitate the growth of a subharmonic sinuous disturbance. For relatively modest amplitudes of the inflow disturbance, the oblique-mode breakdown can lead to transition at the measured location of transition onset during the experiment. On the other hand, the nonlinear development of a planar traveling wave leads to the formation of inclined structures just above the boundary-layer edge and these structures are strongly reminiscent of the transitional events observed during blunt cone experiments by using schlieren flow visualizations.

Entropy↗
NASA NTRS1995

Annual Research Briefs: 1995

This report contains the 1995 annual progress reports of the Research Fellows and students of the Center for Turbulence Research (CTR). In 1995 CTR continued its concentration on the development and application of large-eddy simulation to complex flows, development of novel modeling concepts for engineering computations in the Reynolds averaged framework, and turbulent combustion. In large-eddy simulation, a number of numerical and experimental issues have surfaced which are being addressed. The first group of reports in this volume are on large-eddy simulation. A key finding in this area was the revelation of possibly significant numerical errors that may overwhelm the effects of the subgrid-scale model. We also commissioned a new experiment to support the LES validation studies. The remaining articles in this report are concerned with Reynolds averaged modeling, studies of turbulence physics and flow generated sound, combustion, and simulation techniques. Fundamental studies of turbulent combustion using direct numerical simulations which started at CTR will continue to be emphasized. These studies and their counterparts carried out during the summer programs have had a noticeable impact on combustion research world wide.

Source record↗
DOE OSTI2025

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Numerical Studies of Turbulent Supersonic Plane-Channel Flows

This two-part talk begins with an overview of the `turbulence problem’, the computational strategies used to address it, and the role of direct numerical simulation (DNS). The second part presents DNS results to determine the turbulent Prandtl number Pr t above cold (isothermal) and hot (adiabatic) walls in a family of low-supersonic channel flows. A range of mean temperature/density variations, corresponding to effective/edge Mach numbers between 1.1 to 2.2, and wall-variable-based Reynolds number Reτw from 73 to 3800, is considered. The adiabatic condition is a new feature of special interest. The value of Pr t away from the wall approaches 0.85 above both the isothermal and adiabatic walls. The variations of the near-wall Pr t profiles in both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗wall scaling proposed by Huang et al. (1995), with only a weak dependence on R eτw . This leads to a rather simple proposal for a model of heat transfer, attached to an eddy-viscosity model.

Direct numerical simulation↗
DOE OSTI2022

Physics-informed Machine Learning for Modeling Turbulence in Supernovae

Abstract Turbulence plays an important role in astrophysical phenomena, including core-collapse supernovae (CCSNe), but current simulations must rely on subgrid models, since direct numerical simulation is too expensive. Unfortunately, existing subgrid models are not sufficiently accurate. Recently, machine learning (ML) has shown an impressive predictive capability for calculating turbulence closure. We have developed a physics-informed convolutional neural network to preserve the realizability condition of the Reynolds stress that is necessary for accurate turbulent pressure prediction. The applicability of the ML subgrid model is tested here for magnetohydrodynamic turbulence in both the stationary and dynamic regimes. Our future goal is to utilize this ML methodology (available on GitHub) in the CCSN framework to investigate the effects of accurately modeled turbulence on the explosion of these stars.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Model-free simulations of turbulent reactive flows

The current computational methods for solving transport equations of turbulent reacting single-phase flows are critically reviewed, with primary attention given to those methods that lead to model-free simulations. In particular, consideration is given to direct numerical simulations using spectral (Galerkin) and pseudospectral (collocation) methods, spectral element methods, and Lagrangian methods. The discussion also covers large eddy simulations and turbulence modeling.

Givi, Peyman↗
NASA NTRS2000

CVS Filtering of 3D Turbulent Mixing Layers Using Orthogonal Wavelets

Coherent Vortex Simulation (CVS) filtering has been applied to Direct Numerical Simulation (DNS) data of forced and unforced time-developing turbulent mixing layers. CVS filtering splits the turbulent flow into two orthogonal parts, one corresponding to coherent vortices and the other to incoherent background flow. We have shown that the coherent vortices can be represented by few wavelet modes and that these modes are sufficient to reproduce the vorticity probability distribution function (PDF) and the energy spectrum over the entire inertial range. The remaining incoherent background flow is homogeneous, has small amplitude, and is uncorrelated. These results are compared with those obtained for the same compression rate using large eddy simulation (LES) filtering. In contrast to the incoherent background flow of CVS filtering, the LES subgrid scales have a much larger amplitude and are correlated, which makes their statistical modeling more difficult.

Schneider, Kai↗
NASA NTRS2012

Large Eddy/Reynolds-Averaged Navier-Stokes Simulations of CUBRC Base Heating Experiments

ven with great advances in computational techniques and computing power during recent decades, the modeling of unsteady separated flows, such as those encountered in the wake of a re-entry vehicle, continues to be one of the most challenging problems in CFD. Of most interest to the aerothermodynamics community is accurately predicting transient heating loads on the base of a blunt body, which would result in reduced uncertainties and safety margins when designing a re-entry vehicle. However, the prediction of heat transfer can vary widely depending on the turbulence model employed. Therefore, selecting a turbulence model which realistically captures as much of the flow physics as possible will result in improved results. Reynolds Averaged Navier Stokes (RANS) models have become increasingly popular due to their good performance with attached flows, and the relatively quick turnaround time to obtain results. However, RANS methods cannot accurately simulate unsteady separated wake flows, and running direct numerical simulation (DNS) on such complex flows is currently too computationally expensive. Large Eddy Simulation (LES) techniques allow for the computation of the large eddies, which contain most of the Reynolds stress, while modeling the smaller (subgrid) eddies. This results in models which are more computationally expensive than RANS methods, but not as prohibitive as DNS. By complimenting an LES approach with a RANS model, a hybrid LES/RANS method resolves the larger turbulent scales away from surfaces with LES, and switches to a RANS model inside boundary layers. As pointed out by Bertin et al., this type of hybrid approach has shown a lot of promise for predicting turbulent flows, but work is needed to verify that these models work well in hypersonic flows. The very limited amounts of flight and experimental data available presents an additional challenge for researchers. Recently, a joint study by NASA and CUBRC has focused on collecting heat transfer data on the backshell of a scaled model of the Orion Multi-Purpose Crew Vehicle (MPCV). Heat augmentation effects due to the presence of cavities and RCS jet firings were also investigated. The high quality data produced by this effort presents a new set of data which can be used to assess the performance of CFD methods. In this work, a hybrid LES/RANS model developed at North Carolina State University (NCSU) is used to simulate several runs from these experiments, and evaluate the performance of high fidelity methods as compared to more typical RANS models. .

Salazar, Giovanni↗
NASA NTRSDate not supplied

WMLES of Boundary-Layer Transition

The potential for the accurate prediction of laminar-turbulent transition by Wall-Modeled Large-Eddy Simulation (WMLES) is investigated through analysis of Direct Numerical Simulation (DNS) and WMLES results. The analysis indicates satisfactory prediction of transition by WMLES is possible, provided a) a sufficiently accurate wall model for the initial stage of transition is employed, b) a satisfactory transition is made between this initial-stage wall model and the wall model for the downstream fully-developed turbulence, and c) limitations on temporal and spatial resolution, and on the thickness of the wall-modeled region, imposed by the need to accurately resolve the disturbance waves and the transition region are observed. This is illustrated with a WMLES computation that yields reasonable agreement with DNS results; however, the wall-modeled region thickness, and thus the cost savings, is less than is typical for WMLES of fully developed turbulence.

wall-modeled large-eddy simulation, boundary-layer↗
NASA NTRS2019

Effect of 3D Roughness Patch on Instability Amplification in a Supersonic Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of high-speed vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, sinusoidal shape roughness elements on a Mach 3.5 flat plate boundary layer for flow conditions corresponding to the planned conditions of an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at the NASA Langley Research Center. Analysis of convective instabilities in the wake of the roughness patch was reported in a previous paper and the current work extends that analysis to instability amplification across the length of the roughness patch. Quasi-parallel stability analysis of the modified boundary layer flow over the patch indicates two dominant families of unstable disturbances, namely, a group of high frequency modes that amplify in localized regions along the roughness patch and another group of low frequency modes that have smaller peak amplification rates but amplify steadily both above the roughness patch and in the wake region behind it. The results suggest that the amplification factors associated with the high-frequency modes are sufficiently low, at least for the roughness patches considered in this paper, so that these modes are unlikely to have a major influence on the transition process. The amplification of the low-frequency modes within the region of the roughness patch is further quantified via direct numerical simulations. Results confirm the strongly destabilizing influence of the roughness patch on the first mode instabilities, yielding an N-factor increment of ΔN ≈ 3.6 for a roughness patch length of eight wavelengths.

Laminar turbulent transition↗