SEARCH · Search NASA

Results for “defect evolution”

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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Thin‐Film Ferroelectrics

Abstract Over the last 30 years, the study of ferroelectric oxides has been revolutionized by the implementation of epitaxial‐thin‐film‐based studies, which have driven many advances in the understanding of ferroelectric physics and the realization of novel polar structures and functionalities. New questions have motivated the development of advanced synthesis, characterization, and simulations of epitaxial thin films and, in turn, have provided new insights and applications across the micro‐, meso‐, and macroscopic length scales. This review traces the evolution of ferroelectric thin‐film research through the early days developing understanding of the roles of size and strain on ferroelectrics to the present day, where such understanding is used to create complex hierarchical domain structures, novel polar topologies, and controlled chemical and defect profiles. The extension of epitaxial techniques, coupled with advances in high‐throughput simulations, now stands to accelerate the discovery and study of new ferroelectric materials. Coming hand‐in‐hand with these new materials is new understanding and control of ferroelectric functionalities. Today, researchers are actively working to apply these lessons in a number of applications, including novel memory and logic architectures, as well as a host of energy conversion devices.

Fernandez, Abel↗

DOE OSTI2024

Evolution of dislocation loops in irradiated α-Uranium: An atomistically-informed cluster dynamics investigation

An atomistically-informed mean field cluster dynamics model has been developed to investigate the nucleation and growth of prismatic loops in irradiated a-Uranium. TEM analysis of neutron irradiated a-Uranium shows the evolution of self-interstitial atom and vacancy loops on (010) and (100) crystallographic planes, respectively, resulting in an anisotropic lattice swelling of its face-centered orthorhombic crystal. To provide model parameters, the crystallography of loops and the binding energy of point defects to these loops were studied using an angular dependent EAM potential and classical molecular dynamics (MD) simulations. Furthermore, using the bond-boost hyperdynamics method, the anisotropic diffusion of interstitials and vacancies in a-Uranium was investigated. Here, the mechanisms of point defect diffusion and the associated migration energies were reported and compared with previous DFT studies. The energetics and kinetic quantities mentioned above were used in the cluster dynamics model to predict the cluster density as a function of dose rate, dose and temperature and the results were compared to the reported neutron irradiation experiments. The model predictions reveal an accumulation of small sized vacancy loops along with a population of large and growing self-interstitial loops, which closely corresponds to the TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Non-Equilibrium Synthesis of Highly Active Nanostructured, Oxygen-Incorporated Amorphous Molybdenum Sulfide HER Electrocatalyst

Molybdenum sulfide emerged as promising hydrogen evolution reaction (HER) electrocatalyst thanks to its high intrinsic activity, however its limited active sites exposure and low conductivity hamper its performance. To address these drawbacks, the non-equilibrium nature of pulsed laser deposition (PLD) is exploited to synthesize self-supported hierarchical nanoarchitectures by gas phase nucleation and sequential attachment of defective molybdenum sulfide clusters. Furthermore, the physics of the process are studied by in situ diagnostics and correlated to the properties of the resulting electrocatalyst. The as-synthesized architectures have a disordered nanocrystalline structure, with nanodomains of bent, defective S-Mo-S layers embedded in an amorphous matrix, with excess sulfur and segregated molybdenum particles. Oxygen incorporation in this structure fosters the creation of amorphous oxide/oxysulfide nanophases with high electrical conductivity, enabling fast electron transfer to the active sites. The combined effect of the nanocrystalline pristine structure and the surface oxidation enhances the performance leading to small overpotentials, very fast kinetics (35.1 mV dec -1 Tafel slope) and remarkable long-term stability for continuous operation up to -1 A cm -2 . This work shows possible new avenues in catalytic design arising from a non-equilibrium technique as PLD and the importance of structural and chemical control to improve the HER performance of MoS-based catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

Kurtosis, skewness, and non-Gaussian cosmological density perturbations

Cosmological topological defects as well as some nonstandard inflation models can give rise to non-Gaussian density perturbations. Skewness and kurtosis are the third and fourth moments that measure the deviation of a distribution from a Gaussian. Measurement of these moments for the cosmological density field and for the microwave background temperature anisotropy can provide a test of the Gaussian nature of the primordial fluctuation spectrum. In the case of the density field, the importance of measuring the kurtosis is stressed since it will be preserved through the weakly nonlinear gravitational evolution epoch. Current constraints on skewness and kurtosis of primeval perturbations are obtained from the observed density contrast on small scales and from recent COBE observations of temperature anisotropies on large scales. It is also shown how, in principle, future microwave anisotropy experiments might be able to reveal the initial skewness and kurtosis. It is shown that present data argue that if the initial spectrum is adiabatic, then it is probably Gaussian, but non-Gaussian isocurvature fluctuations are still allowed, and these are what topological defects provide.

