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 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1995

A numerical study of fundamental shock noise mechanisms

The results of this thesis demonstrate that direct numerical simulation can predict sound generation in unsteady aerodynamic flows containing shock waves. Shock waves can be significant sources of sound in high speed jet flows, on helicopter blades, and in supersonic combustion inlets. Direct computation of sound permits the prediction of noise levels in the preliminary design stage and can be used as a tool to focus experimental studies, thereby reducing cost and increasing the probability of a successfully quiet product in less time. This thesis reveals and investigates two mechanisms fundamental to sound generation by shocked flows: shock motion and shock deformation. Shock motion is modeled by the interaction of a sound wave with a shock. During the interaction, the shock wave begins to move and the sound pressure is amplified as the wave passes through the shock. The numerical approach presented in this thesis is validated by the comparison of results obtained in a quasi-one dimensional simulation with linear theory. Analysis of the perturbation energy demonstrated for the first time that acoustic energy is generated by the interaction. Shock deformation is investigated by the numerical simulation of a ring vortex interacting with a shock. This interaction models the passage of turbulent structures through the shock wave. The simulation demonstrates that both acoustic waves and contact surfaces are generated downstream during the interaction. Analysis demonstrates that the acoustic wave spreads cylindrically, that the sound intensity is highly directional, and that the sound pressure level increases significantly with increasing shock strength. The effect of shock strength on sound pressure level is consistent with experimental observations of shock noise, indicating that the interaction of a ring vortex with a shock wave correctly models a dominant mechanism of shock noise generation.

Meadows, Kristine R.↗

NASA NTRS1999

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

Bui, Trong T.↗

NASA NTRS1991

Simulation of flame-turbulence interactions in premixed combustion

Turbulent combustion is a challenging problem of fundamental interest with numerous practical applications. Direct Numerical Simulation is an attractive tool that provides detailed information on flame-flow interactions and, thereby, provides an ideal basis for testing and developing turbulent combustion models. The work being undertaken to investigate the laminar flamelet concepts as applied to premixed combustion is described. More specifically, the influence of non-unity Lewis number with emphasis on the effects of hydrodynamic straining and flame front curvature is examined.

Trouve, Arnaud↗

NASA NTRS2003

Characterization of High-Frequency Excitation of a Wake by Simulation

Insights into the effects of high-frequency forcing on free shear layer evolution are gained through analysis of several direct numerical simulations. High-frequency forcing of a fully turbulent plane wake results in only a weak transient effect. On the other hand, significant changes in the developed turbulent state may result when high-frequency forcing is applied to a transitional wake. The impacts of varying the characteristics of the high-frequency forcing are examined, particularly, the streamwise wavenumber band in which forcing is applied and the initial amplitude of the forcing. The high-frequency excitation is found to increase the dissipation rate of turbulent kinetic energy, to reduce the turbulent kinetic energy production rate, and to reduce the turbulent kinetic energy suppression increases with forcing amplitude once a threshold level has been reached. For a given initial forcing energy, the largest reduction in turbulent kinetic energy density was achieved by forcing wavenumbers that are about two to three times the neutral wavenumber determined from linear stability theory.

Cain, Alan B.↗

NASA NTRS1989

Modeling of high speed chemically reacting flow-fields

The SPARK3D and SPARK3D-PNS computer programs were developed to model 3-D supersonic, chemically reacting flow-fields. The SPARK3D code is a full Navier-Stokes solver, and is suitable for use in scramjet combustors and other regions where recirculation may be present. The SPARK3D-PNS is a parabolized Navier-Stokes solver and provides an efficient means of calculating steady-state combustor far-fields and nozzles. Each code has a generalized chemistry package, making modeling of any chemically reacting flow possible. Research activities by the Langley group range from addressing fundamental theoretical issues to simulating problems of practical importance. Algorithmic development includes work on higher order and upwind spatial difference schemes. Direct numerical simulations employ these algorithms to address the fundamental issues of flow stability and transition, and the chemical reaction of supersonic mixing layers and jets. It is believed that this work will lend greater insight into phenomenological model development for simulating supersonic chemically reacting flows in practical combustors. Currently, the SPARK3D and SPARK3D-PNS codes are used to study problems of engineering interest, including various injector designs and 3-D combustor-nozzle configurations. Examples, which demonstrate the capabilities of each code are presented.

Drummond, J. P.↗

NASA NTRS1997

Time-Domain Filtering for Spatial Large-Eddy Simulation

An approach to large-eddy simulation (LES) is developed whose subgrid-scale model incorporates filtering in the time domain, in contrast to conventional approaches, which exploit spatial filtering. The method is demonstrated in the simulation of a heated, compressible, axisymmetric jet, and results are compared with those obtained from fully resolved direct numerical simulation. The present approach was, in fact, motivated by the jet-flow problem and the desire to manipulate the flow by localized (point) sources for the purposes of noise suppression. Time-domain filtering appears to be more consistent with the modeling of point sources; moreover, time-domain filtering may resolve some fundamental inconsistencies associated with conventional space-filtered LES approaches.

Pruett, C. David↗

NASA NTRS2021

