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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Consistency of Post-Newtonian Waveforms with Numerical Relativity

General relativity predicts the gravitational radiation signatures of mergers of compact binaries,such as coalescing binary black hole systems. Derivations of waveform predictions for such systems are required for optimal scientific analysis of observational gravitational wave data, and have so far been achieved primarily with the aid of the post-Newtonian (PN) approximation. The quaIity of this treatment is unclear, however, for the important late inspiral portion. We derive late-inspiral wave forms via a complementary approach, direct numerical simulation of Einstein's equations, which has recently matured sufficiently for such applications. We compare waveform phasing from simulations covering the last approximately 14 cycles of gravitational radiation from an equal-mass binary system of nonspinning black holes with corresponding 3PN and 3.5PN waveforms. We find phasing agreement consistent with internal error estimates based in either approach, at the level of one radian over approximately 10 cycles. The result suggests that PN waveforms for this system are effective roughly until the system reaches its last stable orbit just prior to the final merger.

Baker, John G.↗
NASA NTRS2006

Consistency of Post-Newtonian Waveforms with Numerical Relativity

General relativity predicts the gravitational radiation signatures of mergers of compact binaries, such as coalescing binary black hole systems. Derivations of waveform predictions for such systems are required for optimal scientific analysis of observational gravitational wave data, and have so far been achieved primarily with the aid of the post-Newtonian (PN) approximation. The quality of this treatment is unclear, however, for the important late inspiral portion. We derive late-inspiral waveforms via a complementary approach, direct numerical simulation of Einstein's equations, which has recently matured sufficiently for such applications. We compare waveform phasing from simulations covering the last approximately 14 cycles of gravitational radiation from an equal-mass binary system of nonspinning black holes with the corresponding 3PN and 3.5PN orbital phasing. We find agreement consistent with internal error estimates based on either approach at the level of one radian over approximately 10 cycles. The result suggests that PN waveforms for this system are effective roughly until the system reaches its last stable orbit just prior to the final merger/

Baker, John G.↗
NASA NTRS2007

Hairpin Vortices and the Final Stages of Transition

A spatially developing direct numerical simulation has been performed for flow over a flat plate that is subjected to a one-time fluid injection through an elongated slit in the wall. The flow parameters have been chosen to closely approximate the experimental conditions of Haidari, Taylow, and Smith (AIAA-89-0964). A hairpin vortex quickly develops near the upstream end of the slit, and a pair of necklace vortices form around the slow-moving injection fluid. As seen in the experiments and reported in Haidari and Smith (in review, JFM), the hairpin vortex spawns both in-line and sidelobe secondary vortices. However, no subdsidiary vortices (those formed by the inviscid deformation of a vortex-line bundle) are observed. At later times, a set of three different types of vortices are identified: hairpin vortex structures with heads that rise away from the wall horseshoe-shaped vortices that do not rise out of the boundary layer, and quasi-streamwise vortices. These structures interact with each other and with the wall layer to generate new vortices that are similar in structure to those mentioned above, although a particular parent vortex may have an offspring that more nearly resembles another member of the set. Perturbation velocity and vertical vorticity contours reveal an arrowhead shape of the highly disturbed region that is reminiscent of a turbulent spot. Spatially averaged velocity profiles in the highly disturbed area are nonlaminar, but as yet do not show typical low-of-the-wall behavior.

Singer, Bart A.↗
NASA NTRS2007

A Research Program for Improving Heat Transfer Prediction Capability for the Laminar to Turbulent Transition Region of Turbine Vanes/Blades

A program sponsored by the National Aeronautics and Space Administration (NASA) for the investigation of the heat transfer in the transition region of turbine vanes and blades with the object of improving the capability for predicting heat transfer is described,. The accurate prediction of gas-side heat transfer is important to the determination of turbine longevity, engine performance and developmental costs. The need for accurate predictions will become greater as the operating temperatures and stage loading levels of advanced turbine engines increase. The present methods for predicting transition shear stress and heat transfer on turbine blades are based on incomplete knowledge and are largely empirical. To meet the objectives of the NASA program, a team approach consisting of researchers from government, universities, a research institute, and a small business is presented. The research is divided into areas of experimentation, direct numerical simulation (DNS) and turbulence modeling. A summary of the results to date is given for the above research areas in a high-disturbance environment (bypass transition) with a discussion of the model development necessary for use in numerical codes.

Simon, Frederick F.↗
NASA NTRS2006

