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.

3,221 records · Page 30

Computational modeling of coupled mechanical damage and electrochemistry in ternary oxide composite electrodes

Performance degradation of ternary layered oxide cathodes largely originates from their loss of structural integrity in cyclic usage. Mechanical damage, such as intergranular fracture of the active particles, is not only a mechanical cleavage process but also interferes with electrochemical kinetics such as infiltration of liquid electrolyte, surface corrosion of the constituent primary particles, and may eventually isolate the primary grains from the electron conducting network. Here, in this work, we develop a computational framework that integrates electrochemistry of a LiNi x Mn y Co 1−x−y O 2 (NMC) composite cathode with mechanical damage of the active particles. To fully examine the intricate chemomechanical behavior of the electrode, we evaluate the effects of the anisotropic material properties, the influence of mechanical potential on Li transport, and the concurrent intergranular fracture and electrolyte penetration along the grain boundaries upon multiple cycles. Electrolyte infiltration benefits capacity retention but aggravates further mechanical damage by corrosion. Structural failure mostly occurs in the first charging due to the anisotropic mechanical strain between the primary grains, while the resulting damage remains stable in the later few cycles. The results are consistent with experimental observations and the integration of electrochemistry and mechanical failure enables a step further understanding of the complex mechanism of battery degradation.

Battery degradation

Spectroscopy in Nanoscopic Cavities: Models and Recent Experiments

The ability of nanophotonic cavities to confine and store light to nanoscale dimensions has important implications for enhancing molecular, excitonic, phononic, and plasmonic optical responses. Spectroscopic signatures of processes that are ordinarily exceedingly weak such as pure absorption and Raman scattering have been brought to the single-particle limit of detection, while new emergent polaritonic states of optical matter have been realized through coupling material and photonic cavity degrees of freedom across a wide range of experimentally accessible interaction strengths. In this review, we discuss both optical and electron beam spectroscopies of cavity-coupled material systems in weak, strong, and ultrastrong coupling regimes, providing a theoretical basis for understanding the physics inherent to each while highlighting recent experimental advances and exciting future directions.

Chemistry

A Parametric Battery Model for the Conceptual Design of Electric Aircraft

Electric aircraft conceptual designs usually have an assumed specific power and specific energy for the aircraft’s energy storage system that dictates the total capacity and the peak power available. The specific energy and power are inherently linked to the state-of-charge and discharge rate from which they were derived; therefore, they will not directly correspond to the diverse circumstances encountered in various aircraft missions and flight segments. These parameters turn the battery system into a black box, disregard potential electrical restrictions, and disallow the aircraft and battery to be optimized as a system. Peering into this box, this study highlights the importance of incorporating a parametric battery model into the conceptual design workflow by splitting high-level terms such as power into voltage and current and investigating their variability during the discharge process. Through the modeling of these more detailed parameters, this methodology shows the feasibility of using low states-of-charge for contingency operations, including the end of the reserve mission, expanding the amount of usable capacity for electric aircraft. Stark differences in aircraft capabilities can arise between varying fidelity battery models due to late-mission, high-power flight operations. This parametric battery model effectively captures these differences by evaluating the limitations that arise within the individual battery cells and the aircraft powertrain. This paper shows that the unusable charge of a battery is set by the balked landing power requirement and can realistically range from 5% to 44% based on assumptions. This sets the analog to unusable fuel capacity in aircraft with liquid fuel systems. Needing only aircraft- and mission-level inputs and only seconds of run time, this model is a prime fit for the fast, accurate exploration of the electric aircraft conceptual design space.

Battery

Test Program

A test program was set up early in the Explorer XIII program to simulate as nearly as possible all expected environmental conditions to which the spacecraft would be subjected during the ascending phase of the trajectory and while in orbit. These tests included spin, acceleration, thermal-vacuum, vibration, and shock. Functional tests were performed on individual modules of the telemeter system. The test procedure started with breadboard construction of telemeter components. Circuits were evaluated in this phase, and components chosen for prototype construction. Prior to the prototype environmental tests, it was found desirable to construct a mock-up structure for preliminary evaluation of the spacecraft hardware.