High-Fidelity Simulation of Turbulent Flow Past a Gaussian Bump

A spanwise-periodic computation of a turbulent flow past a Gaussian bump is performed in the form of a hybrid direct numerical simulation and wall-resolved large-eddy simulation. A fourth-order spatially-accurate flow solver is employed to perform the simulation, using 10.2 billion grid points for a Reynolds number of 170000 based on the bump height. The key findings from the simulation are reported in the acceleration and deceleration flow regions associated with the bump shape. Significant anisotropy in the normal Reynolds stresses, along both the wall-normal and streamwise directions, is observed within the acceleration region. The ratio between the Reynolds shear stress and turbulent kinetic energy in that region also experiences significant deviations from the norms of a zero pressure gradient turbulent boundary layer. The chosen Reynolds number generates strong flow separation in the adverse pressure gradient region, which is in contrast with a previous simulation at half the Reynolds number that only indicated incipient separation. An internal layer generated in the acceleration region evolves into a free shear layer that develops in the deceleration region and separates. Proper modeling of this inner layer appears crucial to predict the flow separation. Surface curvature effects on the attached flow development are also discussed.

Turbulent Boundary Layer↗

NASA NTRS1997

Large Scale Turbulent Structures in Supersonic Jets

Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations (DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those of a spatially evolving jet, a temporal jet problem was solved, using periodicity in the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible appli(,a- tion to active noise suppression. In addition, the data generated can be used to compute, various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved DNS of a turbulent jet flow is produced.

Rao, Ram Mohan↗

NASA NTRS1987

Toward the large-eddy simulations of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model (in terms of Favre-filtered fields) is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The three dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 to the third power grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. Future applications of these compressible subgrid-scale models to the large-eddy simulation of supersonic aerodynamic flows are discussed briefly.

Erlebacher, G.↗

NASA NTRS1992

Features of wavy vortices in a curved channel from experimental and numerical studies

Results are reported from an experimental study obtaining evidence of time-dependent, wavy vortex motions associated with undulating and twisting Dean vortices in a curved channel with 40-to-1 aspect ratio, and mild curvature (radius ratio = 0.979). The results are compared with direct numerical simulations of time-dependent 3D Navier-Stokes equations using boundary conditions in the spanwise and streamwise directions. When viewed in cross section, experimental visualizations of undulating and twisting vortex flows show rocking motion and changes in the direction of the flow between vortices that are like those observed in the simulations. Experimental spectra show that undulating vortices are replaced by the higher-frequency, shorter streamwise wavelength twisting vortices at higher Reynolds numbers. When undulating vortices are present, experimental power spectra and visualizations give frequencies that are somewhat lower than the most unstable frequencies predicted by linear stability analysis.

Ligrani, P. M.↗

NASA NTRS1997

Annual Research Briefs

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Spinks, Debra↗

NASA NTRS1996

Near-Wall Models in Large Eddy Simulations of Flow Behind a Backward-Facing Step

Accurate large eddy simulation (LES) of a wall-bounded flow generally requires a near-wall resolution comparable to that in direct numerical simulation (DNS). As much as 50% of the total grid points and computational costs are expended in the near-wall regions in a typical simulation. This limits LES to fairly low Reynolds numbers on current computers. To perform practical flow applications at realistically high Reynolds numbers, such as flow over an airfoil, it is desirable to replace very thin, near-wall regions in the LES with easily and inexpensively computed wall models to specify the near-wall boundary conditions.

Cabot, W.↗

NASA NTRS1990

One-equation near-wall turbulence modeling with the aid of direct simulation data

The length scales appearing in the relations for the eddy viscosity and dissipation rate in one-equation models were evaluated from direct numerical simulation data for developed channel and boundary-layer flow at two Reynolds numbers each. To prepare the ground for the evaluation, the distribution of the most relevant mean-flow and turbulence quantities is presented and discussed with respect to Reynolds-number influence and to differences between channel and boundary-layer flow. An alternative model is also examined in which bar-(v'(exp 2))(exp 1/2) is used as velocity scale instead of k(exp 1/2). With this velocity scale, the length scales now appearing in the model follow very closely a linear relationship near the wall so that no damping is necessary. For the determination of bar-v'(exp 2) in the context of a one-equation model, a correlation is provided between bar-(v'(exp 2))/k and bar-(u'v')/k.

Rodi, W.↗

NASA NTRS2019

Effect of Spatial Filtering in Implicit Large-Eddy Simulations of Separated Flows

The relatively high Reynolds number of turbulent flows encountered in various applications puts these problems well beyond the reach of direct numerical simulation (DNS) at present. Meanwhile, lower-fidelity Reynolds-averaged Navier-Stokes (RANS) calculations are known to be not accurate enough in complex problems, such as smooth-body flow separation and other flows involving highly-unsteady phenomena. Hence, given the current infeasibility of DNS and the unsatisfactory performance of RANS, intermediate techniques such as large-eddy simulation (LES) and hybrid RANS-LES, whose fidelity lie between RANS and DNS, have received much attention for application to various problems of practical importance. Modeling of the effect of missing scales on resolved scales, also known as subgrid-scale (SGS) modeling, is an important subject for LES. SGS models can be broadly categorized as explicit or implicit approaches. The explicit approach is based on an SGS model that explicitly appears in the governing equations expressed in the form of so-called “filtered Navier-Stokes equations”, which describe the evolution of the turbulence scales resolved by the LES grid. The effect of the scales unresolved by the grid is represented by the SGS model. The implicit modeling approach, on the other hand, does not employ an explicit model but instead treats the intrinsic dissipation of the numerical discretization scheme as an implicit SGS model. An LES without an explicit SGS model is commonly termed as an implicit LES (ILES). The relative merits of one SGS modeling approach over another is a subject of ongoing debate. We have opted to employ an ILES methodology, based on high-order compact finite-difference and spatial filtering schemes, in our recent investigations of separated flow problems [1, 2]. Further discussion of our preference of ILES over explicit LES is provided in Uzun and Malik [2]. The spatial filtering operation, described in the next section, is treated as an implicit SGS model for the ILES. Some observations made during the course of our recent investigations, which pointed out to excessive numerical dissipation in certain parts of the flowfield, prompted us to take a closer look at the potential effect of the spatial filter on ILES predictions. This technical note is therefore devoted to spatial filter effects in the context of a high Reynolds number, transonic shock-induced separated flow.

Ali Uzun↗

NASA NTRS2004

Direct Simulation of Evolution and Control of Nonlinear Instabilities in Attachment-Line Boundary Layers

The unsteady, incompressible Navier-Stokes equations are used for the direct numerical simulation (DNS) of spatially evolving disturbances in a three-dimensional (3-D) attachment-line boundary layer. Two-dimensional (2-D) disturbances are introduced either by forcing at the in ow or by harmonic-source generators at the wall; 3-D disturbances are introduced by harmonic-source generators at the wall. The DNS results are in good agreement with both 2-D non-parallel theory (for small-amplitude disturbances) and weakly nonlinear theory (for finite-amplitude disturbances), which validates the two theories. The 2-D DNS results indicate that nonlinear disturbance growth occurs near branch II of the neutral stability curve; however, steady suction can be used to stabilize this disturbance growth. For 3-D instabilities that are generated o the attachment line, spreading both toward and away from the attachment line causes energy transfer to the attachment-line and downstream instabilities; suction stabilizes these instabilities. Furthermore, 3-D instabilities are more stable than 2-D or quasi-2-D instabilities.

Joslin, Ronald D.↗

NASA NTRS2010

Microscopic Scale Simulation of the Ablation of Fibrous Materials

Ablation by oxidation of carbon-fiber preforms impregnated in carbonized phenolic matrix is modeled at microscopic scale. Direct numerical simulations show that the carbonized phenolic matrix ablates in volume leaving the carbon fibers exposed. This is due to the fact that the reactivity of carbonized phenolic is higher than the reactivity of carbon fibers. After the matrix is depleted, the fibers ablate showing progressive reduction of their diameter. The overall material recession occurs when the fibers are consumed. Two materials with the same carbon-fiber preform, density and chemical composition, but with different matrix distributions are studied. These studies show that at moderate temperatures (< 1000 K) the microstructure of the material influences its recession rate; a fact that is not captured by current models that are based on chemical composition only. Surprisingly, the response of these impregnated-fiber materials is weakly dependent on the microstructure at very high temperatures (e.g., Stardust peak heating conditions: 3360K).

Lachaud, Jean Romain↗

NASA NTRS1990

Large-eddy simulation of transitional channel flow

A large-eddy simulation (LES) of transition in plane channel flow was carried out. The LES results were compared with those of a fine direct numerical simulation (DNS), and with those of a coarse DNS that uses the same mesh as the LES, but does not use a residual stress model. While at the early stages of transition, LES and coarse DNS give the same results: the presence of the residual stress model was found to be necessary to predict accurately mean velocity and Reynolds stress profiles during the late stages of transition (after the second spike stage). The evolution of single Fourier modes is also predicted more accurately by the LES than by the DNS. As small scales are generated, the dissipative character of the residual stress starts to reproduce correctly the energy cascade. As transition progresses, the flow approaches its fully developed turbulent state, the subgrid scales tend towards equilibrium, and the model becomes more accurate.

Piomelli, Ugo↗

NASA NTRS1994

Nonlinear evolution of a large-amplitude circularly polarized Alfven wave: Low beta

The nature of turbulent cascades arising from the parametric instabilities of a monochromatic field-aligned large-amplitude circularly polarized Alfven wave is investigated via direct numerical simulation for the case of low plasma Beta and no wave dispersion. The magnetohydrodynamic code permits nonlinear couplings in the parallel direction to the ambient magnetic field and one perpendicular direction. Compressibility is included in the form of a polytropic equation of state. Anisotropic turbulent cascades, similar to those found in early incompressible two-dimensional simulations, occur after nonlinear saturation of the parallel propagating decay instability. The turbulent spectrum can be divided into three regimes: the lowest wave numbers are dominated by lower sideband remnants of the parametric process, intermediate wave numbers display nearly incompressible dynamics, and the highest wave numbers are dominated by acoustic turbulence.

Ghosh, S.↗