Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TPS material modeling”

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

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

At least 199 records · Page 11

High Velocity Impact Performance of a Dual Layer Thermal Protection System for the Mars Sample Return Earth Entry Vehicle

The Mars Sample Return (MSR) Earth Entry Vehicle (EEV) is currently planned on being released from its Micro-Meteorite/Orbital Debri (MM/OD) shielding housing about two days before the Earth entry phase. This leaves the EEV exposed to incoming MM/OD impacts, potentially damaging the heat shield and compromising its Entry, Decent, and Landing (EDL) integrity. Currently, two materials are proposed to comprise the MSR-EEV heat shield, a dual layer material Heat-shield for Extreme Entry Environment Technology (HEEET) and Phenolic Impregnated Carbon Ablator (PICA). PICA has been well characterized for OD class impacts, ~7km/s high mass impacts, from previous testing done in the Orion program, but hasn't under-gone extensive MM impact testing. HEEET is a rela-tively new material with minimal prior testing in re-gards to High Velocity Impacts (HVI). In order to in-form selectability of a material, it is crucial to under-stand the material performance when faced with an HVI, directly affecting mission success probability.The current measure of a Thermal Protection Sys-tem's (TPS) performance against an HVI is evaluating a thermally sized material against its derived Balistic Limit Equation (BLE). A BLE is generated empirically from multiple shots of HVI testing, and is used as a first order method in evaluating TPS's performance against the expected MM/OD environment. This method has proved useful for previous uniform densi-ty TPS materials, but has never been validated against a dual layer recession material such as HEEET. Testing in the MSR program for FY19 has a re-quirement to assess the effects of a dual layer TPS by testing various thicknesses of HEEET's Recession Layer, seen in Figure 1. From this data, BLE's will be derived for the individual thickness ratio samples, as well as the material as a whole to evaluate if a heritage form of the BLE can capture the complex physics associated with a dual layer system. Crater morpholo-gy will also be assessed with post-test Non-Destructive Evaluation (NDE) methods such as CT scanning to visualize if a BLE can well predict the associated pene-tration depths, since a BLE assumes full disinigration of the impacting particle and is generally only used to size a spherical crater ? disregarding any shrapenel effects from a high density impactor.To inform this analysis, expected MM environ-ments from the Meteoroid Engineering Model (MEM) and analytical equations for mass flux of incoming MM were evaluated against the notional MSR-EEV trajectory [1]. Using those dispersions, a monte-carlo was run to determine the most probable particle pa-rameters, as well as the riskiest in terms of full bondline penetration. From these probabilities, a test matrix was designed to test against bounding cases for the various parameters of the BLE: projectile density, projectile mass, projectile velocity, and the impact angle.This presentation will discuss the performance of the dual layer TPS material HEEET against a wide range of impact kinetic energies and densitites, as well as the comparison of HEEET to PICA in terms of MM/OD performance and selectability criteria.

Libben, B. J.↗

Overview of Micrometeoroid and Orbital Debris Analysis Process for Mars Sample Return Earth Entry System

Introduction: Micrometeoroid and orbital debris (MMOD) risk analyses for the Mars Sample Return (MSR) Earth Entry System (EES) have been significantly more rigorous than previously flown missions because of its categorization as a Class V restricted return mission. This means the returned samples present significant concern for biogenic contamination. These analyses seek to determine if a micrometeoroid or orbital debris strike would result in loss of containment assurance and are summarized in the flowchart in Fig. 1. Methodology: The mission is considered in two MMOD phases: a pre-release phase where the EES is protected by a Micrometeoroid Protection System (MMPS) and a post-release phase called “free-flight” where the EES is exposed directly to the MMOD environment. These phases correspond to interplanetary cruise and imminent re-entry, respectively. To inform the MMPS design, a 30-shot high velocity impact testing (HVIT) series on candidate configurations at NASA White Sands Test Facility was completed in the summer of 2022. Sample post-shot images are shown in Fig 2 [1]. These data are used to baseline the MMPS design and to tune the hydrocode simulations. ALE3D, CTH, and SPHC are the hydrocodes that simulate physics of high-speed impacts [2]. Results generated with these populate a penetration depth versus energy space beyond the testable velocity regime of HVIT (~7 km/s). The penetration depth versus energy space data are used to define a critical projectile diameter function called a Ballistic Limit Equation (BLE), where the projectile “criticality” is determined by zone dependent failure criteria defined a-priori [3]. For example, the nose of the heatshield has a failure criterion of 50% TPS penetration, assigned because the landing loads are concentrated on that region and no substructure damage is permitted. The BLEs for each vehicle material zone are input into the BUMPER 3 code, along with the vehicle surface mesh and the corresponding space environment model, to calculate a probability of penetration or number of penetrations. The environment models, MEM3 for MM and ORDEM 3.2 for OD, simulate the meteoroid environment from 0.2 to 2 au based on the Grün flux equation, and the debris environment up to 40,000 km altitude from Earth surface, respectively [4,5]. Presentation Focus: The presentation or poster will present the results to-date focusing on the full risk analysis process flow seen in Fig. 1. Details on the derivation of the failure criteria will be discussed, along with HVIT results and how these influenced the MMPS configuration baseline decision. Further, results of hydrocode simulations will be presented and the tuning to HVIT outputs will be described. Finally, the strategies that direct the BLE formulation will be reviewed, specifically, for the EES elements that are most exposed to the MMOD environment.

