SEARCH · Engineering Papers

Results for “Direct numerical simulations”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Direct numerical simulation of transitional and turbulent flow over a heated flat plate using finite-difference schemes

This report deals with the direct numerical simulation of transitional and turbulent flow at low Mach numbers using high-order-accurate finite-difference techniques. A computation of transition to turbulence of the spatially-evolving boundary layer on a heated flat plate in the presence of relatively high freestream turbulence was performed. The geometry and flow conditions were chosen to match earlier experiments. The development of the momentum and thermal boundary layers was documented. Velocity and temperature profiles, as well as distributions of skin friction, surface heat transfer rate, Reynolds shear stress, and turbulent heat flux, were shown to compare well with experiment. The results indicate that the essential features of the transition process have been captured. The numerical method used here can be applied to complex geometries in a straightforward manner.

Madavan, Nateri K.↗

NASA NTRS1990

Direct numerical simulation of a three-dimensional turbulent boundary layer

The effects of transverse strain on an initially two-dimensional turbulent boundary layer are studied in a direct numerical simulation of a planar channel flow with impulsively started transverse pressure gradient. Consistent with experiments in three-dimensional boundary layers, the simulation shows a decrease in the Reynolds shear stress with increasing transverse strain. Also, the directions of the Reynolds shear stress vector and the mean velocity gradient vector were found to differ. In addition, the simulation shows a drop in the turbulent kinetic energy. Terms in the Reynolds stress transport equations were computed. The balances indicate that the decrease in turbulent kinetic energy is a result of a decrease in turbulence production, along with an increase in turbulent dissipation. Intuitive reasoning and current turbulence models would predict an increase in kinetic energy along with increases in production and dissipation rates as a result of increased mean-flow strain rate. Later in the evolution of the flow, both turbulence production and dissipation increase.

Moin, P.↗

NASA NTRS1993

Direct numerical simulation of turbulent flow over a backward-facing step

The objectives of this study are as follows: (1) to conduct a direct numerical simulation of turbulent backward facing step flow using inflow and outflow conditions; and (2) to provide data in the form of Reynolds stress budgets for Reynolds averaged modeling. The report presents the basic statistical data and comparisons with the concurrent experiments of Jovic and Driver and budgets of turbulent kinetic energy.

Le, Hung↗

NASA NTRS1985

Direct numerical simulations of a reacting mixing layer with chemical heat release

In order to study the coupling between chemical heat release and fluid dynamics, direct numerical simulations of a chemically reacting mixing layer with heat release are performed. The fully compressible equations as well as an approximate set of equations that is asymptotically valid for low-Mach-number flows are treated. These latter equations have the computational advantage that high-frequency acoustic waves have been filtered out, allowing much larger time steps to be taken in the numerical solution procedure. A detailed derivation of these equations along with an outline of the numerical solution technique is given. Simulation results indicate that the rate of chemical product formed, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release.

Mcmurtry, P. A.↗

NASA NTRS1990

Direct numerical simulation of the transitional zone

Properties of the transitional zone in channel and boundary-layer flow are determined from the results of direct numerical simulations of forced transition. Three channel flow cases produce similar behavior, particularly in the 'transition function' relevant for some algebraic eddy viscosity models of the transitional zone. Moreover, the transitional zone from the computed forced boundary-layer transition is qualitatively similar to experimental data for natural transition and to the computed channel flow transition.

Zang, Thomas A.↗

NASA NTRS2004

Direct Numerical Simulation of Turbulent Condensation in Clouds