Luo, Xiaochun↗

DOE OSTI2025

Mesoscale modeling of restructuring in high burnup UO 2 fuel

Here, this work aims to simulate the restructuring behavior observed in different regions of high burnup fuel, providing a first-of-its-kind restructuring model for the dark zone and rim region of high-burnup UO 2 fuel. We employed a grand-potential-based phase-field model to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of different initial conditions and different modeling parameters were systematically studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to begin around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Restructuring was demonstrated to be influenced by a combination of initial dislocation densities, burnup rate, subgrain formation rate, and temperature. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. Furthermore, bubble size and distribution do not significantly affect the rate of restructuring. The predicted microstructures are consistent with experimental observations of the restructured regions. Finally, we present a correlation demonstrating the evolution of the restructuring volume fraction as a function of local burnup.

UO2↗

DOE OSTI2021

Catalytic materials for lithium-sulfur batteries: mechanisms, design strategies and future perspective

Lithium-sulfur batteries (LSBs) are attractive candidates for post-lithium-ion battery technologies because of their ultrahigh theoretical energy density and low cost of active cathode materials. However, the commercialization of LSBs remains extremely challenging primarily due to poor cycling performance and safety concerns, which are inherently caused by low conductivity of S 8 and Li 2 S, severe polysulfide shuttling, and high polarization by solid Li 2 S 2 /Li 2 S deposition. Catalytic materials could facilitate the large-scale practical application of LSBs by overcoming all these challenges. In this review, we investigate the sulfur species evolution in LSBs and explore the roles of catalytic materials in charge/discharge processes, highlighting the catalysis of solid S 8 to liquid polysulfides and solid Li 2 S 2 to Li 2 S. Furthermore, we offer systematic strategies from atomic to macro levels, including defect engineering, morphology engineering and catalyst compositing, to enhance catalysis efficiency in terms of sulfur supercooling, fast charge transfer, thiosulfate generation, disulfide bond cleavage, tuneable Li 2 S growth and Li 2 S decomposition enhancement. Finally, the design and availability of the proposed catalytic materials will further advance LSB technology from coin cells and pouch cells to the subsequent commercialization scale.

25 ENERGY STORAGE↗

DOE OSTI2020

Evolution of AISI 304L stainless steel part properties due to powder recycling in laser powder-bed fusion

Selective laser melting (SLM) is a powder-bed fusion process in Additive Manufacturing (AM) that uses a laser beam to selectively fuse layers of powder into near net-shape components with little porosity. However, inconsistencies in the part properties due to the presence of defects in as-built components has hindered the widespread adoption of SLM for industrial applications forcing researchers to study the sources of variation for quality control purposes. A critical area suspected of creating variation in the part properties is the feedstock, where batch-to-batch differences as well as changes in the powder properties with reuse have the potential to affect performance. During processing, laser spatter and condensate eject from the melt pool and deposit into the powder-bed surrounding the parts. These particulates, more commonly known as ejecta, differ morphologically and chemically from the virgin powder compromising its reusability. In this study, 304L stainless steel powder was recycled for a total of 5 times through a systematic approach aimed at accelerating powder degradation to reveal the influence on both the tensile and impact toughness properties. Through analysis of variance (ANOVA), it was found that tensile properties did not change with reuse, while the impact toughness showed a steady decline illustrating the differences in static and dynamic part properties due to powder recycling.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Resolving local structural motifs across the phase evolution of zinc titanates with computational x-ray absorption spectroscopy

Resolving the local structure motifs that characterize phase evolution as a function of composition is a key challenge in structure characterization of complex materials. Here, in this study, we combine first-principles simulations and x-ray absorption near-edge structures (XANES) analysis to gain insights into the structure evolution revealed by measurements across a combinatorial zinc titanate thin film, which was grown with smoothly varying composition over a wide range of the Ti:Zn ratio. Specifically, we propose a cluster blind-signal-separation (cBSS) method for XANES spectral analysis based on a library of the structures and spectra of representative local motifs. In addition to motifs from zinc titanate crystals, two types of Ti-defect models constructed in this study are key to the understanding of the structure characteristics in the Zn-rich region. The cBSS method makes use of both spectral clustering of the simulated site-XANES spectra library and the BSS procedure to construct high-fidelity spectral basis functions from an experimental spectral sequence. The method provides a rigorous measure of the spectral sensitivity and basis completeness. The results of the XANES analysis are corroborated with other experimental modalities, including x-ray diffraction and spectroscopic ellipsometry, to validate the cBSS method. The calculated motif weights resulting from fitting the XANES spectra with the cBSS basis probe the atomic structure characteristics of both crystalline and amorphous phases as a function of the Ti/Zn composition. The insights of the local structure motif evolution are pivotal to the understanding of the nonmonotonic trend in the optical gap, which may lead to potential applications through tuning the optical properties of zinc titanate. The workflow of the XANES spectral analysis developed in this work can be generalized to construct the structure-property relationship in a broad material space.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Systematic discovery of new nano-scale metastable intermetallic eutectic phases in laser rapid solidified Aluminum-Germanium alloy

