Results for “multi-region modeling” — Page 4

NASA NTRS2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

mineral dust↗

NASA NTRS2026

Introduction to the Special Issue on Advanced Air Mobility Noise: Predictions, Measurements and Perception

This Special Issue focuses on noise associated with Advanced Air Mobility (AAM), an emerging class of predominantly electric distributed-propulsion aircraft designed for urban and regional transportation. As these vehicles move toward certification and deployment, noise has become a central challenge for regulatory approval and public acceptance, particularly due to operations in densely populated areas and at low altitudes. The 24 contributions in this issue address three key aspects of AAM noise: prediction, measurement, and human perception. Prediction studies span a wide range of modeling fidelities, from high-resolution simulations to improved semi-analytical approaches, and examine complex aeroacoustic mechanisms including rotor interactions, turbulence ingestion, and broadband noise generation. Measurement studies, largely at model scale, provide new insights into tonal and broadband noise characteristics across configurations and operating conditions, while supporting model validation. Perception-focused contributions investigate annoyance, sound quality metrics, and auralization, emphasizing the role of context and operational factors in shaping human response. Together, these works highlight the interdisciplinary nature of AAM noise research and the need for integrated approaches to enable quieter vehicle design.

Perception↗

NASA NTRS2026

Spatially Refined Satellite Gravimetry Captures Human Signatures in Global Terrestrial Water Storage Trends

Human activities have directly altered the water cycle through water management, aquifer pumping, agricultural irrigation, and land use change. Although satellite gravimetry has transformed global hydrological research, its coarse resolution limits attribution of freshwater change to human activities at many management-relevant scales. Here we assessed global terrestrial water storage (TWS) trends from April 2002 to November 2025 using “stacked” regression of Level-1B intersatellite ranging data, which leverages temporal information and variability to dramatically improve effective spatial resolution relative to standard approaches. We combined this refined product with rigorous uncertainty analysis, autocorrelation-robust geostatistical methods, and literature assessment to evaluate TWS trend associations with land and water use, climate variability, and glacial mass loss. We identified TWS trend hotspots exhibiting significant spatial associations with anthropogenic 40 drivers, including groundwater and surface-water irrigation, rainfed agriculture, deforestation, and reservoir impoundment. Across these regions, cumulative TWS losses (3,122 Gt) substantially exceeded gains (2,432 Gt). Compared with traditional regression of monthly mascons, our approach yielded regional trend magnitudes that are on average 33% larger, revealing that global freshwater depletion, particularly from groundwater pumping, is considerably more acute than previously estimated. Multivariate regression models show that humans account for a significant share of the spatial variability in TWS trends on every non-polar continent except Australia. We detected localized TWS gains linked to rainfed agriculture, surface water irrigation, and reservoir filling that were unresolved in earlier gravimetric studies. The methodology provides a foundation for future gravity missions to independently track decadal freshwater change with unprecedented spatial fidelity.

groundwater↗

NASA NTRS2026

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

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

Spacesuit↗

NASA NTRSDate not supplied

The Future of Oaks in the Santa Monica Mountains: A Case Study in Using Remote Sensing Data for Species Distributions Models

The Woolsey Fire began on November 8, 2018, and lasted for almost two weeks, during which it burned almost 100,000 acres of valuable landscape and habitat, including a vast area of woodland. The persistence of key woodland species provides aesthetic, monetary, and ecological value to the landscape through carbon sequestration, air temperature moderation, and erosion mitigation, among other ecosystem services. This study investigated the impact of the Woolsey Fire on native woodland species distributions and identified areas suitable for restoration within the Santa Monica Mountains National Recreation Area. The team partnered with the Resource Conservation District of the Santa Monica Mountains; National Park Service, Santa Monica Mountains National Recreation Area; California Department of Parks and Recreation, Los Angeles County Division; County of Los Angeles Fire Department, Prevention Services Bureau, Forestry Division; County of Los Angeles Department of Regional Planning; and the University of Montana. The Earth observations used include data from Landsat 8 Operational Land Imager, NASA ER-2 Jet Airborne Visible InfraRed Imaging Spectrometer, Shuttle Radar Topography Mission, and RapidEye. The team produced maps of burn severity from the Woolsey Fire, its impact on plant species distributions, and habitat suitability projections for 2050 and 2099 to assist partners in prioritizing areas for restoration. A plant community classification was successfully created using Multiple Endmember Spectral Mixture Analysis (MESMA). Overall accuracy was assessed at 90.54% by comparing the classification to validation pixels derived from ground truth information provided by our partners.

Roger Ly↗

NASA NTRSDate not supplied

Controls for Electrified Aircraft Propulsion

