Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic methods”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Resolving the Coverage Dependence of Surface Reaction Kinetics with Machine Learning and Automated Quantum Chemistry Workflows

Microkinetic models for catalytic systems require estimation of many thermodynamic and kinetic parameters that can be calculated for isolated species and transition states using ab initio methods. However, the presence of nearby coadsorbates on the surface can dramatically alter these thermodynamic and kinetic parameters causing them to be dependent on species coverage fractions. As there are combinatorially many coadsorbed configurations on the surface, computing the coverage dependence of these parameters is far less straightforward. We present a framework for generating and applying machine learning models to predict coverage-dependent parameters for microkinetic models. Our toolkit enables automatic calculation and evaluation of coadsorbed configurations allowing us to sample 2,000 coadsorbed adsorbates and transition states (TSs) for a diverse set of 9 reactions on Cu(111), a challenging surface, with four possible coadsorbates. This dataset was then used to train subgraph isomorphic decision trees (SIDTs) to predict the stability and association energy of configurations. We were able to achieve mean absolute errors (MAEs) of 0.106 eV on adsorbates, 0.172 eV on TSs, and due to natural error cancellation in SIDTs for relative properties, 0.130 eV on reaction energies and 0.180 eV on activation barriers. In conclusion, we describe how to use these models to predict coverage-dependent corrections for adsorbates and TSs and demonstrate on H*, HO*, and O* comparing the generated SIDT model with an iteratively refined version.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗
NASA NTRS1976

Comparison of modal test results - Multipoint sine versus single-point random

The Mariner Jupiter/Saturn (MJS) spacecraft was subjected to the traditional multipoint sine dwell (MPSD) modal test using 111 accelerometer channels, and also to single-point random (SPR) testing using 26 accelerometer channels, and the two methods are compared according to cost, schedule, and technical criteria. A measure of comparison between the systems was devised in terms of the cumulative difference in the kinetic energy distribution of the common accelerometers. The SPR and MPSD method show acceptable agreement with respect to frequencies and mode damping. The merit of the SPR method is that the excitation points are minimized and the test article can be committed to other uses while data analysis is performed. The MPSD approach allows validity of the data to be determined as the test progresses. Costs are about the same for the two methods.

Leppert, E. L.↗
NASA NTRS1971

Wave-wave interactions due to scattering by electrons.

The kinetic wave equation which describes the nonlinear wave-particle interaction in a plasma is considered, and a method which uses the picture of quantized plasmons interacting with particles for the description of nonlinear wave-particle interactions is briefly described. In this method an assumption is made of the Markoffian character of the equation. It is shown that the Markoffian assumption can be justified at least for the case when the plasma is close to a stable stationary state. The diagram method developed by Nishikawa (1966) is used to derive an explicit expression for the kinetic equation. The application of the result to the case of interaction between an electron-wave and an ion-wave is discussed.

Nishikawa, K.↗
NASA NTRS1990

Near-wall turbulence model and its application to fully developed turbulent channel and pipe flows

A near-wall turbulence model and its incorporation into a multiple-timescale turbulence model are presented. The near-wall turbulence model is obtained from a k-equation turbulence model and a near-wall analysis. In the method, the equations for the conservation of mass, momentum, and turbulent kinetic energy are integrated up to the wall, and the energy transfer and the dissipation rates inside the near-wall layer are obtained from algebraic equations. Fully developed turbulent channel and pipe flows are solved using a finite element method. The computational results compare favorably with experimental data. It is also shown that the turbulence model can resolve the overshoot phenomena of the turbulent kinetic energy and the dissipation rate in the region very close to the wall.

Kim, S.-W.↗
DOE OSTI2024

A graph neural network-state predictive information bottleneck (GNN-SPIB) approach for learning molecular thermodynamics and kinetics

Molecular dynamics simulations offer detailed insights into atomic motions but face timescale limitations. Enhanced sampling methods have addressed these challenges but even with machine learning, they often rely on pre-selected expert-based features. Here, in this work, we present a Graph Neural Network-State Predictive Information Bottleneck (GNN-SPIB) framework, which combines graph neural networks and the state predictive information bottleneck to automatically learn low-dimensional representations directly from atomic coordinates. Tested on three benchmark systems, our approach predicts essential structural, thermodynamic and kinetic information for slow processes, demonstrating robustness across diverse systems. The method shows promise for complex systems, enabling effective enhanced sampling without requiring pre-defined reaction coordinates or input features.

Zou, Ziyue↗
DOE OSTI2025

The Role of Kinetic Instabilities and Waves in Collisionless Magnetic Reconnection

