Search NASA⌕ Search

SEARCH · Search NASA

Results for “Advanced Modeling and 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,603 records · Page 89

Development of a Reduced-Order Three-Dimensional Flow Model for Thermal Mixing and Stratification Simulation during Reactor Transients

Mixing, thermal-stratification, and mass transport phenomena in large pools or enclosures play major roles for the safety of reactor systems. Depending on the fidelity requirement and computational resources, various modeling methods, from the 0-D perfect mixing model to 3-D Computational Fluid Dynamics (CFD) models, are available. Each is associated with its own advantages and shortcomings. It is very desirable to develop an advanced and efficient thermal mixing and stratification modeling capability embedded in a modern system analysis code to improve the accuracy of reactor safety analyses and to reduce modeling uncertainties. An advanced system analysis tool, SAM, is being developed at Argonne National Laboratory for advanced non-LWR reactor safety analysis. While SAM is being developed as a system-level modeling and simulation tool, a reduced-order three-dimensional module is under development to model the multi-dimensional flow and thermal mixing and stratification in large enclosures of reactor systems. This paper provides an overview of the three-dimensional finite element flow model in SAM, including the governing equations, stabilization scheme, and solution methods. Additionally, several verification and validation tests are presented, including lid-driven cavity flow, natural convection inside a cavity, laminar flow in a channel of parallel plates. Based on the comparisons with the analytical solutions and experimental results, it is demonstrated that the developed 3-D fluid model can perform very well for a wide range of flow problems.

Hu, Rui↗

Analysis of Radiation Transport Due to Activated Coolant in the ITER Neutral Beam Injection Cell

Detailed spatial distributions of the biological dose rate due to a variety of sources are required for the design of the ITER tokamak facility to ensure that all radiological zoning limits are met. During operation, water in the Integrated loop of Blanket, Edge-localized mode and vertical stabilization coils, and Divertor (IBED) cooling system will be activated by plasma neutrons and will flow out of the bioshield through a complex system of pipes and heat exchangers. This paper discusses the methods used to characterize the biological dose rate outside the tokamak complex due to 16 N gamma radiation emitted by the activated coolant in the Neutral Beam Injection (NBI) cell of the tokamak building. Activated coolant will enter the NBI cell through the IBED Primary Heat Transfer System (PHTS), and the NBI PHTS will also become activated due to radiation streaming through the NBI system. To properly characterize these gamma sources, the production of 16 N, the decay of 16 N, and the flow of activated water through the coolant loops were modeled. The impact of conservative approximations on the solution was also examined. Once the source due to activated coolant was calculated, the resulting biological dose rate outside the north wall of the NBI cell was determined through the use of sophisticated variance reduction techniques. The AutomateD VAriaNce reducTion Generator (ADVANTG) software implements methods developed specifically to provide highly effective variance reduction for complex radiation transport simulations such as those encountered with ITER. Using ADVANTG with the Monte Carlo N-particle (MCNP) radiation transport code, radiation responses were calculated on a fine spatial mesh with a high degree of statistical accuracy. In conclusion, advanced visualization tools were also developed and used to determine pipe cell connectivity, to facilitate model checking, and to post-process the transport simulation results.

59 BASIC BIOLOGICAL SCIENCES↗

Toward Realistic Models of the Electrocatalytic Oxygen Evolution Reaction

