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,405 records · Page 78

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

RELAP-7 User's Guide

The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). The code is based on the INL’s modern scientific software development framework, MOOSE (Multi-Physics Object Oriented Simulation Environment). The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical integration methods, and physical models. The end result will be a reactor systems analysis capability that retains and improves upon RELAP5’s capability and extends the analysis capability for all reactor system simulation scenarios.This document is a user's guide on how to use the RELAP-7 code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Electrical Resistivity Tomography based monitoring of stress perturbations to optimize placement of high-precision strain meters

The Center for Understanding Subsurface Signals and Permeability is a new U.S. Department of Energy Earthshot Center focused on understanding and predicting the long-term evolution of permeability in enhanced geothermal systems. The center will use a highly instrumented testbed within the Sanford Underground Research Facility to conduct field scale experiments that elucidate and test capabilities to simulate geochemical-geomechanical interactions and permeability evolution. Here we demonstrate initial developments using previously collected electrical resistivity tomography (ERT) monitoring data with high-performance multi-physics modelling advancements to inform the optimal location of two new monitoring boreholes. Specifically, ERT monitoring data collected during shear stimulation testing shows marked responses to changes in stress during borehole pressurization. We demonstrate how the same response is being simulated, ultimately to train a machine-learning algorithm to estimate rock properties and enable enhanced prediction of stress and strain responses anticipated during future testing campaigns.

Stress, EGS, CUSSP, 3D Electrical Imaging↗

Efficient and Collision-Free Anticipative Cruise Control in Randomly Mixed Strings

Connected intelligent vehicle following can improve safety and efficiency compared to today's road transport, but real traffic in the near future will not provide an ideal setting for its deployment. Unconnected human-driven vehicles will follow variable behavior patterns that combine with differences in dynamic capability to create heterogeneous scenarios. In this paper, connected automated vehicles employ model predictive control for following in traffic that may include both heavy and passenger vehicles at quasi-random positions. Adding further realism, some quasi-randomly mixed vehicles are completely unconnected and equipped with the reactive intelligent driver model using pseudorandom parameters. A mixed-integer quadratic programming formulation adapts the predictive algorithm to diverse powertrain operating point constraints. The preceding vehicle's control input is estimated from velocity and brake light observations and used to probabilistically generate a preview for the ego vehicle. Here, a terminal constraint designed using particle kinematics prevents collisions due to shortsightedness. Results are simulated at various heavy and predictive vehicle concentrations in the presence of packet loss. Linear fuel economy improvements between 1.4% and 1.9% per 10% increase in predictive vehicle penetration rate are shown.

33 ADVANCED PROPULSION SYSTEMS↗

Preliminary SAM Assessment

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory (Argonne) for advanced non-LWR safety analysis. To assist NRC to assess SAM capabilities for advanced reactor safety analysis and licensing at the NRC, a series of verification and other standard tests (OSTs) are modeled in SAM and code simulation results are compared with available analytical results. This report documents the preliminary SAM assessment using a matrix of six test problems. Each test problem is further examined with different boundary conditions, system configurations, or modeling options. Although relatively simple, these tests cover the basic equation models, basic component models, and basic system level processes and phenomena that must be modeled for advanced reactor safety analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A high-resolution large-eddy simulation framework for wildland fire predictions using TensorFlow