Magnetic reconnection converts magnetic field energy into particle energy by breaking and reconnecting magnetic field lines. Magnetic reconnection is a kinetic process that generates a wide variety of kinetic waves via wave-particle interactions. Kinetic waves have been proposed to play an important role in magnetic reconnection in collisionless plasmas by, for example, contributing to anomalous resistivity and diffusion, particle heating, and transfer of energy between different particle populations. These waves range from below the ion cyclotron frequency to above the electron plasma frequency and from ion kinetic scales down to electron Debye length scales. This review aims to describe the progress made in understanding the relationship between magnetic reconnection and kinetic waves. We focus on the waves in different parts of the reconnection region, namely, the diffusion region, separatrices, outflow regions, and jet fronts. Particular emphasis is placed on the recent observations from the Magnetospheric Multiscale (MMS) spacecraft and numerical simulations, which have substantially increased the understanding of the interplay between kinetic waves and reconnection. Some of the ongoing questions related to waves and reconnection are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Determination of a suitable molar absorption coefficient (ε) for lignin analysis of fibrous plants using the CASA method

In this study, the benchmarking and applicability of the CASA method for the analysis of fibrous plants, such as flax, hemp, and jute, is evaluated. Lignin is a phenolic biopolymer present in plant cell walls and is composed of three primary monomeric units, designated G, S, and H. Various factors, such as the genetic variability, influence the relative proportions of these units in plant samples. Recently, the cysteine-assisted sulfuric acid (CASA) method has been introduced as a rapid method for the quantification of lignin in wood samples. The aim of this study is to establish a suitable molar absorption coefficient (ε, L·g –1 ·cm –1 ) to adapt the CASA method for use with annual plant fibers. This investigation was motivated by the technical advantages of the CASA method, including higher throughput, lower reaction temperatures, and ecological benefits due to the absence of carcinogenic, mutagenic, or reprotoxic (CMR) substances and the need for minimal sample quantities. In this method, lignin solubilization is facilitated by cysteine, which is an amino acid that enhances the reaction kinetics, using a one-hour incubation period. However, as with any spectrophotometric technique, the CASA method depends on a molar absorption coefficient (ε) that varies according to the ratio of the aromatic units within the polymer. To evaluate the suitability of CASA for quantifying lignin in economically significant plant fibers, we investigate the impact of the unit ratios on the accuracy of ε. The results are compared with those of two widely recognized lignin analytical methods: the Klason method, which is a gravimetric reference method, and the acetyl bromide soluble lignin method, which is the most commonly used spectrophotometric approach. The final ε obtained in this study reveals a relative difference in lignin content ranging from –8 % to +9 % based on a comparison between the CASA and Klason methods across different industrial hemp varieties. As a result, using known G:S ratios in annual fibrous plants, ε can be estimated from our results.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2024

Differential roles of kinetic on- and off-rates in T-cell receptor signal integration revealed with a modified Fab’-DNA ligand

Antibody-derived T-cell receptor (TCR) agonists are commonly used to activate T cells. While antibodies can trigger TCRs regardless of clonotype, they bypass native T cell signal integration mechanisms that rely on monovalent, membrane-associated, and relatively weakly binding ligand in the context of cellular adhesion. Commonly used antibodies and their derivatives bind much more strongly than native peptide major histocompatibility complex (pMHC) ligands bind their cognate TCRs. Because ligand dwell time is a critical parameter that tightly correlates with physiological function of the TCR signaling system, there is a general need, both in research and therapeutics, for universal TCR ligands with controlled kinetic binding parameters. To this end, we have introduced point mutations into recombinantly expressed α-TCRβ H57 Fab to modulate the dwell time of monovalent Fab binding to TCR. When tethered to a supported lipid bilayer via DNA complementation, these monovalent Fab’-DNA ligands activate T cells with potencies well-correlated with their TCR binding dwell time. Single-molecule tracking studies in live T cells reveal that individual binding events between Fab'-DNA ligands and TCRs elicit local signaling responses closely resembling native pMHC. The unique combination of high on- and off-rates of the H57 R97L mutant enables direct observations of cooperative interplay between ligand binding and TCR-proximal condensation of the linker for activation of T cells, which is not readily visualized with pMHC. This work provides insights into how T cells integrate kinetic information from TCR ligands and introduces a method to develop affinity panels for polyclonal T cells, such as cells from a human patient.

Science & Technology - Other Topics↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
NASA NTRS1979

Matrix perturbation for analytical model improvement

A matrix perturbation method is proposed to calculate the Jacobian matrix and to compute the new eigendata for the parameter estimation procedure. The advantages of the method are the applicability to large complex structures without knowing the analytical expressions for the mass and stiffness matrices, and a cost effective approach for the re-computation of the eigendata. This method also allows the use of other measurements such as modal forces, kinetic energy distribution, and strain energy distributions in the estimation procedure. A realistic sample problem is presented to demonstrate the effectiveness of the proposed method.

Chen, J. C.↗
NASA NTRS1982

Kinetics and mechanism of the oxidation process of two-component Fe-Al alloys

The oxidation process of two-component Fe-Al alloys containing up to 7.2% Al and from 18 to 30% Al was studied. Kinetic measurements were conducted using the isothermal gravimetric method in the range of 1073-1223 K and 1073-1373 K for 50 hours. The methods used in studies of the mechanism of oxidation included: X-ray microanalysis, X-ray structural analysis, metallographic analysis and marker tests.

Przewlocka, H.↗
NASA NTRS2002

Semiclassical Calculation of Reaction Rate Constants for Homolytical Dissociations

