Search NASA⌕ Search

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.

2,701 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

SimH 2 : an integrated techno-economic modeling framework for hydrogen pipeline infrastructure and network optimization

Large-scale hydrogen (H 2 ) pipeline transport design and network optimization have seldom been reported due to the lack of a cost model accounting for the relationship between transport cost and hydrogen mass flow rate. Here, this work introduced a system-level cost model for hydrogen pipeline transport at supercritical state and integrated it with an existing CO 2 pipeline network tool, SimCCS, for hydrogen-specific pipeline design and optimization. The Intermountain West (I-West) region of the U.S., historically dependent on fossil fuel-based economies, is chosen to demonstrate the capabilities of our H 2 pipeline cost model and transport network optimization platform called SimH 2 . Two scenarios are examined: one where the pipeline is not allowed to pass through disadvantaged communities and the other where it is permitted. The results highlight that incorporating disadvantaged-community constraints lead to longer pipeline routes and increased transport costs, reflecting the trade-offs involved in equitable infrastructure development. It is demonstrated that the newly developed SimH 2 tool not only enables the efficient design of H 2 transportation pipelines but also optimizes the network by accounting for local terrain and the presence of disadvantaged areas.

08 HYDROGEN↗
DOE OSTI2026

The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations

Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2026

Ice in Convective Storm Environments: Assessing Impacts of Microphysics on Simulated Reflectivities

Satellite radar retrievals of clouds and precipitation rely on assumptions about the microphysical properties of hydrometeors such as particle size distributions (PSDs) and massdimensional relationships, which have been shown to vary in distinct environments. A deeper understanding of linkages between mass-dimensional relationships, ice crystal shape, and radar reflectivities across cloud and precipitation regimes is necessary for developing new retrievals for upcoming satellite radar missions, including the INvestigation of Convective UpdraftS (INCUS), as well as future missions that will feature space-borne radars observing clouds, convection, and precipitation. This study investigates ice microphysical characteristics through aircraft in-situ and remote sensing observations from two convective events during the Mid-latitude Continental Convective Clouds Experiment (MC3E) in the Southern Great Plains, and tests implications for ice particle habit in forward modeled radar observations. During MC3E, the Ka/Ku- band High Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) was onboard the NASA ER-2 aircraft, while in-situ measurements were taken onboard the University of North Dakota Citation. In-situ PSDs and particle shape information from optical array probe data are used to forward model radar reflectivities from particle scattering databases and are compared against HIWRAP reflectivities. Various particle combinations can be used to find agreement with HIWRAP observations, and the forward modeled reflectivities were sensitive to particle type. Varying the riming on the aggregate and graupel particles resulted in larger than 10 dB differences in reflectivities. The in-situ and radar observations imply the occurrence of aggregation, size sorting and lofting particles.

Julia A. Shates↗
NASA NTRS1991

Preparation of the Ice Certification of the Dornier 328 Regional Airliner By Numerical Simulation and By Ground Test

The Dornier 328, a new regional airliner, is to be qualified to FAR25/JAR25 requirements for operation into known icing conditions. All lifting surfaces are pneumatically deiced. Impingement limits were determined using a validated computational fluid dynamic (CFD) code. Wind tunnel tests with a model of the whole aircraft and the isolated empennage equipped with artificial ice shapes showed a degradation of handling characteristics due to ice accretion. Extensive two and three dimensional CFD calculations evaluated ice accretion at the deflection elevator horn. Icing tunnel tests were performed with a model of the horizontal tail with icing similitude fulfilled. The Dornier Do328 powerplant air induction system has successfully completed a program of icing tests in support of aircraft certification for operation into known icing conditions. Potential flight icing conditions were extensively analyzed using a CFD potential flow code and 3-D particle tracking routine to identify critical flight icing conditions. The tunnel test program verified the analytically predicted critical ice accretion surfaces of the air induction system, the adequacy of the ice protection provisions and demonstrated compliance with the applicable (JAR25) regulatory material.

D Welte↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science↗
NASA NTRSDate not supplied

Experimental and Numerical Investigation of the NASA High Efficiency Centrifugal Compressor Vaned Stage Geometry and Aerodynamic Performance

