Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Radial Equilibrium”

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.

32 records · Page 2

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method

Extensional Flow Convecting a Reactant Undergoing a First Order Homogeneous Reaction and Diffusional Mass Transfer From a Sphere at Low to Intermediate Peclet and Damkohler Numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 ≤ Pe ≤ 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 ≤ Da ≤ 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r,Θ) = Σ c n (r)P n (cosΘ) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions P n (cosΘ), the variable coefficient PDE can be reduced to a system of N+1 second order ODEs for the radial functions C n (r), n=0,1,2, ... ,N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N+ 1) first order differential equations as a two-point boundary value problem on 1 ≤ r ≤ r •• . Convergence of the expansion for a specific value of N can thus be established and provides "spectral" behavior as well as the full concentration field c(r,Θ).

N Y Shah

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven

Lessons Learned From 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. NASA has undertaken 6-year technical challenge to address this need by designing a 5 MW superconducting motor and cryogenic drive and developing multi-MW prototypes. This poster shares the lessons learned from a completed 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 a motor design code based on analytical calculations (electromagnetic, thermofluid, mechanical, and motor sizing) to determine each configuration’s Pareto front of total efficiency and total specific power. This poster builds upon preliminary results that were previously presented by (1) reporting final results based on a revised design code for each configuration, including a considerable increase in fidelity for a novel radial flux motor and (2) making a down selection of the configuration for future detailed design and demonstration and documenting the rationale for the selection.

Superconductors

Operational Correction for the Temperature Dependence of A Class of Hyperspectral Radiometers

The correction for temperature response of optical radiometers implies access to their working temperature (i.e., either the internal or ambient one with the radiometers in thermal equilibrium). With specific reference to a class of widely used hyperspectral radiometers lacking a thermistor to measure the internal temperature, this work investigated the potential for operational corrections relying on the ambient temperature derived from the radiometer dark digital counts. With the objective to keep within approximately 1% the uncertainty for radiometric corrections in the 400–800-nm spectral range, findings showed the applicability of the method for ambient temperatures higher than approximately 30°C. The alternative use of an ambient temperature proxy, determined from the air temperature incremented by a few degrees Celsius, was shown a valid solution to constrain the uncertainty of radiometric corrections to within 1% over the entire 0°–40°C interval.

Giuseppe Zibordi

Nitrile/Buna N Material Failure Assessment for an O-Ring used on the Gaseous Hydrogen Flow Control Valve (FCV) of the Space Shuttle Main Engine

After the rollout of Space Shuttle Discovery in April 2005 in preparation for return-to-flight, there was a failure of the Orbiter (OV-103) helium signature leak test in the gaseous hydrogen (GH2) system. Leakage was attributed to the Flow Control Valve (FCV) in Main Engine 3. The FCV determined to be the source of the leak for OV-103 is designated as LV-58. The nitrile/Buna N rubber O-ring seal was removed from LV-58, and failure analysis indicated radial cracks providing leak paths in one quadrant. Cracks were eventually found in 6 of 9 FCV O-rings among the three Shuttle Orbiters, though none were as severe as those for LV-58, OV-103. Testing by EM10 at MSFC on all 9 FCV O- rings included: laser dimensional, Shore A hardness and properties from a dynamic mechanical analyzer (DMA) and an Instron tensile machine. The following test data was obtained on the cracked quadrant of the LV-58, OV-103 O-ring: (1) the estimated compression set was only 9.5%, compared to none for the rest of the O-ring; (2) Shore A hardness for the O.D. was higher by almost 4 durometer points than for the rest of the O-ring; and (3) DMA data showed that the storage/elastic modulus E was almost 25% lower than for the rest of the O-ring. Of the 8 FCV O-rings tested on an Instron, 4 yielded tensile strengths that were below the MIL spec requirement of 1350 psi-a likely influence of rubber cracking. Comparisons were made between values of modulus determined by DNA (elastic) and Instron (Young s). Each nitrile/Buna N O-ring used in the FCV conforms to the MIL-P-25732C specification. A number of such O-rings taken from shelf storage at MSFC and Kennedy Space Center (KSC) were used to generate a reference curve of DMA glass transition temperature (Tg) vs. shelf storage time ranging from 8 to 26 years. A similar reference curve of TGA onset temperature (of rubber weight loss) vs. shelf storage time was also generated. The DMA and TGA data for the used FCV O-rings were compared to the reference curves. Correlations were also made between the DMA modulus (at 22 C) and Shore A hardness for all 9 of the FCV O-rings used among the three Shuttle Orbiters. The radial cracking in the FCV O-rings was determined to be due to ozone attack, as nitrile/Buna N rubber is susceptible to such attack. Nitrile/Buna N material under MIL-P25732C should be used in a hydraulic fluid environment to help protect it from cracking. However, the FCV O-rings were used in an air only environment. The FCV design has as much as a 9-mil gap that allows the O.D. of the O-ring to be directly exposed to ozone, pressurized air and some elevated temperatures, accelerating the weathering process that leads to O-ring cracking. Space Shuttle flights will likely not continue past 2010. Therefore, Shuttle management decided to continue using the nitrile/Buna N material for the FCVs, but have each O-ring replaced after 3 years to minimize any chances for crack initiation.