In this brief, we investigate the turbulent condensation of a population of droplets by means of a direct numerical simulation. To that end, a coupled Navier-Stokes/Lagrangian solver is used where each particle is tracked and its growth by water vapor condensation is monitored exactly. The main goals of the study are to find out whether turbulence broadens the droplet size distribution, as observed in in situ measurements. The second issue is to understand if and for how long a correlation between the droplet radius and the local supersaturation exists for the purpose of modeling sub-grid scale microphysics in cloud-resolving codes. This brief is organized as follows. In Section 2 the governing equations are presented, including the droplet condensation model. The implementation of the forcing procedure is described in Section 3. The simulation results are presented in Section 4 together with a sketch of a simple stochastic model for turbulent condensation. Conclusions and the main outcomes of the study are given in Section 5.

Shariff, K.↗

NASA NTRS1987

Analysis of subgrid-scale eddy viscosity with use of results from direct numerical simulations

Without resort to any modeling, subgrid-scale eddy viscosity is computed from the results of high-resolution (64 cubed and 128 cubed grid points) direct numerical simulations of three-dimensional homogeneous isotropic decaying turbulence. In these simulations the eddy viscosity peaks sharply at the cutoff wavenumber, in rough agreement with the results of Kraichnan (1976). In addition, in the low-wavenumber range the eddy viscosity may be negative, contrary to the generally accepted concept of subgrid-scale eddy viscosity. Some possible explanations of this behavior are discussed.

Domaradzki, J. Andrzej↗

NASA NTRS2017

Direct Numerical Simulation of Hypersonic Turbulent Boundary Layer inside an Axisymmetric Nozzle

As a first step toward a study of acoustic disturbance field within a conventional, hypersonic wind tunnel, direct numerical simulations (DNS) of a Mach 6 turbulent boundary layer on the inner wall of a straight axisymmetric nozzle are conducted and the results are compared with those for a flat plate. The DNS results for a nozzle radius to boundary-layer thickness ratio of 5:5 show that the turbulence statistics of the nozzle-wall boundary layer are nearly unaffected by the transverse curvature of the nozzle wall. Before the acoustic waves emanating from different parts of the nozzle surface can interfere with each other and undergo reflections from adjacent portions of the nozzle surface, the rms pressure fluctuation beyond the boundary layer edge increases toward the nozzle axis, apparently due to a focusing effect inside the axisymmetric configuration. Spectral analysis of pressure fluctuations at both the wall and the freestream indicates a similar distribution of energy content for both the nozzle and the flat plate, with the peak of the premultiplied frequency spectrum at a frequency of [(omega)(delta)]/U(sub infinity) approximately 6.0 inside the free stream and at [(omega)(delta)]/U(sub infinity) approximately 2.0 along the wall. The present results provide the basis for follow-on simulations involving reverberation effects inside the nozzle.

Huang, Junji↗

NASA NTRS1993

Direct numerical simulation of hot jets

The ultimate motivation of this work is to investigate the stability of two dimensional heated jets and its implications for aerodynamic sound generation from data obtained with direct numerical simulations (DNS). As pointed out in our last report, these flows undergo two types of instabilities, convective or absolute, depending on their temperature. We also described the limits of earlier experimental and theoretical studies and explained why a numerical investigation could give us new insight into the physics of these instabilities. The aeroacoustical interest of these flows was also underlined. In order to reach this goal, we first need to succeed in the DNS of heated jets. Our past efforts have been focused on this issue which encountered several difficulties. Our numerical difficulties are directly related to the physical problem we want to investigate since these absolutely or almost absolutely unstable flows are by definition very sensitive to the smallest disturbances and are very likely to reach nonlinear saturation through a numerical feedback mechanism. As a result, it is very difficult to compute a steady laminar solution using a spatial DNS. A steady state was reached only for strongly co-flowed jets, but these flows are almost equivalent to two independent mixing layers. Thus they are far from absolute instability and have much lower growth rates.

Jacob, Marc C.↗

NASA NTRS2020

Direct Numerical Simulation of Fog Formation Within the Surface Layer

