Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction Mechanism”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Experimental space weathering of regolith material

Significant advances in the understanding of space weathering processes were recently reported. Submicroscopic iron blebs were produced in lunar simulant glass and natural terrestrial minerals by high-temperature reduction in controlled atmosphere furnaces. These experiments altered the samples' optical properties and microtextures so that they resembled those of extremely mature lunar soil. The results contributed to a revised model for natural reduction in the regolith. Subsequently, supporting results were obtained by reduction of lunar samples. Research to date has focused on reduction of three lunar surface components: basalt, pyroclastic glass, and mare soil. An extensive set of H reduction experiments with simulants has led to a detailed understanding of reaction mechanisms and kinetics. Reduction experiments using lunar basalt were recently conducted by Carbotek. Reduced samples from these test were analyzed. Reduction experiments on lunar glass 74220 were run at temperatures of 900-1100 C. Reduction efficiency of volcanic glass proved to be a function of the sample's FeO abundance and reaction temperature. We also reduced mare soil 75061 at temperatures of 900-1050 C. Partial reduction of FeO in olivine and pyroxene occurred, but was slower and less complete than reduction of ilmenite. Our experiments on simulants and lunar samples have indicated that the most readily reduced phases in the regolith are ilmenite and glass. Based on initial tests with simulants we proposed refinements to the accepted model for space weathering of the regolith. The impact of a micrometeoroid flash heats and melts and ejects from the impact point a small volume of soil that contains trapped solar wind H and C. Reduction occurs rapidly, while the melt volume is still in motion. When a droplet encounters unmelted soil, it envelopes cold mineral grains. The melt is chilled rapidly. Our analyses of experimentally reduced lunar basalt, glass, and mare soil support the proposed model for space weathering.

Mckay, D. S.↗
NASA NTRS2017

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗
NASA NTRS1993

Evaporation kinetics of Mg2SiO4 crystals and melts from molecular dynamics simulations

Computer simulations based on the molecular dynamics (MD) technique were used to study the mechanisms and kinetics of free evaporation from crystalline and molten forsterite (i.e., Mg2SiO4) on an atomic level. The interatomic potential employed for these simulations reproduces the energetics of bonding in forsterite and in gas-phase MgO and SiO2 reasonably accurately. Results of the simulation include predicted evaporation rates, diffusion rates, and reaction mechanisms for Mg2SiO4(s or l) yields 2Mg(g) + 20(g) + SiO2(g).

Kubicki, J. D.↗
NASA NTRS1993

Measurements of soot formation and hydroxyl concentration in near critical equivalence ratio premixed ethylene flame