Explorer XIII Satellite

Elastic-plastic models for multi-site damage

This paper presents recent developments in advanced analysis methods for the computation of stress site damage. The method of solution is based on the p-version of the finite element method. Its implementation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experimental results are available for the panel. The plasticity analysis provided the same limit load as the experimentally determined load. The results of elastic-plastic analysis were compared with the results of linear elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of net-section yielding was determined also. The results show that crack growth rate is accelerated by the presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideally-plastic material model is compared against the Ramberg-Osgood model.

Actis, Ricardo L.

Mathematical Model of a Regenerative Fuel Cell for System Optimization

This thesis developed a system-level optimization model of a regenerative fuel cell (RFC) system for long-duration, off-world energy storage applications. Prior RFC design studies have typically been limited to reduced parameter sets and simplified constraints due to computational limitations relative to the number of relevant degrees of freedom. As a result, important nonlinear interactions between subsystems have not been fully captured. This work began to address that gap by developing a higher-fidelity, nonlinear optimization framework that incorporates a broader set of design variables and coupled constraints, enabling a multidimensional model that captures the coupled behavior of RFC subsystems and demonstrates the feasibility of applying optimization to such systems. An expanded system-level optimization approach was established that captures interactions between electrochemical performance, structural requirements, and storage design. This enabled a more comprehensive evaluation of trade-offs than conventional formulations. The model integrates four coupled subsystems: a fuel cell, an electrolyzer, reactant gas, and high-pressure storage tanks, and was formulated to accommodate a wide range of mission parameters, including operational time and required output power. It incorporates constraints on available solar array power, reactant mass balance between production and consumption, and pressure-dependent storage requirements. To enable reliable convergence, the optimization problem was reformulated to reduce dimensionality and improve numerical stability, with subsystem models organized for efficient evaluation. Problem dimensionality was reduced by consolidating lower-level design variables into higher-level representative quantities, and subsystem behavior was evaluated within the optimization loop. A multi-start initialization strategy was employed to mitigate sensitivity to local minima and improve solution quality, while nonlinear relationships were solved using robust numerical methods. The results showed that convergence was achieved across a range of required output power values. Specific energy reached a maximum at a critical mission power level, where the electrolyzer power matched the available solar input and operated near its voltage and current density limits. Beyond this point, further increases in required power resulted in less mass-efficient operation, increasing total system mass and reducing overall performance. The developed model represents an advancement in RFC system-level optimization by enabling analysis of a broader and more tightly coupled design space than previous considerations. While convergence behavior and computational cost remain challenges, the methods introduced improve solvability and allow inclusion of additional design variables with minimal loss of physical fidelity. However, the numerical results should not be interpreted as definitive design recommendations, as the model includes simplifying assumptions and omits several higher-order effects. Future work should extend this framework by incorporating additional subsystems and loss mechanisms, such as thermal management, parasitic power consumption, and reactant losses, to improve fidelity and ensure more representative design conclusions.

Electrochemistry

Computing the Critical Temperature of the Affine-Transformed $D=3$ Ising Model Using Masked Autoregressive Flow

The simple Ising model provides a rich environment to build and study lattice field theories. As part of an ongoing project to construct a conformal field theory (CFT) on an arbitrarily curved manifold, in this work we develop methods to measure the critical temperature $β_c$ of the affine-transformed Ising model on the face-centered cubic (FCC) lattice. The main challenge in this endeavor is finding a computationally efficient and accurate method of interpolating and extrapolating Monte Carlo observables with respect to coupling coefficients and temperature. Herein, we compare two such methods. A traditional statistical approach uses the multiple histogram (MH) method, while a newer machine learning approach uses a masked autoregressive flow (MAF) to estimate the underlying probability density function of a set of observables. While the MH method is specifically designed to interpolate and extrapolate Monte Carlo observables, we find that MAF is a viable alternative for measuring $β_c$ with a computational cost that scales more favorably. Furthermore, we comment on additional advantages of MAF relevant to our work, such as extrapolating in system volume.

