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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗
NASA NTRS2022

Combined Bluntness and Roughness Effects on Cones at Hypersonic Speeds

This computational study investigates the effects of discrete roughness elements on a blunt cone at zero degrees angle of attack in a Mach 6 flow. Motivation was provided by experiments conducted in the Air Force Research Laboratory Mach 6 High Reynolds Number facility on a 7-degree half-angle cone with a roughness array located at 45 degrees from the apex on two nose tips of different blutness but equivalent roughness Reynolds number. Transition was only affected on the blunter cone, indicating that the transition onset is associated with the combined effects of bluntness and roughness. The present study investigates the 15.24 mm nose radius, 420 azimuthal wavenumber case via Navier-Stokes computations of the laminar base flow and instability analysis. Plane-marching parabolized stability equations (PSE) and inflow-resolvent analysis based on the three-dimensional, harmonic linearized Navier-Stokes equations (HLNSE) are used to calculate the amplification of disturbances along the roughness wake as well as over the roughness nearfield. Results show that the roughness shape can have a great impact on the characteristics of the most amplified wake instabilities. For the experimental configuration with cubic roughness elements of 15 μ m height, the flow is marginally unstable. For prismatic elements of 20 μ m height, the PSE predicts a logarithmic disturbance amplification ratio of N = 5.6 along its wake, but this ratio increases to N = 9.6 when the amplification over the roughness and separation regions is included in the inflow-resolvent analysis.

Boundary-layer transition↗
NASA NTRS1994

Numerical investigation of supersonic turbulent boundary layers with high wall temperature

A direct numerical approach has been developed to simulate supersonic turbulent boundary layers. The mean flow quantities are obtained by solving the parabolized Reynolds-averaged Navier-Stokes equations (globally). Fluctuating quantities are computed locally with a temporal direct numerical simulation approach, in which nonparallel effects of boundary layers are partially modeled. Preliminary numerical results obtained at the free-stream Mach numbers 3, 4.5, and 6 with hot-wall conditions are presented. Approximately 5 million grid points are used in all three cases. The numerical results indicate that compressibility effects on turbulent kinetic energy, in terms of dilatational dissipation and pressure-dilatation correlation, are small. Due to the hot-wall conditions the results show significant low Reynolds number effects and large streamwise streaks. Further simulations with a bigger computational box or a cold-wall condition are desirable.

Guo, Y.↗
NASA NTRS1986

Direct simulation of hypersonic flows over blunt slender bodies

Results of a numerical study of low-density hypersonic flow about cylindrically blunted wedges and spherically blunted cones with body half angles of 0, 5, and 10 deg are presented. Most of the transitional flow regime encountered during entry between the free molecule and continuum regimes is simulated for a reentry velocity of 7.5 km/s by including freestream conditions of 70 to 100 km. The bodies are at zero angle of incidence and have diffuse and finite catalytic surfaces. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method. The numerical simulations show that noncontinuum effects such as surface temperature jump, and velocity slip are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of 96 km for the two-dimensional configurations. Comparisons between the DSMC and continuum viscous shock-layer calculations highlight the significant difference in flowfield structure predicted by the two methods.

Moss, J. N.↗
NASA NTRS2011

Numerical Simulations of Two-Phase Reacting Flow in a Single-Element Lean Direct Injection (LDI) Combustor Using NCC

A series of numerical simulations of Jet-A spray reacting flow in a single-element lean direct injection (LDI) combustor have been conducted by using the National Combustion Code (NCC). The simulations have been carried out using the time filtered Navier-Stokes (TFNS) approach ranging from the steady Reynolds-averaged Navier-Stokes (RANS), unsteady RANS (URANS), to the dynamic flow structure simulation (DFS). The sub-grid model employed for turbulent mixing and combustion includes the well-mixed model, the linear eddy mixing (LEM) model, and the filtered mass density function (FDF/PDF) model. The starting condition of the injected liquid spray is specified via empirical droplet size correlation, and a five-species single-step global reduced mechanism is employed for fuel chemistry. All the calculations use the same grid whose resolution is of the RANS type. Comparisons of results from various models are presented.

Liu, Nan-Suey↗
NASA NTRS2007

On the Origin of Near-Radial Magnetic Fields in the Heliosphere: Numerical Simulations