Validation of High-Speed Turbulent Boundary Layer and Shock-Boundary Layer Interaction Computations with the OVERFLOW Code

The capability of the OVERFLOW code to accurately compute high-speed turbulent boundary layers and turbulent shock-boundary layer interactions is being evaluated. Configurations being investigated include a Mach 2.87 flat plate to compare experimental velocity profiles and boundary layer growth, a Mach 6 flat plate to compare experimental surface heat transfer,a direct numerical simulation (DNS) at Mach 2.25 for turbulent quantities, and several Mach 3 compression ramps to compare computations of shock-boundary layer interactions to experimental laser doppler velocimetry (LDV) data and hot-wire data. The present paper describes outlines the study and presents preliminary results for two of the flat plate cases and two small-angle compression corner test cases.

Oliver, A. B.↗
NASA NTRS2008

Spatial and Temporal Low-Dimensional Models for Fluid Flow

A document discusses work that obtains a low-dimensional model that captures both temporal and spatial flow by constructing spatial and temporal four-mode models for two classic flow problems. The models are based on the proper orthogonal decomposition at two reference Reynolds numbers. Model predictions are made at an intermediate Reynolds number and compared with direct numerical simulation results at the new Reynolds number.

Kalb, Virginia↗
NASA NTRS2008

Testing of Error-Correcting Sparse Permutation Channel Codes

A computer program performs Monte Carlo direct numerical simulations for testing sparse permutation channel codes, which offer strong error-correction capabilities at high code rates and are considered especially suitable for storage of digital data in holographic and volume memories. A word in a code of this type is characterized by, among other things, a sparseness parameter (M) and a fixed number (K) of 1 or "on" bits in a channel block length of N.

Shcheglov, Kirill, V.↗
NASA NTRS2010

Development and Breakdown of Goertler Vortices in High Speed Boundary Layers

The nonlinear development of G rtler instability over a concave surface gives rise to a highly distorted stationary flow in the boundary layer that has strong velocity gradients in both spanwise and wall-normal directions. This distorted flow is susceptible to strong, high frequency secondary instability that leads to the onset of transition. For high Mach number flows, the boundary layer is also subject to the second mode instability. The nonlinear development of G rtler vortices and the ensuing growth and breakdown of secondary instability, the G rtler vortex interactions with second mode instabilities as well as oblique second mode interactions are examined in the context of both internal and external hypersonic configurations using nonlinear parabolized stability equations, 2-D eigenvalue analysis and direct numerical simulation. For G rtler vortex development inside the Purdue Mach 6 Ludwieg tube wind tunnel, multiple families of unstable secondary eigenmodes are identified and their linear and nonlinear evolution is examined. The computation of secondary instability is continued past the onset of transition to elucidate the physical mechanisms underlying the laminar breakdown process. Nonlinear breakdown scenarios associated with transition over a Mach 6 compression cone configuration are also explored.

Li, Fei↗
NASA NTRS2010

Algorithm for Stabilizing a POD-Based Dynamical System

This algorithm provides a new way to improve the accuracy and asymptotic behavior of a low-dimensional system based on the proper orthogonal decomposition (POD). Given a data set representing the evolution of a system of partial differential equations (PDEs), such as the Navier-Stokes equations for incompressible flow, one may obtain a low-dimensional model in the form of ordinary differential equations (ODEs) that should model the dynamics of the flow. Temporal sampling of the direct numerical simulation of the PDEs produces a spatial time series. The POD extracts the temporal and spatial eigenfunctions of this data set. Truncated to retain only the most energetic modes followed by Galerkin projection of these modes onto the PDEs obtains a dynamical system of ordinary differential equations for the time-dependent behavior of the flow. In practice, the steps leading to this system of ODEs entail numerically computing first-order derivatives of the mean data field and the eigenfunctions, and the computation of many inner products. This is far from a perfect process, and often results in the lack of long-term stability of the system and incorrect asymptotic behavior of the model. This algorithm describes a new stabilization method that utilizes the temporal eigenfunctions to derive correction terms for the coefficients of the dynamical system to significantly reduce these errors.

Kalb, Virginia L.↗
NASA NTRS2010

Parameterizing Coefficients of a POD-Based Dynamical System