Svenson, Kai [Texas U.]

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E

Modeling the Potential Impact of Storage on the US Power Sector in a Multisector Dynamic Context

The electric power sector is expected to grow in size, importance, and complexity around the world as economies expand and electric supply technologies and demand patterns evolve in significant ways. At the same time the electric power sector may see substantial increases in VRE generating technologies which can add variability and uncertainty to the diurnal and seasonal profile of electric supply. With these forces in play, the emergence of modular, flexible electricity storage technologies may have profound impacts on the operation of electric power systems. Here we reduce a storage modeling gap in MSD models by incorporating grid-based electricity storage into the electric sector dynamics of the Global Change Analysis Model-USA (GCAM-USA) with improved power sector representation. We find a potentially significant role for storage technologies in the future of the U.S. power system, with storage capacity ranging from 7.7 to 14.7 GW in 2050 and 13.9 to 31.6 GW in 2100 across several techno-economic scenarios. We also find that storage can help to smooth variability in residual load arising from evolving electricity demands and the introduction of high shares of VRE to the electric grid. This reduces reliance on high-cost peaking generators, improves system-wide capacity factors, and limits curtailment of VRE.

Patel, Pralit L.

Earth Observations and Integrative Models in Support of Food and Water Security

Global food production depends upon many factors that Earth observing satellites routinely measure about water, energy, weather, and ecosystems. Increasingly sophisticated, publicly-available satellite data products can improve efficiencies in resource management and provide earlier indication of environmental disruption. Satellite remote sensing provides a consistent, long-term record that can be used effectively to detect large-scale features over time, such as a developing drought. Accuracy and capabilities have increased along with the range of Earth observations and derived products that can support food security decisions with actionable information. This paper highlights major capabilities facilitated by satellite observations and physical models that have been developed and validated using remotely-sensed observations. Although we primarily focus on variables relevant to agriculture, we also include a brief description of the growing use of Earth observations in support of aquaculture and fisheries.

Water Resources

A quality-agnostic combinatoric cost estimation model for large-format directed energy deposition metal additive manufacturing

Directed energy deposition (DED) additive manufacturing (AM) processes are amenable to synergistic combination into multi-process AM systems due to similar requirements for automation and energy sources. This work analyzes the economic performance of such DED AM systems from a quality-agnostic combinatoric standpoint with a model that calculates lowest-cost system combinations based on part geometry and process performance metrics. Common DED AM systems research focuses on a single process and does not consider the process, system, and application in the context of all possible system combinations (e.g., the combined set of process selection(s), motion system(s), and process hardware), leading to limited applicability of the resulting DED AM systems to cost-sensitive components such as those found in energy generation applications. The model developed herein incorporates the capital, material, and energy costs associated with DED AM system combinations into a predictive tool for estimating part and system cost, the output of which is intended to guide deployment of finite research and development resources towards DED AM system combinations with the lowest costs and greatest likelihood of economic impact. The DED AM systems identified by this framework may enable domestic production of the large conventionally cast and forged components necessary for energy generation.

Shanafield, Alexandra [ORNL]

Impact of Crystalline Phases on Low-Activity Waste Glass Durability: Insights from PCT and VHT

During vitrification of nuclear wastes, slow cooling along the container centerline promotes crystalline phase formation, which can alter residual glass composition and reduce chemical durability. This study investigates the effects of crystalline phases on the chemical durability of low-activity waste (LAW) borosilicate glasses using the product consistency test (PCT) and vapor hydration test (VHT) on container centerline cooled (CCC) samples. A preliminary model (R2 = 0.88) was developed to predict CCC PCT responses based on glass composition, PCT data from quenched glasses, and measured crystal fractions. Using the latest LAW glass dataset, the feasibility of predictive modeling is evaluated, limitations in current data and methods are identified, and challenges for improving model accuracy are discussed to guide future data collection and model development.