Deviations from the direction of the "Parker spiral" can be found in in situ measurements of the interplanetary magnetic field on essentially all scales. One intriguing subset is the intervals of near-radial magnetic field, lasting for many hours. Some such intervals are obviously associated with coronal mass ejections, while others appear to be embedded within the ambient solar wind. Most occur on declining speed profiles, such that, when mapped back to the Sun, an entire radial field interval appears to have been launched at approximately the same time. It has been proposed that these events are the result of abrupt, semipermanent speed decreases on these field lines close to the Sun, and that such speed changes might be due to interchange reconnection. In this study, we use a three-dimensional, time-dependent magnetohydrodynamic model to assess to what extent this can account for near-radial magnetic fields observed relatively far out in the heliosphere. We find that sudden speed drops on the trailing portions of high-speed flows can produce strongly underwound (that is near radial) field lines in the heliosphere, although significantly larger speed gradients are required than are typically observed. Moreover, the simulations also reproduce the decreases in density, temperature, and magnetic field strength that are also commonly observed within these events. The question of what produces the abrupt speed drops remains to be answered.

Mikic, Zoran↗
NASA NTRS1988

Numerical investigations in three-dimensional internal flows

An investigation into the use of computational fluid dynamics (CFD) was performed to examine the expected heat transfer rates that will occur within the NASA-Ames 100 megawatt arc heater nozzle. This nozzle was tentatively designed and identified to provide research for a directly connected combustion experiment specifically related to the National Aerospace Plane Program (NASP) aircraft, and is expected to simulate the flow field entering the combustor section. It was found that extremely fine grids, that is very small mesh spacing near the wall, are required to accurately model the heat transfer process and, in fact, must contain a point within the laminar sublayer if results are to be taken directly from a numerical simulation code. In the present study, an alternative to this very fine mesh and its attendant increase in computational time was invoked and is based on a wall-function method. It was shown that solutions could be obtained that give accurate indications of surface heat transfer rate throughout the nozzle in approximately 1/100 of the computer time required to do the simulation directly without the use of the wall-function implementation. Finally, a maximum heating value in the throat region of the proposed slit nozzle for the 100 megawatt arc heater was shown to be approximately 6 MW per square meter.

Rose, William C.↗
NASA NTRS1989

Advances in fluid dynamics

Various papers on advances in fluid dynamics are presented. Individual topics addressed include: nonunique solutions of the Euler equations; numerical solution of the biharmonic equation using a spectral multigrid method; finite-volume shape and placement for efficiency, accuracy, and physics; NND schemes and their application to numerical simulation of complex plume flow problems; symbolic analysis and computational algorithm construction; and velocity and mass flux distribution measurements of spherical glass beads in air flow in a 90-deg vertical-to-horizontal bend. Also discussed are: hyperbolic solution adaptive mesh generation; progress and future directions in CFD; numerical simulation of shock focusing; microbubble formation and splitting in a turbulent boundary layer for turbulence reduction; numerical study of spherical Taylor-Couette flow, trajectories of vortex lines beneath a free surface or above a plane; and a mathematical model for response of the coronary circulation to high sustained gravitational force fields.

Ballhaus, W. F., Jr.↗
NASA NTRS1998

Study of Interesting Solidification Phenomena on the Ground and in Space (MEPHISTO)

Real-time Seebeck voltage variations in a Sn-Bi melt during directional solidification in the MEPHISTO spaceflight experiment flown on the USMP-3 mission, can be correlated with well characterized thruster firings and an Orbiter Main System (OMS) burn. The Seebeck voltage measurement is related to the response of the instantaneous average melt composition at the melt-crystal interface. This allowed us to make a direct comparison of numerical simulations with the experimentally obtained Seebeck signals. Based on the results of preflight and real-time computations, several well-defined thruster firing events were programmed to occur at specific times during the experiment. In particular, we simulated the effects of the thruster firings on melt and crystal composition in a directionally solidifying Sn-Bi alloy. The relative accelerations produced by the firings were simulated by impulsive accelerations of the same magnitude, duration and orientation as the requested firings. A comparison of the simulation results with the Seebeck signal indicates that there is a good agreement between the two. This unique opportunity allows us, for the first time, to quantitatively characterize actual g-jitter effects on an actual crystal growth experiment and to properly calibrate our models of g-jitter effects on crystal growth.

Favier, J.-J.↗
NASA NTRS1999

Study of Interesting Solidification Phenomena on the Ground and in Space (MEPHISTO)

