Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallel in time”

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 325 records · Page 18

The direction of heat flow in an expanding current free stellar atmosphere

The energy budgets of the corona are substantially influenced by losses attributed to thermal conduction losses down the transition region, towards the chromosphere. Although Fourier law for such a plasma regime is now universally criticized as inappropriate, the direction of this heat flow seems to be certain. The purpose of this paper is to point out, based on conservation laws, the direction in which the heat must flow in this regime, provided it is known that there is a time independent expansion devoid of parallel currents. It will be argued that the direction of the heat flux is from the base of the transition region up into the corona, rather than the opposite. Under the circumstances mentioned above, the direction of the heat flow should be opposite to the gradient of the electron pressure gradient. Although the transition region is modeled as isobaric, this is only approximately true. Such a likely scenario makes the heat flow almost certainly to flow from the lower temperature base of the TR up into the incipient corona if a steady state current free expansion is to be maintained.

Scudder, J. D.↗

Temperature and Strain-Rate Effects on Low-Cycle Fatigue Behavior of Alloy 800H

The effects of strain rate (4 x 10(exp -6) to 4 x 10(exp -3)/s) and temperature on the Low-Cycle Fatigue (LCF) behavior of alloy 800H have been evaluated in the range 750 C to 950 C. Total axial strain controlled LCF tests were conducted in air at a strain amplitude of +/- 0.30 pct. LCF life decreased with decreasing strain rate and increasing temperature. The cyclic stress response behavior showed a marked variation with temperature and strain rate. The time- and temperature- dependent processes which influence the cyclic stress response and life have been identified and their relative importance assessed. Dynamic strain aging, time-dependent deformation, precipitation of parallel platelets of M(23)C6 on grain boundaries and incoherent ledges of twins, and oxidation were found to operate depending on the test conditions. The largest effect on life was shown by oxidation processes.

Rao, K. Bhanu Sankara↗

Analysis of the DFP/AFCS Systems for Compensating Gravity Distortions on the 70-Meter Antenna

This paper presents the theoretical computations showing the expected performances for both systems. The basic analysis tool is a Physical Optics reflector analysis code that was ported to a parallel computer for faster execution times. There are several steps involved in computing the RF performance of the various systems. 1 . A model of the RF distortions of the main reflector is required. This model is based upon measured holography maps of the 70-meter antenna obtained at 3 elevation angles. The holography maps are then processed (using an appropriate gravity mechanical model of the dish) to provide surface distortion maps at all elevation angles. 2. From the surface distortion maps, ray optics is used to determine the theoretical shape of the DFP that will exactly phase compensate the distortions. 3. From the theoretical shape and a NASTRAN mechanical model of the plate, the actuator positions that generate a surface that provides the best RMS fit to the theoretical model are selected. Using the actuator positions and the NASTRAN model provides an accurate description of the actual mirror shape. 4. Starting from the mechanical drawings of the feed, a computed RF feed pattern is generated. This pattern is expanded into a set of spherical wave modes so that a complete near field analysis of the reflector system can be obtained. 5. For the array feed, the excitation coefficients that provide the maximum gain are computed using a phase conjugate technique. The basic experimental geometry consisted of a dual shaped 70-meter antenna system; a refocusing ellipse, a DFP and an array feed system. To provide physical insight to the systems performance, focal plane field plots are presented at several elevations. Curves of predicted performance are shown for the DFP system, monopulse tracking system, AFCS and combined DFP/AFCS system. The calculated results show that the combined DFP/AFCS system is capable of recovering the majority of the gain lost due to gravity distortion.

Imbriale, William A.↗

Cell shape, cytoskeletal mechanics, and cell cycle control in angiogenesis

