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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS2002

An Evaluation of Linear Instability Waves as Sources of Sound in a Supersonic Turbulent Jet

Mach wave radiation from supersonic jets is revisited to better quantify the extent to which linearized equations represent the details of the actual mechanism. To this end, we solve the linearized Navier-Stokes equations (LNS) with precisely the same mean flow and inflow disturbances as a previous direct numerical simulation (DNS) of a perfectly expanded turbulent M = 1.92 jet. We restrict our attention to the first two azimuthal modes, n = 0 and n = 1, which constitute most of the acoustic field. The direction of peak radiation and the peak Strouhal number matches the DNS reasonably well, which is in accord with previous experimental justification of the linear theory. However, it is found that the sound pressure level predicted by LNS is significantly lower than that from DNS. In order to investigate the discrepancy, individual frequency components of the solution are examined. These confirm that near the peak Strouhal number, particularly for the first helical mode n = 1, the amplification of disturbances in the LNS closely matches the DNS. However, away from the peak frequency (and generally for the azimuthal mode n = 0), modes in the LNS are damped while those in the DNS grow at rates comparable to those at the peak Strouhal number.

Mohseni, Kamran↗

DOE OSTI2025

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗

NASA NTRS1996

Multiple Scale Simulation for Transitional and Turbulent Flow

A new concept, Multiple Scale Simulation (MSS), is presented in this paper. The basic idea is that the flow is decomposed into several component groups according to spatial and temporal length scales. Each group has its own subdomain, governing system, mesh size, and discretization method. The simulation is then performed groupwise. This approach has been successfully applied in combination with the intergrid dissipation technique for simulation of transitional and turbulent flow in 3-D boundary layers, and it is feasible for 3-D airfoils and other more complex configurations. MSS should prove to ameliorate the scale problems associated with conventional direct numerical simulation.

Liu, Chaoqun↗

NASA NTRS1995

Large-eddy simulation of flow past a circular cylinder

Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Mittal, R.↗

DOE OSTI2026

Neural operator transformers capture bifurcating drift-wave turbulence in fusion plasma simulations

Self-consistent modeling of turbulence-driven transport is critical for optimizing confinement in magnetically confined fusion plasmas, such as tokamaks and stellarators. In particular, capturing the long-term co-evolution of turbulence, flow, and background plasma profiles remains computationally challenging. Direct numerical simulation of these multiscale, highly nonlinear processes is often demanding and impractical for real-time control or design optimization. To address this bottleneck, we investigate transformer-based neural operator partial differential equation surrogates for emulating the dynamics of drift-wave turbulence bifurcation mediated by zonal flows, using the modified Hasegawa–Wakatani (MHW) model as a prototypical system. We find that the finetuned neural operator model has excellent performance in capturing the multi-spatiotemporal-scales of MHW turbulence bifurcation and is robust to testing on rare and out-of-distribution dynamics. Specifically, we demonstrate that a single unified model accurately predicts both quasi-steady-state turbulence and a wide range of dynamical transition processes, such as nonlinear saturation, spontaneous suppression of turbulence, and the emergence of macroscopic zonal flows, over time horizons vastly exceeding the local turbulence correlation time. This computationally efficient approach establishes a strong foundation for fast, AI-based modeling of complex, multiscale phenomena in magnetized fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1996

Hydroacoustic forcing function modeling using DNS database

A wall pressure frequency spectrum model (Blake 1971 ) has been evaluated using databases from Direct Numerical Simulations (DNS) of a turbulent boundary layer (Na & Moin 1996). Good agreement is found for moderate to strong adverse pressure gradient flows in the absence of separation. In the separated flow region, the model underpredicts the directly calculated spectra by an order of magnitude. The discrepancy is attributed to the violation of the model assumptions in that part of the flow domain. DNS computed coherence length scales and the normalized wall pressure cross-spectra are compared with experimental data. The DNS results are consistent with experimental observations.

Zawadzki, I.↗

NASA NTRS1996

Simulation of Nonlinear Instabilities in an Attachment-Line Boundary Layer

The linear and the nonlinear stability of disturbances that propagate along the attachment line of a three-dimensional boundary layer is considered. The spatially evolving disturbances in the boundary layer are computed by direct numerical simulation (DNS) of the unsteady, incompressible Navier-Stokes equations. Disturbances are introduced either by forcing at the in ow or by applying suction and blowing at the wall. Quasi-parallel linear stability theory and a nonparallel theory yield notably different stability characteristics for disturbances near the critical Reynolds number; the DNS results con rm the latter theory. Previously, a weakly nonlinear theory and computations revealed a high wave-number region of subcritical disturbance growth. More recent computations have failed to achieve this subcritical growth. The present computational results indicate the presence of subcritically growing disturbances; the results support the weakly nonlinear theory. Furthermore, an explanation is provided for the previous theoretical and computational discrepancy. In addition, the present results demonstrate that steady suction can be used to stabilize disturbances that otherwise grow subcritically along the attachment line.