The electrocatalytic oxygen evolution reaction (OER) supplies the protons and electrons needed to transform renewable electricity into chemicals and fuels. However, the OER is kinetically sluggish; it operates at significant rates only when the applied potential far exceeds the reversible voltage. The origin of this overpotential is hidden in a complex mechanism involving multiple electron transfers and chemical bond making/breaking steps. Our desire to improve catalytic performance has then made mechanistic studies of the OER an area of major scientific inquiry, though the complexity of the reaction has made understanding difficult. While historically, mechanistic studies have relied solely on experiment and phenomenological models, over the past twenty years ab initio simulation has been playing an increasingly important role in developing our understanding of the electrocatalytic OER and its reaction mechanisms. In this Review we cover advances in our mechanistic understanding of the OER, organized by increasing complexity in the way through which the OER is modeled. We begin with phenomenological models built using experimental data before reviewing early efforts to incorporate ab initio methods into mechanistic studies. We go on to cover how the assumptions in these early ab initio simulations-no electric field, electrolyte, or explicit kinetics-have been relaxed. Through comparison with experimental literature, we explore the veracity of these different assumptions. We summarize by discussing the most critical open challenges in developing models to understand the mechanisms of the OER.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante↗

Conceptual Design of Integrated Energy Systems via Multiscale Market Simulations and Surrogate Models for Market Interactions

This final report describes novel capabilities developed as part of DISPATCHES, Design Integration and Synthesis Platform to Advance Tightly Coupled Hybrid Energy Systems, for designing hybrid energy systems (HES) a.k.a. integrated energy systems (IES) in the context of a larger electricity market. Capabilities are demonstrated on case studies for nuclear and renewable power.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Realistic 3D Radiative Modeling of Turbulent Structure of Moderate-Mass Stars and Sun

Understanding the turbulent dynamics of the Sun and stars is a critical element for interpreting observed processes and phenomena on different scales and for predicting extreme events such as flares and super flares. High-­‐resolution observations of the Sun and high-­‐fidelity radiative MHD numerical simulations have substantially advanced our understanding of solar and stellar local dynamics and magnetism from the upper convection zone to the atmosphere and corona. However, global modeling of the Sun with such a high degree of realism is currently not affordable due to the extremely high computational cost of resolving the scales in the convection zone. The physics of the deep solar dynamics can be effectively addressed through modeling more massive solar-­‐type stars where the convection zone is shallower and the convective overturning time is much shorter than those on the Sun. We present recent 3D realistic simulation results of moderate-­‐mass stars and discuss links between solar and stellar dynamics, such as the multiscale structure of granulation, convective overshoot, and others. In particular, these simulations have provided better understanding of the dynamics of the tachocline (the overshoot layer at the bottom of the convection zone) and have explained long-­‐ standing results from helioseismology.

Realistic↗

Wind Energy Modeling and Simulation, Volume 1: Atmosphere and Plant

In order to optimise the yield of wind power from existing and future wind plants, the entire breadth of the system of a plant, from the wind field to the turbine components, needs to be modelled in the design process. The modelling and simulation approaches used in each subsystem as well as the system-wide solution methods to optimize across subsystem boundaries are described in this reference. Chapters are written by technical experts in each field, describing the current state of the art in modelling and simulation for wind plant design. This comprehensive, two-volume research reference will provide long-lasting insight into the methods that will need to be developed for the technology to advance into its next generation. Volume 1 covers the computing challenges in full turbine modelling, then discusses bridging scales in the atmosphere and turbulence modelling, wind forecasting, wind plant flow, and plant level controller design.

17 WIND ENERGY↗

xEMU Spacesuit Boot Lunar South Pole Thermal Vacuum Testing

The Exploration Extravehicular Mobility Unit (xEMU) spacesuit lunar boot was evaluated in the Jet Propulsion Lab (JPL) CITADEL thermal vacuum chamber at lunar South Pole conditions. The JPL CITADEL chamber provides a radiation environment and contact plate temperature of 48K, simulating the temperature of a lunar Permanently Shadowed Region (PSR) at the South Pole of the Moon. The main objectives of the test series are: to evaluate the risk lunar PSRs pose to advanced spacesuit boots, validate the xEMU lunar boot thermal model, and to develop a test methodology to evaluate lunar spacesuit boots. To simulate the thermal effects of a crewmember’s foot inside the boot, a thermal manikin foot was developed to inject heat into the spacesuit boot. Thermal data was collected through a variety of sensors including thermocouples, thermistors, heat flux sensors, and infrared cameras. A test matrix was developed to investigate the variables of foot contact force, simulated skin temperature, external contact plate temperature, and ventilation flowrate. The JPL CITADEL chamber includes a loadlock (airlock) that afforded the ability to simulate the transition of the boot from a thermally neutral IVA condition to a cold EVA condition within seconds. This paper will document the test design, test objectives, success criteria, test support equipment, spacesuit boot to chamber interfaces, and configuration of the test article. Additionally, this paper will document the high-level results of the testing. More detailed test results will be presented in a different paper.

