Results for “Direct simulation Monte Carlo” — Page 6

NASA NTRSDate not supplied

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

NASA NTRS2016

DSMC Shock Simulation of Saturn Entry Probe Conditions

This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center. To directly compare to the experimental results, the DSMC simulations are post-processed through the NEQAIR line-by-line radiation code. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 kms shock waves are obtained at 0.2 and 0.1 Torr respectively and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.

Saturn↗

NASA NTRS2003

Simulations of Ground and Space-Based Oxygen Atom Experiments

A low-earth orbit (LEO) materials erosion scenario and the ground-based experiment designed to simulate it are compared using the direct-simulation Monte Carlo (DSMC) method. The DSMC model provides a detailed description of the interactions between the hyperthermal gas flow and a normally oriented flat plate for each case. We find that while the general characteristics of the LEO exposure are represented in the ground-based experiment, multi-collision effects can potentially alter the impact energy and directionality of the impinging molecules in the ground-based experiment. Multi-collision phenomena also affect downstream flux measurements.

Finchum, A.↗

NASA NTRS1986

The diffusion of individual molecules within a gas

The Direct Simulation Monte Carlo method is used to study the positional history of the individual molecules in a gas that is homogeneous at the macroscopic level and is in Maxwellian equilibrium at the microscopic level. The behavior at small times is characterized by 'persistence of velocity' effects, and a 'random walk' type of dispersal occurs over a longer timescale. It is shown that the rate of dispersal can be directly related to the self-diffusion coefficient. In addition, the diffusion coefficients are obtained directly from one-dimensional calculations, and the local Knudsen number at which the Chapman-Enskog theory breaks down is determined. Results are presented for both simple gases and gas mixtures.

Bird, G. A.↗

NASA NTRS2016

DSMC Shock Simulation of Saturn Entry Probe Conditions

This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at the NASA Ames Research Center. The DSMC simulations are post-processed through the NEQAIR line-by-line radiation code to compare directly to the experimental results. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 km/s shock waves are obtained at 0.2 and 0.1 Torr, respectively, and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.

Higdon, Kyle J.↗

NASA NTRS1993

Molecular gas dynamics applied to low-thrust propulsion

The Direct Simulation Monte Carlo method is currently being applied to study flowfields of small thrusters, including both the internal nozzle and the external plume flow. The DSMC method is employed because of its inherent ability to capture nonequilibrium effects and proper boundary physics in low-density flow that are not readily obtained by continuum methods. Accurate prediction of both the internal and external nozzle flow is important in determining plume expansion which, in turn, bears directly on impingement and contamination effects.

Zelesnik, Donna↗

NASA NTRS1994

Particle kinetic simulation of high altitude hypervelocity flight

Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration - dissociation - recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.

Boyd, Iain↗

NASA NTRS1994

Particle kinetic simulation of high altitude hypervelocity flight

Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration-dissociation recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.

Boyd, Iain↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

Simulation of Oxidative Etch Pit Formation and Growth on FiberForm

Erosion of carbon surfaces 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 the first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective, and increasing their reactivity, thus leading to localized carbon removal around these “active” sites. In this manner, these highly reactive defects serve as nucleation sites for the formation and growth of etch pits, with 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 the 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, the precursor substrate of the PICA Thermal Protection System (TPS) material. The microstructure of virgin FiberForm obtained directly from X-ray microtomography scans is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported into the Porous Microstructure Analysis (PuMA) software and various material properties such as thermal conductivity, elasticity, and permeability are computed. The variation of these properties because of the complex evolution of the surface topology due to the formation of etch pits is studied and analyzed. Furthermore, the effect of pitting is compared to the case of uniform radial shrinking of fibers, which has been the standard for modelling ablation of carbon structures, and significant differences are observed. Thus, a physically realistic model 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 the growth of etch pits. This will ultimately improve our understanding of the failure modes in these materials due to ablation.

DSMC↗

NASA NTRSDate not supplied

Simulation of Oxidative Etch Pit Formation and Growth on FiberForm

Erosion of carbon surfaces 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 the first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective, and increasing their reactivity, thus leading to localized carbon removal around these “active” sites. In this manner, these highly reactive defects serve as nucleation sites for the formation and growth of etch pits, with 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 the 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, the precursor substrate of the PICA Thermal Protection System (TPS) material. The microstructure of virgin FiberForm obtained directly from X-ray microtomography scans is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported into the Porous Microstructure Analysis (PuMA) software and various material properties such as thermal conductivity, elasticity, and permeability are computed. The variation of these properties because of the complex evolution of the surface topology due to the formation of etch pits is studied and analyzed. Furthermore, the effect of pitting is compared to the case of uniform radial shrinking of fibers, which has been the standard for modelling ablation of carbon structures, and significant differences are observed. Thus, a physically realistic model 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 the growth of etch pits. This will ultimately improve our understanding of the failure modes in these materials due to ablation.