A method of parameterizing the coefficients of a dynamical system based of a proper orthogonal decomposition (POD) representing the flow dynamics of a viscous fluid has been introduced. (A brief description of POD is presented in the immediately preceding article.) The present parameterization method is intended to enable construction of the dynamical system to accurately represent the temporal evolution of the flow dynamics over a range of Reynolds numbers. The need for this or a similar method arises as follows: A procedure that includes direct numerical simulation followed by POD, followed by Galerkin projection to a dynamical system has been proven to enable representation of flow dynamics by a low-dimensional model at the Reynolds number of the simulation. However, a more difficult task is to obtain models that are valid over a range of Reynolds numbers. Extrapolation of low-dimensional models by use of straightforward Reynolds-number-based parameter continuation has proven to be inadequate for successful prediction of flows. A key part of the problem of constructing a dynamical system to accurately represent the temporal evolution of the flow dynamics over a range of Reynolds numbers is the problem of understanding and providing for the variation of the coefficients of the dynamical system with the Reynolds number. Prior methods do not enable capture of temporal dynamics over ranges of Reynolds numbers in low-dimensional models, and are not even satisfactory when large numbers of modes are used. The basic idea of the present method is to solve the problem through a suitable parameterization of the coefficients of the dynamical system. The parameterization computations involve utilization of the transfer of kinetic energy between modes as a function of Reynolds number. The thus-parameterized dynamical system accurately predicts the flow dynamics and is applicable to a range of flow problems in the dynamical regime around the Hopf bifurcation. Parameter-continuation software can be used on the parameterized dynamical system to derive a bifurcation diagram that accurately predicts the temporal flow behavior.

Kalb, Virginia L.↗
NASA NTRS2009

Numerical and Physical Modeling of the Response of Resonator Liners to Intense Sound and Grazing Flow

Two significant advances have been made in the application of computational aeroacoustics methodology to acoustic liner technology. The first is that temperature effects for discrete sound are not the same as for broadband noise. For discrete sound, the normalized resistance appears to be insensitive to temperature except at high SPL. However, reactance is lower, significantly lower in absolute value, at high temperature. The second is the numerical investigation the acoustic performance of a liner by direct numerical simulation. Liner impedance is affected by the non-uniformity of the incident sound waves. This identifies the importance of pressure gradient. Preliminary design one and two-dimensional impedance models were developed to design sound absorbing liners in the presence of intense sound and grazing flow. The two-dimensional model offers the potential to empirically determine incident sound pressure face-plate distance from resonator orifices. This represents an important initial step in improving our understanding of how to effectively use the Dean Two-Microphone impedance measurement method.

Hersh, Alan S.↗
NASA NTRS2004

Part 1 of a Computational Study of a Drop-Laden Mixing Layer

This first of three reports on a computational study of a drop-laden temporal mixing layer presents the results of direct numerical simulations (DNS) of well-resolved flow fields and the derivation of the large-eddy simulation (LES) equations that would govern the larger scales of a turbulent flow field. The mixing layer consisted of two counterflowing gas streams, one of which was initially laden with evaporating liquid drops. The gas phase was composed of two perfect gas species, the carrier gas and the vapor emanating from the drops, and was computed in an Eulerian reference frame, whereas each drop was tracked individually in a Lagrangian manner. The flow perturbations that were initially imposed on the layer caused mixing and eventual transition to turbulence. The DNS database obtained included transitional states for layers with various liquid mass loadings. For the DNS, the gas-phase equations were the compressible Navier-Stokes equations for conservation of momentum and additional conservation equations for total energy and species mass. These equations included source terms representing the effect of the drops on the mass, momentum, and energy of the gas phase. From the DNS equations, the expression for the irreversible entropy production (dissipation) was derived and used to determine the dissipation due to the source terms. The LES equations were derived by spatially filtering the DNS set and the magnitudes of the terms were computed at transitional states, leading to a hierarchy of terms to guide simplification of the LES equations. It was concluded that effort should be devoted to the accurate modeling of both the subgridscale fluxes and the filtered source terms, which were the dominant unclosed terms appearing in the LES equations.

Okong'o, Nora A.↗
NASA NTRS2004

Helicity in Supercritical O2/H2 and C7H16/N2 Mixing Layers

This report describes a study of databases produced by direct numerical simulation of mixing layers developing between opposing flows of two fluids under supercritical conditions, the purpose of the study being to elucidate chemical-species-specific aspects of turbulence, with emphasis on helicity. The simulations were performed for two different fluid pairs -- O2/H2 and C7H16/N2 -- at similar values of reduced pressure.

Okongo, Nora↗
NASA NTRS2003

Irreversible Entropy Production in Two-Phase Mixing Layers