Real-time Seebeck voltage variations in a Sn-Bi melt during directional solidification in the MEPHISTO spaceflight experiment flown on the USMP-3 mission, have been correlated with well-characterized thruster firings and an Orbiter Main System (OMS) burn. The Seebeck voltage measurement is related to the response of the instantaneous average melt composition at the melt-crystal interface. This allowed us to make a direct comparison of numerical simulations with the experimentally obtained Seebeck signals. Based on the results of preflight and real-time computations, several well-defined thruster firing events were programmed to occur at specific times during the experiment. In particular, we simulated the effects of the thruster firings on melt and crystal composition in a directionally solidifying Sn-Bi alloy. The relative accelerations produced by the firings were simulated by impulsive accelerations of the same magnitude, duration and orientation as the requested firings. A comparison of the simulation results with the Seebeck signal indicates that there is a good agreement between the two. This unique opportunity allows us to make the first quantitative characterization of actual g-jitter effects on an actual crystal growth experiment and to calibrate our models of g-jitter effects on crystal growth.

Alexander, J. Iwan D.↗
NASA NTRS1996

Orbiter Aerodynamic Acceleration Flight Measurements in the Rarefied-Flow Transition Regime

Acceleration data taken from the Orbital Acceleration Research Experiment (OARE) during reentry on STS-62 have been analyzed using calibration factors taken on orbit. This is the first Orbiter mission which collected OARE data during the Orbiter reentry phase. The data examined include the flight regime from orbital altitudes down to about 90 km which covers the free-molecule-flow regime and the upper altitude fringes of the rarefied-flow transition into the hypersonic continuum. Ancillary flight data on Orbiter position, orientation, velocity, and rotation rates have been used in models to transform the measured accelerations to the Orbiter center-of-gravity, from which aerodynamic accelerations along the Orbiter body axes have been calculated. Residual offsets introduced in the measurements by unmodeled Orbiter forces are identified and discussed. Direct comparisons are made between the OARE flight data and an independent micro-gravity accelerometer experiment, the High Resolution Accelerometer Package (HiRAP), which also obtained flight data on reentry during the mission down to about 95 km. The resulting OARE aerodynamic acceleration measurements along the Orbiter's body axis, aid the normal to axial acceleration ratio in the free-molecule-flow and transition-flow regimes are presented and compared with numerical simulations from three direct simulation Monte Carlo codes.

Blanchard, Robert C.↗
DOE OSTI2023

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Multi-modal Dataset of a Polycrystalline Metallic Material: 3D Microstructure and Deformation Fields

The development of high-fidelity mechanical property prediction models for the design of polycrystalline materials relies on large volumes of microstructural feature data. Concurrently, at these same scales, the deformation fields that develop during mechanical loading can be highly heterogeneous. Spatially correlated measurements of 3D microstructure and the ensuing deformation fields at the micro-scale would provide highly valuable insight into the relationship between microstructure and macroscopic mechanical response. They would also provide direct validation for numerical simulations that can guide and speed up the design of new materials and microstructures. However, to date, such data have been rare. Here, a one-of-a-kind, multi-modal dataset is presented that combines recent state-of-the-art experimental developments in 3D tomography and high-resolution deformation field measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data-Driven RANS Turbulence Closures for Forced Convection Flow in Reactor Downcomer Geometry

Recent progress in data-driven turbulence modeling has shown its potential to enhance or replace traditional equation-based Reynolds-averaged Navier-Stokes (RANS) turbulence models. Here, this work utilizes invariant neural network (NN) architectures to model Reynolds stresses and turbulent heat fluxes in forced convection flows (when the models can be decoupled). As the considered flow is statistically one dimensional, the invariant NN architecture for the Reynolds stress model reduces to the linear eddy viscosity model. To develop the data-driven models, direct numerical and RANS simulations in vertical planar channel geometry mimicking a part of the reactor downcomer are performed. Different conditions and fluids relevant to advanced reactors (sodium, lead, unitary-Prandtl-number fluid, and molten salt) constitute the training database. The models enabled accurate predictions of velocity and temperature, and compared to the baseline k–τ turbulence model with the simple gradient diffusion hypothesis, do not require tuning of the turbulent Prandtl number. The data-driven framework is implemented in the open-source graphics processing unit–accelerated spectral element solver nekRS and has shown the potential for future developments and consideration of more complex mixed convection flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