Capillary endothelial cells can be switched between growth and differentiation by altering cell-extracellular matrix interactions and thereby, modulating cell shape. Studies were carried out to determine when cell shape exerts its growth-regulatory influence during cell cycle progression and to explore the role of cytoskeletal structure and mechanics in this control mechanism. When G0-synchronized cells were cultured in basic fibroblast growth factor (FGF)-containing defined medium on dishes coated with increasing densities of fibronectin or a synthetic integrin ligand (RGD-containing peptide), cell spreading, nuclear extension, and DNA synthesis all increased in parallel. To determine the minimum time cells must be adherent and spread on extracellular matrix (ECM) to gain entry into S phase, cells were removed with trypsin or induced to retract using cytochalasin D at different times after plating. Both approaches revealed that cells must remain extended for approximately 12-15 h and hence, most of G1, in order to enter S phase. After this restriction point was passed, normally 'anchorage-dependent' endothelial cells turned on DNA synthesis even when round and in suspension. The importance of actin-containing microfilaments in shape-dependent growth control was confirmed by culturing cells in the presence of cytochalasin D (25-1000 ng ml-1): dose-dependent inhibition of cell spreading, nuclear extension, and DNA synthesis resulted. In contrast, induction of microtubule disassembly using nocodazole had little effect on cell or nuclear spreading and only partially inhibited DNA synthesis. Interestingly, combination of nocodazole with a suboptimal dose of cytochalasin D (100 ng ml-1) resulted in potent inhibition of both spreading and growth, suggesting that microtubules are redundant structural elements which can provide critical load-bearing functions when microfilaments are partially compromised. Similar synergism between nocodazole and cytochalasin D was observed when cytoskeletal stiffness was measured directly in living cells using magnetic twisting cytometry. These results emphasize the importance of matrix-dependent changes in cell and nuclear shape as well as higher order structural interactions between different cytoskeletal filament systems for control of capillary cell growth during angiogenesis.

Non-NASA Center↗

A Parallel Incompressible Navier-Stokes Solver With a Parallel Multigrid Elliptic Kernal

The paper describes numerical algorithms and parallel implementations of a time-dependent, incompressible Navier-Stokes flow solver and a multigrid elliptic solver, which is also used as a computation kernal in the flow solver. The implemented solvers are numerically stable and computationally efficient, and they scale well to a large number of processors for problems with moderate granularity.

solver Navier-Stokes solver multigrid elliptic sol↗

Three-Dimensional Engineered High Fidelity Normal Human Lung Tissue-Like Assemblies (TLA) as Targets for Human Respiratory Virus Infections

Unlike traditional two-dimensional (2D) cell cultures, three-dimensional (3D) tissue-like assemblies (TLA) (Goodwin et aI, 1992, 1993, 2000 and Nickerson et aI. , 2001,2002) offer high organ fidelity with the potential to emulate the infective dynamics of viruses and bacteria in vivo. Thus, utilizing NASA micro gravity Rotating Wall Vessel (RWV) technology, in vitro human broncho-epithelial (HBE) TLAs were engineered to mimic in vivo tissue for study of human respiratory viruses. These 3D HBE TLAs were propagated from a human broncho-tracheal cell line with a mesenchymal component (HBTC) as the foundation matrix and either an adult human broncho-epithelial cell (BEAS-2B) or human neonatal epithelial cell (16HBE140-) as the overlying element. Resulting TLAs share several characteristic features with in vivo human respiratory epithelium including tight junctions, desmosomes and cilia (SEM, TEM). The presence of epithelium and specific lung epithelium markers furthers the contention that these HBE cells differentiate into TLAs paralleling in vivo tissues. A time course of infection of these 3D HBE TLAs with human respiratory syncytial virus (hRSV) wild type A2 strain, indicates that virus replication and virus budding are supported and manifested by increasing virus titer and detection of membrane-bound F and G glycoproteins. Infected 3D HBE TLAs remain intact for up to 12 days compared to infected 2D cultures that are destroyed in 2-3 days. Infected cells show an increased vacuolation and cellular destruction (by transmission electron microscopy) by day 9; whereas, uninfected cells remain robust and morphologically intact. Therefore, the 3D HBE TLAs mimic aspects of human respiratory epithelium providing a unique opportunity to analyze, for the first time, simulated in vivo viral infection independent of host immune response.

Goodwin, T. J.↗

Curved Waveguide Based Nuclear Fission for Small, Lightweight Reactors