mmod↗

Technology Readiness Assessment for HEEET TPS

The HEEET project was conceived to develop a heatshield with a high performance ablative thermal protection material that can withstand the extreme entry environment produced as a result of rapid deceleration during high speed entry into Venus, Saturn, Uranus or higher speed entry into Earth's atmosphere. Successful maturation of HEEET supports future New Frontiers and Discovery AO's, as well as Flagship and directed missions in the longer term. In addition, HEEET has the potential to evolve and to support re-entry to Earth, for missions such as Mars Sample Return.The primary goal of the HEEET Project was to develop an ablative TPS heat-shield based on woven TPS technology to Technology Readiness Level (TRL) 6. Key evidence to support the TRL evaluation includes: Demonstration of reproducible manufacturing of a dual layer material over a range of thicknesses and integrated on to a heatshield engineering test unit at a scale that is applicable to near term Discovery as the highest priority and future NF missions as secondary priority set of missions. Demonstration of predictable and stable performance of the dual layer TPS over a range of entry environments that are applicable to near term Discovery and NF missions of interest to SMD.Includes completion of coupon arc jet and laser testing and development of a mid-fidelity thermal response model that correlates with test results. Demonstration of flight heatshield system design for a range of sizes and loads that are relevant to near term Discovery and NF missions of interest to SMD. Includes completion of structural testing to validate analytic thermal/structural models and development of a material property database. Includes structural testing of a ~1m Engineering Test Unit under relevant entry loads.

Gage, Peter↗

Structural Verification of the Redesigned Space Shuttle Bipod Foam Closeout

This document outlines the structural verification approach for the Space Shuttle External Tank Forward Bipod Foam Closeout. Due to the Space Shuttle Columbia accident, debris has become a major concern. The intent of the structural verification is to ensure that any debris shed from the bipod is within acceptable limits. Since cohesive failure due to internal defects was identified as the most likely cause of the STS-107 bipod ramp foam failure, verification for this failure mode receives particular emphasis. However, all failure modes for TPS are considered and appropriate verification rationale is developed for each failure mode. Figure 1 depicts the structural verification of a production design where analysis and test are the primary methods of verification. It can be seen that the successful completion of structural verification is dependent on three main areas: 1. Production process control and quality assurance must ensure that test articles and/or analytical models are representative of (or conservatively envelope) production hardware in terms of geometry, materials and processing. Variability and defects must be considered. 2. Flight environments must be sufficiently characterized to bound driving environments for all failure modes. Applied environments, either test or analytical, must be representative of flight environments and have a load factor that satisfies design requirements. 3. Structural verification must include all failure modes. A comprehensive list of failure modes and the underlying failure mechanisms has been generated based on flight and test experience. Verification tests and / or analyses must address each failure mode. ET TPS Verification is accomplished by a combination of analysis, test, and similarity.

Poole, Eric L.↗

Towards Quantifying and Correlating RTV Intumescence to Entry-Relevant Conditions

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

RTV↗

High-Temperature Strain Sensing for Aerospace Applications