Laser surface remelting of as-cast Al-Ge eutectic alloy is shown to produce ultrafine lamellar eutectic morphology with interlamellar spacing refined up to ∼60 nm and composed of FCC Al solid solution and unusual Al x Ge y intermetallic phases that do not form during near-equilibrium solidification. The microstructures are characterized and analyzed using a combination of selected area electron diffraction, high-resolution scanning transmission electron microscopy, energy dispersive X-ray spectroscopy to obtain high-resolution elemental maps, and atomistic modeling using density functional theory followed by atomic-scale image simulation. Depending on the local solidification conditions, the crystallography of the Al x Ge y intermetallic phases in the eutectic microstructure is either monoclinic (C 2/c) or monoclinic (P 2 1 ), with high densities of defects in both cases. This is in sharp contrast to the as-cast alloys that showed nominally pure Al and Ge phases with significant solute partitioning and equilibrium FCC and diamond cubic crystal structures, respectively. Corresponding kinetic phase diagrams are proposed to interpret the evolution of nano-lamellar eutectic morphologies with equilibrium Al and metastable Al x Ge y phases, and to explain increased solid solubility in the Al phases manifested by precipitation of ultrafine clusters of Ge. Furthermore, the reasons for the formation of these metastable eutectics under laser rapid solidification are discussed from the perspective of the competitive growth criterion.

Al-Ge eutectic↗

DOE OSTI2024

Mechanisms of ion irradiation induced ordering in amorphous TiO 2 nanotubes: Effects of ion mass and energy

Here, amorphous TiO 2 nanotubes were irradiated in-situ in a transmission electron microscope (TEM) with Kr + ions at energies of 46 keV, 150 keV, and 1 MeV and with 46 keV Xe + ions, to investigate the structural and morphological evolution of the nanotubes under irradiation. At all irradiation conditions, amorphous TiO 2 nanotubes exhibited significant morphological instability, and tended to undergo volumetric swelling with increasing ion counts, often until collapse of the original nanotube structure. Molecular dynamics (MD) simulations confirmed that irradiation-induced defects can explain the observed swelling. Structurally, nanotubes remain amorphous following all Kr + irradiation conditions, but irradiation with 46 keV Xe + leads to the formation of anatase nanocrystallites. Importantly, through systematically varying ion energy and ion species, we try to elucidate the influence of nuclear and electronic stopping power on ion irradiation induced changes. By contextualizing these results within the existing literature, we propose that the observed changes in TiO 2 nanotube morphology and structure could be due to a competition between two mechanisms: (1) disorder-induced swelling and (2) irradiation-induced amorphous-to-crystalline transformation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

In Situ Neutron Scattering Studies on the Oxidation and Reduction of CeO 2 and Pt–CeO 2 Nanorods

The oxygen vacancy structure of ceria plays a key role in its performance as a favored material for catalysis applications. Here, in this work, we develop an understanding of the effects of Pt loading on the structural evolution of ceria nanorods under redox gas environments that mimic real automotive catalytic converters. In situ neutron scattering studies under redox flow reveal that both CeO 2 and Pt–CeO 2 nanorods share a bulk fluorite structure with the presence of surface Frenkel-type oxygen defects. However, Pt–CeO 2 nanorods are more easily reducible than CeO 2 rods as evidenced by an increased concentration of Ce 3+ , determined by NAP-XPS. Importantly, this work finds no evidence of oxygen vacancy ordered surface reconstruction which has been reported in earlier ex situ investigations. Thus, this work highlights the discrepancy between ex situ and in situ structural observations and emphasizes the need for robust in situ investigations of catalysts to develop industrially relevant materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Probing and controlling oxygen impurity diffusion in h -BN semi-bulk crystals

