Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Enthalpy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd and associated criterion (KnGLL > 0.05) seems conservative in predicting zones of breakdown in the flow, and associated errors for surface quantities. Updated criteria of KnGLL = 2.0 and 0.5 appear to be more appropriate for surface and flow quantities, respectively. Furthermore, when studying reacting flows, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown is minimal.

DSMC↗

Particle-Based Fiber Models of Woven Materials for Earth Entry Thermal Protection

Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing and deploying woven material heat shields for missions, including Artemis I (3D-MAT for compression pads) and potentially Mars Sample Return - Earth Entry System. These materials are complex, hierarchical and must protect from extreme environments and phenomena, such as deformation, impact and high-enthalpy heating. Woven material performance depends on microstructure, damage and weave geometry. Therefore, fiber-specific models are needed to simulate fiber contacts within the weave hierarchical geometry (fiber to tow to yarn to weave) and the inherent directionality of fibers. Explicit fiber models can simulate how weave microstructure evolution affects thermal and mechanical properties. We parameterize a discrete element bonded particle models (DEM-BPM) of fibers to capture thermal and mechanical behavior within and between fibers, with bonded and contact forces, respectively. We study the proportion of heat transfer and stress via the contact network, fiber bonds and the weave geometry, for example, with respect to yarn warp-weft identity (whether it interlaces weave layers). Our results demonstrate the importance of explicit fiber modeling for connecting microstructure with thermal and mechanical properties.

woven material↗

Development of a Composite Material Aerodynamic Demise Model for the Object Reentry Survival Analysis Tool (ORSAT)

The reentry demise of fiber-reinforced polymer (FRP) composites is an increasing concern for modern spacecraft at end of life. Unlike traditional materials such as metals, shredding of the material by aerodynamic forces appears to be a major component of the reentry demise mechanism. This paper will describe a new mechanical, strength-based, material demise model for the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1. The model is based on laboratory in-situ and residual strength tests of several FRP materials performed during the NASA Orbital Debris Program Office’s Phase II Composite Material Demise test campaign. To test the residual strength of partially charred FRP materials, the authors tested 107 rectangular shaped samples of different thicknesses that were previously exposed to high-enthalpy flow at the University of Texas at Austin’s Inductively Coupled Plasma Torch facility. Tests were performed under normal atmospheric conditions at NASA Johnson Space Center’s Experimental Impact Laboratory using a Chatillon TCD1000 tensile test machine configured with a three-point bending jig. Mean fracture load, residual strength, delamination, and flexural modulus were measured for each test. An engineering model of aeromechanical demise for each material was developed by correlating the measured residual strength of the samples with the duration and magnitude of the applied heat flux and the mass loss and char progression. This model has been implemented in ORSAT 7.1 for the built-in, charring carbon fiber/epoxy and charring glass fiber/epoxy material models. The new model was verified by calculating the demise of hundreds of FRP composite fragments of various shapes and sizes and checking the calculations for consistency, with specific focus on material characterization in one-dimension. Further tests are planned in hypersonic flow facilities to validate the assumptions used in the model.

Priscilla A. Mendoza↗

Irrigated Agriculture Significantly Modifies Seasonal Boundary Layer Atmosphere and Lower Tropospheric Convective Environment

Modification of grasslands into irrigated and nonirrigated agriculture in the Great Plains resulted in significant impacts on weather and climate. However, there has been lack of observational data–based studies solely focused on impacts of irrigation on the PBL and convective conditions. The Great Plains Irrigation Experiment (GRAINEX), conducted during the 2018 growing season, collected data over irrigated and nonirrigated land uses over Nebraska to understand these impacts. Specifically, the objective was to determine whether the impacts of irrigation are sustained throughout the growing season. The data analyzed include latent and sensible heat flux, air temperature, dewpoint temperature, equivalent temperature (moist enthalpy), PBL height, lifting condensation level (LCL), level of free convection (LFC), and PBL mixing ratio. Results show increased partitioning of energy into latent heat relative to sensible heat over irrigated areas while average maximum air temperature was decreased and dewpoint temperature was increased from the early to peak growing season. Radiosonde data suggest reduced planetary boundary layer (PBL) heights at all launch sites from the early to peak growing season. However, reduction of PBL height was much greater over irrigated areas than over nonirrigated croplands. Relative to the early growing period, LCL and LFC heights were also lower during the peak growing period over irrigated areas. Results note, for the first time, that the impacts of irrigation on PBL evolution and convective environment can be sustained throughout the growing season and regardless of background atmospheric conditions. These are important findings and applicable to other irrigated areas in the world.

Emilee Lachenmeier↗

On Characterization of Flow Disturbances in Arc-Jet Testing

