Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Simulations”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Carbon Dioxide Dispersion in the Combustion Integrated Rack Simulated Numerically

When discharged into an International Space Station (ISS) payload rack, a carbon dioxide (CO2) portable fire extinguisher (PFE) must extinguish a fire by decreasing the oxygen in the rack by 50 percent within 60 sec. The length of time needed for this oxygen reduction throughout the rack and the length of time that the CO2 concentration remains high enough to prevent the fire from reigniting is important when determining the effectiveness of the response and postfire procedures. Furthermore, in the absence of gravity, the local flow velocity can make the difference between a fire that spreads rapidly and one that self-extinguishes after ignition. A numerical simulation of the discharge of CO2 from PFE into the Combustion Integrated Rack (CIR) in microgravity was performed to obtain the local velocity and CO2 concentration. The complicated flow field around the PFE nozzle exits was modeled by sources of equivalent mass and momentum flux at a location downstream of the nozzle. The time for the concentration of CO2 to reach a level that would extinguish a fire anywhere in the rack was determined using the Fire Dynamics Simulator (FDS), a computational fluid dynamics code developed by the National Institute of Standards and Technology specifically to evaluate the development of a fire and smoke transport. The simulation shows that CO2, as well as any smoke and combustion gases produced by a fire, would be discharged into the ISS cabin through the resource utility panel at the bottom of the rack. These simulations will be validated by comparing the results with velocity and CO2 concentration measurements obtained during the fire suppression system verification tests conducted on the CIR in March 2003. Once these numerical simulations are validated, portions of the ISS labs and living areas will be modeled to determine the local flow conditions before, during, and after a fire event. These simulations can yield specific information about how long it takes for smoke and combustion gases produced by a fire to reach a detector location, how large the fire would be when the detector alarms, and the behavior of the fire until it has been extinguished. This new capability could then be used to optimize the location of fire detectors and fire-suppression ports as well as to evaluate the effectiveness of fire suppressants and response strategies. Numerical data collected from these simulations could also be used to develop a virtual reality fire event for crew training and fire safety awareness. This work is funded by NASA's Bioastronautics Initiative, which has the objective of ensuring and enhancing the health, safety, and performance of humans in space. As part of this initiative, the Microgravity Combustion Science Branch at the NASA Glenn Research Center is conducting spacecraft fire safety research to significantly improve fire safety on inhabited spacecraft.

Wu, Ming-Shin↗
DOE OSTI2020

Toward Exascale: Overview of Large Eddy Simulations and Direct Numerical Simulations of Nuclear Reactor Flows with the Spectral Element Method in Nek5000

At the beginning of the last decade, Petascale supercomputers (i.e., computers capable of more than 1 petaFLOP) emerged. Now, at the dawn of exascale supercomputing, we provide a review of recent landmark simulations of portions of reactor components with turbulence-resolving techniques that this computational power has made possible. In fact, these simulations have provided invaluable insight into flow dynamics, which is difficult or often impossible to obtain with experiments alone. We focus on simulations performed with the spectral element method, as this method has emerged as a powerful tool to deliver massively parallel calculations at high fidelity by using large eddy simulation or direct numerical simulation. We also limit this paper to constant-property incompressible flow of a Newtonian fluid in the absence of other body or external forces, although the method is by no means limited to this class of flows. We briefly review the fundamentals of the method and the reasons it is compelling for the simulation of nuclear engineering flows. We review in detail a series of Petascale simulations, including the simulations of helical coil steam generators, fuel assemblies, and pebble beds. Even with Petascale computing, however, limitations for nuclear modeling and simulation tools remain. In particular, the size and scope of turbulence-resolving simulations are still limited by computing power and resolution requirements, which scale with the Reynolds number. In the final part of this paper, we discuss the future of the field, including recent advancements in emerging architectures such as GPUbased supercomputers, which are expected to power the next generation of high-performance computers.

computational fluid dynamics↗
NASA NTRS1990

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

This research is involved with the implementation of advanced computational schemes based on large eddy simulations (LES) and direct numerical simulations (DNS) to study the phenomenon of mixing and its coupling with chemical reactions in compressible turbulent flows. In the efforts related to LES, a research program to extend the present capabilities of this method was initiated for the treatment of chemically reacting flows. In the DNS efforts, the focus is on detailed investigations of the effects of compressibility, heat release, and non-equilibrium kinetics modelings in high speed reacting flows. Emphasis was on the simulations of simple flows, namely homogeneous compressible flows, and temporally developing high speed mixing layers.