Joslin, Ronald D.↗

NASA NTRS2009

Criteria for Modeling in LES of Multicomponent Fuel Flow

A report presents a study addressing the question of which large-eddy simulation (LES) equations are appropriate for modeling the flow of evaporating drops of a multicomponent liquid in a gas (e.g., a spray of kerosene or diesel fuel in air). The LES equations are obtained from the direct numerical simulation (DNS) equations in which the solution is computed at all flow length scales, by applying a spatial low-pass filter. Thus, in LES the small scales are removed and replaced by terms that cannot be computed from the LES solution and instead must be modeled to retain the effect of the small scales into the equations. The mathematical form of these models is a subject of contemporary research. For a single-component liquid, there is only one LES formulation, but this study revealed that for a multicomponent liquid, there are two non-equivalent LES formulations for the conservation equations describing the composition of the vapor. Criteria were proposed for selecting the multicomponent LES formulation that gives the best accuracy and increased computational efficiency. These criteria were applied in examination of filtered DNS databases to compute the terms in the LES equations. The DNS databases are from mixing layers of diesel and kerosene fuels. The comparisons resulted in the selection of one of the multicomponent LES formulations as the most promising with respect to all criteria.

Bellan, Josette↗

NASA NTRS1990

Pre-mixed flame simulations for non-unity Lewis numbers

A principal effect of turbulence on premixed flames in the flamelet region is to wrinkle the flame fronts. For non-unity Lewis numbers (Le), the local flame structure is altered in curved regions. This effect is examined using direct numerical simulations of the three dimensional, constant density, decaying isotropic turbulence with a single step, finite rate chemical reaction. Simulations of Lewis numbers 0.8, 1.0, and 1.2 are compared. The turbulent flame speed, S(sub T), increases as Le decreases. The correlation between S(sub T) and u prime found in previous Le = 1 simulations has a strong Lewis number dependency. The variance of the pdf of the flame curvature increases as Le decreases, indicating that the flames become more wrinkled. A strong correlation between local flame speed and curvature was found. For Le greater than 1, the flame speed increases in regions concave towards the products and decreases in convex regions. The opposite correlation was found for Le less than 1. The mean temperature of the products was also found to vary with Lewis number. For Le = 0.8, it is less than the adiabatic flame temperature and for Le = 1.2 it is greater.

Rutland, C. J.↗

NASA NTRS1990

Subgrid-scale backscatter in transitional and turbulent flows

Most subgrid-scale (SGS) models for large-eddy simulations are absolutely dissipative (that is, they remove energy from the large scales at each point in the physical space). The actual SGS stresses, however, may transfer energy to the large scales (backscatter) at a given location. Direct numerical simulations of turbulent channel flow and compressible isotropic turbulence are used to study the backscatter phenomena. In all flows considered roughly 50 percent of the grid points were experiencing backscatter when a Fourier cutoff filter was used. The backscatter fraction was less with a Gaussian filter, and intermediate with a box filter in physical space. Moreover, the backscatter and forward scatter contributions to the SGS dissipation were comparable, and each was often much larger than the total SGS dissipation. The SGS dissipation (normalized by total dissipation) increased with filter width almost independently of filter type and Reynolds number. The amount of backscatter showed an increasing trend with Reynolds numbers. In the near-wall region of the channel, events characterized by strong Reynolds shear stress correlated fairly well with areas of high SGS dissipation (both forward and backward). In compressible isotropic turbulence similar results were obtained, independent of fluctuation Mach number.

Piomelli, Ugo↗

NASA NTRS1992

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗

NASA NTRS2000

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.↗

NASA NTRS2010

Detailed Comparison of DNS to PSE for Oblique Breakdown at Mach 3

A pair of oblique waves at low amplitudes is introduced in a supersonic flat-plate boundary layer. Their downstream development and the concomitant process of laminar to turbulent transition is then investigated numerically using Direct Numerical Simulations (DNS) and Parabolized Stability Equations (PSE). This abstract is the last part of an extensive study of the complete transition process initiated by oblique breakdown at Mach 3. In contrast to the previous simulations, the symmetry condition in the spanwise direction is removed for the simulation presented in this abstract. By removing the symmetry condition, we are able to confirm that the flow is indeed symmetric over the entire computational domain. Asymmetric modes grow in the streamwise direction but reach only small amplitude values at the outflow. Furthermore, this abstract discusses new time-averaged data from our previous simulation CASE 3 and compares PSE data obtained from NASA's LASTRAC code to DNS results.