DSMC↗

NASA NTRSDate not supplied

A Parallelized Oxidation-Driven Surface Recession Framework in DSMC Code, SPARTA

Spacecrafts rely on ablative thermal protection systems (TPS) made of composites consisting of a carbon-based reinforcement and a polymeric matrix. These materials are designed to withstand high-temperature oxidation and surface recession during re-entry into the Earth's atmosphere. However, ablation occurs due to a complex interplay of thermal, mechanical, and chemical factors, making it challenging to determine the individual impact of each on the TPS's overall degradation. In this study, we have developed an ablation model that can leverage a finite rate carbon oxidation model to predict material recession and surface states more accurately. 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 can perform both gas-surface and pure-surface reactions and is now adapted to tally 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 briefly discuss some unique challenges associated with parallelizing this dynamic surface state and geometry. Finally, we analyze the performance of this parallelized implicit chemistry model with simple 2D and 3D benchmark cases by producing surface state statistics, area changes over time, and visualization across a range of surface temperatures and processors with and without load-balancing.

DSMC↗

NASA NTRS1988

Direct simulation of three-dimensional hypersonic flow about intersecting blunt wedges

The general three-dimensional direct simulation Monte Carlo method is used to study the hypersonic flow around two blunt wedges that intersect at a 90 deg angle. Results are obtained for the transitional flow regimes found at 85 and 100 km altitude with a reentry velocity of 7.5 km/s. The disturbance field in front of the double-wedge body is found to be larger than that produced by a single wedge. Surface pressures and flow densities are higher near the wedge intersection, whereas surface heating and shear streses are greater at locations removed from the corner. Results also show that three-dimensional flow structure occurs only near the wedge corner, and that the flow monotonically approaches the limiting two-dimensional wedge flow case in the spanwise direction.

Celenligil, M. Cevdet↗

NASA NTRS2004

Remote Sensing of Icy Galilean Moon Surface and Atmospheric Composition Using Low Energy (1 eV-4 keV) Neutral Atom Imaging

We describe a low energy neutral atom imager suitable for composition measurements Europa and other icy Galilean moons in the Jovian magnetosphere. This instrument employs conversion surface technology and is sensitive to either neutrals converted to negative ions, neutrals converted to positive ions and the positive ions themselves depending on the power supply. On a mission such as the Jupiter Icy Moons Orbiter (JIMO), two back-to-back sensors would be flown with separate power supplies fitted to the neutral atom and iodneutral atom sides. This will allow both remote imaging of 1 eV < E < 4 keV neutrals from icy moon surfaces and atmospheres, and in situ measurements of ions at similar energies in the moon ionospheres and Jovian magnetospheric plasma. The instrument provides composition measurements of the neutrals and ions that enter the spectrometer with a mass resolution dependent on the time-of-flight subsystem and capable of resolving molecules. The lower energy neutrals, up to tens of eV, arise from atoms and molecules sputtered off the moon surfaces and out of the moon atmospheres by impacts of more energetic (keV to MeV) ions from the magnetosphere. Direct Simulation Monte Carlo (DSMC) models are used to convert measured neutral abundances to compositional distributions of primary and trace species in the sputtered surfaces and atmospheres. The escaping neutrals can also be detected as ions after photo- or plasma-ionization and pickup. Higher energy, keV neutrals come from charge exchange of magnetospheric ions in the moon atmospheres and provide information on atmospheric structure. At the jovicentric orbits of the icy moons the presence of toroidal gas clouds, as detected at Europa's orbit, provide M e r opportunities to analyze both the composition of neutrals and ions originating from the moon surfaces, and the characteristics of magnetospheric ions interacting with neutral cloud material. Charge exchange of low energy ions near the moons, and directional distributions of the resultant neutrals, allow indirect global mapping of magnetic field structures around the moons. Temporal variation of the magnetic structures can be linked to induced magnetic fields associated with subsurface oceans.

Collier, M. R.↗

NASA NTRS2001

Chemical and Flowfield Modeling for Enhanced Analysis of Contamination Experiments

This paper describes the application of a new Direct Simulation Monte Carlo (DSMC) code, the Molecular Beam Simulator (MBS), which is designed to analyze laboratory scale molecular beam-surface (and crossed-beam) experiments. The MBS is primarily intended to model experiments associated with spacecraft contamination effects, but it can also be used to simulate a variety of surface chemistry and reactive flow measurements. The MBS code is fully three-dimensional, includes a wide-range of chemical processes, and can model one or multiple pulsed (non-steady) sources. As an example application of the MBS code, a fast, pulsed, oxygen atom-surface experiment which examines the chemistry behind erosion of graphite by oxygen atoms is analyzed. Unsteady DSMC simulations show that experimental observations of excited molecular states after the pulse has hit the surface are consistent with two distinct chemical mechanisms: a direct one where the excited molecules are formed on the surface, and a two-step mechanism where ground state molecules formed on the surface are collisionally excited after they leave the surface by trailing oxygen atoms in the pulse. Further DSMC calculations suggest experiments which can distinguish between these mechanisms.