RKKY coupled local-moment magnetism in NaFe 1 – x Cu x As

A central question in a large class of strongly correlated electron systems, including heavy fermion compounds and iron pnictides, is the identification of different phases and their origins. It has been shown that the antiferromagnetic (AFM) phase in some heavy fermion compounds is induced by Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between localized moments, and that the competition between this interaction and Kondo effect is responsible for quantum criticality. However, conclusive experimental evidence of the RKKY interaction in pnictides is lacking. Here, using high resolution 23 Na NMR measurements on lightly Cu-doped metallic single crystals of NaFe 1–x Cu x As (x ≈ 0.01) and numerical simulation, we show direct evidence of the RKKY interaction in this pnictide system. Aided by computer simulation, we identify the 23 Na NMR satellite resonances with the RKKY oscillations of spin polarization at Fe sites. Our 23 Na spin-lattice and spin-spin relaxation data exhibit a signature of an itinerant and inhomogeneous AFM phase in this system, accompanied by a simulation of Cu-induced perturbation to the ordered moments on the Fe sites. Our NMR results indicate coexistence of local and itinerant magnetism in lightly Cu-doped NaFe 1–x Cux As.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A numerical approach for detecting switch-like bistability in mass action chemical reaction networks with conservation laws

Theoretical analysis of signaling pathways can provide a substantial amount of insight into their function. One particular area of research considers signaling pathways capable of assuming two or more stable states given the same amount of signaling ligand. This phenomenon of bistability can give rise to switch-like behavior, a mechanism that governs cellular decision making. Investigation of whether or not a signaling pathway can confer bistability and switch-like behavior, without knowledge of specific kinetic rate constant values, is a mathematically challenging problem. Recently a technique based on optimization has been introduced, which is capable of finding example parameter values that confer switch-like behavior for a given pathway. Although this approach has made it possible to analyze moderately sized pathways, it is limited to reaction networks that presume a uniterminal structure. It is this limited structure we address by developing a general technique that applies to any mass action reaction network with conservation laws. In this paper we developed a generalized method for detecting switch-like bistable behavior in any mass action reaction network with conservation laws. The method involves (1) construction of a constrained optimization problem using the determinant of the Jacobian of the underlying rate equations, (2) minimization of the objective function to search for conditions resulting in a zero eigenvalue, (3) computation of a confidence level that describes if the global minimum has been found and (4) evaluation of optimization values, using either numerical continuation or directly simulating the ODE system, to verify that a bistability region exists. The generalized method has been tested on three motifs known to be capable of bistability. We have developed a variation of an optimization-based method for the discovery of bistability, which is not limited to uniterminal chemical reaction networks. Successful completion of the method provides an S-shaped bifurcation diagram, which indicates that the network acts as a bistable switch for the given optimization parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Investigation of Automated Corrective Actions for Special Protection Schemes

The constantly evolving nature of the grid is compelling the design process of Remedial Action Schemes (RAS) to keep up with the changes. This document proposes a flexible and computationally efficient approach to automatically determine RAS corrective actions that alleviate line overloading violations. Statistical and functional characteristics summarized from RAS implemented in real power systems are used to guide the design parameters. This report also leverages sensitivity-based techniques to determine corrective actions for specific contingencies quickly without repeated numerical simulations. Finally, future directions for implementing this approach for a fully automated, online RAS are discussed.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1979

Retarding potential analyzer for the Pioneer-Venus Orbiter Mission

The retarding potential analyzer on the Pioneer-Venus Orbiter Mission has been designed to measure most of the thermal plasma parameters within and near the Venusian ionosphere. Parameters include total ion concentration, concentrations of the more abundant ions, ion temperatures, ion drift velocity, electron temperature, and low-energy (0-50 eV) electron distribution function. To accomplish these measurements on a spinning vehicle with a small telemetry bit rate, several functions, including decision functions not previously used in RPA's, have been developed and incorporated into this instrument. The more significant functions include automatic electrometer ranging with background current compensation; digital, quadratic retarding potential step generation for the ion and low-energy electron scans; a current sampling interval of 2 ms throughout all scans; digital logic inflection point detection and data selection; and automatic ram direction detection. Extensive numerical simulation and plasma chamber tests have been conducted to verify adequacy of the design for the Pioneer Mission.

Knudsen, W. C.↗