Since its inception in the early 2010s, the NASA High Efficiency Centrifugal Compressor (HECC) has been enigmatic for the propulsion research community: experimental data and numerical simulations of the stage have generally not aligned in their quantifications of performance metrics. Typically, the fault for these disagreements is assigned to the numerical simulations as the simulations are models, and models are inherently incomplete representations of the experiment. This “incompleteness” may manifest in assumptions regarding roughness or heat transfer, simplifications of the flow path (i.e., neglecting bleed flows), or the oft-scapegoated turbulence model. In the case of HECC, recent work showed unexpected discrepancies between the intended impeller geometry defined in the design report (termed Design-Intent) and the manufactured impeller used in the experimental campaigns (termed As-Manufactured). That work used numerical simulations to establish that the geometric differences between the Design-Intent and As-Manufactured impellers were significant enough to result in drastically different performance predictions for the HECC vaneless diffuser configuration. The Design-Intent impeller simulations over predicted the performance relative to the experiment, whereas the As-Manufactured simulations better represented the experimental data, both in terms of one-dimensional performance metrics and spanwise flow profiles. This effort expands on that work by examining in detail the geometry and aerodynamic performance of the HECC vaned diffuser configuration. Further differences between the Design-Intent and As-Manufactured geometries have been discovered in the vaned diffuser and exit guide vanes, and these differences are documented herein. The summations of the geometric differences for all of the components were used to create two numerical models of HECC vaned diffuser configuration: the Design-Intent simulations which are generated from the original geometry definitions given in the design report and the As-Manufactured simulations which are the best available representation of the manufactured compressor hardware used in the experimental test campaigns. In congruence with the earlier vaneless diffuser work, the numerical predictions of the Design-Intent choked mass flow rate, total pressure ratio, and efficiency were notably greater than that of the As-Manufactured simulations. To increase confidence in the experimental dataset, measurements from a recent test campaign conducted in 2024 are used to validate the original experimental data acquired from 2012 to 2014 with good repeatability overall, especially considering the passage of time and differences in the data acquisition systems between the test campaigns. Both numerical simulations were then extensively evaluated against the experimental data. The As-Manufactured simulations provided better estimates of the stage performance than the Design-Intent cases in terms of most performance metrics. Nonetheless, more detailed results still show opportunities for improvement. Despite a more accurate representation of the physical hardware, characterization of the impeller work input remains a challenge even for rigorously developed numerical models.

centrifugal compressor↗
DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗
DOE OSTI2026

SOLPS-ITER modeling of a dedicated divertor for negative triangularity operation on DIII-D

The design of a new dedicated divertor for negative triangularity (NT) operation on DIII-D with neutral baffles and pumping is informed by SOLPS-ITER transport modeling. This dedicated NT divertor is the latest step in a progression of NT shapes with various divertor characteristics explored on DIII-D, including NT shapes at reduced triangularity and a campaign with stronger shaping that included new armored components on the outboard side. SOLPS simulations played a key role in these divertor designs. Interpretive simulations, using cross-field diffusivities constrained by experimental data in the NT Shelf shape were used to inform the design of the 2023 armor campaign components. A similar procedure used armor campaign data to predict conditions for the dedicated NT divertor. The predictive simulations were used to assess the divertor fluxes, detachment threshold, pumped flux, and neutral leakage. For the dedicated NT divertor, SOLPS simulations and two-point-modeling were used to show the relative impact of magnetic topology (mainly longer connection length) and divertor closure on the divertor conditions relative to the armor campaign. It is predicted that the dedicated NT divertor reaches detachment (measured by target ion flux rollover) at a lower upstream density (≈(1.75−−2.4)×1⁢019m−3) as compared to the armor campaign shape. For the preliminary design geometry, divertor closure reduces the neutral leakage by ≈10%. Parametric optimization indicating further ≈20%–60% improvement in the leakage flux and recycled flux crossing the pump entrance is possible for relatively minor changes to the divertor and baffle layout.

Lore, Jeremy [ORNL] (ORCID:000000029192465X)↗
NASA NTRS1985

Battery-Charge-State Model

Charge-state model for lead/acid batteries proposed as part of effort to make equivalent of fuel gage for battery-powered vehicles. Models based on equations that approximate observable characteristics of battery electrochemistry. Uses linear equations, easier to simulate on computer, and gives smooth transitions between charge, discharge, and recuperation.

Vivian, H. C.↗
NASA NTRSDate not supplied

Thermal Data-driven Model Reduction for Enhanced Battery Health Monitoring

