Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct 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 1,765 records · Page 98

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗

Direct Simulation Monte Carlo Studies of the Gas Sampling for the VATMOS-SR Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept which would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere, and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus, with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a full picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of atmospheric volatile signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process in order to determine how measured isotopic ratios of noble gases in the sample compare to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept, and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed at Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield a comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters as well as the modeling parameters will be performed.

direct simulation Monte Carlo↗

Interaction Between Tropical Convection and its Embedding Environment: An Energetics Analysis of a 2-D Cloud Resolving Simulation

The phase relation between the perturbation kinetic energy (K') associated with the tropical convection and the horizontal-mean moist available potential energy (bar-P) associated with environmental conditions is investigated by an energetics analysis of a numerical experiment. This experiment is performed using a 2-D cloud resolving model forced by the TOGA-COARE derived vertical velocity. The imposed upward motion leads to a decrease of bar-P directly through the associated vertical advective cooling, and to an increase of K' directly through cloud related processes, feeding the convection. The maximum K' and its maximum growth rate lags and leads, respectively, the maximum imposed large-scale upward motion by about 1-2 hours, indicating that convection is phase locked with large-scale forcing. The dominant life cycle of the simulated convection is about 9 hours, whereas the time scales of the imposed large-scale forcing are longer than the diurnal cycle. In the convective events, maximum growth of K' leads maximum decay of the perturbation moist available potential energy (P') by about 3 hours through vertical heat transport by perturbation circulation, and perturbation cloud heating. Maximum decay of P' leads maximum decay of bar-P by about one hour through the perturbation radiative, processes, the horizontal-mean cloud heating, and the large-scale vertical advective cooling. Therefore, maximum gain of K' occurs about 4-5 hours before maximum decay of bar-P.

Li, Xiaofan↗

CFD 2030 Grand Challenge: CFD-in-the-Loop Monte Carlo Flight Simulation for Space Vehicle Design

Flight qualification of space vehicles is markedly different from those typically employed for aircraft. The concept of an extensive flight test campaign for a space vehicle does not exist, and vehicle designers must look to alternative techniques for demonstrating robust and reliable performance of their vehicles prior to operational flight. A space vehicle may undergo only a handful of flight tests in its development cycle, with each flight representing a drastically different flight phase or flight configuration. For instance, NASA’s Space Launch System (SLS) launch vehicle and Orion spacecraft will only see a total of four flight demonstrations before flying a crew on its first operational mission, and each flight demonstrates a unique vehicle configuration and/or set of flight conditions. The SLS will be flown only one time before it becomes operational (Artemis 1). The Orion spacecraft Crew Module (CM) will have been tested twice, once on a Delta IV launch vehicle (Exploration Flight Test 1) and once as a fully integrated system with the SLS launch vehicle (Artemis 1). The Orion Launch abort system will have been tested twice, once in a pad abort scenario (Pad Abort 1) and once in an inflight abort scenario (Ascent Abort 2) on a modified Peacekeeper booster. Both of these latter tests involve only a boiler plate CM, not a functional Orion spacecraft. Thus, unlike aircraft, there is very little opportunity for engineers to assess and evaluate their preflight predictions. Instead, space vehicle designers rely on Monte Carlo flight simulations with detailed dispersions of predicted nominal flight behavior to determine how robust their design is to errors and uncertainties in the flight conditions their vehicle may encounter. These Monte Carlo analyses entail thousands of trajectory simulations to demonstrate that the vehicle can meet design requirements at a specified level of reliability. From an aerodynamics and aerothermodynamics perspective, these trajectory simulations are fueled by an extensive aerodynamic database that covers the complete range of expected flight conditions, vehicle configurations, and flight attitudes expected in a given mission. Today, these databases amount to a table of engineering parameters that can be quickly interrogated by the trajectory simulator. The aerodynamic and aerothermodynamic databases are assembled via a series of ground tests, empirical and analytical analysis, physics-based computational analysis, applicable past flight performance data, and in some cases, engineering judgment. These databases generally take years to assemble for a new space vehicle system and in the case of SLS/Orion, over a decade of test and analysis have been expended to develop the extensive databases required to cover the myriad of configurations and potential flight conditions required for the system. Recently, it has been proposed that Computational Fluid Dynamic (CFD) and computing capability may be reaching a point where it is foreseeable that CFD could be integrated directly into the production trajectory simulation tools used to design NASA’s space vehicles. To demonstrate this, NASA has embarked on two demonstrations of this type of capability, one where six degree of freedom flight trajectory simulation equations are embedded in an existing CFD solver and another where a production CFD solver is loosely coupled with a production trajectory simulation tool. These efforts represent an initial demonstration of a future approach to flight trajectory simulation, but they are a far cry from the capability required to perform a full-up CFD-in-the-loop Monte Carlo trajectory simulation. Therefore, this represents a viable grand challenge for computational methods addressing space vehicle design and development. The final paper/presentation will discuss the many hurdles, beyond simply raw computational power, to realizing this grand challenge and how they map directly to the CFD Vision 2030 ojectives. Among these are the wide range of flight conditions, including accelerating/decelerating flight, encountered by a space vehicle during launch and/or entry. The vehicle can also encounter numerous configuration changes, some of which can be quite drastic, during the course of its flight, so robust, automated geometry modeling, grid generation, and adaptation will play a huge role in reaching this goal. Multiply this by 1000’s of trajectory simulations occurring simultaneously in a given Monte Carlo analysis, and the problem readily scales to absorb virtually any size of supercomputer envisioned today. The concept of CFD-in-the-loop Monte Carlo trajectory simulation poses a formidable challenge for emerging and future computing systems, and it has the potential to shave years off the development cycle for aerodynamic and aerothermodynamic performance predictions as compared to today’s space vehicle design approach.

CFD 2030↗

Nonlinear multi-magnon scattering in artificial spin ice

Abstract Magnons, the quantum-mechanical fundamental excitations of magnetic solids, are bosons whose number does not need to be conserved in scattering processes. Microwave-induced parametric magnon processes, often called Suhl instabilities, have been believed to occur in magnetic thin films only, where quasi-continuous magnon bands exist. Here, we reveal the existence of such nonlinear magnon-magnon scattering processes and their coherence in ensembles of magnetic nanostructures known as artificial spin ice. We find that these systems exhibit effective scattering processes akin to those observed in continuous magnetic thin films. We utilize a combined microwave and microfocused Brillouin light scattering measurement approach to investigate the evolution of their modes. Scattering events occur between resonance frequencies that are determined by each nanomagnet’s mode volume and profile. Comparison with numerical simulations reveals that frequency doubling is enabled by exciting a subset of nanomagnets that, in turn, act as nanosized antennas, an effect that is akin to scattering in continuous films. Moreover, our results suggest that tunable directional scattering is possible in these structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Amplified spontaneous emission in the extreme ultraviolet by expanding xenon clusters

Focused short-wavelength free-electron laser (FEL) pulses interacting with gas phase samples can induce by inner-shell ionization a short-lived population inversion, followed by coherent collective emission of directed, short, and strong radiation bursts. We extend our studies into the warm-dense matter (WDM) regime by investigating the nanoplasmas produced in an ensemble of nanometer-sized clusters by FEL irradiation. Here, additional pathways can also lead to strong, laserlike emission: Electron-ion collisions can yield a long-lived population inversion, and subsequent amplified spontaneous emission. We observe amplified spontaneous emission (ASE) in the extreme ultraviolet in xenon clusters excited by soft x-ray FEL pulses, we diagnose the generated nanoplasmas by fluorescence spectroscopy, and we study under various cluster and FEL parameters the directed ASE from the Xe 2+ 65 nm line. We show its exponential increase as a function of FEL irradiation power, and an accompanying collisional broadening of the emission spectra. These findings are corroborated by extensive numerical simulations based on theory, combining detailed hydrodynamic and kinetic simulations with time-dependent calculations of radiation transport, amplification, and collective emission in the WDM nanoplasma. Furthermore, our theoretical findings underline that population inversion is due to electron-ion collisions and that the observed decoherence processes can be empirically characterized by a phenomenological decoherence time in the range of 100–200 fs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A spacecraft-borne gradiometer mission analysis

Numerical simulations were performed to obtain the orbit- and attitude-determination requirements of a spacecraft-borne gradiometer mission. Results demonstrated that position determination of 300 meters in the along-track and cross-track directions and 50 meters in the radial direction are mission requirements. The optimal orientation of the gradiometer sensing plane is achieved when the spin vector elevation is 0 degrees. The attitude-determination requirements are 5 degrees resolution for spin-vector azimuth and 0.2 degree resolution for spin-vector elevation. When these requirements are met, 3-degree gravity anomalies can be recovered globally with an accuracy of 0.025/mm/sq s (2.5 mgals). The Appendix documents the mathematical procedures for estimating detailed gravity fields from gradiometer data.

Argentiero, P.↗

Efficient numerical simulation of an electrothermal de-icer pad

In this paper, a new approach to calculate the transient thermal behavior of an iced electrothermal de-icer pad was developed. The method of splines was used to obtain the temperature distribution within the layered pad. Splines were used in order to create a tridiagonal system of equations that could be directly solved by Gauss elimination. The Stefan problem was solved using the enthalpy method along with a recent implicit technique. Only one to three iterations were needed to locate the melt front during any time step. Computational times were shown to be greatly reduced over those of an existing one dimensional procedure without any reduction in accuracy; the curent technique was more than 10 times faster.

Roelke, R. J.↗

Analysis and visualization of complex unsteady three-dimensional flows

Flow field animation is the natural choice as a tool in the analysis of the numerical simulations of complex unsteady three-dimensional flows. The PLOT4D extension of the widely used PLOT3D code to allow the interactive animation of a broad range of flow variables was developed and is presented. To allow direct comparison with unsteady experimental smoke and dye flow visualization, the code STREAKER was developed to produce time accurate streaklines. Considerations regarding the development of PLOT4D and STREAKER, and example results are presented.

Van Dalsem, William R.↗

NASA's supercomputing experience

A brief overview of NASA's recent experience in supercomputing is presented from two perspectives: early systems development and advanced supercomputing applications. NASA's role in supercomputing systems development is illustrated by discussion of activities carried out by the Numerical Aerodynamical Simulation Program. Current capabilities in advanced technology applications are illustrated with examples in turbulence physics, aerodynamics, aerothermodynamics, chemistry, and structural mechanics. Capabilities in science applications are illustrated by examples in astrophysics and atmospheric modeling. Future directions and NASA's new High Performance Computing Program are briefly discussed.

Bailey, F. Ron↗

Mean Line Pump Flow Model in Rocket Engine System Simulation

A mean line pump flow modeling method has been developed to provide a fast capability for modeling turbopumps of rocket engines. Based on this method, a mean line pump flow code PUMPA has been written that can predict the performance of pumps at off-design operating conditions, given the loss of the diffusion system at the design point. The pump code can model axial flow inducers, mixed-flow and centrifugal pumps. The code can model multistage pumps in series. The code features rapid input setup and computer run time, and is an effective analysis and conceptual design tool. The map generation capability of the code provides the map information needed for interfacing with a rocket engine system modeling code. The off-design and multistage modeling capabilities of the code permit parametric design space exploration of candidate pump configurations and provide pump performance data for engine system evaluation. The PUMPA code has been integrated with the Numerical Propulsion System Simulation (NPSS) code and an expander rocket engine system has been simulated. The mean line pump flow code runs as an integral part of the NPSS rocket engine system simulation and provides key pump performance information directly to the system model at all operating conditions.

Veres, Joseph P.↗

A statistical model of the human core-temperature circadian rhythm

We formulate a statistical model of the human core-temperature circadian rhythm in which the circadian signal is modeled as a van der Pol oscillator, the thermoregulatory response is represented as a first-order autoregressive process, and the evoked effect of activity is modeled with a function specific for each circadian protocol. The new model directly links differential equation-based simulation models and harmonic regression analysis methods and permits statistical analysis of both static and dynamical properties of the circadian pacemaker from experimental data. We estimate the model parameters by using numerically efficient maximum likelihood algorithms and analyze human core-temperature data from forced desynchrony, free-run, and constant-routine protocols. By representing explicitly the dynamical effects of ambient light input to the human circadian pacemaker, the new model can estimate with high precision the correct intrinsic period of this oscillator ( approximately 24 h) from both free-run and forced desynchrony studies. Although the van der Pol model approximates well the dynamical features of the circadian pacemaker, the optimal dynamical model of the human biological clock may have a harmonic structure different from that of the van der Pol oscillator.

NASA Discipline Regulatory Physiology↗

Doubly-Shadowed Regions in Lunar Polar Craters: Hydrogen Accumulation in the Presence of Recursive Plasma Wakes

Permanently shadowed regions (PSRs) of the Moon have been identified as unique environments of extreme cold and comprise a natural cold trap for sequestering volatiles [Paige et al. 2010]. The diverse chemical composition of the LCROSS impact plume provided evidence for a volatile-rich and chemically-complex PSR environment [Cola prete et al. 2010, Schultz et al. 2010]. Additionally, the polar electrostatic environment is highly complex, with the possibility of strong, localized electric fields that divert solar wind ions directly into polar cold traps [Farrell et al. 2010, Zimmerman et al. 2011]. Thus, regional plasma physics processes couple directly with volatile sequestration. In the present work, kinetic simulations show that recursive plasma wake structure arises in the presence of step-like topographic features (Le. doubly-shadowed craters). Combining the plasma code with a numerical sputtering model demonstrates that solar wind protons can be either a hydrogen source via implantation or a volatile loss mechanism via sputtering, depending on properties of the regolith and solar wind. The present model provides a novel theoretical pathway toward understanding the lunar surface/solar wind physical and chemical interactions for complex topography near the poles.

Zimmerman, Michael I.↗

Updates to Simulation of a Single-Element Lean-Direct Injection Combustor Using Arbitary Polyhedral Meshes

This paper summarizes the procedures of (1) generating control volumes anchored at the nodes of a mesh; and (2) generating staggered control volumes via mesh reconstructions, in terms of either mesh realignment or mesh refinement, as well as presents sample results from their applications to the numerical solution of a single-element LDI combustor using a releasable edition of the National Combustion Code (NCC).

Power↗

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗

An Object Oriented Extensible Architecture for Affordable Aerospace Propulsion Systems

Driven by a need to explore and develop propulsion systems that exceeded current computing capabilities, NASA Glenn embarked on a novel strategy leading to the development of an architecture that enables propulsion simulations never thought possible before. Full engine 3 Dimensional Computational Fluid Dynamic propulsion system simulations were deemed impossible due to the impracticality of the hardware and software computing systems required. However, with a software paradigm shift and an embracing of parallel and distributed processing, an architecture was designed to meet the needs of future propulsion system modeling. The author suggests that the architecture designed at the NASA Glenn Research Center for propulsion system modeling has potential for impacting the direction of development of affordable weapons systems currently under consideration by the Applied Vehicle Technology Panel (AVT). This paper discusses the salient features of the NPSS Architecture including its interface layer, object layer, implementation for accessing legacy codes, numerical zooming infrastructure and its computing layer. The computing layer focuses on the use and deployment of these propulsion simulations on parallel and distributed computing platforms which has been the focus of NASA Ames. Additional features of the object oriented architecture that support MultiDisciplinary (MD) Coupling, computer aided design (CAD) access and MD coupling objects will be discussed. Included will be a discussion of the successes, challenges and benefits of implementing this architecture.

Follen, Gregory J.↗

An Overview of the Waste-to-Energy System Simulation (WESyS) Model

The leveraging of waste streams for energy and chemical production could add revenue to waste disposal operations, and it presents opportunities for addressing a variety of economic and environmental objectives at the local, state, and national levels. The Waste-to-Energy System Simulation (WESyS) model is a system dynamics model that was created to simulate the development of the U.S. waste-to-energy industry over time. For each of the three primary waste resources modeled (landfills, concentrated animal feeding operations, and publicly owned treatment works), WESyS simulates technically feasible scenarios for use of the waste, including direct conversion to fuels, and anaerobic digestion followed by flaring, electricity generation, combined heat and power, cleanup and compression to compressed natural gas, and cleanup and injection into an existing pipeline. The model allows users to explore numerous plausible future scenarios for the development of the U.S. waste-to-energy industry. This report provides an overview of the WESyS model and documents the key assumptions, equations, and data sources used to create the model.

09 BIOMASS FUELS↗