This paper reports computational simulations and analysis of flow characterization tests in high enthalpy arc-jet facilities at NASA Ames Research Center. Flow disturbances and their characterization are presented through case studies, and their implications for testing are discussed. The case studies include free-jet test configurations from three different arc heaters and nozzles: the 60-MW IHF 30-inch conical nozzle, the 10-MW TP3 15-inch conical nozzle, and the 20-MW AHF 12-inch conical nozzle. For all cases, test articles are placed in the jet exiting the conical nozzle, and the existence of a flow disturbance is confirmed through flow survey data using pitot pressure and heat flux probes and accompanying analysis. The paper focuses on cases where the effects of the disturbances on test article surface quantities are important. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzles and test box, including the models tested. Comparisons of computations with the experimental measurements are presented. The computations that reproduce the probe sweep data approximately are essential to interpret the arc-jet test data accurately, while providing insights into several observed calorimeter anomalies caused by flow disturbances.

Arc-jet testing↗

Active Region Simulation with EBTEL

Active regions (ARs) are areas in the sun’s upper atmosphere that are bright in the extreme ultraviolet (EUV) and X-ray spectrum. Coronal loops are one key feature of an AR. They can be described as arcs of plasma resulting from loop-shaped magnetic fields on the sun. The average temperature of these loops is over one million kelvin, which is significantly hotter than the sun’s surface. The mechanism for the heat transfer of these high temperatures remains unknown and is referred to as the coronal heating problem. One leading hypothesis regarding this mystery is that small, impulsive heating events called nanoflares are a major contributor. Here, we explore the contribution of nanoflares in the heating of AR plasma. We utilize the Enthalpy-Based Thermal Evolution of Loops (EBTEL) program to simulate NOAA 12846, as observed on July 25, 2021. We then use that simulation to analyze how the frequency of nanoflare heating events affects the EUV and X-ray observations in existing instruments, including SDO’s AIA and Hinode’s XRT. We also predict how this AR would appear in recently developed X-ray instruments, such as NASA’s MaGIXS.

Active region↗

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗

On Characterization of Flow Disturbances in Arc-Jet Testing

This paper reports computational simulations and analysis of flow characterization tests in high enthalpy arc-jet facilities at NASA Ames Research Center. Flow disturbances and their characterization are presented through case studies, and their implications for testing are discussed. The case studies include free-jet test configurations from three different arc heaters and nozzles: the 60-MW IHF 30-inch conical nozzle, the 10-MW TP3 15-inch conical nozzle, and the 20-MW AHF 12-inch conical nozzle. For all cases, test articles are placed in the jet exiting the conical nozzle, and the existence of a flow disturbance is confirmed through flow survey data using pitot pressure and heat flux probes and accompanying analysis. The paper focuses on cases where the effects of the disturbances on test article surface quantities are important. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzles and test box, including the models tested. Comparisons of computations with the experimental measurements are presented. The computations that reproduce the probe sweep data approximately are essential to interpret the arc-jet test data accurately, while providing insights into several observed calorimeter anomalies caused by flow disturbances.

Arc-jet testing↗

A Summary of the 2023 Joint High-Speed Boundary Layer Transition Workshop

A joint workshop between the Air Force Office of Scientific Research, the Office of Naval Research, and the NASA Hypersonic Technology Project was held in June 2023 on the topic of high-speed boundary-layer transition. The purpose of this workshop was to gather subject matter experts in the areas of supersonic and hypersonic boundary-layer transition in the areas of geometries of interest, computational methods, facilities and diagnostics, flight experiments, flow control, and new approaches and to discuss the current state of the art in the field. As gaps and accomplishments were identified, roadmaps were also created to guide future investments in these areas. The top recommendations are (1) to define relevant common geometries of interest that can be worked by multiple research groups; (2) improve modeling of high-enthalpy flows involving chemistry and ablation, studies of receptivity of freestream disturbances and atmospheric particulates into the boundary layer, surface feature characterization and modeling of changes to geometry; (3) design and determine viable flow control methods; and (4) move towards incorporating boundary layer transition prediction into a standard practice for design of high-speed vehicles. To execute this last recommendation successfully, investments beyond fundamental (basic and early applied) research are needed.

boundary layer transition↗

BLAYER User Guide

Estimation of boundary layer edge or thickness is an essential aspect of utilizing high fidelity aerothermodynamic flow solutions. For viscous flow over a flat plate, the boundary layer thickness is usually defined as the distance from the wall at which the flow velocity (zero at the wall) has increased through some profile to 99% or 99.5% of the free stream velocity. In real gas flows over aerospace vehicles, a total enthalpy profile is employed for edge calculations. The importance of knowing boundary layer thickness and its relation to possible transition from laminar to turbulent flow is illustrated by Space Shuttle mission STS-114, when the appearance of tile gap filler protrusions on the nose was deemed serious enough that a crew member was called upon to exit the vehicle and remove the potential threat prior to its descent from orbit. The BLAYER utility used as part of post-processing Orion/MPCV flow solutions is a generalization of the earlier BLAYER_RESULTS utility developed following the STS-107 Columbia accident. BLAYER handles a variable number of gas species (not just 5), allows for more than one temperature and for optional extra flow quantities, and omits optional handling of Orbiter tile datasets.

momentum thicknes↗