Mayer, Christian S. J.↗

NASA NTRS1993

Direct simulations of premixed turbulent flames with nonunity Lewis numbers

A principal effect of turbulence on premixed flames in the flamelet regime is to wrinkle the flame fronts. For nonunity Lewis numbers, Le is not equal to 1, the local flame structure is altered in curved regions. This effect is examined using direct numerical simulations of 3D isotropic turbulence with constant density, single-step Arrhenius kinetics chemistry. Simulations of Lewis numbers 0.8, 1.0, and 1.2 are compared. At the local level, curvature effects dominated changes to the flame structure while strain effects were insignificant. A strong Lewis-number-dependent correlation was found between surface curvature and the local flame speed. The correlation was positive for Le less than 1 and negative for Le greater than 1. At the global level, strain-related effects were more significant than curvature effects. The turbulent flame speed changed significantly with Lewis number, increasing as Le decreased. This was found to be due to strain effects that have a nonzero mean over the flame surface, rather than to curvature effects that have a nearly zero mean. The mean product temperature was also found to vary with Lewis number, being higher for Le greater than 1 and lower for Le less than 1.

Rutland, C. J.↗

NASA NTRS1994

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

The Reynolds numbers that characterize geophysical and astrophysical turbulence (Re approximately equals 10(exp 8) for the planetary boundary layer and Re approximately equals 10(exp 14) for the Sun's interior) are too large to allow a direct numerical simulation (DNS) of the fundamental Navier-Stokes and temperature equations. In fact, the spatial number of grid points N approximately Re(exp 9/4) exceeds the computational capability of today's supercomputers. Alternative treatments are the ensemble-time average approach, and/or the volume average approach. Since the first method (Reynolds stress approach) is largely analytical, the resulting turbulence equations entail manageable computational requirements and can thus be linked to a stellar evolutionary code or, in the geophysical case, to general circulation models. In the volume average approach, one carries out a large eddy simulation (LES) which resolves numerically the largest scales, while the unresolved scales must be treated theoretically with a subgrid scale model (SGS). Contrary to the ensemble average approach, the LES+SGS approach has considerable computational requirements. Even if this prevents (for the time being) a LES+SGS model to be linked to stellar or geophysical codes, it is still of the greatest relevance as an 'experimental tool' to be used, inter alia, to improve the parameterizations needed in the ensemble average approach. Such a methodology has been successfully adopted in studies of the convective planetary boundary layer. Experienc e with the LES+SGS approach from different fields has shown that its reliability depends on the healthiness of the SGS model for numerical stability as well as for physical completeness. At present, the most widely used SGS model, the Smagorinsky model, accounts for the effect of the shear induced by the large resolved scales on the unresolved scales but does not account for the effects of buoyancy, anisotropy, rotation, and stable stratification. The latter phenomenon, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Canuto, V. M.↗

NASA NTRS2013

Towards Bridging the Gaps in Holistic Transition Prediction via Numerical Simulations

The economic and environmental benefits of laminar flow technology via reduced fuel burn of subsonic and supersonic aircraft cannot be realized without minimizing the uncertainty in drag prediction in general and transition prediction in particular. Transition research under NASA's Aeronautical Sciences Project seeks to develop a validated set of variable fidelity prediction tools with known strengths and limitations, so as to enable "sufficiently" accurate transition prediction and practical transition control for future vehicle concepts. This paper provides a summary of selected research activities targeting the current gaps in high-fidelity transition prediction, specifically those related to the receptivity and laminar breakdown phases of crossflow induced transition in a subsonic swept-wing boundary layer. The results of direct numerical simulations are used to obtain an enhanced understanding of the laminar breakdown region as well as to validate reduced order prediction methods.

Choudhari, Meelan M.↗

NASA NTRS1987

On the subgrid-scale modeling of compressible turbulence

A subgrid-scale model recently derived for use in the large-eddy simulation of compressible turbulent flows is examined from a fundamental theoretical and computational standpoint. It is demonstrated that this model, which is applicable only to compressible turbulent flows in the limit of small density fluctuations, correlates somewhat poorly with the results of direct numerical simulations of compressible isotropic turbulence at low Mach numbers. An alternative model, based on Favre-filtered fields, is suggested which appears to reduce these limitations.

Speziale, C. G.↗

NASA NTRS1988

The subgrid-scale modeling of compressible turbulence

Speziale, C. G.↗