We investigate the effects of shear-driven turbulence on fog formation in the stably stratified surface layer, using recent direct numerical simulations (DNS) of dry and moist open-channel flows. A constant cooling rate applied at the ground mimics longwave radiative cooling, leading to either turbulent (weakly stable) or laminar (very stable) flows. The transition of the Obukhov, Ozmidov and buoyancy length scales are investigated for different cooling rates and Reynolds numbers. Compared to the dry case, the condensation of liquid water in the moist case supports slightly higher cooling rates before leading to turbulence collapse. In the very stable moist cases, runaway cooling leads to substantial condensation of liquid water close to the ground and fog (visibility less than 1 km) results over much of the domain. In the weakly stable cases, turbulent mixing narrowly yields visibilities of 1 km close to the ground over a similar time period. Despite the idealised nature of the system, the present results suggest that turbulence impedes, although will not necessarily inhibit, fog formation.

Teixeira, Joao↗

NASA NTRS1991

Spatial Direct Numerical Simulation of Boundary-Layer Transition Mechanisms: Validation of PSE Theory

A study of instabilities in incompressible boundary-layer flow on a flat plate is conducted by spatial direct numerical simulation (DNS) of the Navier-Stokes equations. Here, the DNS results are used to critically evaluate the results obtained using parabolized stability equations (PSE) theory and to study mechanisms associated with breakdown from laminar to turbulent flow. Three test cases are considered: two-dimensional Tollmien-Schlichting wave propagation, subharmonic instability breakdown, and oblique-wave break-down. The instability modes predicted by PSE theory are in good quantitative agreement with the DNS results, except a small discrepancy is evident in the mean-flow distortion component of the 2-D test problem. This discrepancy is attributed to far-field boundary- condition differences. Both DNS and PSE theory results show several modal discrepancies when compared with the experiments of subharmonic breakdown. Computations that allow for a small adverse pressure gradient in the basic flow and a variation of the disturbance frequency result in better agreement with the experiments.

Joslin, R. D.↗

NASA NTRS1995

Applications of direct numerical simulation of turbulence in second order closures

This paper discusses two methods of developing models for the rapid pressure-strain correlation term in the Reynolds stress transport equation using direct numerical simulation (DNS) data. One is a perturbation about isotropic turbulence, the other is a perturbation about two-component turbulence -- an extremely anisotropic turbulence. A model based on the latter method is proposed and is found to be very promising when compared with DNS data and other models.

Shih, Tsan-Hsing↗

NASA NTRS1992

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗

NASA NTRS1991

Direct numerical simulation and analysis of shock turbulence interaction

Two kinds of linear analysis, rapid distortion theory (RDT) and linear interaction analysis (LIA), were used to investigate the effects of a shock wave on turbulence. Direct numerical simulations of two-dimensional isotropic turbulence interaction with a normal shock were also performed. The results from RDT and LIA are in good agreement for weak shock waves, where the effects of shock front curvature and shock front unsteadiness are not significant in producing vorticity. The linear analyses predict wavenumber-dependent amplification of the upstream one-dimensional energy spectrum, leading to turbulence scale length scale decrease through the interaction. Instantaneous vorticity fields show that vortical structures are enhanced while they are compressed in the shock normal direction. Entrophy amplfication through the shock wave compares favorably with the results of linear analyses.

Lee, Sangsan↗

NASA NTRS1987

A statistical investigation of the single-point pdf of velocity and vorticity based on direct numerical simulations

Vorticity plays a fundamental role in turbulent flows. The dynamics of vorticity in turbulent flows and the effect on single-point closure models were investigated. The approach was to use direct numerical simulations of turbulent flows to investigate the pdf of velocity and vorticity. The preliminary study of homogeneous shear flow has shown that the expectation of the fluctuating pressure gradient, conditioned with a velocity component, is linear in the velocity component, and that the coefficient is independent of velocity and vorticity. In addition, the work shows that the expectation of the pressure gradient, conditioned with a vorticity component, is essentially zero.

Mortazavi, M.↗

NASA NTRS2016

