Engineering PapersSearch

SEARCH · Engineering Papers

Results for “model validation”

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.

453 records · Page 11

A Probabilistic Model of a Porous Heat Exchanger

This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

O P Agrawal

The Contaminant Footprint of Landed Spacecraft: Toward an Inventory and Modelling Framework

All spacecraft generate and carry contaminants, i.e., unwanted and potentially harmful material. When a spacecraft lands and operates in vacuum, as onto Earth’s Moon, it introduces contaminants into its environment that may compromise mission science objectives and engineering performance. Contamination may degrade sites of unique value to planetary science or in situ resource utilization. This presentation will identify and compare source terms and transport vectors for contaminants – in particular, organic material – generated by landed spacecraft. An integrated modeling framework for the organic contamination footprint of spacecraft missions will be described and presented.

Gas Dynamics

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

Advances in Entry Systems Modeling: Modeling Summer Visit 2025 Technical Reports

This NASA Technical Memorandum presents the collective work carried out during the Modeling Summer Visit (MSV) 2025, a collaborative research program between NASA Ames Research Center and three leading European research institutions: the University of Bordeaux, the von Karman Institute (VKI), and Ecole Polytechnique. Held from July 7th to August 1st, 2025, the NASA ARC MSV program brought together 42 researchers from Europe to work alongside NASA staff on 22 projects in entry systems modeling. The program was organized by Bruno Dias, Sergio Fraile-Izquierdo, and Jeremie Meurisse from Analytical Mechanics Associates, Inc., in collaboration with Marc Massot (Ecole Polytechnique), Jean Lachaud (University of Bordeaux), and Thierry Magin (VKI). The MSV program traces its origins to an initiative started by Nagi N. Mansour in 2018, with this 2025 edition representing a consolidation and expansion of that original vision. The projects were organized into six thematic groups: Porous Media, Numerical Methods, Plasma/Rarefied Methods, Experimental Reconstruction, Flow and Material Coupling Methods, and Transition and High-Speed Flows. Throughout the four-week program, participants presented their work in three public presentations (July 7th, July 21st, and August 1st), allowing for continuous feedback and discussion of their evolving research.

Plasmas

Technology Assessment and Modelling of Three Terminal Tandem Solar Cells for In-Space Utilization

NASA seeks a sustained human presence on the lunar surface. High-efficiency and resilient energy systems are critical to supporting habitats, scientific outposts, and lunar surface operations. Photovoltaics have long been the backbone of space power systems, transitioning from silicon solar cells to multi-junction III-V solar cells, which now dominate state-of-the-art (SOA) technology. Multi-junction cells achieve high efficiency by layering semiconductors, each absorbing a specific portion of the solar spectrum. However, efficiency gains are increasingly constrained by device physics and manufacturing complexity. Additionally, higher-order junctions pose challenges for accurate modelling and performance characterization due to difficulties in replicating the AM0 spectrum.

solar cells

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

Cryogenic Propellant Tank Drain Simulation with a Stratified Ullage Nodal Model

Cryogenic propellant tank pressurization systems maintain tank pressure within specified limits to comply with propellant thermodynamic and tank structural design requirements. Many pump-fed cryogenic liquid propulsion systems control tank pressure using an autogenous pressurization system. Autogenous pressurization systems vaporize and heat propellant to produce the pressurant gas delivered to the tank. While autogenous pressurization systems eliminate the need for a separate pressurant gas storage and delivery system, they couple pressurant mass requirements to propellant mass requirements. Accurate prediction of autogenous pressurant mass during propellant drain is essential in optimizing cryogenic propellant tank designs.

Richard G Hibbs

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG

Exploring the Model Design Space for Battery Health Management

Battery Health Management (BHM) is a core enabling technology for the success and widespread adoption of the emerging electric vehicles of today. Although battery chemistries have been studied in detail in literature, an accurate run-time battery life prediction algorithm has eluded us. Current reliability-based techniques are insufficient to manage the use of such batteries when they are an active power source with frequently varying loads in uncertain environments. The amount of usable charge of a battery for a given discharge profile is not only dependent on the starting state-of-charge (SOC), but also other factors like battery health and the discharge or load profile imposed. This paper presents a Particle Filter (PF) based BHM framework with plug-and-play modules for battery models and uncertainty management. The batteries are modeled at three different levels of granularity with associated uncertainty distributions, encoding the basic electrochemical processes of a Lithium-polymer battery. The effects of different choices in the model design space are explored in the context of prediction performance in an electric unmanned aerial vehicle (UAV) application with emulated flight profiles.

Saha, Bhaskar

High-Fidelity Modeling Approach for Coupling Thermal, Fluid, and Neutronic Analysis in an NTP Engine