Thermal protection systems (TPS) and hot structures are utilizing advanced materials that operate at temperatures that exceed abilities to measure structural performance. Robust strain sensors that operate accurately and reliably beyond 1800 F are needed but do not exist. These shortcomings hinder the ability to validate analysis and modeling techniques and hinders the ability to optimize structural designs. This presentation examines high-temperature strain sensing for aerospace applications and, more specifically, seeks to provide strain data for validating finite element models and thermal-structural analyses. Efforts have been made to develop sensor attachment techniques for relevant structural materials at the small test specimen level and to perform laboratory tests to characterize sensor and generate corrections to apply to indicated strains. Areas highlighted in this presentation include sensors, sensor attachment techniques, laboratory evaluation/characterization of strain measurement, and sensor use in large-scale structures.

Piazza, Anthony↗

The Mysteries of Real Materials

The presentation will consist of showing arc jet data mysterious to the modelers. It will show pictures from an arc jet test where a material (unidentified) exhibited a failure mode that nobody understands followed by thermocouple data from arc jet tests on another (unidentified) material of interest in which the T/Cs exhibit repeatable, consistent, fascinating yet frustrating response characteristics that have the modelers stumped. This all happens between RT and 200 F. Unless we figure out what the responsible phenomenology is and can model it, we can't size the TPS with any confidence.

Laub, Bernard↗

A Combined CFD/Material Response Analysis of 3MDCP Arcjet Experiments

The Mars Sample Return Earth Entry System (MSR-EES) project has selected 3-D Woven Mid-Density Carbon Phenolic (3MDCP) as the baseline thermal protection system (TPS) material for the capsule that will return Martian soil samples to Earth sometime in the 2030’s. A series of experiments in the NASA Ames AHF and IHF arcjets will be used to characterize the performance of 3MDCP and to develop and refine material response models for it. The test objectives for these experiments include obtaining in-depth and surface temperature data, mass loss and surface recession measurements, and char depth measurements. To provided pre- and post-test support for the arcjet experiments, a combined CFD/material response was performed using the DPLR and Icarus code. Pre-and post-test analysis on the 3MDCP arcjet experiments was performed using the DPLR CFD and Icarus material response codes. The codes were not tightly coupled, but boundary condition data required by Icarus (e.g., heat transfer coefficients and surface pressures), were extracted from the CFD solutions. This poster shows comparisons between DPLR/Icarus and the experimental data taken during arcjet experiment AHF-348 where 4-inch diameter iso-q models made from 3MDCP were tested in the 12-inch nozzle of the AHF arcjet. This poster also includes a top-level overview of the CFD/material response arcjet simulation process. Best-practices for running CFD solutions of arcjet experiments including how to determine the inflow conditions for the CFD solution from the arc heater settings are discussed. The required inputs to the material response solver are presented, as well as a discussion of how they can be extracted from the CFD solution.

thermal protection systems↗

CHARACTERIZATION ON ANISOTROPIC THERMAL CONDUCTIVITY FOR BIG AREA ADDITIVE MANUFACTURING WITH POLYMERS

Additive manufacturing with polymers has been used mainly for prototyping. A recent development of Big Area Additive Manufacturing (BAAM) at Oak Ridge National Laboratory has opened its applications in the mold and die industry. A numerical simulation and prediction for a mold heating performance requires accurate anisotropic thermal properties of the printed material, which are challenging to obtain, and often requires the use of multiple techniques. The transient plane source (TPS) technique has been widely used due to its ability to measure the thermal properties of an extensive range of materials (solids, liquids, powder). Despite the capability to characterize thermal conductivity k of isotropic and anisotropic materials, the measurements of latter materials are limited to the cases, where the samples have the same thermal conductivity k along x- and y-axis that form the radial plane. In this work, the method for a characterization of k in all three dimensions is developed, and the application of TPS is extended to the determination of thermal properties along the x-, y-, and z-axis individually. The materials are represented by additively manufactured polymers including polylactic acid (PLA) and styrene maleic anhydride (SMA). The developed method consists of (1) a determination of the heat capacity of the polymers by means of TPS in combination with the developed in this work data analysis procedure, (2) a machining three types of cylindrical samples from the same material, with the height corresponding either to x-, y-, or z-direction of printing, and (3) a determination of axial thermal conductivity employing anisotropic model and using previously determined heat capacity

Trofimov, Artem↗

Evaluation of SRB phenolic TPS material made by an alternate vendor