Direct Numerical Simulation of an Airfoil with Sand Grain Roughness on the Leading Edge

As part of a computational study of acoustic radiation due to the passage of turbulent boundary layer eddies over the trailing edge of an airfoil, the Lattice-Boltzmann method is used to perform direct numerical simulations of compressible, low Mach number flow past an NACA 0012 airfoil at zero degrees angle of attack. The chord Reynolds number of approximately 0.657 million models one of the test conditions from a previous experiment by Brooks, Pope, and Marcolini at NASA Langley Research Center. A unique feature of these simulations involves direct modeling of the sand grain roughness on the leading edge, which was used in the abovementioned experiment to trip the boundary layer to fully turbulent flow. This report documents the findings of preliminary, proof-of-concept simulations based on a narrow spanwise domain and a limited time interval. The inclusion of fully-resolved leading edge roughness in this simulation leads to significantly earlier transition than that in the absence of any roughness. The simulation data is used in conjunction with both the Ffowcs Williams-Hawkings acoustic analogy and a semi-analytical model by Roger and Moreau to predict the farfield noise. The encouraging agreement between the computed noise spectrum and that measured in the experiment indicates the potential payoff from a full-fledged numerical investigation based on the current approach. Analysis of the computed data is used to identify the required improvements to the preliminary simulations described herein.

Ribeiro, Andre F. P.↗

NASA NTRS1992

Direct numerical simulations of a spatially developing plane wake

In the present paper, direct numerical methods by which to simulate the spatially developing free shear flows in the transitional region are described and the numerical results of a spatially developing plane wake are presented. The incompressible time-dependent Navier-Stokes equations were solved using Pade finite difference approximations in the streamwise direction, a mapped pseudospectral Fourier method in the cross-stream direction, and a third-order compact Runge-Kutta scheme for time advancement. The unstable modes of the Orr-Sommerfeld equations were used to perturb the inlet of the wake. Statistical analyses were performed and some numerical results were compared with experimental measurements. When only the fundamental mode is forced, the energy spectra show amplification of the fundamental and its higher harmonics. In this case, unperturbed alternate vortices develop in the saturation region of the wake. The phase jitter around the fundamental frequency plays a critical role in generating vortices of random shape and spacing. Large- and small-scale distortions of the fundamental structure are observed. Pairing of vortices of the same sign is observed, as well as vortex coupling of vortices of the opposite sign.

Maekawa, Hiroshi↗

DOE OSTI2024

Direct numerical simulation of open-channel flow over a heterogeneous particle bed at low relative submergence

In this article, this study investigates turbulent open-channel flows over beds of irregularly arranged particles, using direct numerical simulations at a friction Reynolds number of Re T = 300⁠. Two distinct cases are examined: a polydisperse bed (P800) composed of multiple layers of randomly distributed spheres of varying sizes, and a monodisperse bed (M1015) formed by a random distribution of uniform sized spheres, with a bottommost single layer of varied-sized particles to introduce realistic randomness. Our investigation unveils a rich network of low- and high-speed streaks within the flow field, exhibiting distinctive behaviors in different bed configurations. The P800 case presents a poorly organized flow pattern induced by the varied particle sizes and arrangements, while the M1015 case shows a more regular flow pattern, marked by larger streaks. We also observe that total wall shear stress is substantially influenced by surface roughness-induced drag, extending beyond the effects documented in existing studies of open-channel flows. The present study reveals intricate secondary flow patterns over irregular particle beds. Large-scale circulations are discerned around particle crests in the P800 case and localized circulations with increased turbulence in the M1015 case. Furthermore, analysis of Reynolds stress tensor components indicates that roughness disrupts coherent turbulent eddies, consequently mitigating peak stress. We quantify correlations between drag force and local fluid velocity fluctuations. Notably, a larger deviation in drag is observed in the P800 case compared to M1015, accentuating the influence of particle size and distribution on fluid–particle interactions.

42 ENGINEERING↗