Nuclear thermal propulsion (NTP), with its increased efficiency compared to traditional rocket engines, is a potentially enabling technology for space missions of higher complexity, such as delivering heavy machinery or crew to the moon or other planets. Thus, exploring modeling approaches that can effectively couple the thermal, fluid, and neutronic behaviors of NTP reactors is essential to mitigate cost and schedule challenges with hardware development. This paper details one such approach which utilizes Cardinal to unify reactor physics analysis performed using OpenMC, heat transfer analysis performed using the MOOSE framework, and thermal hydraulic analysis performed using the MOOSE Thermal Hydraulic Module. This approach is demonstrated with a proof of concept that focuses on a legacy reactor geometry from the historic NERVA program. This preliminary model found an overprediction in, when compared to the reference data, axial power peaking and the temperatures of both solids and fluids by roughly 20%. It is clear, however, that the model conserved a number of broad trends, and a few model refinements should reduce the observed error. This approach has definite merit and should enable the production of informative results for more modern NTP systems. This work is a first step in creating a transient reactor model that can be updated across a sequence of time steps based on continually fluctuating, time-dependent boundary conditions.

Nuclear Thermal Propulsion

Dynamic Regulation of Sub-Atmospheric Pressure for Constant and Cyclic Gas Loads During Testing of Spacesuit Components

New iterations of the various subsystems within the spacesuit will benefit from a new sorbent technology. For example, sorbents used in the Trace Contaminant Control (TCC) and the Rapid Cycle Amine (RCA) systems within the Exploration Portable Life Support System (xPLSS), part of the Extra-Vehicular Mobility Unit (xEMU). For proper validation, the sub-atmospheric pressure needs to be maintained within a simulated 2 ft 3 spacesuit volume. The traditional approach, which involves placing the entire vent loop within a hypobaric chamber, is not practical for the widespread testing of components and prototypes, as these specialized chambers are costly and not widely available. A two-stage regulator based sub-atmospheric pressure system is presented in this work. The method regulates the pressure based on pressure differentials between the pressure regulation pump, the test system, and the ambient environment. The performance of this system is demonstrated using data collected during a 150+ hour non-regenerative TCC sorbent evaluation, several 24-hour regenerative TCC sorbent evaluations with different regeneration cycles, and an 8-hour RCA sorbent evaluation, all under xEMU operating conditions. For regulation, the inlet regulator was set to propagate a small leak to increase the stability of the system as the gas load from the testing system changes. It was also found that controlling the input flow rate that replaces the material lost during regeneration is critical for ensuring the stability of the system. This system provided excellent pressure regulation, without adjustments, for constant loads during shorter time evaluations (hours), while longer time evaluations (days) are easily obtainable with periodic regulator adjustments. A second method is under development to address longer-term stability and to automate pressure regulation. It will incorporate flow and pressure measurements to dynamically control a set of electronically controlled proportional valves for adjusting the pumping speed and air bleed supply.

Nicholas F Materer

A review of nickel hydrogen battery technology

This paper on nickel hydrogen batteries is an overview of the various nickel hydrogen battery design options, technical accomplishments, validation test results and trends. There is more than one nickel hydrogen battery design, each having its advantage for specific applications. The major battery designs are individual pressure vessel (IPV), common pressure vessel (CPV), bipolar and low pressure metal hydride. State-of-the-art (SOA) nickel hydrogen batteries are replacing nickel cadmium batteries in almost all geosynchronous orbit (GEO) applications requiring power above 1 kW. However, for the more severe low earth orbit (LEO) applications (greater than 30,000 cycles), the current cycle life of 4000 to 10,000 cycles at 60 percent DOD should be improved. A NASA Lewis Research Center innovative advanced design IPV nickel hydrogen cell led to a breakthrough in cycle life enabling LEO applications at deep depths of discharge (DOD). A trend for some future satellites is to increase the power level to greater than 6 kW. Another trend is to decrease the power to less than 1 kW for small low cost satellites. Hence, the challenge is to reduce battery mass, volume and cost. A key is to develop a light weight nickel electrode and alternate battery designs. A common pressure vessel (CPV) nickel hydrogen battery is emerging as a viable alternative to the IPV design. It has the advantage of reduced mass, volume and manufacturing costs. A 10 Ah CPV battery has successfully provided power on the relatively short lived Clementine Spacecraft. A bipolar nickel hydrogen battery design has been demonstrated (15,000 LEO cycles, 40 percent DOD). The advantage is also a significant reduction in volume, a modest reduction in mass, and like most bipolar designs, features a high pulse power capability. A low pressure aerospace nickel metal hydride battery cell has been developed and is on the market. It is a prismatic design which has the advantage of a significant reduction in volume and a reduction in manufacturing cost.

Smithrick, John J.

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES

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

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.