Doug Wingard

Combustor Technology for Future Small Gas Turbine Aircraft

To enhance fuel efficiency, future advanced small gas turbine engines will utilize engine cycles calling for overall engine pressure ratios, leading to higher combustor inlet pressures and temperatures. Further, the temperature rise through the combustor and the corresponding exit temperature are also expected to increase. This report describes future combustor technology needs for small gas turbine engines. New fuel injectors with large turndown ratios which produce uniform circumferential and radial temperature patterns will be required. Uniform burning will be of greater importance because hot gas temperatures will approach turbine material limits. The higher combustion temperatures and increased radiation at high pressures will put a greater heat load on the combustor liners. At the same time, less cooling air will be available as more of the air will be used for combustion. Thus, improved cooling concepts and/or materials requiring little or no direct cooling will be required. Although presently there are no requirements for emissions levels from small gas turbine engines, regulation is anticipated in the near future. This will require the development of low emission combustors. In particular, nitrogen oxides will increase substantially if new technologies limiting their formation are not evolved and implemented. For example, staged combustion employing lean, premixed/prevaporized, lean direct injection, or rich burn-quick quench-lean burn concepts could replace conventional single stage combustors. Due to combustor size considerations, staged combustion is more easily accommodated in large engines. The inclusion of staged combustion in small engines will pose greater combustor design challenges.

Valerie J Lyons

Helmholtz Motor: A Novel Electric Machine That Enables High Temperature Superconducting Armatures

Electrified aircraft are being developed to increase the efficiency and reduce the cost of operating subsonic transport aircraft. Achieving a substantial impact necessitates focusing on single- and twin-aisle aircraft which use propulsion systems with >20 MW ratings. For these aircraft, multi-MW superconducting electric machines are being developed due to their high specific power and high efficiency. Cryogenic electric machines are a promising technology for multi megawatt electric aircraft drivetrains. Second generation high temperature superconductors (HTS) are the ideal conductor for the field winding (typically the rotor) of superconducting machines, and they have attractive features for the armature winding (typically the stator) of these machines. However, their use in armature windings has been severely limited due to excessive AC losses in conventional electric machines. This presentation will present the working principle, design, and a performance trade study of a novel electric machine – termed a Helmholtz machine – that enables HTS armatures by significantly reducing AC losses. The working principle of the patent pending Helmholtz motor is to optimize the rotor-produced magnetic field (i.e., the field winding) so that it produces a magnetic field almost purely in the plane of the HTS. This is accomplished using 2 rotors that contain a set of matched Helmholtz coil pairs (or matched permanent magnet Halbach arrays). The armature (stator) HTS coils are positioned between the rotors and oriented to minimize the out of plane magnetic field. The armature HTS coils are also designed to minimize the out of plane component of their self field. The presentation will include the lessons learned in the sizing and design of this motor, including equations for how to match the magnet arrays. A performance trade study will be presented for a partially superconducting, radial flux version of the motor. The trade study involved an optimization of a motor design code based on analytical calculations (electromagnetic, AC loss, thermofluid, mechanical, and motor sizing) that determine the motor’s total efficiency and total specific power. The design code will be described in the presentation, with an emphasis on the calculation of magnetic fields, AC loss, and temperature distribution.

Superconducting electric machines

The Origin of Palimpsests

Palimpsests are the most enigmatic impact features observed on icy worlds of the Solar System. Palimpsests are large impact features which are bright, old, mostly flat, and roughly circular. They range from >150 to ~350 km in diameter and sometimes exhibit recognizable associated radial secondary ejecta crater strings beyond their rims. Curiously, among all the Solar System's icy worlds, palimpsests are exclusively seen on the icy Galilean moons. Palimpsests are obviously different than similarly sized impact features seen on other icy satellites. For example, Evander, a 350 km diameter peak-ring impact basin on Saturn’s moon Dione, has a very low depth/diameter ratio of 0.008, indicating significant viscous relaxation. Nevertheless, it still retains short wavelength topography including a prominent central peak, peak ring, and a continuous rim, none of which are seen at palimpsests. Palimpsests display an approximately circular, and dominating, expanse of undulating plains, which is level with surrounding terrain and shows no recognizable rim. At the center of palimpsests are smooth central plains that are enclosed by a low scarp. Central plains material superposes the surrounding undulating plains and is littered with small blocks reminiscent of “small chaos” on Europa. Palimpsests have yet to be explored using modern impact simulation tools. We have begun to investigate two competing hypotheses: (A) Palimpsests form on targets which have a preexisting near-subsurface layer that is liquid or solid but very weak, and is excavated as proximal ejecta; or (B) They form in a solid ice target underlain by a steep thermal gradient that strongly influences the resulting impact feature, including the production of impact-induced melting at depth that is then a major component of proximal ejecta. Our results will have profound implications for the early thermal history of the icy Galilean satellites. We will test our two hypotheses for the formation of palimpsests on Ganymede (where they are both best observed and preserved) by means of 2D iSALE hydrodynamic simulations of impacts into ice targets that will be structured according to the scenarios for each hypothesis, with select runs in 3D. We have already conducted a pilot study using 2D iSALE that has yielded promising results, which we will report out at the meeting.