Advanced electrified aircraft propulsion (EAP) concepts with integrated power, propulsion, and thermal systems require the development of equally sophisticated controllers to fly safe and efficient missions. Hybrid electric aircraft utilize electric machines mechanically coupled to the engine shafts to extract and insert power for a variety of purposes. Additionally, electric machines may drive propulsive fans pulling from a combination of on-board energy storage devices and engine extracted power. NASA’s investment in hybrid electric aircraft hardware-in-the-loop testing enables controls research on representative engine models using a novel emulation and scaling methodology. High-voltage, high-power electrical powertrain presents several technical risks at altitude. High-power requires high efficiency to minimize losses. Superconducting electric machines and power distribution research at NASA has identified key challenges with potential controls solutions. These risks necessitate the development of controllers robust to model uncertainty, disturbances, and fault conditions. System health management and fault detection schemes play an important role in a multi-layer controls approach that utilizes a supervisor which oversees inner loops for the highly coupled power, propulsion, and thermal systems. Such a supervisory controller enables coordinates energy transfer between engines, energy storage devices, electric machines, propulsors, and heat exchangers. Existing methods have shown an improvement in engine operability using the hybrid electric powertrain.

Halle E Buescher↗

NASA NTRS2026

Towards an Improved Understanding of the Antarctic Coastal Zone and Its Contribution to Future Global Sea Level

Understanding the coastal zone of the Antarctic Ice Sheet (AIS), where it interacts with the Southern Ocean and warmer air masses, is crucial for predicting Antarctica's influence on the global climate and sea level. This region has multiple tipping mechanisms that could trigger large, rapid, and potentially irreversible changes in the AIS, the Southern Ocean and their global connections in the coming centuries. The AIS remains the largest source of uncertainty in future sea-level projections. Bed topography beneath the ice shelves and the coastal ice sheet is not yet well documented, and is a major source of this uncertainty. This review assesses current knowledge of the coastal zone and highlights methods to investigate it, including aerogeophysical surveys, ground- and ship-based measurements, satellite observations, and computer modeling. An ensemble analysis of published bed topography data sets identifies significant data gaps and their regional distribution, framed in the context of current ice-sheet behavior and potential instability. We propose scientific priorities and guidelines for future aerogeophysical surveys, advocating for a comprehensive, coordinated international effort to build a next-generation data set of Antarctic bed properties. Such an initiative would significantly advance understanding of the role of coastal processes in ice-sheet dynamics, reducing uncertainties in sea-level rise projections and improving predictions of future ocean and climate changes.

Kenichi Matsuoka↗

NASA NTRSDate not supplied

Advancing Wildfire Monitoring with TEMPO and ML tools: Hourly Smoke and Fire‑Front Mapping and Near‑Surface NO₂ Predictions

Wildfires impose substantial impacts on communities and regions downwind of wildfire smoke. We present a TEMPO‑enabled workflow that generates value‑added Level‑3 smoke‑plume masks and fire‑front maps for large wildfires, such as 2024 Park Fire, using the self‑supervised deep learning system SIT‑FUSE, along with near‑surface NO₂ predictions produced by a foundation model (Microsoft Aurora). We conclude by outlining a roadmap for expanding these capabilities to additional Western U.S. wildfire events and for delivering actionable tools to stakeholders. This open-source, reproducible workflow provides a scalable framework for cross-agency wildfire monitoring to overcome traditional limitations in smoke-cloud discrimination and air-quality forecasting by incorporating TEMPO data and beyond.

Xiaohua Pan↗

NASA NTRSDate not supplied

Conclusions of a Cryogenic Motor Configuration Study for Aircraft Propulsion

Cryogenic electric machines are a promising technology for multi megawatt electric aircraft drivetrains. Superconducting stators that make use of low AC loss superconducting wires are considered to be a key component for these cryogenic machines. This paper shares the lessons learned from a design study that evaluated several configurations of a 5 MW cryogenic motor for this application. The lessons are drawn from a qualitative study of 16 configurations and a quantitative study of 8 down selected configurations. The quantitative study involved using analytical calculations (electromagnetic, thermofluid, mechanical, and motor sizing) to determine each configuration’s Pareto front of total efficiency and total specific power. The conclusions from a preliminary assessment of different stator conductors (superconductors and cryogenic metals) are also presented. One motor configuration is down selected for future detailed design and demonstration at the multi megawatt level.

Electrified aircraft propulsion↗

NASA NTRS2026

Thermal Inspection for the Assessment of Adhesively Bonded Metal Adherents

Certification of adhesively bonded structures is currently a challenge for aircraft manufacturers. The ability to certify bonds in primary structure can reduce dependence on fasteners and this will enable more efficient manufacturing. For example, drilling holes and installing fasteners in bonded joints can be a potentially significant bottleneck in airframe manufacturing. In addition, fasteners increase the airframe weight and can add stress concentration areas. Reducing fastener count can accelerate manufacturing and improve aircraft performance. NASA is currently investigating nondestructive evaluation (NDE) techniques to assess bond integrity. For example, bond thickness influences bond strength and therefore NDE techniques are being investigated to determine bond thickness. For bonded metal structures, there is a large difference in the thermal diffusivity between the 7075 aluminum alloy metal adherends and aircraft grade adhesive. Multi-layered thermal models show a large variation in the thermal response for bondline adhesive thicknesses that vary from 100 to 300 microns for adherend thicknesses of 0.163 cm. Experimental through-transmission thermal measurements reveal promise to quantitatively characterize the bondline thickness. Results were validated with X-ray computed tomography and optical microscopy measurements, and influence of porosity on the thermal model and measurements are investigated.