Ben Swartout↗

An improved pressure drop correlation for modeling localized effects in a pebble bed reactor

Advances in the development of pebble bed reactors (PBRs) has created a desire for accurate and cost-effective simulation tools for design scoping studies and safety analysis. The current state-of-the-art for these simulations is the use of porous media models, although these models rely on correlations to capture the effects of flow features that are not explicitly modeled. One of the areas where correlation accuracy is currently lacking is in the near-wall region of the bed. In this region, the presence of the wall causes the pebbles to pack more orderly, drastically changing the geometry and flow behavior in this region. This work presents a new generalized pressure drop correlation for PBRs based on the KTA equation. A high-to-low methodology is applied, where large eddy simulation (LES) is performed on two beds of 1568 and 1700 pebbles to generate a high-fidelity dataset. The flow fields are then averaged in time and separated into concentric rings of 0.05 D peb width. Average porosity, velocity, and pressure drop are extracted for each ring and the friction and form losses are calculated. The Reynolds number range for this study is 625–10,000, and thus the form losses are dominant over the friction losses. The form losses across the rings are investigated, and a correction term for the form loss calculation is determined and applied to the KTA equation to drastically improve the capability of modeling localized porosity effects in a porous media code. Finally, the improved correlation reduces near-wall velocity prediction error from over 50% with the KTA correlation to around 5%. Agreement in pressure drop prediction between LES and porous media simulations is also improved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

Simulated Polarimetric Fields of Ice Vapor Growth Using the Adaptive Habit Model. Part II: A Case Study from the FROST Experiment

Abstract A new adaptive habit model (AHM) grows ice crystals through vapor deposition while evolving ice particle properties, including shape and effective density. The AHM provides an opportunity to investigate observed microphysical processes through the computation of polarimetric variables and corroboration with microphysical model output. This study is unique because the polarimetric scattering calculations are computed using predicted microphysical parameters rather than a priori assumptions that are imposed within the scattering calculations in the forward simulator, allowing for a more effective comparison to radar observations. Through the simulation of a case in the Front Range of the Rocky Mountains in Colorado using the Advanced Research version of the Weather Research and Forecasting Model, it is found that the AHM approximates ice mass, shape, cloud vertical structure, and temporal evolution as reflected through polarimetric quantities compared to observations. AHM reflectivity magnitudes are similar to those observed with radar and are an improvement over spherical ice crystal assumptions. Further analyses are completed to examine the effect of microphysical processes on the evolution of the differential reflectivity and specific differential phase, both of which are simulated using the AHM. Simulations reveal a polarimetric response to ice crystal mass, number, size, density, and aspect ratio. While results reveal the need for model improvements (e.g., parameterizations for aggregation rate), testing forward-simulated radar fields against observations is a first step in the validation of model microphysical and precipitation processes.

54 ENVIRONMENTAL SCIENCES↗

Compilation of Abstracts for SC12 Conference Proceedings