Givi, Peyman↗
NASA NTRS1991

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

This research is involved with the implementations of advanced computational schemes based on large eddy simulations (LES) and direct numerical simulations (DNS) to study the phenomenon of mixing and its coupling with chemical reactions in compressible turbulent flows. In the efforts related to LES, a research program was initiated to extend the present capabilities of this method for the treatment of chemically reacting flows, whereas in the DNS efforts, focus was on detailed investigations of the effects of compressibility, heat release, and nonequilibrium kinetics modeling in high speed reacting flows. The efforts to date were primarily focussed on simulations of simple flows, namely, homogeneous compressible flows and temporally developing hign speed mixing layers. A summary of the accomplishments is provided.

Givi, Peyman↗
NASA NTRS1991

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The main objective is to extend the boundaries within which large eddy simulations (LES) and direct numerical simulations (DNS) can be applied in computational analyses of high speed reacting flows. In the efforts related to LES, we were concerned with developing reliable subgrid closures for modeling of the fluctuation correlations of scalar quantities in reacting turbulent flows. In the work on DNS, we focused our attention to further investigation of the effects of exothermicity in compressible turbulent flows. In our previous work, in the first year of this research, we have considered only 'simple' flows. Currently, we are in the process of extending our analyses for the purpose of modeling more practical flows of current interest at LaRC. A summary of our accomplishments during the third six months of the research is presented.

Givi, P.↗
NASA NTRS1992

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The basic objective of this research is to extend the capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. In the efforts related to LES, we were primarily involved with assessing the performance of the various modern methods based on the Probability Density Function (PDF) methods for providing closures for treating the subgrid fluctuation correlations of scalar quantities in reacting turbulent flows. In the work on DNS, we concentrated on understanding some of the relevant physics of compressible reacting flows by means of statistical analysis of the data generated by DNS of such flows. In the research conducted in the second year of this program, our efforts focused on the modeling of homogeneous compressible turbulent flows by PDF methods, and on DNS of non-equilibrium reacting high speed mixing layers. Some preliminary work is also in progress on PDF modeling of shear flows, and also on LES of such flows.

Givi, P.↗
NASA NTRS1993

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The primary objective of this research is to extend current capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first two years of this research have been concentrated on a priori investigations of single-point Probability Density Function (PDF) methods for providing subgrid closures in reacting turbulent flows. In the efforts initiated in the third year, our primary focus has been on performing actual LES by means of PDF methods. The approach is based on assumed PDF methods and we have performed extensive analysis of turbulent reacting flows by means of LES. This includes simulations of both three-dimensional (3D) isotropic compressible flows and two-dimensional reacting planar mixing layers. In addition to these LES analyses, some work is in progress to assess the extent of validity of our assumed PDF methods. This assessment is done by making detailed companions with recent laboratory data in predicting the rate of reactant conversion in parallel reacting shear flows. This report provides a summary of our achievements for the first six months of the third year of this program.

Givi, P.↗
NASA NTRS1993

Large eddy simulation and direct numerical simulation of high speed turbulent reacting flows

The objective of this research is to make use of Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first phase of this research conducted within the past three years have been directed in several issues pertaining to intricate physics of turbulent reacting flows. In our previous 5 semi-annual reports submitted to NASA LaRC, as well as several technical papers in archival journals, the results of our investigations have been fully described. In this progress report which is different in format as compared to our previous documents, we focus only on the issue of LES. The reason for doing so is that LES is the primary issue of interest to our Technical Monitor and that our other findings were needed to support the activities conducted under this prime issue. The outcomes of our related investigations, nevertheless, are included in the appendices accompanying this report. The relevance of the materials in these appendices are, therefore, discussed only briefly within the body of the report. Here, results are presented of a priori and a posterior analyses for validity assessments of assumed Probability Density Function (PDF) methods as potential subgrid scale (SGS) closures for LES of turbulent reacting flows. Simple non-premixed reacting systems involving an isothermal reaction of the type A + B yields Products under both chemical equilibrium and non-equilibrium conditions are considered. A priori analyses are conducted of a homogeneous box flow, and a spatially developing planar mixing layer to investigate the performance of the Pearson Family of PDF's as SGS models. A posteriori analyses are conducted of the mixing layer using a hybrid one-equation Smagorinsky/PDF SGS closure. The Smagorinsky closure augmented by the solution of the subgrid turbulent kinetic energy (TKE) equation is employed to account for hydrodynamic fluctuations, and the PDF is employed for modeling the effects of scalar fluctuations. The implementation of the model requires the knowledge of the local values of the first two SGS moments. These are provided by additional modeled transport equations. In both a priori and a posteriori analyses, the predicted results are appraised by comparison with subgrid averaged results generated by DNS. Based on these results, the paths to be followed in future investigations are identified.

