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,549 records · Page 86

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

Stochastic weighted particle methods for population balance equations with coagulation, fragmentation and spatial inhomogeneity

Graphical abstract: -- Highlights: •Problems concerning multi-compartment population balance equations are studied. •A class of fragmentation weight transfer functions is presented. •Three stochastic weighted algorithms are compared against the direct simulation algorithm. •The numerical errors of the stochastic solutions are assessed as a function of fragmentation rate. •The algorithms are applied to a multi-dimensional granulation model. -- Abstract: This paper introduces stochastic weighted particle algorithms for the solution of multi-compartment population balance equations. In particular, it presents a class of fragmentation weight transfer functions which are constructed such that the number of computational particles stays constant during fragmentation events. The weight transfer functions are constructed based on systems of weighted computational particles and each of it leads to a stochastic particle algorithm for the numerical treatment of population balance equations. Besides fragmentation, the algorithms also consider physical processes such as coagulation and the exchange of mass with the surroundings. The numerical properties of the algorithms are compared to the direct simulation algorithm and an existing method for the fragmentation of weighted particles. It is found that the new algorithms show better numerical performance over the two existing methods especially for systems with significant amount of large particles and high fragmentation rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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

Displacement Threshold Energy and Recovery in an Al-Ti Nanolayered System with Intrinsic Point Defect Partitioning

A method is established and validated using molecular dynamics (MD) to determine the displacement threshold energies as Ed in nanolayered, multilayered systems of dissimilar metals. The method is applied to specifically oriented nanolayered films of Al-Ti where the crystal structure and interface orientations are varied in atomic models and Ed is calculated. Methods for defect detection are developed and discussed based on prior research in the literature and based on specific crystallographic directions available in the nanolayered systems. These are compared and contrasted to similar calculations in corresponding bulk materials, including fcc Al, fcc Ti, hcp Al, and hcp Ti. In all cases, the calculated Ed in the multilayers are intermediate to the corresponding bulk values but exhibit some important directionality. In the nanolayer, defect detection demonstrated systematic differences in the behavior of Ed in each layer. Importantly, collision cascade damage exhibits significant defect partitioning within the Al and Ti layers that is hypothesized to be an intrinsic property of dissimilar nanolayered systems. This type of partitioning could be partly responsible for observed asymmetric radiation damage responses in many multilayered systems. In addition, a pseudo-random direction was introduced to approximate the average Ed without performing numerous simulations with random directions.

Al-Ti↗

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↗

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↗

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

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

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

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

Guest Molecule Exchange Kinetics for the 2012 Ignik Sikumi Gas Hydrate Field Trial

A commercially viable technology for producing methane from natural gas hydrate reservoirs remains elusive. Short-term depressurization field tests have demonstrated the potential for producing natural gas via dissociation of the clathrate structure, but the long-term performance of the depressurization technology ultimately requires a heat source to sustain the dissociation. A decade of laboratory experiments and theoretical studies have demonstrated the exchange of pure CO2 and N2-CO2 mixtures with CH4 in sI gas hydrates, yielding critical information about molecular mechanisms, recoveries, and exchange kinetics. Findings indicated the potential for producing natural gas with little to no production of water and rapid exchange kinetics, generating sufficient interest in the guest-molecule exchange technology for a field test. In 2012 the U.S. DOE/NETL, ConocoPhillips Company, and Japan Oil, Gas and Metals National Corporation jointly sponsored the first field trial of injecting a mixture of N2-CO2 into a CH4-hydrate bearing formation beneath the permafrost on the Alaska North Slope. Known as the Ignik Sikumi #1 Gas Hydrate Field Trial, this experiment involved three stages: 1) the injection of a N2-CO2 mixture into a targeted hydrate-bearing layer, 2) a 4-day pressurized soaking period, and 3) a sustained depressurization and fluid production period. Data collected during the three stages of the field trial were made available after an extensive quality check. These data included continuous temperature and pressure logs, injected and recovered fluid compositions and volumes. The Ignik Sikumi #1 data set is extensive, but contains no direct evidence of the guest-molecule exchange process. This investigation is directed at using numerical simulation to provide an interpretation of the collected data. A numerical simulator, STOMP-HYDT-KE, was recently completed that solves conservation equations for energy, water, mobile fluid guest molecules, and hydrate guest molecules, for up to three gas hydrate guest molecules: CH4, CO2, and N2. The independent tracking of mobile fluid and hydrate guest molecules allows for the kinetic exchange of guest molecules between the mobile fluids and hydrate. The particular interest of this numerical investigation is to determine whether kinetic exchange parameters, determined from laboratory-scale experiments, are directly applicable to interpreting the Ignik Sikumi #1 data.

Alaska North Slope↗

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

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↗

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↗

Analysis of InGaN light-emitting diodes with GaN-AlGaN and AlGaN-GaN composition-graded barriers

The effects of InGaN-based light-emitting diodes (LEDs) with Al composition increasing and decreasing GaN-AlGaN barriers along the growth direction are studied numerically. Simulation results suggest that the LEDs with GaN-AlGaN composition-decreased barriers show more significant enhancement of light-output power and internal quantum efficiency than LEDs with composition-increasing GaN-AlGaN barriers when compared with the conventional LED with GaN barriers, due to the improvement in hole injection efficiency and electron blocking capability. Moreover, the optical performance is further improved by replacing GaN-AlGaN barriers with AlGaN-GaN barriers of the same Al composition-decreasing range, which are mainly attributed to the modified band diagrams. In addition, the major causes of the different efficiency droop behaviors for all the designed structures are explained by the electron leakage current and the different increase rates of hole concentration with injection current.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ion sound instability driven by the ion flows

Ion sound instabilities driven by the ion flow in a system of a finite length are considered by analytical and numerical methods. The ion sound waves are modified by the presence of stationary ion flow resulting in negative and positive energy modes. The instability develops due to coupling of negative and positive energy modes mediated by reflections from the boundary. It is shown that the wave dispersion due to deviation from quasineutrality is crucial for the stability. In finite length system, the dispersion is characterized by the length of the system measured in units of the Debye length. The instability is studied analytically and the results are compared with direct, initial value numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