1 A Breakthrough in Rotorcraft Prediction Accuracy Using Detached Eddy Simulation; 2 Adjoint-Based Design for Complex Aerospace Configurations; 3 Simulating Hypersonic Turbulent Combustion for Future Aircraft; 4 From a Roar to a Whisper: Making Modern Aircraft Quieter; 5 Modeling of Extended Formation Flight on High-Performance Computers; 6 Supersonic Retropropulsion for Mars Entry; 7 Validating Water Spray Simulation Models for the SLS Launch Environment; 8 Simulating Moving Valves for Space Launch System Liquid Engines; 9 Innovative Simulations for Modeling the SLS Solid Rocket Booster Ignition; 10 Solid Rocket Booster Ignition Overpressure Simulations for the Space Launch System; 11 CFD Simulations to Support the Next Generation of Launch Pads; 12 Modeling and Simulation Support for NASA's Next-Generation Space Launch System; 13 Simulating Planetary Entry Environments for Space Exploration Vehicles; 14 NASA Center for Climate Simulation Highlights; 15 Ultrascale Climate Data Visualization and Analysis; 16 NASA Climate Simulations and Observations for the IPCC and Beyond; 17 Next-Generation Climate Data Services: MERRA Analytics; 18 Recent Advances in High-Resolution Global Atmospheric Modeling; 19 Causes and Consequences of Turbulence in the Earths Protective Shield; 20 NASA Earth Exchange (NEX): A Collaborative Supercomputing Platform; 21 Powering Deep Space Missions: Thermoelectric Properties of Complex Materials; 22 Meeting NASA's High-End Computing Goals Through Innovation; 23 Continuous Enhancements to the Pleiades Supercomputer for Maximum Uptime; 24 Live Demonstrations of 100-Gbps File Transfers Across LANs and WANs; 25 Untangling the Computing Landscape for Climate Simulations; 26 Simulating Galaxies and the Universe; 27 The Mysterious Origin of Stellar Masses; 28 Hot-Plasma Geysers on the Sun; 29 Turbulent Life of Kepler Stars; 30 Modeling Weather on the Sun; 31 Weather on Mars: The Meteorology of Gale Crater; 32 Enhancing Performance of NASAs High-End Computing Applications; 33 Designing Curiosity's Perfect Landing on Mars; 34 The Search Continues: Kepler's Quest for Habitable Earth-Sized Planets.

HPC↗

Comparison of the Effect of 2 at. % Additions of Nb and Ta on the 1100 °C Oxidation Behavior of Ni-6Al-(4,6,8) Cr Model Alloys

To continue improving alloy performance in harsh service environments, the development of alumina-forming nickel-based superalloys is essential. Current generations of these alloys heavily rely on the addition of refractory elements to enhance their mechanical properties at high temperatures; however, a systematic understanding of how such additions affect the overall oxidation behavior is still not well established, particularly from the standpoint of predicting the transition from internal to external alumina formation. The present work seeks to better understand the intrinsic effects that common minor additions of Ta and Nb have on the oxidation behavior of alumina-scale-forming γ-Ni model alloys. By combining a novel simulation approach with high-temperature oxidation experiments and advanced characterization techniques, the present study provides insightful details on the differing effects that 2 at. % addition of Ta and Nb have on the alumina scale formation of Ni-based alloys during 1100 °C oxidation.

Rodriguez, Rafael↗

New preconditioning strategy for Jacobian-free solvers for variably saturated flows with Richards’ equation

We develop a new approach for solving the nonlinear Richards’ equation arising in variably saturated flow modeling. The growing complexity of geometric models for simulation of subsurface flows leads to the necessity of using unstructured meshes and advanced discretization methods. Typically, a numerical solution is obtained by first discretizing PDEs and then solving the resulting system of nonlinear discrete equations with a Newton-Raphson-type method. Efficiency and robustness of the existing solvers rely on many factors, including an empiric quality control of intermediate iterates, complexity of the employed discretization method and a customized preconditioner. We propose and analyze a new preconditioning strategy that is based on a stable discretization of the continuum Jacobian. We will show with numerical experiments for challenging problems in subsurface hydrology that this new preconditioner improves convergence of the existing Jacobian-free solvers 3-20 times. Furthermore, we show that the Picard method with this preconditioner becomes a more efficient nonlinear solver than a few widely used Jacobian-free solvers.