Adumitroaie, V.↗
DOE OSTI2023

Investigation of Farm-to-Farm Interactions and Blockage Effects from AWAKEN using Large-Scale Numerical Simulations

A large-scale numerical computation of five wind farms was performed as a part of the American Wake Experiment (AWAKEN). This high-fidelity computation used the ExaWind/AMR-Wind LES solver to simulate a 100km × 100km domain containing 541 turbines under unstable atmospheric conditions matching previous measurements. The turbines were represented by Joukowski and OpenFAST coupled actuator disk models. Results of this qualitative comparison illustrate the interactions with wind farms with large scale ABL structures in the flow, as well as the extent of downstream wake penetration in the flow and blockage effects around wind farms.

Cheung, Lawrence C.↗
DOE OSTI2022

Modeling of Cube Array Roughness: RANS, Large Eddy Simulation, and Direct Numerical Simulation

Abstract Flow over arrays of cubes is an extensively studied model problem for rough wall turbulent boundary layers. While considerable research has been performed in computationally investigating these topologies using direct numerical simulation (DNS) and large eddy simulation (LES), the ability of sublayer-resolved Reynolds-averaged Navier–Stokes (RANS) to predict the bulk flow phenomena of these systems is relatively unexplored, especially at low and high packing densities. Here, RANS simulations are conducted on six different packing densities of cubes in aligned and staggered configurations. The packing densities investigated span from what would classically be defined as isolated, up to those in the d-type roughness regime, filling in the gap in the present literature. Three different sublayer-resolved turbulence closure models were tested for each case: a low Reynolds number k–ϵ model, the Menter k–ω SST model, and a full Reynolds stress model. Comparisons of the velocity fields, secondary flow features, and drag coefficients are made between the RANS results and existing LES and DNS results. There is a significant degree of variability in the performance of the various RANS models across all comparison metrics. However, the Reynolds stress model demonstrated the best accuracy in terms of the mean velocity profile as well as drag partition across the range of packing densities.

Engineering↗
NASA NTRS1982

Calculation of pressure statistics in turbulent free shear flows by direct numerical simulation

Accurate numerical simulations of turbulent free shear flows were performed in order to study the behavior of certain critical turbulence quantities which have been difficult to measure in the laboratory. The three-dimensional time-dependent nonlinear Navier-Stokes equations are solved by the direct numerical simulation technique without turbulence modeling or ad hoc closure assumptions. The pressure interaction terms were factored into diffusion and return-to-isotropy components according to two models. Simulations of the evolution of a turbulent wake with and without swirl present were performed at sufficiently high Reynolds numbers to insure that the flows were strongly nonlinear and turbulent, but low enough so that all scales of motion containing significant energy are adequately resolved. It is noted that the presence of swirl tends to enhance the amplitude of the rms pressure levels.

Metcalfe, R. W.↗
DOE OSTI2024

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗
DOE OSTI2021

Investigation of high-pressure turbulent jets using direct numerical simulation

