Search NASASearch

SEARCH · Search NASA

Results for “boundary tracking”

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.

207 records · Page 10

Thermochemical Stability of Ca2Yb8(SiO4)6O2 Apatite in Presence of Molten Calcium-Magnesium-Aluminosilicate (CMAS)

Thermochemical stability of ytterbium silicon oxyapatite Ca 2 Yb 8 (SiO 4 ) 6 O 2 (CYbS) in the presence of molten calcium-magnesium aluminosilicate (CMAS) has been investigated at elevated temperatures for consideration as a thermal and environmental barrier coating (T/EBC) material. CYbS apatite powder was synthesized from the constituent oxides via a solid-state reaction method. Hot-pressed apatite substrates were exposed to molten CMAS at 1200, 1300, and 1400 °C for 1, 10, and 50 h. Development of phases in the interaction region of the heat-treated specimens was monitored using scanning electron microscopy, transmission electron microscopy, high-angle annular dark-field imaging, selected area electron diffraction, and energy dispersive x-ray spectroscopy. Monoclinic cyclosilicate Ca 3 Yb 2 (Si 3 O 9 ) 2 formed from interaction of CYbS apatite with CaO in the CMAS melt at the apatite-CMAS reaction front and continued to nucleate and grow within the residual CMAS in diffusion couples annealed for 1 to 50 h at 1200 °C as well as in those heat treated at 1300 °C for 1 h. Residual CMAS was depleted of Ca when cyclosilicate was present. Dendritic wollastonite (CaSiO3) was observed within the residual CMAS in couples annealed at 1200 and 1300 °C. Ingress of molten CMAS, because of its exponential decrease in viscosity, occurred through open pores and along the grain boundaries of the apatite substrates without any detectable chemical reaction at 1300 and 1400 °C. Results of this study indicate that Ca 2 Yb 8 (SiO 4 ) 6 O 2 apatite has the potential to mitigate the CMAS corrosion up to about 1200 °C but not at higher temperatures.

Narottam P Bansal

Investigating two-dimensional adjoint QCD on the lattice

We present our investigations of SU(N) adjoint QCD in two dimensions with one Majorana fermion on the lattice. We determine the relevant parameter range for the simulations with Wilson fermions and present results for Polyakov loop, chiral condensate, and string tension. In the theory with massive fermions, all observables we checked show qualitative agreement between numerical lattice data and theory, while the massless limit is more subtle since chiral and non-invertible symmetry of the continuum theory are explicitly broken by lattice regularization. In thermal compactification, we observe N perturbative vacua for the holonomy potential at high-T with instanton events connecting them, and a unique vacuum at low-T. At finite-N, this is a cross-over and it turns to a phase transition at large-N thermodynamic limit. In circle compactification with periodic boundary conditions, we observe a unique center-symmetric minimum at any radius. In continuum, the instantons in the thermal case carry zero modes (for even N) and indeed, in the lattice simulations, we observe that chiral condensate is dominated by instanton centers, where zero modes are localized. We present lattice results on the issue of confinement vs. screening in the theory and comment on the roles of chiral symmetry and non-invertible symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Three-dimensional nucleation and growth of deformation twins in magnesium

At two-thirds the weight of aluminum, magnesium alloys have the potential to reduce the fuel consumption of transportation vehicles. These advancements depend on our ability to optimize the desirable versus undesirable effects of deformation twins, which are three-dimensional (3D) microstructural domains that form under mechanical stresses. Previously only characterized through surface or thin-film measurements, we here present 3D in situ characterization of deformation twinning inside an embedded grain over mesoscopic fields of view using dark-field x-ray microscopy supported by crystal plasticity finite element analysis. The results revealed the role of triple junctions on twin nucleation and the sequence and irregularity of twin growth and showed that twin-grain junctions, twin-twin junctions, and twin boundaries were the sites of localized dislocation accumulation.