There is growing interest in extending organometallic chemical vapor deposition (OMCVD) to III-V materials that exhibit large thermal decomposition at their optimum growth temperature, such as indium nitride. The group III nitrides are candidate materials for light-emitting diodes and semiconductor lasers operating into the blue and ultraviolet regions. To overcome decomposition of the deposited compound, the reaction must be conducted at high pressures, which causes problems of uniformity. Microgravity may provide the venue for maintaining conditions of laminar flow under high pressure. Since the selection of optimized parameters becomes crucial when performing experiments in microgravity, efforts are presently geared to the development of computational OMCVD models that will couple the reactor fluid dynamics with its chemical kinetics. In the present study, we developed a method to calculate reaction rate constants for the homolytic dissociation of III-V compounds for modeling OMCVD. The method is validated by comparing calculations with experimental reaction rate constants.

Beatriz H. Cardelino↗
NASA NTRS2006

BLIMPK/Streamline Surface Catalytic Heating Predictions on the Space Shuttle Orbiter

This paper describes the results of an analysis of localized catalytic heating effects to the U.S. Space Shuttle Orbiter Thermal Protection System (TPS). The analysis applies to the High-temperature Reusable Surface Insulation (HRSI) on the lower fuselage and wing acreage, as well as the critical Reinforced Carbon-Carbon on the nose cap, chin panel and the wing leading edge. The object of the analysis was to use a modified two-layer approach to predict the catalytic heating effects on the Orbiter windward HRSI tile acreage, nose cap, and wing leading edge assuming localized highly catalytic or fully catalytic surfaces. The method incorporated the Boundary Layer Integral Matrix Procedure Kinetic (BLIMPK) code with streamline inputs from viscous Navier-Stokes solutions to produce heating rates for localized fully catalytic and highly catalytic surfaces as well as for nominal partially catalytic surfaces (either Reinforced Carbon-Carbon or Reaction Cured Glass) with temperature-dependent recombination coefficients. The highly catalytic heating results showed very good correlation with Orbiter Experiments STS-2, -3, and -5 centerline and STS-5 wing flight data for the HRSI tiles. Recommended catalytic heating factors were generated for use in future Shuttle missions in the event of quick-time analysis of damaged or repaired TPS areas during atmospheric reentry. The catalytic factors are presented along the streamlines as well as a function of stagnation enthalpy so they can be used for arbitrary trajectories.

Marichalar, Jeremiah J.↗
NASA NTRS2014

Evaluation of Human and Anthropomorphic Test Device Finite Element Models under Spaceflight Loading Conditions

In an effort to develop occupant protection standards for future multipurpose crew vehicles, the National Aeronautics and Space Administration (NASA) has looked to evaluate the test device for human occupant restraint with the modification kit (THOR‐K) anthropomorphic test device (ATD) in relevant impact test scenarios. With the allowance and support of the National Highway Traffic Safety Administration, NASA has performed a series of sled impact tests on the latest developed THOR‐K ATD. These tests were performed to match test conditions from human volunteer data previously collected by the U.S. Air Force. The objective of this study was to evaluate the THOR‐K finite element (FE) model and the Total HUman Model for Safety (THUMS) FE model with respect to the tests performed. These models were evaluated in spinal and frontal impacts against kinematic and kinetic data recorded in ATD and human testing. Methods: The FE simulations were developed based on recorded pretest ATD/human position and sled acceleration pulses measured during testing. Predicted responses by both human and ATD models were compared to test data recorded under the same impact conditions. The kinematic responses of the models were quantitatively evaluated using the ISO‐metric curve rating system. In addition, ATD injury criteria and human stress/strain data were calculated to evaluate the risk of injury predicted by the ATD and human model, respectively. Results: Preliminary results show well‐correlated response between both FE models and their physical counterparts. In addition, predicted ATD injury criteria and human model stress/strain values are shown to positively relate. Kinematic comparison between human and ATD models indicates promising biofidelic response, although a slightly stiffer response is observed within the ATD. Conclusion: As a compliment to ATD testing, numerical simulation provides efficient means to assess vehicle safety throughout the design process and further improve the design of physical ATDs. The assessment of the THOR‐K and THUMS FE models in a spaceflight testing condition is an essential first step to implementing these models in the computational evaluation of spacecraft occupant safety. Promising results suggest future use of these models in the aerospace field.

Putnam, Jacob P.↗
NASA NTRSDate not supplied

Direct Simulation Monte Carlo Studies of the Gas Sampling for the VATMOS-SR Mission Concept

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept which would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere, and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus, with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a full picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of atmospheric volatile signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process in order to determine how measured isotopic ratios of noble gases in the sample compare to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept, and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed at Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield a comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters as well as the modeling parameters will be performed.

direct simulation Monte Carlo↗
NASA NTRSDate not supplied

The VATMOS-SR Mission Concept: DSMC Studies of the Gas Sampling

VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept that would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a complete picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O, and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of volatile atmospheric signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process to determine how measured isotopic ratios of noble gases in the sample compared to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open-source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed on Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters, as well as the modeling parameters, will be performed.

direct simulation Monte Carlo↗