Direct numerical simulations of free round jets at a Reynolds number (Re D ) of 5000, based on jet diameter (D) and jet-exit bulk velocity (U e ), are performed to study jet turbulence characteristics at supercritical pressures. The jet consists of nitrogen (N 2 ) that is injected into N 2 at the same temperature. To understand turbulent mixing, a passive scalar is transported with the flow at unity Schmidt number. Two sets of inflow conditions that model jets issuing from either a smooth contraction nozzle (laminar inflow) or a long pipe nozzle (turbulent inflow) are considered. By changing one parameter at a time, the simulations examine the jet-flow sensitivity to the thermodynamic condition (characterized in terms of the compressibility factor ($Z$) and the normalized isothermal compressibility), inflow condition and ambient pressure (p ∞ ) spanning perfect- to real-gas conditions. The inflow affects flow statistics in the near field (containing the potential core closure and the transition region) as well as further downstream (containing fully developed flow with self-similar statistics) at both atmospheric and supercritical p ∞ . The sensitivity to inflow is larger in the transition region, where the laminar-inflow jets exhibit dominant coherent structures that produce higher mean strain rates and higher turbulent kinetic energy than in turbulent-inflow jets. Decreasing $Z$ at a fixed supercritical p ∞ enhances pressure and density fluctuations (non-dimensionalized by local mean pressure and density, respectively), but the effect on velocity fluctuations depends also on the local flow dynamics. When $Z$ is reduced, large mean strain rates in the transition region of laminar-inflow jets significantly enhance velocity fluctuations (non-dimensionalized by local mean velocity) and scalar mixing, whereas the effects are minimal in jets from turbulent inflow.

42 ENGINEERING↗
NASA NTRS1987

Numerical simulations of aerodynamic contribution of flows about a space-plane-type configuration

The slightly supersonic viscous flow about the space-plane under development at the National Aerospace Laboratory (NAL) in Japan was simulated numerically using the LU-ADI algorithm. The wind-tunnel testing for the same plane also was conducted with the computations in parallel. The main purpose of the simulation is to capture the phenomena which have a great deal of influence to the aerodynamic force and efficiency but is difficult to capture by experiments. It includes more accurate representation of vortical flows with high angles of attack of an aircraft. The space-plane shape geometry simulated is the simplified model of the real space-plane, which is a combination of a flat and slender body and a double-delta type wing. The comparison between experimental results and numerical ones will be done in the near future. It could be said that numerical results show the qualitatively reliable phenomena.

Matsushima, Kisa↗
NASA NTRS1974

Thermal relaxation of a two dimensional plasma in a dc magnetic field. Part 2: Numerical simulation

The thermal relaxation process for a spatially uniform two dimensional plasma in a uniform dc magnetic field is simulated numerically. Thermal relaxation times are defined in terms of the time necessary for the numerically computer Boltzman H-function to decrease through a given part of the distance to its minimum value. Dependence of relaxation time on two parameters is studied: number of particles per Debye square and ratio of gyrofrequency to plasma frequency.

Hsu, J. Y.↗
NASA NTRS1983

Three-dimensional flow over a conical afterbody containing a centered propulsive jet - A numerical simulation

The supersonic flow field over a body of revolution incident to the free stream is simulated numerically on a large, array processor (the CDC Cyber 205). The configuration is composed of a cone-cylinder forebody followed by a conical afterbody from which emanates a centered, supersonic propulsive jet. The free-stream Mach number is 2, the jet-exit Mach number is 2.5, and the jet-to-free-stream static pressure ratio is 3. Both the external flow and the exhaust are ideal air at a common total temperature. The thin-layer approximation to the time-dependent, compressible, Reynolds-averaged Navier-Stokes equations are solved using an implicit finite-difference algorithm. The data base, of 5 million words, is structured in a 'pencil' format so that efficient use of the array processor can be realized. The computer code is completely vectorized to take advantage of the data structure. Turbulence closure is achieved using an empirical algebraic eddy-viscosity model. The configuration and flow conditions correspond to published experimental tests and the computed solutions are consistent with the experimental data.

Deiwert, G. S.↗
NASA NTRS1984

Three-dimensional flow over a conical afterbody containing a centered propulsive jet: A numerical simulation

The supersonic flow field over a body of revolution incident to the free stream is simulated numerically on a large, array processor (the CDC CYBER 205). The configuration is composed of a cone-cylinder forebody followed by a conical afterbody from which emanates a centered, supersonic propulsive jet. The free-stream Mach number is 2, the jet-exist Mach number is 2.5, and the jet-to-free-stream static pressure ratio is 3. Both the external flow and the exhaust are ideal air at a common total temperature.

Deiwert, G. S.↗