Combining its unique features of ultrawide bandgap (UWBG) and two-dimensional nature, h-BN has been explored for emerging applications such as deep ultraviolet optoelectronic devices and single photon emitters. One of the unusual applications of h-BN is for solid-state neutron detectors by utilizing the property of high thermal neutron capture cross section of B-10 as well as its UWBG properties. Although a record high detection efficiency of 59% has been attained by h-BN detectors, the understanding/minimization of defects and impurities is still needed to further advance the h-BN material and detector technologies. We report metal organic chemical vapor deposition growth and oxygen (O) impurity diffusion in thick h-BN. The diffusion coefficient (D) of O impurities has been measured via the evolution of an oxygen related emission with the etching depth, providing a value of D of ∼ 2 × 10−13 cm2/s at 1450 °C and supporting the interpretation that oxygen in h-BN is a substitutional donor. A multiple-buffer-layer approach was employed to mitigate to a certain degree the issue of oxygen diffusion from sapphire substrate during growth. It was demonstrated that the performance of h-BN neutron detectors fabricated from the wafer incorporating multiple buffer layers was significantly improved, as manifested by the enhanced thermal neutron detection efficiency. The advancement of the crystal growth technology of h-BN semi-bulk crystals creates applications in optoelectronic and power electronic devices utilizing the UWBG semiconductor properties of h-BN, while high efficiency h-BN neutron detectors have the potential to supplant the traditional He-3 gas detectors in various application areas by offering the obvious advantages of UWBG semiconductor technologies.

Physics↗

DOE OSTI2020

Topological structure and dynamics of three-dimensional active nematics

Topological structures are effective descriptors of the nonequilibrium dynamics of diverse many-body systems. For example, motile, point-like topological defects capture the salient features of two-dimensional active liquid crystals composed of energy-consuming anisotropic units. We dispersed force-generating microtubule bundles in a passive colloidal liquid crystal to form a three-dimensional active nematic. Light-sheet microscopy revealed the temporal evolution of the millimeter-scale structure of these active nematics with single-bundle resolution. The primary topological excitations are extended, charge-neutral disclination loops that undergo complex dynamics and recombination events. Our work suggests a framework for analyzing the nonequilibrium dynamics of bulk anisotropic systems as diverse as driven complex fluids, active metamaterials, biological tissues, and collections of robots or organisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Strain-Driven Evolution of Structural and Carrier Recombination Dynamics in Germanium

Strain engineering in group-IV semiconductors, in particular germanium (Ge), offers a promising route for advancing next-generation optoelectronic device performance. High minority carrier lifetime serves as a direct indicator of superior material quality and device efficiency. Several factors govern the minority carrier lifetime in both strained Ge and strained GeSn, such as strain level, epilayer thickness, defects, surface roughness, and Sn alloying. In this study, the impact of biaxial strain on effective minority carrier lifetimes of tensile-strained Ge (0% ..epsilon..-Ge/Al?As, 0.78% ..epsilon..-Ge/In0.115?Ga0.885?As, and 0.93% ..epsilon..-Ge/In0.16?Ga0.84?As) and compressively strained GeSn (-0.53% c-Ge0.97?Sn0.03/Al?As) epilayers, grown by molecular beam epitaxy, was experimentally investigated via the microwave-reflection photoconductive decay (..mu..-PCD) technique as a simultaneous function of thickness and strain. The structural and strain-relaxation properties were analyzed by high-resolution x-ray diffraction and Raman spectroscopy, while surface morphology analysis using atomic force microscopy provided insights into the surface roughness of strained Ge-based epilayers. Defect characterization by cross-section transmission electron microscopy revealed superior crystalline quality of ..epsilon..-Ge in one epilayer, while twin boundary formation was observed in the other. We investigated the carrier recombination dynamics in epitaxially strained Ge by extracting the bulk recombination lifetime and surface recombination velocity using the thickness-dependent minority carrier lifetime. High effective minority carrier lifetimes of approximately 60-340 ns were achieved at high strain levels with minimal dependence on the thickness of the strained Ge-based epilayers by means of the ..mu..-PCD technique. Minority carrier lifetimes in Ge-based epilayers exhibit a quadratic dependence on the biaxial strain, ranging from -0.88% (compressive) to +1.6% (tensile). The analysis of the carrier recombination dynamics in strained Ge-based heterostructures offers a pathway to designing group-IV based optoelectronic devices.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Structural damage response of lanthanum and yttrium aluminate crystals to nuclear collisions and electronic excitation: Threshold assessment of irradiation damage