palimpsests

Orbital Processing of High-Quality CdTe Compound Semiconductors

CdZnTe crystals were grown in one-g and in micro-g for comparative analysis. The two micro-g crystals were grown in the Crystal Growth Furnace during the First United States Microgravity Laboratory mission (USML-1). The samples were analyzed for chemical homogeneity, structural perfection, and optoelectronic performance (infrared transmission). Fourier Transform Infrared (FTIR) transmission of both ground and flight materials showed that the infrared transmission was close to theoretical, 63% versus 66%, suggesting that the material was close to the stochiometric composition during both the ground and flight experiments. Infrared microscopy confirmed that the principal precipitates were Te and their size (1-10 microns) and density suggested that the primary flight and ground base samples experienced similar cooling rates. Macrosegregation was predicted, using scaling analysis, to be low even in one-g crystals and this was confirmed experimentally, with nearly diffusion controlled growth achieved even in the partial mixing regime on the ground. Radial segregation was monitored in the flight samples and was found to vary with fraction solidified, but was disturbed due to the asymmetric grvitational and thermal fields experienced by the flight samples. The flight samples, however, were found to be much higher in structural perfection than the ground samples produced in the same furnace under identical growth conditions except for the gravitational level. Rocking curve widths were found to be substantially reduced, from 20/35 (one-g) to 9/20 (micro-g) for the best regions of the crystals. The full width at half maximum (FWHM) of 9 arc seconds is as good as the best reported terrestrially for this material. The ground samples were found to have a fully developed mosaic structure consisting of subgrains, whereas the flight sample dislocations were discrete and no mosaic substructure was evident. The defect density was reduced from 50-100,000 (one-g) to 500-25000 EPD (micro-g). These results were confirmed using rocking curve analysis, synchrotron topography, and etch pit analysis. The low dislocation density is thought to have resulted from the near-absence of hydrostatic pressure which allowed the melt to solidify with minimum or no wall contact, resulting in very low stress being exerted on the crystal during growth or during post-solidification cooling.

D J Larson, Jr.

Coma Physics of an Interstellar Object: JWST Spatial-Spectral Mapping of 3I/ATLAS

We report a survey of molecular emission from cometary volatiles using the James Webb Space Telescope (JWST) toward interstellar object 3I/ATLAS carried out on UT 2025 December 22 and 23 at a heliocentric distance (\rh{}) of $2.37-2.41$ au. These measurements of CO, \ce{CO2}, \ce{H2O}, \ce{CH3OH}, and \ce{CH4} sampled molecular chemistry in 3I/ATLAS as it receded from its encounter with our Sun and entered the vicinity of the \ce{H2O} ice line --- the region between \rh{} = $2-3$ au where the temperature becomes too low for H$_2$O to vigorously sublime and CO and \ce{CO2} begin to control the overall activity. CO was the most abundant molecule, followed by \ce{H2O} and \ce{CO2}, whose molecular abundances with respect to CO were $(40.5\pm3.1)\%$ and ($41.6\pm0.3)\%$, respectively. This work presents spatial-spectral maps of column density and rotational temperature as a function of distance from the nucleus for all detected species. The spatial distributions of both quantities were highly anisotropic for the apolar species in the coma of 3I/ATLAS, yet were more nearly symmetric for the polar molecules. These results demonstrate how volatiles were segregated in the nucleus ices of 3I/ATLAS and reveal heating and cooling mechanisms in its coma. Derived maps of the ortho-to-para ratio (OPR) for \ce{H2O} were flat with increasing distance from the nucleus and consistent with a coma-averaged value $\mathrm{OPR}=2.7\pm0.2$, slightly less than the expected equilibrium value of three.

Nathan X Roth

A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations

Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.

Radiation

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

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

Bohdan O Wesely

An Overview of Experiments and Modeling of Polysiloxane-Coated Thermal Protection Systems for Missions to Mars, Titan, and Beyond.