Background: Wildfires are becoming more severe, so we need improved tools to predict them over a wide range of conditions and scales. One approach towards this goal entails the use of coupled fire/atmosphere modelling tools. Although significant progress has been made in advancing their physical fidelity, existing tools have not taken full advantage of emerging programming paradigms and computing architectures to enable high-resolution wildfire simulations. Aims: The aim of this study was to present a new framework that enables landscape-scale wildfire simulations with physical representation of combustion at an affordable cost. Methods: We developed a coupled fire/atmosphere simulation framework using TensorFlow, which enables efficient and scalable computations on Tensor Processing Units. Key Results: Simulation results for a prescribed fire were compared with experimental data. Predicted fire behavior and statistical analysis for fire spread rate, scar area, and intermittency showed overall reasonable agreement. Scalability analysis was performed, showing close to linear scaling. Conclusions: While mesh refinement was shown to have less impact on global quantities, such as fire scar area and spread rate, it benefits predictions of intermittent fire behavior, buoyancy-driven dynamics, and small-scale turbulent motion. Implications: This new simulation framework is efficient in capturing both global quantities and unsteady dynamics of wildfires at high spatial resolutions.

54 ENVIRONMENTAL SCIENCES↗

DiaMonD: An Integrated Multifaceted Approach to Mathematics at the Interfaces of Data, Models, and Decisions

For many problems in computational science and engineering, the so-called forward problem-solution of the underlying mathematical model to yield output quantities of interest, given input parameters--is difficult enough for frontier complex models, which are often characterized by coupled multiphysics and possibly stochastic behavior over a wide range of length and time scales. However, we face the need to move beyond the forward problem, to address what is often the ultimate goal of CS&E: decision-making. This requires us to tackle a spectrum of mathematical problems that subsume, and are thus even more difficult than, the forward problem. First, given experimental data, we wish to estimate unknown parameters characterizing a model of a natural or engineered system by solving an inverse problem. Second, we seek the optimal con figuration of the system (or experiments) by solving an optimal design problem. Third, the optimal operation of the system must be determined by solving an optimal control problem. And fourth, we must quantify uncertainties as they propagate through all of the preceding problems, from data to model inference to prediction and finally to optimal design and control. Optimal mathematical methods and algorithms for such end-to-end, data-to-decisions modeling and simulation of complex problems--requiring integrated solution of forward, inverse, optimization, and UQ problems for large, multiphysics, multiscale models--entails challenges of the highest order, while also presenting great opportunities for applied mathematics research. Unfortunately, research on forward, inverse, optimization, and UQ problems has all-too-often progressed in isolation. This has led to mathematical methods that perform well when considering the forward simulation in isolation, but are prohibitive or suboptimal or unstable when combined with other methods within the framework of inversion, optimization, or UQ. The reverse is also true: general-purpose methods developed within the inversion, optimization, and UQ fields often become prohibitive when applied to complex CS&E problems since they do not exploit their structure. DiaMonD is a multi-institutional DOE MMICC effort that aims to address the challenges of end-to-end, data-to-decisions modeling and simulation for complex CS&E problems in a unifi ed, systematic, and integrated fashion. Institutions involved are Colorado State, Florida State, Los Alamos, MIT, Oak Ridge, Stanford, and UT-Austin. The goals of DiaMonD are (1) to develop advanced mathematical methods and analysis for multimodel, multiphysics, multiscale model problems driven by frontier DOE applications, including those in subsurface energy and environmental flows, materials for energy storage and conversion, and ice sheet dynamics; (2) to create theory and algorithms for integrated inversion, optimization, and uncertainty quantification for these complex problems; and (3) to disseminate the philosophy of an integrated end-to-end, data-to-decisions approach to modeling and simulation of complex problems to the broader applied math and computational science communities through workshops and other outreach.

97 MATHEMATICS AND COMPUTING↗

Advanced Simulation Tools for Muon-Based Accelerators

The project developed new software tools incorporating the most accurate theoretical calculations and experimental data available for crucial and not-yet-considered physics processes specific to muon accelerators. These new tools were applied to modeling muon ionization cooling channels for a muon collider or a neutrino factory, the facilities that have the potential to discover and explore new exciting fundamental physics but require the development of demanding technologies and innovative concepts. G4beamline, COSY Infinity, and WARP computer codes are particularly well-suited for addressing muon accelerator modeling problems. G4beamline was chosen for its flexibility and modularity, and for being one of the de facto standard codes used by the muon accelerator community. COSY was selected for speed, ability to take into account high-order nonlinear effects essential for the large beams under consideration, and availability of built-in optimization tools. WARP was chosen for its space charge calculation capabilities, speed and flexibility of implementing new lattice elements such as RF cavities. This final technical report summarizes the efforts on implementing advanced new tools in these computer codes addressing the shortcomings in modeling muon physics processes and using these tools to model muon ionization cooling channels.