The focus of the presented work is on the creation of a system of grazing incidence, supermirror waveguides for the capture and reuse of fission sourced neutrons. Within research reactors, neutron guides are a well known tool for directing neutrons from the confined and hazardous central core to a more accessible testing or measurement location. Typical neutron guides have rectangular, hollow cross sections, which are crafted as thin, mirrored waveguides plated with metal (commonly nickel). Under glancing angles with incoming neutrons, these waveguides can achieve nearly lossless transport of neutrons to distant instruments. Furthermore, recent developments have created supermirror surfaces which can accommodate neutron grazing angles up to four times as steep as nickel. A completed system will form an enclosing ring or spherical resonator system to a coupled neutron source for the purpose of capturing and reusing free neutrons to sustain and/or accelerate fission. While grazing incidence mirrors are a known method of directing and safely using neutrons, no method has been disclosed for capture and reuse of neutrons or sustainment of fission using a circular waveguide structure. The presented work is in the process of fabricating a functional, highly curved, neutron supermirror using known methods of Ni-Ti layering capable of achieving incident reflection angles up to four times steeper than nickel alone. Parallel work is analytically investigating future geometries, mirror compositions, and sources for enabling sustained fission with applicability to the propulsion and energy goals of NASA and other agencies. Should research into this concept prove feasible, it would lead to development of a high energy density, low mass power source potentially capable of sustaining fission with a fraction of the standard critical mass for a given material and a broadening of feasible materials due to reduced rates of release, absorption, and non-fission for neutrons. This advance could be applied to direct propulsion through guided fission products or as a secondary energy source for high impulse electric propulsion. It would help meet national needs for highly efficient energy sources with limited dependence on fossil fuels or conflict materials, and it would improve the use of low grade fissile materials which would help reduce national stockpiles and waste.

Coker, Robert↗

Delta-Doped Back-Illuminated CMOS Imaging Arrays: Progress and Prospects

In this paper, we report the latest results on our development of delta-doped, thinned, back-illuminated CMOS imaging arrays. As with charge-coupled devices, thinning and back-illumination are essential to the development of high performance CMOS imaging arrays. Problems with back surface passivation have emerged as critical to the prospects for incorporating CMOS imaging arrays into high performance scientific instruments, just as they did for CCDs over twenty years ago. In the early 1990's, JPL developed delta-doped CCDs, in which low temperature molecular beam epitaxy was used to form an ideal passivation layer on the silicon back surface. Comprising only a few nanometers of highly-doped epitaxial silicon, delta-doping achieves the stability and uniformity that are essential for high performance imaging and spectroscopy. Delta-doped CCDs were shown to have high, stable, and uniform quantum efficiency across the entire spectral range from the extreme ultraviolet through the near infrared. JPL has recently bump-bonded thinned, delta-doped CMOS imaging arrays to a CMOS readout, and demonstrated imaging. Delta-doped CMOS devices exhibit the high quantum efficiency that has become the standard for scientific-grade CCDs. Together with new circuit designs for low-noise readout currently under development, delta-doping expands the potential scientific applications of CMOS imaging arrays, and brings within reach important new capabilities, such as fast, high-sensitivity imaging with parallel readout and real-time signal processing. It remains to demonstrate manufacturability of delta-doped CMOS imaging arrays. To that end, JPL has acquired a new silicon MBE and ancillary equipment for delta-doping wafers up to 200mm in diameter, and is now developing processes for high-throughput, high yield delta-doping of fully-processed wafers with CCD and CMOS imaging devices.

imagers↗

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗

Development of a Detailed Surface Chemistry Framework in DSMC

A generalized finite-rate surface chemistry framework incorporating a comprehensive list of reaction mechanisms is developed and implemented into the Direct Simulation Monte Carlo (DSMC) solver SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer). The various mechanisms include adsorption, desorption, Eley-Rideal (ER), and several types of Langmuir-Hinshelwood (LH) mechanisms. The approach is to stochastically model the various competing reactions occurring on a set of active sites. Both gas-surface (e.g., adsorption, ER) and pure-surface (e.g., desorption) reaction mechanisms are incorporated, and the framework also includes catalytic or surface altering mechanisms involving the participation of the bulk-phase species (e.g., bulk carbon atoms). Marschall and MacLean developed a general formulation in which multiple phases and surface sites are used and a similar convention is adopted in the current work. Expressions for the microscopic parameters of reaction probabilities (for gas-surface reactions) and frequencies (for pure-surface reactions) that are required for DSMC are derived from the surface properties and macroscopic parameters such as rate constants, sticking coefficients, etc. The energy and angular distributions of the products are specified according to the reaction type and input parameters. This framework also presents physically consistent procedures to accurately compute the reaction probabilities and frequencies in the case of multiple reactions. The result is a modeling tool with a wide variety of surface reactions characterized via user-specified reaction rate constants, surface properties and parameters.

Surface Chemistry↗

One Giant Leap for Womankind: An Aerodynamic Study of the SLS Rocket ft. Pressure-Sensitive Paint

Combining computer vision techniques, high-speed cameras, pressure-sensitive paint and a transonic wind tunnel, NASA studies unsteady aerodynamic forces on the Space Launch System rocket with unprecedented temporal and spatial resolution. NASA’s most powerful supercomputer, Pleiades, enables parallel processing and real-time visualization to investigate buffet forces and aeroacoustic physics.