Tests conducted to evaluate the adequacy of solid rocket booster (SRB) phenolic thermal protection system (TPS) material supplied by an alternate vendor chosen by United Space Boosters, Inc. (USBI), to replace the current phenolic TPS sections used thus far on the first four Shuttle flights. The phenolic TPS is applied mainly to the attach and kick rings of the solid rocket booster (SRB). Full-scale sectional models of both the attach and kick ring structure were made up with 0.0265 in. thick stainless steel thin skin covers with thermocouples on them to determine the heating rates. Such models were made up for both the forward and rear faces of the kick ring which has a different configuration on each side. The thin skins were replaced with the alternate phenolic TPS sections cut from flight hardware configuration phenolic parts as supplied by the new vendor. Two tests were performed for each configuration of the attach and kick rings and the samples were exposed to the flow for a duration that gave a heat load equivalent to that obtained in the series of runs made for the current line of phenolic TPS. The samples performed very well with no loss of any phenolic layers. The post-test samples looked better than those used to verify the current phenolic TPS.

Karu, Z. S.↗

Orion EM-2 TPS Radiometer & Spectrometer

Radiometers and Spectrometers embedded as flight instrumentation in a spacecraft's thermal protection system (TPS) can generate data on shock layer radiation to a spacecraft's surface for high-speed trajectories such as lunar return. Robust and accurate flight measurements of shock layer radiation can validate models and inform TPS selection and sizing, to reduce risk and TPS launch mass. Contamination from soot deposition and ablation products in the boundary layer condensing on a cold sapphire radiometer window or blocking open ports interferes with measurements. This poster describes characterization of the optical properties and performance of relevant materials, components, and contaminated windows, as well as investigation of contamination control strategies for the optics.

Ablation↗

Experimental Investigations into the Material Response of Heatshields, High-Temperature Coatings, and Meteors.

Phenolic Impregnated Carbon Ablator (PICA) gained heritage during the Stardust mission and is the baseline thermal protection system (TPS) material for missions to Mars. However, PICA is friable, so a thin layer of siloxane resin (NuSil) is applied to the surface of the heatshield to suppress the spread of particulate matter during cleanroom activities. Until recently, the thermochemical response of NuSil had not been accounted for in material response models. Therefore, a test campaign was executed at the Hypersonic Materials Environmental Test System (HyMETS) facility to investigate the material response of NuSil and provide crucial insight for model development. Finally, salient results from an experiment supported by the Asteroid Threat Assessment Project (ATAP) will be discussed, which probed the material response of the Tamdakht meteorite and notable analogs. Bio: Dr. Bessire is a senior research scientist in the thermal protection system materials branch at the NASA Ames Research Center. His research focuses on developing analytical techniques designed to inform and validate the development of material response models.

Brody K. Bessire↗

Probabalistic Risk Analysis and Thermal Margin Process for an Inflatable Aeroshell

Uncertainties always exist in atmospheric entry aeroheating environments and the thermal response of thermal protection system (TPS) material. These uncertainties are mitigated in the design by ap-plying margin and factors of safety to the TPS. Entry vehicle TPS is often conservatively over-sized for the heat loads that are experienced along the entry trajectory by stacking worst-case scenarios together. Additionally, the current TPS design and margin process used by NASA offers very little insight into the risk of over-temperature during flight and the reliability of the heat shield performance [1,3]. A probabilistic margin process can be used to calculate the amount of TPS margin necessary to survive a given entry heat load at a specified level of risk [2,3,4]. The vehicle’s initial entry state (entry velocity, flight path angle, and entry mass) determines the expected atmospheric entry environmental conditions and resulting heat load that the entry vehicle will experience. If there is flexibility in the entry state, then this process can be used to select an appropriate combination of entry state parameters and TPS size to target a desired reentry reliability. This probabilistic margin process allows engineers to make informed aeroshell design, entry-trajectory design, and TPS performance risk trades while preventing excessive TPS margin from being applied. The probabilistic TPS margin process has been performed to determine TPS thickness and entry heating constraints given an acceptable risk level for the Low Earth Orbit Flight Experiment of an Inflatable Decelerator (LOFTID) flight project. The process is used in a manner to size the entry heat load for a given flexible TPS (FTPS) thickness so that it meets project reliability standards while allowing the FTPS and the underlying inflatable structure (IS) to be pushed to adequately high temperatures. Since the LOFTID project is an experimental flight demonstration, it is de-sired to drive the FTPS and IS to temperatures that cover a large range of their thermal response models’ applicability. This will allow the thermal response models to be better improved and validated post-flight using LOFTID’s extensive instrumentation embedded within the aeroshell. The presentation demonstrates how uncertainty analysis is carried out using an end-to-end Monte Carlo process where three separate Monte Carlo simulations are run in sequence. The first Monte Carlo simulation operates on the entry trajectory model to generate trajectory parameter dispersions that are fed into the second Monte Carlo simulation. The second Monte Carlo simulation operates on the aerothermodynamics model to generate aeroheating parameter dispersions that are fed into the third Monte Carlo simulation. The third Monte Carlo simulation operates on the FTPS material thermal response model to generate the final FTPS/IS thermal response dispersions. The end-to-end Monte Carlo simulation propagates the uncertainties of each model into the next to quantify the resulting uncertainty of the FTPS/IS thermal response. The fractional contributions of the uncertain parameters in the trajectory, aerothermal, and thermal response models to the variance in the FTPS/IS thermal response is determined as a byproduct of the Monte Carlo analysis. The structural uncertainty of the FTPS thermal response model is evaluated by flight relevant ground testing and model error analysis using test measurements. This probabilistic TPS margin process had never been applied to an entry vehicle and it is one of the LOFTID project’s goals to demonstrate its merits.