borosilicate glass

Predicting Operational Performance of xEMU Boot at Lunar South Pole Temperatures using Thermal Desktop ®

The spacesuit boots that will be used on Artemis lunar south pole surface missions will be exposed to extremely cold temperatures (down to ~50 K). To assess the performance of the government’s Exploration Extravehicular Mobility Unit (xEMU) lunar boot in these permanently shadowed regions, testing was performed at the Jet Propulsion Lab (JPL) in the Cryogenic Ice Transfer, Acquisition Development, and Excavation Laboratory (CITADEL) thermal vacuum (TVAC) chamber. This paper documents the data analysis, thermal boot model correlation, and operational predictions conducted using data from the xEMU CITADEL TVAC test. Expected thermal conductances within the boot and between the boot and environment were calculated from test data, which was then used as an initial guess for conductances within a Thermal Desktop (TD) model. Correlation of the TD model using the internal SOLVER feature was performed across 10 different test points which varied external temperature, internal boot ventilation flowrate, and contact pressure. Operational performance at the lunar south pole was then predicted using results from the correlated model. While the predictions provide evidence for acceptable performance of the boots at the 100K environment test point, there is still substantial uncertainty in performance, especially at the 48K test point. This uncertainty is due in part to testing limitations such as contacting the foot to a hard metal plate rather than granular regolith, and model limitations such as the lack of a realistic foot model. These limitations and their impacts are addressed in detail in this paper. The results of this test series and model correlation underscore the importance of additional improved testing and modeling for characterizing the expected thermal resistance between the outside of the boot and the lunar surface.

Spacesuit

High-resolution modeling of indoor radon exposure with uncertainty quantification in Utah

Indoor radon accounts for 37% of population-level exposure to ionizing radiation in the United States. However, radon metrics are typically reported at coarse spatial scales, potentially obscuring meaningful local variation. We developed a high-resolution modeling framework to estimate indoor radon concentrations across Utah while explicitly quantifying predictive uncertainty. A total of 19,497 residential radon measurements collected between 2006 and 2017 were combined with environmental and housing characteristics and analyzed using a geospatial neural network that accommodates spatial dependence and nonlinear associations. Predictions were generated on a uniform hexagonal grid at 0.73 km2 resolution (H3 level 8). Out-of-sample predictions aggregated to the H3 level 8 grid showed good agreement with observed concentrations (Pearson r=0.64), while household-level predictions exhibited more moderate agreement (r=0.45). The model produced well-calibrated uncertainty estimates, with 24.1% of held-out observations exceeding the predicted 75th-percentile threshold. Maps of predicted radon concentrations and the probability of exceeding the U.S. EPA action level of 148 Bq/m3 (4 pCi/L) revealed substantial fine-scale spatial heterogeneity that was not apparent in conventional coarse-resolution summaries, with greater local variability observed in densely monitored urban counties than in sparsely sampled regions. High-resolution radon models that explicitly quantify uncertainty provide a useful framework for characterizing the spatial distribution of indoor radon and identifying areas of elevated exceedance risk. These findings highlight the value of fine-scale monitoring data and uncertainty-aware modeling approaches for radon exposure assessment, environmental risk characterization, and radon-related health research.

Wu, Yunhan [ORNL] (ORCID:0000000178842994)

Shakeup and shakeoff spectra in the electron capture decay of atomic 7 Be