Phenolic Impregnated Carbon Ablator (PICA) gained heritage during the historic Stardust mission, where it successfully returned samples from a comet’s tail and has since been instrumental in delivering payloads to the surface of Mars [1-3]. Most recently, PICA enabled the safe return of samples collected from asteroid Bennu as part of the OSIRIS-REx mission. This rich legacy underscores PICA’s critical role in allowing NASA’s most ambitious exploration missions. However, the friable nature of its phenolic phase presents challenges during handling and pre-launch activities. To mitigate this issue, PICA is coated with a polysiloxane resin system, which serves to suppress particulate dispersion and thereby safeguard spacecraft components. A comprehensive understanding of the polysiloxane resin’s behavior is imperative, as it profoundly shapes the material response of PICA during atmospheric entry by influencing its thermal and oxidative stability. This influence extends to thermocouple plugs embedded within thermal protection systems. These plugs have demonstrated their significance in missions such as Mars Science Laboratory (MSL) and Mars 2020, where the MEDLI and MEDLI2 instrumentation suites delivered in-valuable insights into the performance of thermal protection systems during entry into the Martian atmosphere [4]. Looking ahead, missions such as Dragonfly, set to descend into Titan’s dense atmosphere, aim to leverage advanced sensor technologies to further refine our understanding of thermal protection response [5]. Moreover, thermocouple plugs play an essential role in validating cutting-edge material response models, such as those pioneered under NASA’s Entry Systems Modeling Project (ESM), designed, in-part, to predict the operational integrity of thermal protection systems under the extreme stresses of atmospheric entry. To achieve these modeling goals, ground-based experiments are crucial to provide the foundational data necessary for developing and refining these predictive tools. To this end, an extensive test campaign was conducted at the Hypersonic Materials Environmental Test System (HyMETS) to investigate the high-temperature behavior of the polysiloxane resin in an air environment [6]. These experiments revealed critical phenomena, including the formation of a silicon oxycarbide layer that enhances oxidation resistance, moderates surface temperatures, and alters in-depth thermal response. Building on these findings, subsequent tests were designed to simulate atmospheric entry conditions in reactive gases, such as CO2 and N2, to mimic the environments of Mars and Titan, respectively, as well as non-reactive gases representing the atmospheres of the Ice Giants (Neptune and Uranus). A heating rate dependent decomposition mechanism has been identified for the polysiloxane resin under oxidizing conditions (Fig. 1). In the initial stage, the resin and the underlying thermal protection system undergo pyrolysis, rapidly generating a thin amorphous silicon oxycarbide interwoven with carbonaceous char and residual fibers from PICA. During the second stage, the nascent oxide layer establishes a robust, oxidation-resistant thermal barrier coating, which significantly impedes heat transfer to the underlying carbonaceous char, resulting in a stagnation of the surface temperature. A key factor contributing to this thermal resistance is the low recombination efficiency of atomic oxygen (γ), which further diminishes the heat load on the material’s interior layers [7]. Moreover, as the surface temperature stagnates, the silicon oxycarbide phase separates into distinct regions of silica and free graphite. Ultimately, when the heat flux reaches a critical threshold, a third stage is triggered, leading to the breakdown of the coating through carbothermal reduction, exposing the underlying char layer. This exposure leads to a dramatic surface temperature spike, driven by highly exothermic reactions between atomic oxygen and the char layer, further accelerating material degradation. A detailed mass and heat transfer model of PICA coated with polysiloxane resin was implemented in the Porous material Analysis Toolbox based on OpenFOAM, PATO [8]. The initial stage was considered negligible in this model because the resin decomposition occurs rapidly within a thin surface layer. Instead, the coating was directly considered as an oxygen-resistant thermal barrier coating. For the second stage, the thin amorphous silicon oxycarbide was treated as a pure silica surface to simplify the thermochemical behavior. The model ac-counts for surface equilibrium processes using representative elements of the coating-environment system. For the third stage, specific boundary conditions were developed to estimate the onset and progression of the coating removal. Two-dimensional material response simulations were conducted to compare uncoated and coated PICA using boundary conditions calibrated with HyMETS data. Fig. 2 illustrates that the simulations closely align with experimental data, successfully reproducing measured temperature profiles. This work will include the latest advancements in the coating model, including the calibration of recombination of atomic oxygen at the surface during the second phase. These simulated results will be further validated against additional CO2 data points from HyMETS, reinforcing the models’ predictive capabilities. These mechanisms and their effects on thermal protection systems, including thermochemical behavior and thermocouple probe performance in extreme environments, provide crucial insights for optimizing spacecraft designs that safeguard scientific payload and ensure mission success in future planetary exploration endeavors.

Active Oxidation