A comparative analysis is performed on the structural damage response and associated mechanisms in lanthanum aluminate and yttrium aluminate crystals under various irradiation conditions by a combination of experimental and theoretical approaches. Under low-energy Au+ irradiation, the damage accumulation curve shows a higher damage rate for LaAlO 3 crystals than YAlO 3 crystals. The relatively low irradiation tolerance of LaAlO 3 to the action of nuclear collisions is ascribed to the large amorphization cross-section and effective cross-section for defect-stimulated amorphization. Under swift Ar 12+ , Ni 19+ and Kr 17+ irradiation with different ion energies and velocities, the formed highly-disordered/amorphous latent tracks with different morphologies in pristine and predamaged crystals are discussed, and the corresponding electronic energy loss and lattice temperature thresholds are quantitatively determined. Compared to YAlO 3 , LaAlO 3 exhibits lower sensitivity and higher damage tolerance to the electronic energy loss process, attributing to its relatively high recrystallization efficiency during the rapid quenching process. Furthermore, the introduction of lattice defects into LaAlO 3 and YAlO 3 crystals considerably enhances the sensitivity and intensity of thermal spike response to the electronic energy loss, and the induced effective modification of track morphologies demonstrates the synergistic effect between the electronic energy loss and pre-existing defects created by nuclear collisions. In this case, even under the action of electronic energy loss below the threshold, the lattice temperature in the nuclear-collision damaged crystalline system could still meet the criterion for track production. The irradiation energy deposited to atoms and induced lattice temperature evolution discussed in this work provide a deeper insight into the complex processes involved in irradiation-induced latent track behaviors.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Predicting Failure Using Deep Learning SAND Report

Accurate prediction of ductile failure is critical to Sandia’s NW mission, but the models are computationally heavy. The costs of including high-fidelity physics and mechanics that are germane to the failure mechanisms are often too burdensome for analysts either because of the person-hours it requires to input them or because of the additional computational time, or both. In an effort to deliver analysts a tool for representing these phenomena with minimal impact to their existing workflow, our project sought to develop modern data-driven methods that would add microstructural information to business-as-usual calculations and expedite failure predictions. The goal is a tool that receives as input a structural model with stress and strain fields, as well as a machine-learned model, and output predictions of structural response in time, including failure. As such, our project spent substantial time performing high-fidelity, three-dimensional experiments to elucidate materials mechanisms of void nucleation and evolution. We developed crystal-plasticity finite-element models from the experimental observations to enrich the findings with fields not readily measured. We developed engineering length-scale simulations of replicated test specimens to understand how the engineering fields evolve in the presence of fine-scale defects. Finally, we developed deep learning convolutional neural networks, and graph-based neural networks to encode the findings of the experiments and simulations and make forward predictions in time for structural performance. This project demonstrated the power of data-driven methods for model development, which have the potential to vastly increase both the accuracy and speed of failure predictions. These benefits and the methods necessary to develop them are highlighted in this report. However, many challenges remain to implementing these in real applications, and these are discussed along with potential methods for overcoming them.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2019

Myco-Architecture off Planet: Growing Surface Structures at Destination

Our work focused on filling major select key technical knowledge gaps, and technical aspects to be addressed in a Phase II proposal were identified such as the secretion of glues or plastics by the fungi or bacteria to form biocomposites. Specifically, we raised the technical TRL by assessing the growth of in-house mycelial-producing strains on potential food substrates and analyzed the advantages and disadvantages to their use off planet. Growth of the mycelia on sawdust and powdered nutrients including dried, powdered algae, was used as a baseline. We measured growth as a function of temperature, tested the relevant material properties of the mycelia products produced, and measured the mass of input and output volume for production per volume of material. Mycotecture was produced in a bag simulating the proposed mission implementation. Post-production such as heat treatment (as is done with terrestrial mycotecture)was assessed.The material properties of the dried, frozen and baked mycelial outputs included density, strength,thermal insulation and flame retardation. Tensile/compression testing on biocomposite material specimens was performed to obtain three-dimensional elastic constants, as well as to examine deformation and fracture behavior. These mechanical tests were performed on both dormant and activated samples to understand the evolution of the structural material. In addition, experimental data obtained from the mechanical testing was used to build a failure prediction model that accounts for material anisotropy. This aided in future structural design through a quantitative understanding of the mechanical limits of the material. To examine the texture, defects, fracture surfaces of the material, high-resolution microstructural imaging was utilized before and after the mechanical tests. This microstructural analysis informed us about the macro-structural behavior and influenced structural design. Thermal decomposition analysis was also performed tounderstand the thermal limits.

Biocomposites↗

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↗