Steven A. Tobin↗

Thermo Physics Facilities Branch Brochure ARC Jet Complex Fact Sheets, Hypervelocity Free-Flight Aerodynamic Facility Fact Sheets, Ames Vertical Gun Range Fact Sheets

The Ames Arc Jet Complex has a rich heritage of over 40 years in Thermal Protection System (TPS) development for every NASA Space Transportation and Planetary program, including Apollo, Space Shuttle, Viking, Pioneer-Venus, Galileo, Mars Pathfinder,Stardust, NASP,X-33,X-34,SHARP-B1 and B2,X-37 and Mars Exploration Rovers. With this early TPS history came a long heritage in the development of the arc jet facilities. These are used to simulate the aerodynamic heating that occurs on the nose cap, wing leading edges and on other areas of the spacecraft requiring thermal protection. TPS samples have been run in the arc jets from a few minutes to over an hour,from one exposure to multiple exposures of the same sample, in order t o understand the TPS materials response to a hot gas flow environment (representative of real hyperthermal environments experienced in flight). The Ames Arc l e t Complex is a key enabler for customers involved in the three major areas of TPS development: selection, validation, and qualification. The arc jet data are critical for validating TPS thermal models, heat shield designs and repairs, and ultimately for flight qualification.

Fretter, E. F.↗

Data Parallel Line Relaxation (DPLR) Code User Manual: Acadia - Version 4.01.1

Data-Parallel Line Relaxation (DPLR) code is a computational fluid dynamic (CFD) solver that was developed at NASA Ames Research Center to help mission support teams generate high-value predictive solutions for hypersonic flow field problems. The DPLR Code Package is an MPI-based, parallel, full three-dimensional Navier-Stokes CFD solver with generalized models for finite-rate reaction kinetics, thermal and chemical non-equilibrium, accurate high-temperature transport coefficients, and ionized flow physics incorporated into the code. DPLR also includes a large selection of generalized realistic surface boundary conditions and links to enable loose coupling with external thermal protection system (TPS) material response and shock layer radiation codes.

Wright, Michael J.↗

Analysis of PICA-NuSil at the Hypersonic Materials Environmental Test System (HyMETS)

Phenolic Impregnated Carbon Ablator (PICA) gained heritage during theStardust mission and is the baseline thermal protection system (TPS) material for missions to Mars. However, PICA is friable, so a thin layer of siloxane resin (NuSil) is applied to the surface of the heat shield to suppress the spread of particulate matter during cleanroom activities. Until recently, the thermochemical response of NuSil had not been accounted for in material response models. Therefore, a test campaign was executed at the Hypersonic Materials Environmental Test System (HyMETS) facility to investigate the material response of NuSil and provide crucial insight for model development. Finally, salient results from an experiment supported by the AsteroidThreat Assessment Project (ATAP) will be discussed, which probed the materialr esponse of the Tamdakht meteorite and notable analogs.

Brody K. Bessire↗

Simulation of etch pit formation through active sites in carbon fiber micro-structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation as a result of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gasses such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

Simulation of Etch Pit Formation in DSMC Through Active Sites in Carbon Fiber Micro-Structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation because of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gases such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various detailed surface reaction mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