Search NASA⌕ Search

SEARCH · Search NASA

Results for “General Circulation Model”

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.

1,378 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Solar Heat for Industrial Processes: Integration with Chemical Reactors

The integration of solar thermal systems with chemical reactors has been proposed as part of a larger effort to develop and deploy solar heat for industrial processes (SHIP) technologies. A strong motivation for SHIP processes and technologies is the potential for high thermal efficiency coupled with low-cost thermal energy storage (TES) which can enable commercial deployment of such systems. While there are different ways to categorize SHIP technologies, one important such distinction is between directly irradiated systems and indirect off-sun process driven by a SHIP system. While directly irradiated systems can provide high thermal efficiencies and high fluxes, they usually require complex engineering solutions due to the need for redesigning the established processes and unit operations. In most cases, it is also more challenging to couple such a process to a TES system, losing some of the benefits of SHIP. On the other hand, using a SHIP system to drive an industrial process off-sun can allow better integration with existing process chains, easier TES capabilities, and potential for more applications fitting a specific SHIP technology. However, the integration of SHIP systems with the industrial processes is not fully explored in detail, especially in the case of high-temperature processes such as reforming, cracking, cement manufacturing, and iron/steelmaking. Many of these systems require heating fluxes of >50 kW/m^2, supplied via combustion of hydrocarbons in a fire box and benefitting from radiative heat transfer between the flue gases and the reaction zones. As such, using SHIP systems for such processes is more complex than providing the same thermal input in the form of a heat transfer medium (HTM) entering the reactor, kiln, or furnace. Moreover, in case convective heat transfer using SHIP is envisioned, for example using supercritical CO2 as the HTM from a particle receiver, the thermal integration might be more challenging than initially envisioned: lower heat transfer coefficients and limited approach temperature might require large flow rates, causing a mismatch between the process thermal requirements and the thermal capacity of the SHIP system. In addition, even if the heat exchange between SHIP and reactor is effective, there is still a cold leg HTM at the reaction temperature or slightly below it. Chemical plants usually include a set of heat exchangers, heat recovery steam generators, and even power generation units - in a tightly integrated design - to recover the flue gases which are eventually vented. With SHIP systems mostly operating on a closed HTM loop, bottoming the cold leg is crucial. In this talk, we will present different modeling results for a variety of syngas production reactions, using catalytic and chemical looping processes, and discuss some of the challenges and design considerations for off-sun chemical reactors using SHIP systems.

14 SOLAR ENERGY↗
NASA NTRSDate not supplied

Integrated Process-Structure-Property Simulations for Additive Manufacturing Using the Open-Source Materialite Package

The microstructure and properties of additively manufactured (AM) metals are strongly dependent on process conditions. Therefore, process-structure-property (PSP) simulations are a useful tool for exploring process parameter space, studying process variations, and quantifying uncertainty in material properties. However, integrating process-structure and structure-property simulations often involves connecting multiple software packages. Each package may use unique data structures and require substantial domain knowledge. This presentation demonstrates PSP simulation capabilities of Materialite, an open-source package developed at NASA Langley Research Center. Materialite simplifies model linkages by using a common data structure and model interface, enabling straightforward simulation across a PSP model chain. Physics-based models, including kinetic Monte Carlo and crystal plasticity, are implemented within the package. The model interface is also intended to simplify implementation of new models and enable integration with external simulation tools. Example use cases include uncertainty quantification with PSP models and GPU-accelerated powder bed fusion AM process models.

additive manufacturing↗
DOE OSTI2026

Interfacial Adhesion Mechanism of Anionic Polyelectrolyte Brushes Induced by Oppositely Charged Macromolecular Counterions