bondline thickness↗

NASA NTRS2026

Towards an improved understanding of the Antarctic coastal zone and its contribution to future global sea level

Understanding the coastal zone of the Antarctic Ice Sheet, where it interacts with the Southern Ocean and warmer air masses, is crucial for predicting Antarctica's influence on the global climate. This region has multiple tipping mechanisms that could trigger large, rapid, and potentially irreversible changes in the coming centuries. The Antarctic Ice Sheet remains the largest source of uncertainty in future sea-level projections. Insufficient knowledge of bed topography beneath the ice shelves and the coastal ice sheet is not yet well documented, but is a major source of this uncertainty. This review assesses current knowledge of the coastal zone and highlights methods to investigate it, including aerogeophysical surveys, ground- and ship-based measurements, satellite observations, and computer modeling. An ensemble analysis of published bed topography datasets identifies significant data gaps and their regional distribution, framed in the context of current ice-sheet behavior and potential instability. We propose scientific priorities and guidelines for future aerogeophysical surveys, advocating for a comprehensive, coordinated international effort to build a next-generation dataset of Antarctic bed properties. Such an initiative would significantly advance understanding of the role of coastal processes in ice-sheet dynamics, reducing uncertainties in sea-level rise projections and enhancing predictions of future ocean and climate changes.

Kenichi Matsuoka↗

NASA NTRS2026

Thermal Inspection for the Assessment of Adhesively Bonded Metal Adherends

Certification of adhesively bonded structures is currently a challenge for aircraft manufacturers. The ability to certify bonds in primary structure can reduce dependence on fasteners and this will enable more efficient manufacturing. For example, drilling holes and installing fasteners in bonded joints can be a potentially significant bottleneck in airframe manufacturing. In addition, fasteners increase the airframe weight and can add stress concentration areas. Reducing fastener count can accelerate manufacturing and improve aircraft performance. NASA is currently investigating nondestructive evaluation (NDE) techniques to assess bond integrity. For example, bond thickness influences bond strength and therefore NDE techniques are being investigated to determine bond thickness. For bonded metal structures, there is a large difference in the thermal diffusivity between the 7075 aluminum alloy metal adherends and aircraft grade adhesive. Multi-layered thermal models show a large variation in the thermal response for bondline adhesive thicknesses that vary from 100 to 300 microns for adherend thicknesses of 0.163 cm. Experimental through-transmission thermal measurements reveal promise to quantitatively characterize the bondline thickness. Results were validated with X-ray computed tomography and optical microscopy measurements, and influence of porosity on the thermal model and measurements are investigated.

X-ray computed tomography↗

NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

NASA NTRS2025

JWST Reveals Spectral Tracers of Recent Surface Modification on Europa

Europa has been modified by a variety of geologic processes, exposing internally-derived materials that are heavily irradiated by charged particles trapped in Jupiter’s magnetosphere. Prior spectral analysis of H2O ice on Europa relied on low signal-to-noise data at wavelengths >2.5 µm, limiting assessment of a 3.1 µm Fresnel peak that is diagnostic of exposed crystalline ice. We report new measurements of H2O ice spectral features using high signal-to-noise data collected by the NIRSpec spectrograph (1.48 – 5.35 µm) on the James Webb Space Telescope. These data reveal a narrow 3.1 µm crystalline H2O ice Fresnel peak, which is primarily located at southern latitudes in Tara and Powys Regiones. Our analysis indicates that crystalline ice exposed in these low-latitude regiones is likely sustained by ongoing thermal (re)crystallization outpacing charged particle amorphization of the top ∼10 µm of Europa’s regolith over short timescales (<15 days). We also measured H2O ice features centered near 1.5 µm, 1.65 µm, and 2.0 µm, and a broad 3.6 µm H2O continuum peak, which are all stronger at northern latitudes, in contrast to the 3.1 µm Fresnel peak identified at southern latitudes. These results support the hypothesis that H2O ice in Europa’s regolith is vertically stratified, with amorphous ice grains dominating its exposed surface, except in Tara and Powys Regiones. We also find that a previously detected 4.38 µm 13CO2 feature is present almost exclusively at southern latitudes in Tara and Powys Regiones, likely derived from an internal source of carbon-bearing material.

Surface processes (2116)↗

NASA NTRS1977

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↗

NASA NTRS1993

Experimental and analytical ion thruster research

The results of further spectroscopic studies on the plume from a 3 cm ion source operated on an argon propellant is reported on. In particular, it is shown that it should be possible to use the spectroscopic technique to measure the plasma density of the ion plume close to the grids, where it is difficult to use electrical probe measurements. How the technique, along with electrical probe measurements in the far downstream region of the plume, can be used to characterize the operation of a three-grid, 15 cm diameter thruster from NASA JPL is outlined. Pumping speed measurements on the Vacuum Research Facility have shown that this facility should be adequate for testing the JPL thruster at pressures in the low 10(exp -5) Torr range. Finally, we describe a simple analytical model which can be used to calculate the grid impingement current which results from charge-exchange collisions in the ion plume.

Wilhelmus M. Ruyten↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

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

silicone intumescence↗