Search NASASearch

SEARCH · Search NASA

Results for “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.

Mechanical, Electrochemical & Thermal Modeling of Electric Vehicle Batteries for Crash Simulation (CRADA CRD-19-00811 Final Report)

Under the proposed effort in partnership with Hyundai Motor Company (HMC), the National Laboratory of the Rockies (NLR) will cooperate with HMC to develop mathematical models for battery cells and modules for simulating abuse response in batteries subject to the type of mechanical crushing that can occur in a full motor vehicle crash.

33 ADVANCED PROPULSION SYSTEMS

Thermomechanical Modeling of Woven Materials With Particle-Based, Explicit-Fiber Simulations

Fiber-based materials are extensively used to protect spacecraft during entry. Insulative fibers, often in a fiber network or woven, provide rigidity, strength, and control of material anisotropy and density. Woven thermal protection materials, such as ADEPT (Adaptable, Deployable Entry and Placement Technology), 3D-MAT (3-Dimensional Multifunctional Ablative Thermal Protection), and 3MDCP (3D Woven Mid-Density Carbon Phenolic), enable missions with stronger and denser materials for entry profiles with high shear and heat flux. Vulnerabilities to woven thermal protection materials include manufacturing-induced material property variation, and impact from micrometeoroids. Simulating woven materials under these conditions require models that can resolve hierarchal structures, thermomechanical behavior, and failure. To address this, we simulate weave thermal conduction and mechanical deformation. We simulate the full weave with a coarse-grained yarn model is presented. The model combines a validated, high-resolution single 3MDCP yarn model and phenolic resin model. Instead of modeling every fiber, each yarn ply with order 10, instead of order 1000, fibers. The discrete element bonded particle model (DEM-BPM) of fibers captures the thermal and mechanical behavior within and between fibers. We study the proportion of heat transfer and stress via the contact network, fiber bonds, and overall weave geometry.

bonded particle

Advanced Model Development for Large Eddy Simulation of Oxy-Combustion and Supercritical Carbon Dioxide Power Cycles

A joint experimental and numerical study is performed to observe the characteristics of a supercritical carbon dioxide turbulent mixing layer in the presence of strong nonlinearities in the thermodynamic and transport properties. A bespoke experimental setup is designed and employed for this purpose and provides insight into macroscopic mixing behavior. The mixing is experimentally observed using two techniques: shadowgraphy and spontaneous Raman scattering. Qualitative and quantitative intensity fields obtained via these techniques yield instantaneous and mean density data. Spanwise temperature data is also collected using analogue resistance temperature detectors. These measurements are used to quantify the level of mixed material within the field. The experimental data are supplemented by a companion high-fidelity numerical study. The numerical results are obtained through fully resolved, three-dimensional direct numerical simulation. The numerical dataset permits observation of the near-field mixing characteristics, which are difficult to measure experimentally due to the rapid dynamics and sharp thermophysical gradients in this area. Qualitative field visualizations are presented, followed by quantitative mixed material results and observations regarding thermodynamic property trends at select locations within the field. One-dimensional spectra of the turbulent kinetic energy and solenoidal dissipation are provided to observe the spectral characteristics of the flow. Reynolds stress anisotropy is analyzed graphically through anisotropy invariance maps (Lumley triangles). The mixing quantification, spectral data and anisotropy analysis of a flow at these thermodynamic conditions represent the main outcomes of the work.

20 FOSSIL-FUELED POWER PLANTS

Spacecraft Fire Safety Predictions using Verified Saffire Model

A model developed using Fire Dynamics Simulator (FDS) that aimed to determine the effect of a fire in a spacecraft was validated by data collected during the Saffire campaign. The model used inlet and outlet temperatures and CO 2 concentrations of the Saffire payload where fire spread was taking place to determine the amount of heat and combustion products that made it into Northrop Grumman’s Cygnus vehicle. The model was then validated using six remote sensors in various places in the vehicle, as well as a far field device (FFD) in the open zenith section that was representative of average vehicle values. The current work focuses on using the model to predict fire safety scenarios. One simulation aimed to determine the fate of HCl, which sticks to surfaces. The model prediction showed that the HCl was removed from the atmosphere rapidly. This compared well against the FFD data in the Saffire VI campaign event where a bottle of 5% HCl was released into the vehicle. An additional simulation where the Environmental Control and Life Support System (ECLSS) was shut off once the FFD reached 5 ppm of HCl showed that HCl stayed in the atmosphere considerably longer. Continuing the simulation with the ECLSS activated and after temperatures returned to their initial conditions, resulted in a rapid removal of HCL similar to that observed in the original HCl release scenario model. Finally, a simulation that used the heat release rate from a lithium-ion battery test to determine the effect it would have on a spacecraft was performed. This simulation used a heat addition rate that is considerably higher than what was determined from the burning of solid fuels in the Saffire campaign and hence produced a non-trivial temperature increase in more locations within the vehicle.