Braunstein, Matthew↗

NASA NTRS1993

Simulated rarefied entry of the Galileo probe into the atmosphere of Jupiter

Flow properties and aerodynamics are computed with a direct simulation Monte Carlo (DSMC) method for rarefied entry of the Galileo Probe into the atmosphere of Jupiter. Accurate predictions of vehicle drag coefficients are needed in order to assess atmospheric properties from the onboard Atmosphere Structure Experiment where highly-sensitive accelerometers will measure the drag force to within 10-6 barr during the initial entry phase at high altitudes. The corresponding flow rarefraction extends from the free molecule limit to the near continuum transition regime (Re less than 1000). Simulation results indicate that C(sub D) varies from 2.1 at the free molecule limit down to 1.6 at Re(infinity) = 1,000. Temperatures, densities, and internal energies throughout the flow field were also computed at each altitude ranging from 735 km to 353 km above the 1 barr level in the Jovian atmosphere. Surface heating and temperatures of the probe were computed directly in the DSMC code by assuming radiative equilibrium. Material response was re-asssessed accurately during entry by accounting for conductivity, heat capacity, and pyrolysis which led to surface material mass efflux several times that of the freestream mass influx. The simulation also accounted for the quantum nature of the rotational energy mode of the dominant atmospheric species H2 through partial internal excitation in the freestream gas.

Hass, Brian L.↗

NASA NTRS2002

Simulations of Ground and Space-Based Oxygen Atom Experiments

Fast, pulsed atomic oxygen sources are a key tool in ground-based investigations of spacecraft contamination and surface erosion effects. These technically challenging ground-based studies provide a before and after picture of materials under low-earth-orbit (LEO) conditions. It would be of great interest to track in real time the pulsed flux from the source to the surface sample target and beyond in order to characterize the population of atoms and molecules that actually impact the surface and those that make it downstream to any coincident detectors. We have performed simulations in order to provide such detailed descriptions of these ground-based measurements and to provide an assessment of their correspondence to the actual LEO environment. Where possible we also make comparisons to measured fluxes and erosion yields. To perform the calculations we use a detailed description of a measurement beam and surface geometry based on the W, pulsed apparatus at Montana State University. In this system, a short pulse (on the order of 10 microseconds) of an O/O2 beam impacts a flat sample about 40 cm downstream and slightly displaced &om the beam s central axis. Past this target, at the end of the beam axis is a quadrupole mass spectrometer that measures the relative in situ flux of 0102 to give an overall normalized erosion yield. In our simulations we use the Direct Simulation Monte Carlo (DSMC) method, and track individual atoms within the atomic oxygen pulse. DSMC techniques are typically used to model rarefied (few collision) gas-flows which occur at altitudes above approximately 110 kilometers. These techniques are well suited for the conditions here, and multi-collision effects that can only be treated by this or a similar technique are included. This simulation includes collisions with the surface and among gas atoms that have scattered from the surface. The simulation also includes descriptions of the velocity spread and spatial profiles of the O/O2 beam obtained from separate measurements. These computations use basic engineering models for the gas-gas and gas-surface scattering and focus on the influence of multi-collision effects. These simulations characterize many important quantities of interest including the actual flux of atoms that reach the surface, the energy distribution of this flux, as well as the direction of the velocity of the flux that strikes the surface. These quantities are important in characterizing the conditions which give rise to measured surface erosion. The calculations also yield time- snapshots of the pulse as it impacts and flows around the surface. These snapshots reveal the local environment of gas near the surface for the duration of the pulse. We are also able to compute the flux of molecules that travel downstream and reach the spectrometer, and we characterize their velocity distribution. The number of atoms that reach the spectrometer can in fact be influenced by the presence of the surface due to gas-gas collisions from atoms scattered h m the surface, and it will generally be less than that with the surface absent. This amounts to an overall normalization factor in computing erosion yields. We discuss these quantities and their relationship to the gas-surf$ce interaction parameters. We have also performed similar calculations corresponding to conditions (number densities, temperatures, and velocities) of low-earth orbit. The steady-state nature and lower overall flux of the actual space environment give rise to differences in the nature of the gas-impacts on the surface from those of the ground-based measurements using a pulsed source.

Minton, T. K.↗