43 PARTICLE ACCELERATORS↗

Accelerating Our Understanding of Supernova Explosion Mechanism via Simulations and Visualizations with GenASiS

Core-collapse supernovae are among the most powerful explosions in the Universe, releasing about 1053 erg of energy on timescales of a few tens of seconds. These explosion events are also responsible for the production and dissemination of most of the heavy elements, making life as we know it possible. Yet exactly how they work is still unresolved. One reason for this is the sheer complexity and cost of a self-consistent, multi-physics, and multi-dimensional core-collapse supernova simulation, which is impractical, and often impossible, even on the largest supercomputers we have available today. To advance our understanding we instead must often use simplified models, teasing out the most important ingredients for successful explosions, while helping us to interpret results from higher fidelity multi-physics models. In this paper we investigate the role of instabilities in the core-collapse supernova environment. We present here simulation and visualization results produced by our code GenASiS.

Budiardja, R. D.↗

Fuel Performance Modeling Internship Final Presentation

This study examines the performance of U-Zr and U-Pd-Zr annular metallic fuel rodlets and details the current status of modeling efforts regarding U-Pu-Zr solid metallic fuel rodlets which were experimentally evaluated as part of the Advanced Fuels Campaign (AFC). The AFC mission is to develop novel fuel technologies and facilitate the implementation of those technologies by industry partners. A key objective is to improve steady-state and transient performance over current fuel types. The experiments of interest in this study included metallic fuel rodlets within HT-9 cladding which were placed in SS-316 capsules and inserted in the Advanced Test Reactor (ATR). Certain mechanical and thermal conditions cannot be directly evaluated through experiments and fuel performance modeling is used to shed light on this evolution over time. In this study, BISON Multiphysics simulations are leveraged to investigate the state of the fuel system throughout and after the experimental conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Using simulation to study solid oxide cell degradation

Slides of the presentation given by Harry Abernathy from the National Energy Technology Laboratory at the 2023 International Conference and Expo on Advanced Ceramics and Composites. The presentation was an overview of NETL's solid oxide cell performance degradation modeling, including a demonstration of the impact of microstructural analysis methods on the resulting degradation simulations.

Abernathy, Harry↗

Advancing the science of explosive fragmentation and afterburn fireballs though experiments and simulations at the benchtop scale

Detonation of explosive devices produces extremely hazardous fragments and hot, luminous fireballs. Prior experimental investigations of these post-detonation environments have primarily considered devices containing hundreds of grams of explosives. While relevant to many applications, such large- scale testing also significantly restricts experimental diagnostics and provides limited data for model validation. As an alternative, the current work proposes experiments and simulations of the fragmentation and fireballs from commercial detonators with less than a gram of high explosive. As demonstrated here, reduced experimental hazards and increased optical access significantly expand the viability of advanced imaging and laser diagnostics. Notable developments include the first known validation of MHz-rate optical fragment tracking and the first ever Coherent Anti-Stokes Raman Scattering (CARS) measures of post-detonation fireball temperatures. While certainly not replacing the need for full-scale verification testing, this work demonstrates new opportunities to accelerate developments of diagnostics and predictive models of post-detonation environments.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Design Analysis Kit for Optimization and Terascale Applications 6.0