Abstract Polyelectrolyte brushes are widely used as model systems for investigating electrostatic interactions at soft interfaces and to achieve exceptional lubrication properties. Previous studies have primarily focused on their behavior in the presence of multivalent counterions, which induce brush collapse, ionic crosslinking, and pronounced changes in interfacial structures. However, interactions between polyelectrolyte brushes and oppositely charged macromolecular counterions, such as polycations, remain poorly understood. Here, we prepare well-defined polystyrene sulfonate (PSS) brushes fabricated via surface-initiated grafting and employ surface forces apparatus measurements to investigate their interactions with oppositely charged polycations. In contrast to multivalent ions, polycations do not induce noticeable brush collapse, but instead generate significant adhesion between symmetric PSS brush layers. This adhesion increases with both contact time and applied load, eventually reaching a steady plateau, indicating that the interaction is governed not simply by electrostatic screening but by the gradual formation of polycation-mediated bridging under confinement. Furthermore, the introduction of monovalent Na+ ions disrupts the adhesive interaction even at very low concentrations, suggesting that the bridging function of the adsorbed polycation is highly sensitive to competitive ionic screening. In comparison, measurements with trivalent counterions reveal the expected brush collapse behavior but minimal dependence of adhesion on contact time, highlighting a clear mechanistic distinction from the polymeric counterion case. Collectively, these results demonstrate that molecular size, configurational restriction, and confinement-induced rearrangement, rather than charge valency alone, govern the interaction behavior of macromolecular counterions at brush interfaces. This work provides new insight into the molecular origins of adhesion and the regulation of interfacial interactions in charged polymer brush systems.