This report presents a study of dissipation (irreversible production of entropy) in three-dimensional, temporal mixing layers laden with evaporating liquid drops. The purpose of the study is to examine the effects of evaporating drops on the development of turbulent features in flows. Direct numerical simulations were performed to analyze transitional states of three mixing layers: one without drops, and two that included drops at different initial mass loadings. Without drops, the dissipation is essentially due to viscous effects. It was found that in the presence of drops, the largest contribution to dissipation was made by heating and evaporation of the drops, and that at large length scales, this contribution is positive (signifying that the drops reduce turbulence), while at small scales, this contribution is negative (the drops increase turbulence). The second largest contribution to dissipation was found to be associated with the chemical potential, which leads to an increase in turbulence at large scales and a decrease in turbulence at small scales. The next smaller contribution was found to be that of viscosity. The fact that viscosity effects are only third in order of magnitude in the dissipation is in sharp contrast to the situation for the mixing layer without the drops. The next smaller contribution - that of the drag and momentum of the vapor from the drops - was found to be negative at lower mass loading but to become positive at higher mass loading.

Okongo, Nora↗
NASA NTRS2003

Update on Area Production in Mixing of Supercritical Fluids

The focus of this research is on supercritical C7H16/N2 and O2/H2 mixing layers undergoing transitions to turbulence. The C7H16/N2 system serves as a simplified model of hydrocarbon/air systems in gas-turbine and diesel engines; the O2/H2 system is representative of liquid rocket engines. One goal of this research is to identify ways of controlling area production to increase disintegration of fluids and enhance combustion in such engines. As used in this research, "area production" signifies the fractional rate of change of surface area oriented perpendicular to the mass-fraction gradient of a mixing layer. In the study, a database of transitional states obtained from direct numerical simulations of the aforementioned mixing layers was analyzed to investigate global layer characteristics, phenomena in regions of high density-gradient magnitude (HDGM), irreversible entropy production and its relationship to the HDGM regions, and mechanisms leading to area production.

Okongo, Nora↗
NASA NTRS2003

Noncircular Cross Sections Could Enhance Mixing in Sprays

A computational study has shown that by injecting drops in jets of gas having square, elliptical, triangular, or other noncircular injection cross sections, it should be possible to increase (relative to comparable situations having circular cross section) the entrainment and dispersion of liquid drops. This finding has practical significance for a variety of applications in which it is desirable to increase dispersion of drops. For example, in chemical-process sprays, increased dispersion leads to increases in chemical- reaction rates; in diesel engines, increasing the dispersion of drops of sprayed fuel reduces the production of soot; and in household and paint sprays, increasing the dispersion of drops makes it possible to cover larger surfaces. It has been known for some years that single-phase fluid jets that enter flow fields through noncircular inlets entrain more fluid than do comparable jets entering through circular inlets. The computational study reported here was directed in part toward determining whether and how this superior mixing characteristic of noncircular single phase jets translates to a similar benefit in cases of two-phase jets (that is, sprays). The study involved direct numerical simulations of single- and two-phase free jets with circular, elliptical, rectangular, square, and triangular inlet cross sections. The two-phase jets consisted of gas laden with liquid drops randomly injected at the inlets. To address the more interesting case of evaporating drops, the carrier gas in the jets was specified to be initially unvitiated by the vapor of the liquid chemical species and the initial temperature of the drops was chosen to be smaller than that of the gas. The mathematical model used in the study was constructed from the conservation equations for the two-phase flow and included complete couplings of mass, momentum, and energy based on thermodynamically self-consistent specification of the enthalpy, internal energy, and latent heat of vaporization of the vapor.

Bellan, Josette↗
NASA NTRS2012

Effects of Riblets on Skin Friction in High-Speed Turbulent Boundary Layers

Direct numerical simulations of spatially developing turbulent boundary layers over riblets are conducted to examine the effects of riblets on skin friction at supersonic speeds. Zero-pressure gradient boundary layers with an adiabatic wall, a Mach number of M1 = 2.5, and a Reynolds number based on momentum thickness of Re = 1720 are considered. Simulations are conducted for boundary-layer flows over a clean surface and symmetric V- groove riblets with nominal spacings of 20 and 40 wall units. The DNS results confirm the few existing experimental observations and show that a drag reduction of approximately 7% is achieved for riblets with proper spacing. The influence of riblets on turbulence statistics is analyzed in detail with an emphasis on identifying the differences, if any, between the drag reduction mechanisms for incompressible and high-speed boundary layers.

Duan, Lian↗