Fire Safety

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

Accurate prediction of aerothermodynamic loads in thermal non-equilibrium flows requires precise modeling of internal energy exchange. While previous direct simulation Monte Carlo frameworks have successfully implemented discrete rotational energy models for diatomic species, the treatment of polyatomic molecules has traditionally relied on continuous energy assumptions that break down at low temperatures and neglect critical high-temperature corrections. This study extends the discrete rotational energy models of Boyd and Gimelshein to fully encompass polyatomic molecules. The proposed framework implements quantized rotational energy level sampling for linear, spherical, and symmetric/asymmetric top rotors. Crucially, the model incorporates centrifugal distortion to address the limitations of the rigid-rotor assumption at hypersonic temperatures, and accounts for nuclear spin parity, which dictates the permissible rotational states and macroscopic specific heats at low temperatures. The model is verified through equilibrium sampling procedures, demonstrating agreement with theoretical quantum Boltzmann distributions and accurately reproducing thermophysical properties across a wide range of temperatures.

DSMC

Preliminary Dynamic Modeling of the Quarter-Scale Distributed Electric Propulsion Aircraft

This presentation describes the early-stage modeling of a quarter-scale distributed electric propulsion aircraft, based on the SUbsonic Single Aft eNgine (SUSAN) Electrofan, a transformative concept aircraft for which a model exists. The full-scale 180 passenger SUSAN concept has a single turbofan engine in the tail that both produces thrust and provides electrical power to 16 electric fans distributed across the wings utilizing a series/parallel hybrid architecture. The quarter-scale version would have an internal combustion piston engine to provide power to 17 electric fans–16 on the wings, one in the tail–in a series hybrid configuration, and no vertical or horizontal stabilizers. It would also have a reduced flight envelope in terms of both altitude and speed. The initial modeling approach is to scale down the original airframe model to capture the dynamic behavior of a much smaller aircraft, albeit with an empennage. The original powertrain model is then replaced with one representing the physical components of that of the quarter-scale vehicle. This powertrain model is suitable for control design and analysis, and the fully integrated, although preliminary, aircraft model allows flight simulator testing and evaluation. Results from simulations are presented.

electrified aircraft propulsion

Preliminary Dynamic Modeling of the Quarter-Scale Distributed Electric Propulsion Aircraft

This paper describes the early-stage modeling of a quarter-scale distributed electric propulsion aircraft, based on the SUbsonic Single Aft eNgine (SUSAN) Electrofan, a transformative concept aircraft for which a model exists. The full-scale 180 passenger SUSAN concept has a single turbofan engine in the tail that both produces thrust and provides electrical power to 16 electric fans distributed across the wings utilizing a series/parallel hybrid architecture. The quarter-scale version would have an internal combustion piston engine to provide power to 17 electric fans–16 on the wings, one in the tail–in a series hybrid configuration, and no vertical or horizontal stabilizers. It would also have a reduced flight envelope in terms of both altitude and speed. The initial modeling approach is to scale down the original airframe model to capture the dynamic behavior of a much smaller aircraft, albeit with an empennage. The original powertrain model is then replaced with one representing the physical components of that of the quarter-scale vehicle. This powertrain model is suitable for control design and analysis, and the fully integrated, although preliminary, aircraft model allows flight simulator testing and evaluation. Results from simulations are presented.

Electrified powertrain

Probabilistic Resource Adequacy Suite (PRAS) v0.8 Model Documentation

The Probabilistic Resource Adequacy Suite, or PRAS, is a software package for studying power system resource adequacy. It allows the user to simulate power system operations under a wide range of operating conditions, in order to study the system's risk of failing to meet demand due to a resource shortfall, and identify the time periods and regions in which that risk occurs. This reports documents version 0.8 of the tool.

29 ENERGY PLANNING, POLICY, AND ECONOMY

An Investigation Regarding the Effects of Volume Dynamics on Gas Turbine Engine Performance and Operability

This work explores the impact of adding volume dynamics to an electrified turbofan engine model. The zero-dimensional single-input single-output volume elements capture the physics of fluid compressibility of flow in the engine core gas path turbomachinery during transient operation with greater fidelity. Some background regarding the motivation for this effort is provided. A model of the volume element is developed from first principles. The baseline engine model is described and an explanation is given of the modifications made to incorporate the volume elements into the new model. Details of the simulated inputs, designed to elicit exaggerated dynamic responses and to facilitate comparison between the baseline and modified engine models, are explicated. The results of these simulations are presented and discussed, demonstrating that the transient performance and operability of the modified engine model, measured by parameters such as shaft speeds, thrust, and high-pressure compressor (HPC) stall margin, are impacted by the additional dynamics. Commentary is provided regarding potential explanations for the differences in behavior exhibited by the two models. Finally, suggestions are offered for future work that might improve the volume models developed in this effort and leverage the improved fidelity enabled by the modeling of turbomachinery volume dynamics.

Volume dynamics

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science