Sandia's Dakota software (available at http://dakota.sandia.gov) supports science and engineering transformation through advanced exploration of simulations. Specifically it manages and analyzes ensembles of simulations to provide broader and deeper perspective for analysts and decision makers. This enables them to: (1) enhance understanding of risk, (2) improve products, and (3) assess simulation credibility. In its simplest mode, Dakota can automate typical parameter variation studies through a generic interface to a computational model. However, Dakota also delivers advanced parametric analysis techniques enabling design exploration, optimization, model calibration, risk analysis, and quantification of margins and uncertainty with such models. It directly supports verification and validation activities. The algorithms implemented in Dakota aim to address challenges in performing these analyses with complex science and engineering models from desktop to high performance computers.

Source record↗

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗

DiaMonD: An Integrated Multifaceted Approach to Mathematics at the Interfaces of Data, Models, and Decisions

DiaMonD is a multi-institutional DOE MMICC effort that aims to address the challenges of end-to-end, data-to-decisions modeling and simulation for complex computational science and engineering problems in a unified, systematic, and integrated fashion. Institutions involved are Colorado State, Florida State, Los Alamos, MIT, Oak Ridge, Stanford, and UT-Austin. The goals of DiaMonD are (1) to develop advanced mathematical methods and analysis for multimodel, multiphysics, multiscale model problems driven by frontier DOE applications, including those in subsurface energy and environmental flows, materials for energy storage and conversion, and ice sheet dynamics; (2) to create theory and algorithms for integrated inversion, optimization, and uncertainty quantification for these complex problems; and (3) to disseminate the philosophy of an integrated end-to-end, data-to-decisions approach to modeling and simulation of complex problems to the broader applied math and computational science communities through workshops and other outreach.

42 ENGINEERING↗

Recent advancements and remaining challenges associated with inner magnetosphere cross-energy/population interactions (IMCEPI)

The geospace inner magnetosphere, within about ten Earth radii, contains various plasma populations with energy from a few eV to MeV and plays important roles in regulating the energy density of the magnetosphere, the magnetic field configuration, and wave dynamics. As an integrated part of the magnetosphere, the inner magnetosphere region also ties to other regions and can change the global geospace circulation. Therefore understanding both internal and external cross–energy/population interactions can help further our knowledge of the inner magnetosphere dynamics and non–linear feedback processes. In view of this, in the past five years (2014–2018), the GEM Focus Group (FG) “Inner magnetosphere cross–energy/population interactions (IMCEPI)” has gathered and boosted community–wide interactions among observation, simulation and modeling studies. As a result, this commentary reports some major accomplishments of the interactive inner magnetosphere community that were advanced by the IMCEPI FG discussions and layouts remaining challenges that need to be carried on.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Sepsis reconsidered: Identifying novel metrics for behavioral landscape characterization with a high-performance computing implementation of an agent-based model