Lee, Sangwon [Univ. of Michigan, Ann Arbor, MI (Un

A topological superconductor tuned by electronic correlations

A topological superconductor, characterized by either a chiral order parameter or a topological surface state in proximity to bulk superconductivity, is foundational to topological quantum computing. A key open challenge is whether electron-electron interactions can tune such topological superconducting phases. Here, we provide experimental signatures of a unique topological superconducting phase in competition with electronic correlations in 10-unit-cell thick FeTe x Se 1-x films grown on SrTiO 3 substrates. When the Te content x exceeds 0.7, we observe a topological transition marked by the emergence of a superconducting surface state. Near the FeTe limit, the system undergoes another transition where the surface state disappears, and superconductivity is suppressed. Theory suggests that electron-electron interactions in the odd-parity xy− band drives this second topological transition. The flattening and eventual decoherence of d xy -derived bands track the superconducting dome, linking correlation effects directly to superconducting coherent transport. Our work establishes many-body electronic correlations as a sensitive knob for tuning topology and superconductivity, offering a pathway to engineer new topological phases in correlated materials.

36 MATERIALS SCIENCE

The Effect of Gravity on the Combustion of Bulk Metals

In recent years, metal combustion studies at the University of Colorado have focused on the effects of gravity (g) on the ignition and burning behavior of bulk metals. The impetus behind this effort is the understanding of the ignition conditions and flammability properties of structural metals found in oxygen (O,sub>2 ) systems for space applications. Since spacecraft are subjected to higher-than-1g loads during launch and reentry and to a zero-gravity environment while in orbit, the study of ignition and combustion of bulk metals at different gravitational accelerations is of great practical concern. From the scientific standpoint, studies conducted under low gravity conditions provide simplified boundary conditions, since buoyancy is removed, and make possible the identification of fundamental ignition and combustion mechanisms. This investigation is intended to provide experimental verification of the influence of natural convection on the burning behavior of metals. In addition, the study offers the first findings of the influence of gravity on ignition of bulk metals and on the combustion mechanism and structure of metal-oxygen, vapor-phase diffusion flames in a buoyancy-free environment. Titanium (Ti) and magnesium (Mg) metals were chosen because of their importance as elements of structural materials and their simple chemical composition-pure metals instead of multicomponent alloys to simplify chemical and spectroscopic analyses. In addition, these elements present the two different combustion modes observed in metals: heterogeneous or surface burning (for Ti) and homogeneous or gas-phase reaction (for Mg). Finally, Mg, Ti, and their oxides exhibit a wide range of thermophysical and chemical properties. Metal surface temperature profiles, critical and ignition temperatures, propagation rates, burning times, and spectroscopic measurements are obtained under normal and reduced gravity. Visual evidence of all phenomena is provided by high-speed photography.

Melvyn C Branch

Gas permeation properties of amorphous zeolitic imidazolate framework membranes made by atomic/molecular layer deposition

Amorphous metal-organic framework (MOF) membranes are desirable because they may retain some of the molecular sieving properties of their crystalline counterparts while being free of grain boundary defects, which often hinder the consistent achievement of high membrane performance. However, current methods, like melting and compression, for fabricating amorphous MOF membranes involve multi-step processes that require the formation of a crystalline membrane first, that is then amorphized, and therefore, could be challenging to scale. Here, we utilize atomic/molecular layer deposition (ALD/MLD) of diethylzinc (DEZ) and 2-methylimidazole (2mIm) to directly synthesize ultrathin amorphous zeolitic imidazolate framework (aZIF) deposits on γ-alumina-coated α-alumina supports. As the number of ALD/MLD cycles increased from 10 to 300, gas permeances decreased while ideal selectivities increased. Mixture separation factors for C3H6/C3H8, CO2/N2, and H2/C3H8 as high as 4, 37, and 194, respectively, were obtained. At 200 °C and 2.5 bar equimolar feed of H2 and C3H8, a H2/C3H8 mixture separation factor of 185 is obtained with an H2 permeance of ca. 4.15x10-8 mol/m2-s-Pa (124 GPU). Additionally, the membrane achieves a CO2/N2 mixture separation factor of 37 at 25 °C and 1 bar with a CO2 permeance of ca. 3.96x10-8 mol/m2-s-Pa (118 GPU). Considering the vast array of compositionally distinct aZIFs that can be potentially deposited by this approach, an enormously large parameter space for membrane design is emerging to be explored.

36 MATERIALS SCIENCE

Overview of IMPACT Data Acquisition System and Data Reduction Process

This report documents the development of the data acquisition system (DAS) and data reduction methodologies for the Irradiated Material Property Accelerated Characterization Test (IMPACT) experiment at the Advanced Test Reactor (ATR). The IMPACT experiment is designed to enable in-pile measurement of thermal conductivity in metallic nuclear fuels, specifically U-10Zr, using an instrumented thermal conductivity probe. The DAS supports both passive temperature monitoring and active thermal interrogation of the probe through controlled AC and DC excitation. Significant modifications to laboratory-scale systems were required to accommodate the higher resistance paths associated with the in-pile application. Custom electronics and relay-controlled measurement sequencing were developed to enable the measurement and sufficient power delivery to the sensing region. A reduced-order, axisymmetric thermal model based on the thermal quadrupoles method is presented to support data interpretation. This model enables efficient evaluation of transient heat transfer behavior and facilitates solution of the inverse problem required to extract thermal properties from measured signals. Multiple boundary condition formulations are discussed to address varying experimental time scales and geometries. Additionally, machine learning techniques are introduced to support data reduction and improve confidence in inverse solutions. Convolutional neural networks are applied to identify the presence of gas gaps and other evolving geometric features that significantly impact thermal response during irradiation. These efforts contribute to the broader integration of digital twin frameworks and real-time modeling capabilities within the Advanced Fuels Campaign.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN

Identifying microstructures susceptible to pulverization in commercially irradiated high burnup UO2 under LOCA conditions

High burnup fuel fragmentation (HBFF) has been a concern in the nuclear industry for many years. When UO2 reaches and exceeds pellet average burnups around 55 GWd/tU, microstructural changes in the fuel pellet begin to occur that render the pellet susceptible to fine fragmentation under loss-of-coolant accident (LOCA) conditions. During a LOCA event, the cladding balloons and bursts, potentially releasing part of these fine fuel fragments into the reactor pressure vessel. This process is known as fuel fragmentation, relocation, and dispersal (FFRD). Presently, the US nuclear industry is developing a safety basis for FFRD with the goal of increasing burnups and pressurized water reactor cycle lengths. Increased cycle lengths would push the burnup limits of the fuel rods past the known threshold for HBFF susceptibility. In this work, the microstructural impact on HBFF was investigated by comparing the microstructures of as-irradiated rod segments to their post-LOCA tested counterparts. This investigation found that varied operational histories result in different microstructural evolutions across the pellet radius, which impacts the fragmentation behavior and radial location of the fragmentation. It is noted that regions with a high bubble density and a high density of grain boundaries fragment under LOCA-relevant conditions. In addition to the fragmentation that has been previously reported at the periphery of the fuel, fragmentation was also noted in the dark zone toward the fuel center. Analysis of the power histories of the fuel samples suggests that the dark zone forms in the central regions of the pellet between temperatures of approximately 740 and 960 °C.

McKinney, Casey [ORNL] (ORCID:0000000335383614)

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions to the finite element problem, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization

Hydrology Copilot: A Cloud-Native Ai System for Hydrological Data Analysis

The emergence of AI-driven Earth observation systems promises to broaden access to petabyte-scale geospatial data beyond domain specialists. However, translating this vision into operational scientific infrastructure requires addressing fundamental challenges in data virtualization, code transparency, and domain-specific reasoning. We present Hydrology Copilot, a cloud-native AI framework for natural-language-driven analysis of Earth observation data. To demonstrate operational capabilities at scale, we implement the system using NASA's North American Land Data Assimilation System version 3 (NLDAS-3), which provides surface meteorological forcing and land-surface model output across North and Central America at 1-km resolution, from which drought diagnostics are derived. The system integrates five core contributions: (1) scalable data virtualization using Kerchunk-based cloud optimized access, achieving a 1.5 to 4.6 times improvement in I/O latency across benchmark queries spanning regional single-day extractions (4.6 times speedup) to continental monthly aggregations (1.5 times speedup); (2) transparent code generation through Microsoft Azure AI Foundry agents that expose executable Python workflows for scientific verification; (3) persistent conversational memory enabling multi-turn analytical discourse across sessions; (4) intelligent query validation that enforces dataset boundaries and resolves ambiguous requests before execution; and (5) a multi-agent architecture coordinating query parsing, code generation, and visualization. We evaluate the system through drought-monitoring workflows, demonstrating reliable code generation, accurate results validated against reference computations and the operational U.S. Drought Monitor, and efficient operation across increasingly complex tasks. By bridging natural-language interfaces with rigorous hydrological analysis, Hydrology Copilot advances beyond proof-of-concept demonstrations to provide a deployable framework for operational Earth science applications.

Data virtualization

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys

Self-Healing Lithium Dendrites through Spontaneous Passivating Layer Formation for Stable Solid-State Lithium–Metal Batteries

All-solid-state lithium–metal batteries have attracted significant attention, owing to their high energy density and superior safety. However, lithium–metal penetration through the solid electrolyte, leading to short-circuiting, remains a critical failure mode that demands comprehensive mitigation strategies. Most existing strategies are effective only prior to the initiation of lithium-dendrite formation and fail once dendrites begin to propagate through the electrolyte. In this study, we propose a self-healing mechanism in which the penetrated lithium reacts with a self-healing agent to form a passivating layer along the particle boundaries. Lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) was incorporated into a Li 6 PS 5 Cl solid electrolyte as the self-healing agent to suppress lithium-dendrite propagation even after dendrite formation initiated under high current densities. The self-healing induced by LiTFSI was verified through comprehensive experimental analyses and was further demonstrated in a full-cell configuration. Moreover, LiTFSI incorporation plays an important role in increasing the critical current density by reducing the overall electronic conductivity of the solid electrolyte and facilitating the formation of a robust LiF-containing solid-electrolyte interphase.

25 ENERGY STORAGE

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization

Observations of Offshore Low‐Level Jets Off the U.S. East Coast Reveal Systematic Biases in ERA5 and HRRR

Low-level jets (LLJs)—wind speed maxima typically occurring a few hundred meters above the surface—are common off the U.S. East Coast and influence many atmospheric processes with societal importance, including cloud formation, aviation safety, and search-and-rescue. However, their vertical structure and frequency remain poorly quantified due to limited offshore observations. This study presents new scanning Doppler LiDAR and infrared spectroradiometer data from the 2024 summer deployment of an offshore barge during the Wind Forecast Improvement Project 3. These coupled wind and temperature profiles provide unprecedented resolution to assess LLJ behavior and model performance. LLJs occurred in over 21% of observed profiles, with a weak diurnal preference for nighttime and early morning hours and maximum winds typically near 300 m. Both ERA5 and High-Resolution Rapid Refresh analysis underestimate jet wind speeds and misrepresent the boundary layer thermal structure. These results highlight persistent model biases and the critical need for high-resolution offshore observations.

17 WIND ENERGY

Analysis of Gas Absorption to A Thin Liquid Film in the Presence of A Zero-Order Chemical Reaction

The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

S Rajagopalan

Lunar Elemental Analysis Obtained From the Apollo Gamma-Ray and X-Ray Remote Sensing Experiment

Gamma-ray and X-ray spectrometers carried in the Service Modules of the Apollo 15 and Apollo 16 spacecraft were employed for compositional mapping of the lunar surface. The measurements involved the observation of the intensity and characteristic energy distribution of gamma rays and X-rays emitted from the lunar surface. A large-scale compositional map of over 10 percent of the lunar surface was obtained from an analysis of the observed spectra. The Apollo 15 flight was at a lunar orbital inclination of 29 degrees as compared with a 9-degree inclination of the Apollo 16 flight; thus, the projected ground track of the Apollo 15 flight covered a larger projected surface area than that of the Apollo 16 flight. The objective of the X-ray experiment was to measure the K spectral lines from Mg, Al, and Si. Spectra were obtained and the data were reduced to Al/Si and Mg/Si intensity ratios and ultimately to chemical ratios. Analyses of the results have indicated (1) that the Al/Si ratios are highest in the lunar highlands and considerably lower in the maria, and (2) that the Mg/Si concentrations generally show the opposite relationship. There is a tendency for the Al/Si values to increase from the western mare areas to the eastern limb highlands. There are distinct chemical contrasts between such features: as the small mare basins and the highland rims. The objective of the gamma-ray experiment was to measure the natural and cosmic-ray-induced activity emission spectrum. At this time, the elemental abundances for Th, U, K, Fe, Ti, Si, and O have been determined over a number of major lunar regions. Regions of relatively high natural radioactivity were found in the Mare Imbrium and Oceanus Procellarum regions. High spots of natural radioactivity were also found south of Fra Mauro; somewhat southwest of Archimedes, and south of Aristarchus. An enhanced region of natural radioactivity was found around Van de Graaff on the far side of tile Moon. In regions other than Mare Imbrium and Oceanus Procellarum, an anticorrelation between natural radioactivity and lunar elevation, as determined from the Apollo laser, altimeter, has been found. From the combined results of the gamma-ray and X-ray spectrometer experiments almost complete information concerning the major element composition of over 10 percent of the lunar surface has been obtained. Distributions have also been mapped for K, Th, and Ti. Interesting correlations between lunar topography and magnetic and gravitational properties have been found.

J I Trombka

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