The testing and development of existing global and detailed chemical kinetic models for soot formation requires measurements of soot and radical concentrations in flames. A clearer understanding of soot particle inception relies upon the evaluation and refinement of these models in comparison with such measurements. We present measurements of soot formation and hydroxyl (OH) concentration in sequences of flat premixed atmospheric-pressure C2H4/O2/N2 flames and 80-torr C2H4/O2 flames for a unique range of equivalence ratios bracketting the critical equivalence ratio (phi(sub c)) and extending to more heavily sooting conditions. Soot volume fraction and number density profiles are measured using a laser scattering-extinction apparatus capable of resolving a 0.1 percent absorption. Hydroxyl number density profiles are measured using laser-induced fluorescence (LIF) with broadband detection. Temperature profiles are obtained from Rayleigh scattering measurements. The relative volume fraction and number density profiles of the richer sooting flames exhibit the expected trends in soot formation. In near-phi(sub c) visibility sooting flames, particle scattering and extinction are not detected, but an LIF signal due to polycyclic aromatic hydrocarbons (PAH's) can be detected upon excitation with an argon-ion laser. A linear correlation between the argon-ion LIF and the soot volume fraction implies a common mechanistic source for the growth of PAH's and soot particles. The peak OH number density in both the atmospheric and 80-torr flames declines with increasing equivalence ratio, but the profile shape remains unchanged in the transition to sooting, implying that the primary reaction pathways for OH remain unchanged over this transition. Chemical kinetic modeling is demonstrated by comparing predictions using two current reaction mechanisms with the atmospheric flame data. The measured and predicted OH number density profiles show good agreement. The predicted benzene number density profiles correlate with the measured trends in soot formation, although anomalies in the benzene profiles for the richer and cooler sooting flames suggest a need for the inclusion of benzene oxidation reactions.

Inbody, Michael Andrew↗
NASA NTRSDate not supplied

Modeling Electrolytic O2 Recovery from Metabolic CO2 for Advanced Closed Loop Life Support Systems in Extraterrestrial Human Missions

The International Space Station (ISS) is currently equipped with a complex, heavy, and power consuming system that recovers approximately 50% of O2 from metabolic CO2. Future long duration human missions to the Moon and Mars will necessitate a sustainable and highly efficient metabolic oxygen recovery system capable of yielding a minimum of 75% O2 recovery. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to significantly increase current metabolic O2 recovery efficiency, expand mission sustainability, and reduce complexity of the system. The novel design combines CO2 conversion to O2 along with C2H4 as byproduct and water electrolysis (currently conducted in two separate units) into a single compact unit that runs at standard conditions and is theoretically capable of generating O2 with a theoretical maximum metabolic CO2 conversion of 73% while consuming less than metabolic water. This paper presents a comprehensive multi-physic 3D model developed at MSFC on CO2 conversion to O2 and C2H4 at standard conditions via MFECR. The 3D spatial domain of the model is a replica of the actual MFECR’s 3D drawing generated for the MFECR fabrication and operated to recover O2 from CO2 yielding C2H4 as byproduct. Electrochemical (EC) physics that includes EC multicomponent reaction mechanisms, mass transport, and electrical current density distributions is coupled in the model with all the other physics phenomena involved in the MFECR’s process, such as two-phase flow, free and porous fluid regimes, multicomponent mass transfer, heat transfer, and DC electrical current generation along with Joule heating effect. The EC reaction sections of the MFECR consists of two porous gas diffusion electrodes (GDE) and an electrolyte serpentine channel sandwiched in the middle. The CO2 feeds the cathode serpentine channel and part of the O2 product is fed back to the anode serpentine chamber. An alkaline solution feeds the electrolyte serpentine chamber wetting the GDEs of both, the anode and cathode allowing the OH- ionic transport between them. The EC reactions in the cathode’s GDE yield C2H4 from CO2 and H2 from water while the EC reaction in the anode’s GDE yields O2. The authors will present in this paper the validation of the model using experimental data and the utilization of the validated model in building a reliable simulator that will not only assist the authors on the MFECR design but also the optimization of its operation in the ISS and future spatial human missions.

Jesus A Dominguez↗
NASA NTRS1986

Shock-tube pyrolysis of chlorinated hydrocarbons - Formation of soot

Soot formation in pyrolysis of chlorinated methanes, their mixtures with methane, and chlorinated ethylenes were studied behind reflected shock waves by monitoring the attenuation of an He-Ne laser beam. An additional single-pulse shock-tube study was conducted for the pyrolysis of methane, methyl chloride, and dichloromethane. The experiments were performed at temperatures 1300-3000 K, pressures of 0.4-3.6 bar, and total carbon atom concentrations of 1-5 x 10 to the 17th atoms cu cm. The amounts of soot produced in the pyrolysis of chlorinated hydrocarbons are larger than that of their nonchlorinated counterparts. The sooting behavior and product distribution can be generally explained in terms of chlorine-catalyzed chemical reaction mechanisms. The pathway to soot from chlorinated methanes and ethylenes with high H:Cl ratio proceeds via the formation of C2H, C2H2, and C2H3 species. For chlorinated hydrocarbons with low H:Cl ratio, the formation of C2 and its contribution to soot formation at high temperatures becomes significant. There is evidence for the importance of CHCl radical and its reactions in the pyrolysis of dichloromethane.

Frenklach, M.↗
NASA NTRS1977

High performance N2O4/amine elements: Blowapart

The mechanisms controlling hypergolic propellant reactive stream separation (RRS) were studied and used to develop design criteria for injectors free from both steady state RSS and cyclic propellant stream separation. This was accomplished through the analysis of single element injectors using N204/MMH propellants; the injectors were representative of the space shuttle orbit maneuvering engine and space tug applications. A gas phase/surface reaction mechanism which controls RSS was identified. Injector design criteria were developed, which defined a critical chamber pressure for those operating conditions above which RSS occurs. It was found that the amount of interfacial surface area at impingement is controlled by injector hydraulics.

Lawver, B. R.↗
NASA NTRSDate not supplied

INCORPORATION OF NITROGEN INTO ORGANICS PRODUCED BY FISCHER-TROPSCH TYPE CHEMISTRY

Laboratory simulations have demonstrated that hydrothermal systems have the potential to produce a range of organic compounds through Fischer-Tropsch type (FTT) chemistry. The distribution of products depends on several factors, including the abundance and composition of feedstock molecules, reaction temperature, and the physical and chemical characteristics of catalytic materials included in the reactions. The majority of studies per-formed to date have focused solely on inclusion of CO2 or CO and H2 as the carbon, oxygen and hydrogen sources, which limits the possible products to hydrocarbons, alcohols and carboxylic acids. A few studies have included nitrogen in the form of ammonia, which led to the production of amino acids and nitrogenous bases. Although these demonstrations provide compelling evidence that FTT reactions can produce compounds of interest for the origins of life, such reactions have been conducted under a very limited range of conditions and the synthetic reaction mechanisms have generally not been well-characterized. We have begun a series of laboratory experiments that will incorporate a range of precursor molecules in varying compositions to determine how these variables affect the relative amounts and speciation of life-essential elements in organic molecules produced under FTT conditions. In the present work, we focus on systems containing C, H, O and N. Reaction mixtures contained CO (~60-80 mM), 150 mM ammonium chloride, and water with a native iron catalyst, which also served as a source of H2 through the reduction of water. Three initial experiments were conducted at 110, 154 and 185 °C and heated for several days. Following the reaction, samples were analyzed by gas chromatography-mass spectrometry, direct analysis in real time mass spectrometry and liquid chromatography-mass spectrometry.

A. S. Burton↗
NASA NTRS1989

Modeling of carbon monoxide oxidation kinetics over NASA carbon dioxide laser catalysts

The recombination of CO and O2 formed by the dissociation of CO2 in a sealed CO2 laser discharge zone is examined. Conventional base-metal-oxide catalysts and conventional noble-metal catalysts are not effective in recombining the low O2/CO ratio at the low temperatures used by the lasers. The use of Pt/SnO2 as the noble-metal reducible-oxide (NMRO), or other related materials from Group VIIIA and IB and SnO2 interact synergistically to produce a catalytic activity that is substantially higher than either componet separately. The Pt/SnO2 and Pd/SnO2 were reported to have significant reaction rates at temperatures as low as -27 C, conditions under which conventional catalysts are inactive. The gas temperature range of lasers is 0 + or - 40 C. There are three general ways in which the NMRO composite materials can interact synergistically: one component altering the properties of another component; the two components each providing independent catalytic functions in a complex reaction mechanism; and the formation of catalytic sites through the combination of two components at the atomic level. All three of these interactions may be important in low temperature CO oxidation over NMRO catalysts. The effect of the noble metal on the oxide is discussed first, followed by the effect of the oxide on the noble metal, the interaction of the noble metal and oxide to form catalytic sites, and the possible ways in which the CO oxidation reaction is catalyzed by the NMRO materials.

Herz, Richard K.↗
NASA NTRS2005

Updated Chemical Kinetics and Sensitivity Analysis Code

An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗
NASA NTRSDate not supplied

Overview of NASA Research Activities in Shock Layer Kinetics and Radiation

Shock layer radiation is an important heating mechanism for entry probes to most planetary destinations and re-entry to Earth from beyond low Earth orbit. Our understanding of shock layer radiation phenomena has improved tremendously over the last decade thanks to NASA investment in fundamental radiation research through Entry Systems Modeling’s (ESM) Shock Layer Kinetics and Radiation (SLKR) task. This talk will overview the recent activities within SLKR, including validation of models through ground testing, flight instrumentation and remote observations; development of first principles ab initio calculations for reaction mechanisms and spectroscopic databases and advancing numerical and computational methods for prediction of radiation in reacting hypersonic flows. 1 Sr. Research Scientist, Aerothermodynamics Branch, and AIAA Associate Fellow.

Brett A. Cruden↗
NASA NTRS1992

Preliminary results of an experimental study of the interactions of basalt glass and a water vapor atmosphere: Implications for weathering on Mars

Models of weathering processes on the surface of Mars invoke hydrothermal alteration as the primary mechanism responsible for clay formation. Previous experimental studies of basalt glass interactions with water under hydrothermal conditions demonstrate that phyllosilicates and zeolites are the primary alteration minerals. Gas-solid weathering is thought to be less thermodynamically favorable and relatively unimportant, however, the experimental alteration of basalt glass under vapor ion conditions (large SA/V) can result in the formation of clay minerals, zeolites, and hydrated calcium silicates. We have undertaken a study of the reacted layers formed on basalt glasses experimentally altered under vapor hydration conditions to resolve this issue. High SA/V ratios promote reaction product buildup in solution and promote alteration mineral formation. It was previously shown that these reaction conditions promote weathering processes similar to those found in nature for tektite glasses, rhyolitic glasses, and basalt glasses. The hydration of basalt glass is described for 100 percent relative humidity experiments at temperatures of 150, 175, and 200 C for up to 400 days. Preliminary characterization of the alteration layers with analytical electron and scanning microscopy suggest that the reaction mechanism includes precipitation of a smectite clay on the outermost surface of the glass. Between the clay and the unreacted glass is an amorphous gel-like phase (palagonite). Our results provide an experimental basis for proposing that basalt glass interactions with vapor atmospheres can be an important source of clay minerals.

Mazer, J. J.↗
NASA NTRS1999

Laminar Dust Flames: A Program of Microgravity and Ground Based Studies at McGill

Fundamental knowledge of heterogeneous combustion mechanisms is required to improve utilization of solid fuels (e.g. coal), safe handling of combustible dusts in industry, and solid propulsion systems. The objective of the McGill University research program on dust combustion is to obtain a reliable set of data on basic combustion parameters for dust suspensions (i.e. laminar burning velocity, flame structure, quenching distance, flammability limits, etc.) over a range of particle sizes, dust concentrations, and types of fuel. This set of data then permits theoretical models to be validated and, when necessary, new models to be developed to describe the detailed reaction mechanisms and transport processes. Microgravity is essential to the generation of a uniform dust suspension of arbitrary particle size and concentration. When particles with a characteristic size on the order of tens of microns are suspended, they rapidly settle in a gravitational field. To maintain a particulate in suspension for time duration adequate to carry out combustion experiments invariably requires continuous convective flow in excess of the gravitational settling velocity (which is comparable with and can even exceed the dust laminar burning velocity). This makes the experiments turbulent in nature and thus renders it impossible to study laminar dust flames. Even for small particle sizes on the order of microns, a stable laminar dust flow can be maintained only for relatively low dust concentrations at normal gravity conditions. High dust loading leads to gravitational instability of the dust cloud and to the formation of recirculation cells in the dust suspension in a confined volume, or to the rapid sedimentation of the dense dust cloud, as a whole, in an unconfined volume. Many important solid fuels such as carbon and boron also have low laminar flame speeds (of the order of several centimeters per second). Convection that occurs in combustion products due to buoyancy disrupts the low speed dust flames and makes observation of such flames at normal gravity difficult.

Goroshin, Sam↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗
NASA NTRSDate not supplied

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↗