Towards Quantifying and Correlating RTV Intumescence to Entry-Relevant Conditions

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also adopted to bond instrumentation plugs of temperature and pressure sensors into heatshield’s tiles. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Experimental data has also shown heating rate-dependent swelling (or intumescence) and shrinking of RTV, these volume changes may cause roughness-induced boundary layer transition. Outgassing and surface oxide formation of RTV upon decomposition also occur, which can be sources of contamination of heat shield sensors. This combination of detrimental effects motivates a critical need to develop a high-fidelity model for RTV ablation. As a first critical step in RTV modeling, a comprehensive material properties database for RTV was collected that account for the lack of data in properties such as pyrolysis char yield, shape change, virgin and char porosity. The next step is to quantify the volume change of RTV as a function of temperature. While this was done in previous campaigns, limited data was collected which showed a possible heating rate dependent intumescent behavior. To further understand the intumescence phenomenon of RTV, numerous dedicated experiments were performed at Beamline 8.3.2 of the Advanced Light Source, where RTV samples were heated while X-ray scans were taken in situ. Micro-Computed Tomography scans were taken in situ for low heating rates (≤ 60 °C/min), whereas continuous radiography scans were taken for higher heating rates (≤ 1500 ºC/min). Unconstrained samples and samples constrained in different holders (quartz, graphite, FiberForm and PICA) were tested to determine the change in intumescence due to confinement. Using deep learning techniques, the tomographies and radiographies are segmented to obtain volume change, porosity, pore size distributions and connectivity. From the results obtained, modifications to the previous RTV intumescence model are proposed, along with preliminary comparisons to testing RTV at high-enthalpy facilities.

RTV↗

An Overview of Nuclear Thermal Rocket Element Environmental Simulator (NTREES) Capabilities, Upgrades, and Overlaps in Advanced Material Testing Areas

The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) was designed to perform non-nuclear testing of nuclear thermal rocket (NTR) fuel elements and materials. NTREES can simulate the thermal hydraulic environment within the reactor of an NTR and has been leveraged to provide data on the thermochemical, thermomechanical and thermohydraulic performance of components and materials. Active upgrades to the NTREES facility include installing a DC power supply for DC powered joule heating, a supply gas chiller for a cryogenically cool gas supply, gas pre-heaters with sample holders for testing samples too small to be inductively heated by bathing them in hot gas, and a high temperature digital image correlation (DIC) capability for measuring high temperature strains real time. In addition to NTR testing, given the commonality of high enthalpy flow environments between the reactor of an NTP engine and the aerothermal heating of high-speed flight, the unique capabilities of NTREES to provide high temperature gas flows at high pressure without combustion or electrode spallation byproducts may lend itself to additional usage of the facility beyond the original intent.

Michael P Schoenfeld↗

Simulations of Large-Scale Zero Boiloff, Densification, and Solidification of Hydrogen

An Integrated Refrigeration and Storage (IRAS) experimental system called the Ground Operations and Demonstration Unit for Liquid Hydrogen (GODU-LH2) demonstrated the ability store cryogens in a zero boiloff (ZBO) process, and to densify liquid hydrogen by reducing the temperature and pressure down to the triple point resulting in solidification. The incompressible Navier-Stokes was utilized to simulate ZBO, densification, and solidification of liquid hydrogen (LH2) in the IRAS tank. The simulations were performed using a commercially available Computational Fluid Dynamics (CFD) pressure-based mass and momentum flow model and an enthalpy-porosity energy model. Results demonstrated the simulation’s ability to predict time-dependent flow and temperature fields and solid-liquid phase locations for hydrogen during ZBO, densification, and solidification. The simulations showed good agreement with experimental data, with errors within acceptable ranges for temperature and pressure predictions, provided detailed insights into natural convection and solid-liquid phase change dynamics. These findings are critical for the design of future cryogenic fluid management systems.

Computational Fluid Dynamics↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials↗

Thermodynamic Properties of Uranium-Zirconium Carbide Solid Solutions with Oxygen

Uranium-zirconium carbide solid solutions are potential candidates for fuel in space nuclear reactors. To support development of advanced robust fuels for space applications, (U y Zr 1-y )C x O z solid solutions were synthesized via carbothermic reduction with uranium, carbon, and oxygen contents raging from y=0.05 – 0.3, x=0.95 – 0.98, and z = 0.02 – 0.07. These solid solutions were subsequently characterized by X-ray diffraction and elemental combustion analysis. High-temperature oxidative drop solution calorimetry in molten sodium molybdate solvent at 1073 K was employed to determine the formation, oxidation, and mixing enthalpies of (U y Zr 1-y )C x O z solid solutions, and thereby for evaluating their thermodynamic stability. The solid solutions are enthalpically less stable than the UC and ZrC end members. They also become less stable and more reactive towards oxidation with increasing uranium content. Incorporation of oxygen into solid solutions increases their solid solution-stability and improves their oxidation resistance.

Calorimetry↗

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.

sting arm↗