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

Collisional evolution of mass-distribution spectrum of planetesimals

The early growth of planetesimals by mutual direct collisions is numerically simulated with a Monte Carlo technique to show how planetesimals with specific mass and velocity distributions evolve into a full-size planet. Four types of collisions are taken into account: rebound, erosion, catastrophic break-up, and coagulation. It is shown that evolution of mass-distribution spectrum is expressed by an inverse power relation. It is also suggested that whether or not a planetesimal can survive catastrophic collision is primarily dependent on mean relative velocity and mechanical properties of planetesimals. It is necessary for the early growth of rocky (basaltic) materials that mean relative velocity be much smaller than 0.1 km/s. Otherwise it is necessary to introduce something like nucleating agents (such as iron bodies which have plastic properties at temperatures higher than 200 K) for the formation of terrestrial planets.

Matsui, T.↗
NASA NTRS1991

Numerical investigation on Benard instability in a finite liquid layer

A numerical procedure for directly simulating the Benard-Marangoni instabilities (B-M-I) in a bounded liquid layer is presented in this paper. The procedure consists of applying a finite amplitude disturbance to the basic static state, and then integrating the Navier-Stokes equations to determine whether the disturbance will die down or will reach a state of finite-strength steady convection. The critical Marangoni number (Mac) for the onset of B-M-I can thus be determined and can be correlated as a function of the aspect ratio (Ar), Prandtl number (Pr), and Rayleigh number (Ra). The Biot number (B) between the liquid and the air is analyzed to approximate the heat transfer condition along the free surface. A 2D calculation is performed to investigate the effect of various initial disturbances, and the Mac is determined for Ar = 2, Ra = 0, and Pr = 0.7. Current results show that disturbances of different nature and amplitude have little effect on Mac. The Mac determined in this study also clearly demonstrates the dominant effect of the sidewalls.

Duh, J. C.↗
NASA NTRS2003

Numerical Simulations of Wake/Boundary Layer Interactions

Direct and large-eddy simulations of the interaction between the wake of a circular cylinder and a flat-plate boundary layer are conducted. Two Reynolds numbers are examined. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re(sub theta) achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established.

Piomelli, Ugo↗