The most stringent laboratory-based experimental limits on the existence of submegaelectronvolt sterile neutrinos are currently set by decay spectroscopy of radioactive 7 Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shakeup and shakeoff spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wave functions ab initio and compute all single and double shake processes in the electron capture decay of atomic 7 Be . The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the 7 Be is embedded into. The new models also show that the 𝐿/𝐾 electron capture ratio of 7 Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Developing a Pyrolysis Gas Thermal Blocking Model for Reentry Demise

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These heating coefficients were derived for typical metallic materials that melt and do not have a strong gas-phase contribution to the flow in the boundary layer. Modern satellites typically feature fiber-reinforced polymer (FRP) components, such as solar array booms, facesheets of sandwich panels, or overwraps for composite-overwrapped pressure vessels (COPV). These FRP materials do not behave the same as metals in the reentry environment, but instead will pyrolyze and develop significant volumes of gas into the boundary layer. Accurately predicting the reentry demise of FRP components is critical to assessing the reentry casualty risk for modern spacecraft. Research in recent years has shown that this demisability can depend heavily on how the expulsion of gaseous pyrolysis products through the outer surface of the material affects the heat flux at the surface. The ODPO has been developing a reduced-order model of the effect of pyrolysis gas blowing on the heat flux based on correlations between a blowing factor and a non-dimensional heat flux to be incorporated in the upcoming version 7.3 of the Object Reentry Survivability Analysis Tool (ORSAT). This presentation discusses the progress of this development project and the challenges remaining for generalizing the model across families of FRP materials.

Benton Greene

SERFE Ground Unit EVA Series After Three Year Spacesuit Stowage Period

NASA’s spacesuit government reference design for returning to the Moon is called the Exploration Extravehicular Mobility Unit (xEMU). The xEMU subassembly that provides life support, such as oxygen and thermal control, is the Portable Life Support System (PLSS). Inside the PLSS is a new technology that NASA wanted to test to provide cooling to the crew during EVAs (ExtraVehicular Activity). This technology is called the Spacesuit Water Membrane Evaporator (SWME). In order to test SWME in an improved thermal control loop (TCL) both on Earth and in Space, the Spacesuit Evaporation Rejection Flight Experiment (SERFE) was created. The Ground unit, or testbed at Johnson Space Center (JSC), tested the cooling technology in Earth’s gravity, while the Flight unit, or payload on the International Space Station (ISS), tested the cooling technology in micro-gravity. Since fluids flow differently in micro-gravity, testing in both environments would provide important data for improving cooling performance during EVAs. Both units completed 25 simulated EVAs with the same settings so SWME performance on the ground could be compared to the ISS. The Flight unit was completed first and performed EVAs on the ISS between 2020 and 2022. The Ground unit performed EVAs between 2021 and 2022. When the Flight unit came back from the ISS, it was taken apart for analysis. This analysis looked at how well SWME was able to maintain its heat rejection capability after various dwell times, such as a 90 day Airlock Coolant Loop Recovery (ALCLR) cycle, when the Extravehicular Mobility Unit (EMU) currently on the ISS would normally need maintenance. After a three year simulated spacesuit dwell, the Ground unit performed another EVA series in 2025 to test SWME’s shelf life. The results from this test series will inform mission planning as NASA plans to go back to the Moon and beyond.

Michael Lewandowski

Prime Time for Model-Predictive Control? Assessing the Technical and Market Readiness of Advanced Controls in Buildings

Despite three decades of extensive research and field testing that have consistently validated the benefits of Model Predictive Control (MPC) in building applications, the technology has seen limited market adoption. This paper evaluates the readiness of MPC for widespread deployment, showcases recent demonstrations and field tests across diverse building types, including residential, small commercial, large commercial, and campus settings. Our results demonstrate that MPC can optimize system operations to achieve load shifting, minimize curtailment of on-site generation, and reduce energy costs by up to 80 %, while maintaining or improving occupant comfort. We also show that MPC can effectively control large assets, such as MW-sized thermal storage systems, and respond to dynamic pricing signals. However, achieving scale remains difficult due to labor-intensive workflows, reliance on a “PhD-in-the-loop” for MPC design and maintenance, susceptibility to fragile data infrastructure, and persistent workforce education and acceptance barriers. To bridge this gap, we outline a transition from bespoke, labor intensive prototypes toward streamlined, segment-targeted deployment strategies that leverage model templates, semantic tools, and generative AI. By automating control configuration and reducing engineering effort, these recommendations provide a pathway for transforming successful research demonstrations into scalable, market ready solutions for MPC-based controls.

Pritoni, Marco