Extending mARC II Arc-Jet Test Duration via Design and Implementation of Model System Cooling Sleeve

The mARC II is a 30 kW arc-jet facility at NASA Ames Research Center used to generate high-enthalpy flows for low-cost thermal protection system (TPS) technology development. Sustained operation of downstream instrumentation and material samples is constrained by thermal loading transmitted through the arc-jet test environment, limiting achievable run times and experimental throughput. This work presents the design, integration, and validation of a cooling sleeve implemented on the sweep arm drive motor feedthrough to mitigate thermal accumulation during testing. The addition of the cooling sleeve is a simple, robust upgrade that translates directly into enhanced facility capability by supporting longer run durations, reduced turnaround time, and higher throughput.

numerical simulations

AutoBEM: A scalable framework for nationwide building energy simulation and retrofit evaluation in the United States

This paper presents AutoBEM, an integrated, automated framework for nationwide building energy modeling and retrofit evaluation in the United States. Unlike prior UBEM platforms that either rely primarily on representative stock sampling or operate at city scale, AutoBEM automates the generation of building-resolved, physics-based EnergyPlus/OpenStudio simulation models at national scale using GIS-derived geometry, prototype-based assumptions, and standardized scalable workflows. Leveraging the Model America dataset and high-performance computing, AutoBEM generates and simulates energy models for 122.9 million buildings, representing 97.8% of the U.S. building stock. These models are being made publicly and freely available as the Model America v1.0 (MAv1) dataset. AutoBEM supports detailed, building-level assessments of energy consumption, CO2 emissions, and post-processed anthropogenic heat emissions (AHE), and evaluates 151 energy conservation measures (ECMs) using localized utility pricing and building characteristics. In addition, AutoBEM incorporates both typical and future climate conditions through integration with Typical Meteorological Year (TMY) and Future TMY (fTMY) weather data derived from IPCC scenarios. In a case study of Phoenix, Arizona, AutoBEM identified several high-efficiency HVAC upgrades and selected envelope measures with short modeled payback periods (1.5 years) for certain building types and standards. Simulations under future climate scenarios (SSP5–RCP8.5) project an 11.3% increase in electricity use and a 32% reduction in natural gas demand by 2100, underscoring the need for climate-adaptive retrofit planning. By enabling reproducible, bottom-up, and location-specific analysis at scale, AutoBEM provides a step toward a national digital twin of the built environment and supports data-driven screening and planning for decarbonization, resilience, and energy equity.

Li, Hang [ORNL] (ORCID:0000000306001920)

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics

Interactions Between Climate Policy and Technology-influenced Travel Behavior: Mitigating Induced Demand from CACC

Advances in vehicle technology have influenced the development of automated vehicle systems, where vehicles that do not require human intervention are already deployed in the roadway networks. While these advances are proved to increase roadway safety and highway capacity, more research is needed to understand the long-term and regional-level impacts on mobility, land use, energy consumption, and emissions. This study proposes a multi-model approach to analyze the effect of vehicle automation and deep decarbonization policies over a period from 2020 to 2040 in Austin, Texas. We use the Global Change Analysis Model (GCAM) to develop internally the scenarios that are then passed to the SMART Mobility modeling workflow, a large-scale simulation framework combining the POLARIS activity-based travel demand model and mesoscopic traffic simulator with the Autonomie vehicle energy consumption model and the UrbanSim land use simulator. Results suggest that the introduction of vehicles with advanced automation could increase fuel consumption when no decarbonization policies are implemented. Also, advances in vehicle technology research and development could lead to a decline in energy use in the long-term. Energy pricing and vehicle electrification incentives could help reduce the impact of vehicle automation. Finally, our analysis indicates the relevance of introducing land use processes in longterm vehicle automation studies.

land use

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven

Hierarchical Multi-agent Large Language Model Reasoning for Autonomous Heterogeneous Catalyst Discovery

Artificial intelligence is reshaping scientific exploration, but most methods automate procedural tasks without engaging in scientific reasoning, limiting autonomy in discovery. We demonstrate that hierarchical agentic large language model reasoning can efficiently drive simulation and scientific exploration. Across two chemical applications, CO adsorption on Cu surface transition metal adatoms and on M–N–C catalysts, reasoning-guided exploration reduces required atomistic simulations by up to 90% relative to heuristic or random selection. Comparisons across single-agent, multi-agent, and stochastic baselines show that hierarchical strategies yield more coherent and information-efficient search trajectories. Reasoning traces reveal chemically grounded decisions that cannot be explained by semantic bias or stochastic sampling. We realize these agentic reasoning strategies in Materials Agents for Simulation and Theory in Electronic-structure Reasoning (MASTER), a multimodal system that translates natural language into density functional theory workflows. Altogether, multi-agent collaboration accelerates heterogeneous catalyst discovery and marks a step toward more autonomous, reasoning-guided scientific exploration.

30 DIRECT ENERGY CONVERSION