58 GEOSCIENCES↗

Simulation of supersonic molecular beam injection fueling into H-mode plasmas on EAST using BOUT++

Transport dynamic 2D simulations for supersonic molecular beam injection (SMBI) fueling into H-mode deuterium plasmas on the Experimental Advanced Superconducting Tokamak (EAST) is first simulated using a seven-field two-fluid model in the BOUT++ framework. The SMB is assumed to be injected into plasma from the midplane at the low field side with a fixed width and constant molecular flux. The different densities and injection velocities of SMBI are investigated within the upper single-null geometry of EAST. The simulations indicate that the SMBI has a self-shielding effect on molecules' inward transport into the plasma, and the deposition of SMBI leads to a large increase in plasma density and decrease in plasma temperature. There is a velocity threshold for SMB penetrating the pedestal and depositing at the top of density pedestal. However, the deposition point would be back toward the plasma boundary after SMB arriving at the deepest penetration position. Comparing the different molecular injection velocities and densities, the outcomes demonstrate that the depth of deposition is closely related to the injection velocity rather than the molecular density. The simulated results show good agreement with the EAST experiments by comparing the electron density profiles obtained by simulation and experiment, respectively. These results will be helpful for guiding future experiments, as well as the design of the SMBI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Perspectives on Research in Computational Plasma Physics With Applications to Experiments

Computational plasma physics was born some 50 years ago with the development of high-speed computers. Computer simulation is well suited to the challenging nature of fundamental and applied plasma science, which exhibits enormous ranges of time and space scales, and whose underlying mathematical framework is comprised of nonlinear partial differential equations and nonlinear kinetic equations. Computational plasma physics has advanced from relatively simple calculations with limited dimensionality and scope to sophisticated and comprehensive simulation models. Moreover, computational plasma physics has matured as a significant scientific discipline with a rigorous mathematical foundation and voluminous literature. While computational plasma physics has greatly benefitted from the growth in computing capability by many orders of magnitude, the contributions from innovation in methods and algorithms have been no less important. This talk presents a personal perspective on the development and application of computational plasma physics to plasmas in nature and in the laboratory. The examples described will be based on experience over the course of a career in computational plasma physics and illustrate both the fundamental plasma phenomena and the behavior of laboratory plasmas. Some of the specific examples include simulations of microinstabilities in the magnetic mirror, tokamak and spheromak plasmas, laser-plasma interactions, and the Knudsen-layer phenomena, as it affects fusion performance. The examples also illustrate the development of appropriate models and algorithms that are well suited to simulating the phenomena of interest efficiently. Finally, a particular interest of the author has been the development of multiple time-scale algorithms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FY2018 Analysis Annual Progress Report

During fiscal year 2018 (FY 2018), the U.S. Department of Energy Vehicle Technologies Office (VTO) funded analysis projects supportive of VTO’s goals to pursue early stage research in vehicle and mobility system technologies to improve energy efficiency, increase energy reliability and security, improve transportation affordability, and promote economic growth. VTO analysis projects result in a foundation of data, analytical models, and applied analyses that provide insights into critical transportation energy problems and assist in research investment prioritization and portfolio planning.

33 ADVANCED PROPULSION SYSTEMS↗

Commercialization of a Comprehensive Spark-Ignition Model for Automotive Engine Applications in CONVERGE CFD

The work completed under this CRADA commercializes an advanced hybrid Lagrangian-Eulerian spark-ignition (LESI) model for multi-dimensional engine simulations. In this report, LESI is coupled with combustion models (WSR, TFM, G-equation) to model spark channel elongation and kernel and flame development in a constant volume combustion vessel, and a four-valve direct-injection spark-ignition (4VDISI) engine in two operating conditions: stoichiometric E30 with no dilution, and stoichiometric PACE-20 with dilution using CONVERGE CFD solver.

42 ENGINEERING↗