Park, Jinwoo [Argonne National Laboratory , , , ,;↗
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
NASA NTRS2023

CFD Modeling of Bi-Directional PMD inside Cryogenic Propellant Tanks Onboard Parabolic Flights

Future cryogenic propulsion systems will require efficient methods with which to transfer cryogenic propellants from a depot storage tank to a customer receiver tank to minimize cost and maximize reusability. The Reduced Gravity Cryogenic Transfer project is currently developing advanced cryogenic fluid management technology and developing and validating new numerical models for three phases of transfer: line chilldown, tank chilldown, and tank fill. Additionally, multiple liquid nitrogen (LN 2 ) parabolic flight transfer rigs are being designed by universities and NASA to investigate the gravitational sensitivities that exist in these three technologies. In order to maximize the collection of low-g data during flights, it is required to extract as much (LN 2 as possible from the supply tank, despite variable gravity levels. The purpose of this paper is to present computational fluid dynamics (CFD) volume of fluid simulations of (LN 2 behavior in the supply tank onboard parabolic flights to validate the optimal design of a bi-directional propellant management device (PMD) using the commercial software FLOW-3D. A parametric study is conducted on the effects of gravity level, fill level, pore size, open area, thickness, and type of baffle on PMD performance. Based on results, the PMD as designed exceeds the targeted expulsion efficiency.

Jason Hartwig↗
NASA NTRSDate not supplied

Pinned in Microgravity: How Colloidal Suspensions Inhibit Droplet Motion

Understanding the fundamental behavior of colloidal fluids in microgravity is critical for advancing material science, manufacturing, and biological processes. In this study we investigate the dynamics of a droplet containing colloidal suspension when ejected in short term microgravity. These conditions are generated using a 2.2 second drop tower. A custom droplet generator was designed to produce consistent size-controlled droplets containing passive colloidal particles. The particles are suspended in water and subsequently ejected out of a hydrophobic wedge under microgravity conditions, where the droplet ejection time can be predicted [Torres & Weislogel, 2021]. The ejection times will be compared with de-ionized water ejected from the same system to investigate the effects of change in density and surface tension caused by the colloids. Successful passive colloidal ejection will lay the foundation for active colloid research in the future.

colloids↗
DOE OSTI2026

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Space Power

Appropriate directions for the applied research and technology programs that will develop space power systems for U.S. future space missions beyond 1995 are explored. Spacecraft power supplies; space stations, space power reactors, solar arrays, thermoelectric generators, energy storage, and communication satellites are among the topics discussed.

Source record↗
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
DOE OSTI2025

Characterizing and modeling the mechanical behavior of an anion exchange membrane for carbon capture applications

A new direct air capture (DAC) technology uses a moisture swing (MS) process with anion exchange membranes, potentially offering a more energy-efficient way to remove CO₂ from the air. In this MS process, the membrane absorbs CO₂ as it dries and releases it when water is added. Understanding the mechanical behavior of these membranes is essential for improving the design and efficiency of DAC systems and prolonging sorbent lifetime. This study tested one anion exchange membrane, Fumasep’s FAA-3, under mechanical loading and various temperature and humidity conditions to measure its swelling, stiffness, strength, plastic deformation, and stress relaxation. Experimental results were used to identify a mechanical model for FAA-3 that can be used to predict the material’s nonlinear viscous behavior under various loads and environments. Unlike prior studies that assumed linear elastic behavior, this work incorporates humidity-dependent swelling, thermal expansion, and nonlinear viscoelasticity of FAA-3 in both the experiments and the model.

Direct air capture↗
DOE OSTI2023

Four Years of Type Ia Supernovae Observed by TESS: Early-time Light-curve Shapes and Constraints on Companion Interaction Models

Abstract We present 307 type Ia supernova (SN) light curves from the first 4 yr of the Transiting Exoplanet Survey Satellite mission. We use this sample to characterize the shapes of the early-time light curves, measure the rise times from first light to peak, and search for companion star interactions. Using simulations, we show that light curves must have noise <10% of the peak flux to avoid biases in the early-time light-curve shape, restricting our quantitative analysis to 74 light curves. We find that the mean power-law index t β 1 of the early-time light curves isβ 1 = 1.93 ± 0.57, and the mean rise time to peak is 15.7 ± 3.5 days. The underlying population distribution forβ 1 may instead consist of a Gaussian component with mean 2.29, width 0.34, and a long tail extending to values less than 1.0. We find that the data can rarely distinguish between models with and without companion interaction models. Nevertheless, we find three high-quality light curves that tentatively prefer the addition of a companion interaction model, but the statistical evidence for the companion interactions is not robust. We also find two SNe that disfavor the addition of a companion interaction model to a curved power-law model. Taking the 74 SNe together, we calculate 3σupper limits on the presence of companion signatures to control for orientation effects that can hide companions in individual light curves. Our results rule out common progenitor systems with companions having Roche lobe radii >31R ⊙ (separations >5.7 × 10 12 cm, 99.9% confidence level) and disfavor companions having Roche lobe radii >10R ⊙ (separations >1.9 × 10 12 cm, 95% confidence level). Lastly, we discuss the implications of our results for the intrinsic fraction of single degenerate progenitor systems.

Astronomy & Astrophysics↗
DOE OSTI2025

Grain boundary segregation and chemical ordering in CoCrFeMnNi multi-principal element alloy

Owing to their far-from-dilute compositions, multi-principal element alloys (MPEAs) can exhibit unique combinations of engineering properties. As nearly all MPEAs are polycrystalline aggregates, it is necessary to understand the interactions of various elemental species with grain boundaries (GBs). This is of particular importance in extreme environments, such as radiation and elevated temperatures, where such interactions have implications on the properties of MPEAs. Herein, we employ atomistic simulations to generate a series of [001] asymmetric tilt GBs in a model CoCrFeMnNi MPEA and quantify solute interactions and segregation to these boundaries. We employ the Warren-Cowley order parameters to investigate the interplay between GB segregation and chemical short-range order (SRO). At temperatures above 800 K, simulation results reveal the segregation of Cr and Mn to CoCrFeMnNi GBs and show weak dependence of boundary solute excess on GB geometry, at least for the boundaries explored in this work. At temperatures in the range of 673–800 K, formation of domains rich in Cr is observed at GBs in agreement with experimental observations. Quantitative analysis shows that solute excess of various alloying elements decreases rapidly with the increase in temperature in the range of 1000–1200 K. Furthermore, we show that GB regions exhibit SRO characteristics that are distinct from the bulk crystals, leading to spatial variations in SRO. In broad terms, our study highlights the need to account for GB interactions with alloying elements when designing advanced MPEAs with novel chemistries.

Wang, Yitao [Lehigh University, Bethlehem, PA (Uni↗
DOE OSTI2026

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Fast Growth of Ga2O3 on Highly Offcut (100)-Oriented Substrates

Beta-Ga2O3 has emerged as a leading candidate for next-generation power electronics, radio frequency (RF) switches, and extreme environment electronics due to a wide band gap (4.6 - 4.9 eV), high dopability (approximately 40 meV activation energy for an isolated silicon donor), and melt growth characteristics resulting in commercially available 4-inch substrates and commercial demonstrations of 6-inch substrates by multiple techniques. The (100) surface of Ga2O3 is highly desirable from a device and epitaxy standpoint - bulk growth of (100) material is more scalable than (010), the surface is nearly lattice-matched to p-type partner NiO, and Al2O3 incorporates at higher concentrations without phase separation. More importantly, the impact ionization coefficients along the [100] direction are low, leading to the highest possible critical fields. This is advantageous compared to the current state of the art, (001), due to reduced surface defects and increased possible breakdown voltage. However, the epitaxial growth rate on (100) surfaces is less than 10% of other faces due to weak bonding and favorable desorption, and on-axis (100) growth easily forms twin domains. Recent demonstrations have shown growth rate improvements from 0.4 nm/min to 1.5 nm/min by growing on (100) wafers that are offcut 6 degrees in the -c direction. These films show step-flow growth from (20-1) step-edges and high electron mobility due to suppressed twins. Despite these exciting results, offcuts greater than 6 degrees have not been explored due to the waste associated with grinding and polishing large offcuts. In this talk we will discuss the molecular beam epitaxy (MBE) growth and properties of Beta-Ga2O3 grown on (100) substrates offcut in the -c direction up to 13.4 degrees. These large offcuts are enabled by edge-fed film-defined growth (EFG) where the offcut is grown into the surface by pulling the crystal through the EFG die with the seed crystal rotated by the desired offcut angle. We will demonstrate that 13.4 degrees offcut substrates still exhibit a terraced (100) surface, and that a >10x increase (>5 nm/min) in growth rate is achieved. As previously reported on lower offcuts, we observe reversal of substrate twin domains around the (001) direction at the substrate-epilayer interface. We will discuss electrical properties including record-low (by MBE) unintentional doping densities of < 5E15 cm-3 and critical breakdown field in Schottky barrier diodes comparable with state-of-the-art (001) Ga2O3 without edge termination.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analytical Model for Atomic Relaxation in Twisted Moiré Materials

By virtue of being atomically thin, the electronic properties of heterostructures built from two-dimensional materials are strongly influenced by atomic relaxation. The atomic layers behave as flexible membranes rather than rigid crystals. Here we develop an analytical theory of lattice relaxation in twisted moiré materials. We obtain analytical results for the lattice displacements and corresponding pseudo gauge fields, as a function of twist angle. We benchmark our results for twisted bilayer graphene and twisted WSe 2 bilayers using large-scale molecular dynamics simulations. Our single-parameter theory is valid in graphene bilayers for twist angles 𝜃 ≳ 0.7°, and in twisted WSe 2 for 𝜃 ≳ 1.6°. Furthermore, we also investigate how relaxation alters the electronic structure in twisted bilayer graphene, providing a simple extension to the continuum model to account for lattice relaxation.

36 MATERIALS SCIENCE↗
NASA NTRS2026

GVIPS Prediction of Cyclic Ratchetting Behavior and Cycles with Interrupted Relaxation Periods of TIMETAL 21S at 650 °C

The predictive performance of a multimechanism GVIPS model with saturating hardening function is assessed against a set of high temperature cyclic experiments under both strain-control and stress-control at varying rates of loading as well as several complex interrupted cyclic responses. The material specimens were composed of the titanium alloy, TIMETAL 21S, and tested at 650 °C, see Lissenden et al., 2007. The model did a very good job of predicting the variety of tests, particularly given the fact that only monotonic loading tests (tensile, creep, relaxation) and a single fully reversed cycle was used for characterization. It was particularly interesting that the model was capable of reasonably capturing the rate of accumulated strain during ratchetting under tensile mean stress. Given complex interrupted cyclic response, the model was able to both qualitatively and even quantitatively predict reasonably well both the cyclic and associated relaxation periods (in all quadrants of the stress-strain space) thus confirming the validity of the chosen functional forms for both hardening and thermal recovery. The quantitative inaccuracy is associated with the fact that the material specimens tested by Lissenden et al., 2007 was significantly “softer” than those tested by Castelli in Saleeb et al., 1994 which were used for characterization of the GVIPS model parameters.

Viscoplaticity↗
DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗
NASA NTRSDate not supplied

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. 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. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO↗