Objectives: Sepsis affects nearly 1 million people in the United States per year, has a mortality rate of 28–50% and requires more than $20 billion a year in hospital costs. Over a quarter century of research has not yielded a single reliable diagnostic test or a directed therapeutic agent for sepsis. Central to this insufficiency is the fact that sepsis remains a clinical/physiological diagnosis representing a multitude of molecularly heterogeneous pathological trajectories. Advances in computational capabilities offered by High Peformance Computing (HPC) platforms call for an evolution in the investigation of sepsis to attempt to define the boundaries of traditional research (bench, clinical and computational) through the use of computational proxy models. We present a novel investigatory and analytical approach, derived from how HPC resources and simulation are used in the physical sciences, to identify the epistemic boundary con- ditions of the study of clinical sepsis via the use of a proxy agent-based model of systemic inflammation. Design: Current predictive models for sepsis use correlative methods that are limited by patient heterogeneity and data sparseness. We address this issue by using an HPC version of a system-level validated agent-based model of sepsis, the Innate Immune Response ABM (IIRBM), as a proxy system in order to identify boundary conditions for the possible behavioral space for sepsis. We then apply advanced analysis derived from the study of Random Dynamical Systems (RDS) to identify novel means for characterizing system behavior and providing insight into the tractability of traditional investigatory methods. Results: The behavior space of the IIRABM was examined by simulating over 70 million sepsis patients for up to 90 days in a sweep across the following parameters: cardio-respiratory-metabolic resilience; microbial invasiveness; microbial toxigenesis; and degree of nosocomial exposure. In addition to using established methods for describing parameter space, we developed two novel methods for characterizing the behavior of a RDS: Probabilistic Basins of Attraction (PBoA) and Stochastic Trajectory Analysis (STA). Computationally generated behavioral landscapes demonstrated attractor structures around stochastic re- gions of behavior that could be described in a complementary fashion through use of PBoA and STA. The stochasticity of the boundaries of the attractors highlights the challenge for correlative attempts to characterize and classify clinical sepsis. Conclusions: HPC simulations of models like the IIRABM can be used to generate approximations of the behavior space of sepsis to both establish “boundaries of futility” with respect to existing investi- gatory approaches and apply system engineering principles to investigate the general dynamic properties of sepsis to provide a pathway for developing control strategies. Here, the issues that bedevil the study and treatment of sepsis, namely clinical data sparseness and inadequate experimental sampling of system behavior space, are fundamental to nearly all biomedical research, manifesting in the “Crisis of Reproducibility” at all levels. HPC-augmented simulation-based research offers an investigatory strategy more consistent with that seen in the physical sciences (which combine experiment, theory and simulation), and an opportunity to utilize the leading advances in HPC, namely deep machine learning and evolutionary computing, to form the basis of an iterative scientific process to meet the full promise of Precision Medicine (right drug, right patient, right time).

59 BASIC BIOLOGICAL SCIENCES↗

LES and RANS Spray Combustion Analysis of OME 3-5 and n-Dodecane

Clean-burning oxygenated and synthetic fuels derived from renewable power, so-called e-fuels, are a promising pathway to decarbonize compression–ignition engines. Polyoxymethylene dimethyl ethers (PODEs or OMEs) are one candidate of such fuels with good prospects. Their lack of carbon-to-carbon bonds and high concentration of chemically bound oxygen effectively negate the emergence of polycyclic aromatic hydrocarbons (PAHs) and even their precursors like acetylene (C 2 H 2 ), enabling soot-free combustion without the soot-NO x trade-off common for diesel engines. The differences in the spray combustion process for OMEs and diesel-like reference fuels like n-dodecane and their potential implications on engine applications include discrepancies in the observed ignition delay, the stabilized flame lift-off location, and significant deviations in high-temperature flame morphology. For CFD simulations, the accurate modeling and prediction of these differences between OMEs and n-dodecane proved challenging. This study investigates the spray combustion process of an OME 3 − 5 mixture and n-dodecane with advanced optical diagnostics, Reynolds-Averaged Navier–Stokes (RANS), and Large-Eddy Simulations (LESs) within a constant-volume vessel. Cool-flame and high-temperature combustion were measured simultaneously via high-speed (50 kHz) imaging with formaldehyde (CH 2 O) planar laser-induced fluorescence (PLIF) representing the former and line-of-sight OH* chemiluminescence the latter. Both RANS and LES simulations accurately describe the cool-flame development process with the formation of CH 2 O. However, CH 2 O consumption and the onset of high-temperature reactions, signaled by the rise of OH* levels, show significant deviations between RANS, LES, and experiments as well as between n-dodecane and OME. A focus is set on the quality of the simulated results compared to the experimentally observed spatial distribution of OH*, especially in OME fuel-rich regions. The influence of the turbulence modeling is investigated for the two distinct ambient temperatures of 900 K and 1200 K within the Engine Combustion Network Spray A setup. The capabilities and limitations of the RANS simulations are demonstrated with the initial cool-flame propagation and periodic oscillations of CH 2 O formation/consumption during the quasi-steady combustion period captured by the LES.

CFD↗