Electric aviation faces a major challenge of avoiding potentially catastrophic consequences of the battery’s thermal runaway while keeping the weight of the battery low. Detection of early warning signals of battery failures requires accurate monitoring of the battery’s health throughout its lifespan. However, identifying the parameters of the battery from field data is notoriously difficult. We investigate this problem within the framework of modeling the temperature dynamics of a Li-ion cell during tests simulating loading in electric aircraft flights. It is found that the parameters of a higher-fidelity physics-based thermal model cannot be identified from the simulated flight data. To resolve this issue, we reduce the higher-fidelity thermal model to a model with fewer parameters. The resulting reduced-order model can predict temperature dynamics accurately and is identifiable throughout the cell’s lifespan which allows using the model’s parameters to monitor the state-of-health of the aging cell and detect anomalies in thermal behavior.

Li ion batteries↗
NASA NTRS1992

Monte Carlo modeling of atomic oxygen attack of polymers with protective coatings on LDEF

Characterization of the behavior of atomic oxygen interaction with materials on the Long Duration Exposure Facility (LDEF) will assist in understanding the mechanisms involved, and will lead to improved reliability in predicting in-space durability of materials based on ground laboratory testing. A computational simulation of atomic oxygen interaction with protected polymers was developed using Monte Carlo techniques. Through the use of assumed mechanistic behavior of atomic oxygen and results of both ground laboratory and LDEF data, a predictive Monte Carlo model was developed which simulates the oxidation processes that occur on polymers with applied protective coatings that have defects. The use of high atomic oxygen fluence-directed ram LDEF results has enabled mechanistic implications to be made by adjusting Monte Carlo modeling assumptions to match observed results based on scanning electron microscopy. Modeling assumptions, implications, and predictions are presented, along with comparison of observed ground laboratory and LDEF results.

Bruce A. Banks↗
DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗
DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Multiscale and multidimensional modeling of particle acceleration and transport in solar flares

Multi-messenger, multi-viewpoint, and time-resolved observations of solar flares are now providing unprecedented constraints on particle acceleration sites, energy conversion, and energy transport. The interpretation of current observations, including microwave imaging spectroscopy from EOVSA, hard x-ray (HXR) imaging from Solar Orbiter/STIX, gamma-ray diagnostics from Fermi, and in situ measurements from Parker Solar Probe and Solar Orbiter, collectively demands modeling frameworks that go beyond traditional spatially unresolved, one-zone models or single-mechanism descriptions. This review surveys multiscale and multidimensional modeling approaches, including kinetic, magnetohydrodynamic (MHD), and macroscopic particle models, that are being developed to meet the need. Kinetic simulations reveal that three-dimensional (3D) effects, including field-line chaos and self-generated turbulence, are essential for sustained power-law particle acceleration. MHD simulations now capture flux-rope eruptions, plasmoid-unstable current sheets, and turbulent flare regions in realistic magnetic topologies. Macroscopic models coupling MHD with energetic-particle models produce spatially resolved electron distributions and synthetic HXR and microwave emissions for direct comparison with observations. Despite these advances, outstanding challenges remain in bridging kinetic and global scales, improving MHD simulations and macroscopic particle models, and achieving quantitative model-observation closure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Modeling Occupant Core Temperatures Across the Boston Building Stock to Advance Public Health

In recent history, extreme heat has been the cause of most deaths from a natural disaster. Exposure to extreme heat can aggravate preexisting conditions, increase hospitalization, and even cause death. Thus, modeling the thermal resilience of households across the United States will allow for a quantitative assessment of the health and safety risks posed by extreme heat. For the first time, we simulate occupant comfort in representative households across Boston by combining the granular results of the ResStock(TM) model with a two-node heat strain model. This model calculates occupants' core body temperature in each simulated household over a year. We compare the simulated core temperatures to two public health metrics: hyperthermia and heat stroke. Our results show for households in Boston that do not have or use aid conditioning, thousands potentially experience many dangerous heat events each summer and these events can last for more than a day at a time. These heat events peak during the late afternoon and evening, just as residents are coming home, cooking meals, and trying to go to sleep. We have found that multi-family buildings, renters, and low-income homes in more airtight and insulated homes are the most at risk for these events. These results demonstrate the scale and urgency of exposure to extreme heat.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRSDate not supplied

Multiscale Modeling of Fracture Strength in Fibrous Thermal Protection System Materials

This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.

Jaehyun Cho↗
DOE OSTI2026

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