Lucy Tang↗

Coupling Carbon Oxidation and Surface Recession in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that are often made out of carbon-based reinforcement and a polymeric matrix. They endure high-temperature oxidation and surface recession when re-entering Earth’s atmosphere. Ablation is the result of many coupled and competing thermal, mechanical, and chemical phenomena, and it is difficult to isolate the role of each on the overall degradation of the TPS. Here we develop an ablation model for material recession coupled explicitly to finite rate carbon oxidation in complex microstructures. In this work, Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al, was adapted to tally surface reactions and other surface data on a per-grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also develop robust algorithms which handle the evolution of the flow cells and solid material regions, including split cells (flow cell divided in two by a solid surface). Finally, we demonstrate our implicit chemistry model for 2D and 3D geometries by producing reaction statistics and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗

Implementation of active sites to capture pitting of oxidizing carbon materials in DSMC.

In this work we demonstrate a newly developed capability to capture pitting of carbon fibers in DSMC simulations, specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code [1]. State-of-the-art reactive surface models in DSMC compute collision dependent carbon consumption rates (usually through desorption of CO) based on a set of surface reactions that has been derived from molecular beam experiments [2]. The reactivity on each carbon surface element is constant in those models, such that the carbon surface recedes uniformly as a result of ablation. However, it is well known that in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale [3]. These defective sites have a much higher reactivity than the average sites (2-3 orders of magnitude) and are first to react during ablation leading to its removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these ”active” sites (as shown in Fig. 1). In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. 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 [4]. Within this framework, we have implemented the capability of a single surface having multiple site sets with different reactivities. Using this feature, we can simulate the presence of active sites on carbon surfaces, whose reactivity is much greater than an average site as a result of defects. We have implemented the active site fraction as a property of surface elements within SPARTA, which 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 are able to capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

Implementation of Active Sites in DSMC to Capture Pitting of Oxidizing Carbon Materials

In this work we demonstrate a newly developed capability to capture pitting of carbon fibers in DSMC simulations, specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. State-of-the-art reactive surface models in DSMC compute collision dependent carbon consumption rates (usually through desorption of CO) based on a set of surface reactions that has been derived from molecular beam experiments. The reactivity on each carbon surface element is constant in those models, such that the carbon surface recedes uniformly as a result of ablation. However, it is well known that in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. These defective sites have a much higher reactivity than the average sites (2-3 orders of magnitude) and are first to react during ablation leading to its removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these ”active” sites. In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. 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. Within this framework, we have implemented the capability of a single surface having multiple site sets with different reactivities. Using this feature, we can simulate the presence of active sites on carbon surfaces, whose reactivity is much greater than an average site as a result of defects. We have implemented the active site fraction as a property of surface elements within SPARTA, which 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 are able to capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗

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↗

DSMC 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 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↗

Computational Study of Oxidative Etch Pitting in FiberForm and the Effect on Its Material Properties

Erosion of carbon due to oxidation does not occur uniformly but through the formation of localized etch pits because of active surface sites. These active sites are formed due to the presence of atomic defects on the carbon surface, and have much higher reactivity compared to average non-defective sites. Thus, these active sites are first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these “active” sites. In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. In order to understand the influence of these etch pits on the material properties of carbon fiber microstructures, we have developed a new capability within direct simulation Monte Carlo (DSMC) to capture the etch pit formation process. This capability is developed within the DSMC code SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) and can model the material removal in presence of active sites leading to the formation of etch pits. The focus of the current work will be to study the effect of etch pits on the material properties of FiberForm, a commonly used base material within many thermal protection system materials (TPS). The microstructure of virgin FiberForm obtained directly from X-ray microtomography experiments is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported within the Porous Microstructure Analysis (PuMA) software and various material properties such as elasticity, thermal conductivity, and permeability are computed. The variation of these properties as a result of the complex evolution of the surface topology due to etch pit formation is studied and analyzed. Furthermore, the effect of pitting is compared to the case of shrinking fibers, which has been the standard for modelling ablation of carbon structures; and significant differences are observed. Thus, such a physically realistic modeling of material removal through the formation of etch pits will be helpful in predicting the degradation of carbon-based TPS more accurately during oxidation; as well as other mechanisms such as spallation, which involves the removal of chunks of material into the flow due to etch pit growth. This will ultimately improve our understanding of the failure modes in these materials due to ablation.

Carbon Ablators↗