SEARCH · Engineering Papers

Results for “Point Kinetics”

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 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Pressure and Temperature Dependence of the Reaction of Vinyl Radical with Ethylene

This work reports measurements of absolute rate coefficients and Rice-Ramsperger-Kassel-Marcus (RRKM) master equation simulations of the C2H3 + C2H4 reaction. Direct kinetic studies were performed over a temperature range of 300-700 K and pressures of 20 and 133 mbar. Vinyl radicals (H2C=CH) were generated by laser photolysis of vinyl iodide (C2H31) at 266 nm, and time-resolved absorption spectroscopy was used to probe vinyl radicals through absorption at 423.2 nm. Measurements at 20 mbar are in good agreement with previous determinations at higher temperature. A weighted three-parameter Arrhenius fit to the experimental rate constant at 133 mbar, with the temperature exponent fixed, gives k = (7 +/- 1) x 10(exp -l4) cu cm/molecule/s (T/298 K)(exp 2) exp[-(1430 +/- 70) K/T]. RRKM master equation simulations, based on G3 calculations of stationary points on the C4H7 potential energy surface, were carried out to predict rate coefficients and product branching fractions. The predicted branching to 1-methylallyl product is relatively small under the conditions of the present experiments but increases as the pressure is lowered. Analysis of end products of 248 nm photolysis of vinyl iodide/ethylene mixtures at total pressures between 27 and 933 mbar provides no direct evidence for participation of I -methylallyl.

Ismail, Huzeifa↗

DOE OSTI2024

Spatial decay of electromagnetic waves from hypervelocity impact plasmas

Due to the high kinetic energy of impactors in space, hypervelocity impacts on spacecraft can ablate spacecraft material and produce plasmas. These plasmas have been observed in ground-based experiments to produce electromagnetic pulses (EMPs) and are suspected to be a cause of spacecraft anomalies. The initial EMP that is produced after plasma formation is replicated using a 2D electromagnetic plasma simulation, consisting of a discontinuous Galerkin, particle-in-cell (DG-PIC) framework. The spontaneous formation of EMPs from an initial quasineutral plasma is consistent with previous works and show behavior similar to that of waves emitted by antennas. The EMP emits from the plasma with a frequency less than the plasma frequency associated with the point of maximum density, showing the importance of the density gradient, and the far-field region of the “plasma antenna” is estimated. From the proportionality between plasma frequency and EMP frequency, there is a direct correlation between plasma density and radiation. Finally, as the plasma density decreases, the EMP frequency decreases, which suggests that several observable phenomena, such as optical flashes, radio frequency emission, and microwave emission, all result from the same plasma-related mechanism within a single hypervelocity impact.

42 ENGINEERING↗

DOE OSTI2023

Molecular dynamics study of grain boundary and radiation effects on tritium population and diffusion in zirconium

Tritium population thermodynamics and transport kinetics critically define the tritium storage performance of zirconium tritides that can be used for a variety of nuclear applications including tritium-producing burnable absorber rods. Both thermodynamic and kinetic properties can be sensitive to grain sizes of materials and can be significantly altered by irradiated defects during operation under the reactor environments. A thorough experimental characterization of how these properties evolve under different reactor conditions and different initial grain structures is extremely challenging. Here molecular dynamics simulations are used to investigate tritium population and diffusion in zirconium with and without different planar symmetric and asymmetric tilt grain boundaries and irradiated defects. Here, we found that in addition to trapping tritium, the most significant effect of planar grain boundaries is to increase tritium diffusivity on the boundary plane. Furthermore, fine grain structures are found to mitigate the change of tritium diffusivity due to irradiated point defects as these point defects are likely to migrate to and sink at grain boundaries.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

A Simulation of the Photoionization of H- Together With the Subsequent Tracking of the Liberated Electrons

The Proton Improvement Plan - II (PIP-II) is a new linear accelerator (LINAC) complex being built at Fermilab. It is based on superconducting radiofrequency cavities and will accelerate H- ions to 800 MeV kinetic energy before injection into the existing Booster ring. Measurements of the profile of the beam along the LINAC must be done by non-intercepting methods due to the superconducting cavities. The method chosen is photoionization of a small number of H- by a focused infrared laser, aka laserwire. The number of ionized electrons is measured as a function of laser position within the H- beam. To aid in the design of the collection mechanism, a simulation was written in MATLAB with input from the commercial electromagnetic simulation, CST. This simulation calculates the number and positions of the liberated electrons and tracks them through the magnetic collection and H- beam fields to the collection point. Results from this simulation for various points along the LINAC will be shown.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2017

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗

NASA NTRSDate not supplied

Supercritical Water Oxidation: Testing of Ersatz Wastewater

Supercritical Water Oxidation (SCWO) technology is under consideration by NASA for treatment of wastewater and other wet waste streams. The dramatic changes in water's thermophysical properties near its critical point (374◦C and 22 MPa) result in nearly complete solubility of organics and gases. The absence of inter–phase transport processes and phase separation results in dramatically reduced reaction timescales and the elimination of complications of two–phase transport and processing in reduced gravity. Other advantages are that the product streams are microbially inert,benign, amenable to resource recovery and organic waste. This paper first described the design and construction of tubular SCWO reactor built at NASA Glenn. It consists of separate wastewater and oxidizer inlets, which are independently heated and pressurized to conditions approaching the critical point of water. These flows enter a reactor chamber having a volume of 62 ml. The reactor section is heated to the desired temperature in the range of 450 ◦C – 650 ◦C. A pressure regulator maintains the operating pressure nominally between 27.6 MPa and 28.3 MPa. The reacted fluid is captured downstream of the pressure regulator during an experiment run. A separate vent line may also be used for collection of the effluent. Next, the results of SCWO experiments using both a representative organic (ethanol) as well as a wastewater simulant provided by NASA Ames are presented. Raman analysis of the pre–reacted and post–reacted solutions are presented. Finally, initial efforts on supplementary chemical kinetics modeling of ethanol oxidation with and without inclusion of a nitrogenous compounds (ammonia) are presented.

Michael C Hicks↗

DOE OSTI2023

A Simulation of the Photoionization of H- Together with the Subsequent Tracking of the Liberated Electrons

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

Nonequilibrium information entropy approach to ternary fission of actinides

Ternary fission of actinides probes the state of the nucleus at scission. Light clusters are produced in space and time very close to the scission point. Within the nonequilibrium statistical operator method, a generalized Gibbs distribution is constructed from the information given by the observed yields of isotopes. Using this relevant statistical operator, yields are calculated taking excited states and continuum correlations into account, in accordance with the virial expansion of the equation of state. Furthermore, clusters with mass number A≤10 are well described using the nonequilibrium generalizations of temperature and chemical potentials. Improving the virial expansion, in-medium effects may become of importance in determining the contribution of weakly bound states and continuum correlations to the intrinsic partition function. Yields of larger clusters, which fail to reach this quasiequilibrium form of the relevant distribution, are described by nucleation kinetics, and a saddle-to-scission relaxation time of about 7000 fm/c is inferred. Light-charged particle emission, described by reaction kinetics and virial expansions, may therefore be regarded as a very important tool to probe the nonequilibrium time evolution of actinide nuclei during fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1993

Electronic excitation of ground state atoms by collision with heavy gas particles

Most of the important chemical reactions which occur in the very high temperature air produced around space vehicles as they enter the atmosphere were investigated both experimentally and theoretically, to some extent at least. One remaining reaction about which little is known, and which could be quite important at the extremely high temperatures that will be produced by the class of space vehicles now contemplated - such as the AOTV - is the excitation of bound electron states due to collisions between heavy gas particles. Rates of electronic excitation due to free electron collisions are known to be very rapid, but because these collisions quickly equilibrate the free and bound electron energy, the approach to full equilibrium with the heavy particle kinetic energy will depend primarily on the much slower process of bound electron excitation in heavy particle collisions and the subsequent rapid transfer to free electron energy. This may be the dominant mechanism leading to full equilibrium in the gas once the dissociation process has depleted the molecular states so the transfer between molecular vibrational energy and free electron energy is no longer available as a channel for equilibration of free electron and heavy particle kinetic energies. Two mechanisms seem probable in electronic excitation by heavy particle impact. One of these is the collision excitation and deexcitation of higher electronic states which are Rydberg like. A report, entitled 'Semi-Classical Theory of Electronic Excitation Rates', was submitted previously. This presented analytic expressions for the transition probabilities, assuming that the interaction potential is an exponential repulsion with a perturbation ripple due to the dipole-induced dipole effect in the case of neutral-neutral collisions, and to the ion-dipole interaction in the case of ion-neutral collisions. However the above may be, there is little doubt that excitation of ground state species by collision occurs at the point where the initial and final potentials cross, or at least come very close. Therefore, this mechanism would be applicable to the case where a gas is initially at very low temperature suddenly subjected to high energy heavy particle bombardment. This situation would model the measurement of excitation cross section by molecular beam techniques, for example. The purpose is to report values of cross sections and rate coefficients for collision excitation of ground state atoms estimated with the Landau-Zener transition theory and to compare results with measurement of excitation cross sections for a beam of Hydrogen atoms impacting Argon atom targets. Some very dubious approximations are used, and the comparison with measurement is found less than ideal, but results are at least consistent within order of magnitude. The same model is then applied to the case of N-N atom collisions, even though the approximations then become even more doubtful. Still the rate coefficients obtained are at least plausible in both magnitude and functional form, and as far as I am aware these are the only estimates available for such rate coefficients.

Hansen, C. Frederick↗

NASA NTRS1998

Rejoining and misrejoining of radiation-induced chromatin breaks. IV. Charged particles

We have recently reported the kinetics of chromosome rejoining and exchange formation in human lymphocytes exposed to gamma rays using the techniques of fluorescence in situ hybridization (FISH) and premature chromosome condensation (PCC). In this paper, we have extended previous measurements to cells exposed to charged particles. Our goal was to determine differences in chromatin break rejoining and misrejoining after exposure to low- and high-linear energy transfer (LET) radiation. Cells were irradiated with hydrogen, neon, carbon or iron ions in the LET range 0.3-140 keV/microm and were incubated at 37 degrees C for various times after exposure. Little difference was observed in the yield of early prematurely condensed chromosome breaks for the different ions. The kinetics of break rejoining was exponential for all ions and had similar time constants, but the residual level of unrejoined breaks after prolonged incubation was higher for high-LET radiation. The kinetics of exchange formation was also similar for the different ions, but the yield of chromosome interchanges measured soon after exposure was higher for high-LET particles, suggesting that a higher fraction of DNA breaks are misrejoined quickly. On the other hand, the rate of formation of complete exchanges was slightly lower for densely ionizing radiation. The ratios between the yields of different types of aberrations observed at 10 h postirradiation in prematurely condensed chromosome preparations were dependent on LET. We found significant differences between the yields of aberrations measured in interphase (after repair) and metaphase for densely ionizing radiation. This difference might be caused by prolonged mitotic delay and/or interphase death. Overall, the results point out significant differences between low- and high-LET radiation for the formation of chromosome aberrations.

NASA Center JSC↗

NASA NTRS2016

A Numerical Approach to Estimate the Ballistic Coefficient of Space Debris from TLE Orbital Data

Low Earth Orbit (LEO) is full of space debris, which consist of spent rocket stages, old satellites and fragments from explosions and collisions. As of 2009, more than 21,000 orbital debris larger than 10 cm are known to exist], and while it is hard to track anything smaller than that, the estimated population of particles between 1 and 10 cm in diameter is approximately 500,000, whereas small as 1 cm exceeds 100 million. These objects orbit Earth with huge kinetic energies speeds usually exceed 7 kms. The shape of their orbit varies from almost circular to highly elliptical and covers all LEO, a region in space between 160 and 2,000 km above sea level. Unfortunately, LEO is also the place where most of our active satellites are situated, as well as, International Space Station (ISS) and Hubble Space Telescope, whose orbits are around 400 and 550 km above sea level, respectively.This poses a real threat as debris can collide with satellites and deal substantial damage or even destroy them.Collisions between two or more debris create clouds of smaller debris, which are harder to track and increase overall object density and collision probability. At some point, the debris density couldthen reach a critical value, which would start a chain reaction and the number of space debris would grow exponentially. This phenomenon was first described by Kessler in 1978 and he concluded that it would lead to creation of debris belt, which would vastly complicate satellite operations in LEO. The debris density is already relatively high, as seen from several necessary debris avoidance maneuvers done by Shuttle, before it was discontinued, and ISS. But not all satellites have a propulsion system to avoid collision, hence different methods need to be applied. One of the proposed collision avoidance concepts is called LightForce and it suggests using photon pressure to induce small orbital corrections to deflect debris from colliding. This method is very efficient as seen from theoretical simulations, even few continuous mode 10 kW ground-based lasers, focused by 1.5 m telescopes with adaptive optics, were enough to prevent significant amount of the debris collisions. Simulations were done by propagating all space objects in LEO by 1 year into the future and checking whether the probability of collision was high. For those space objects different ground-based lasers were used to divert them, afterwards collision probabilities were reevaluated. However, the actual accuracy of the LightForce software, which has been developed at NASA AmesResearch Center, depends on the veracity of the input parameters, one of which is the objects ballistic coefficient. It is a measure of bodys ability to overcome air resistance, which has a significant impact on the debris in LEO, and thus it is responsible for the shape of the trajectory of the debris. Having the exact values of the ballistic coefficient would make significantly better collision predictions, unfortunately, we do not know what are the values for most of the objects.In this research, we were working with part of LightForce code, which estimates the ballistic coefficient from ephemerides. Previously used method gave highly inaccurate values, when compared to known objects, and it needed to be changed. The goal of this work was to try out a different method of estimating the ballistic coefficient and to check whether or not it gives noticeable improvements.

Collisions↗

DOE OSTI2024

Low-temperature formation of pyridine and (iso)quinoline via neutral–neutral reactions

Aromatic molecules represent fundamental building blocks in prebiotic chemistry and are contemplated as vital precursors to DNA and RNA nitrogen bases. However, despite the identification of some 300 molecules in extraterrestrial environments, the pathways to pyridine (C 5 H 5 N), pyridinyl (C 5 H 4 N•), and (iso)quinoline (C 9 H 7 N) – the simplest represen-tative of mono and bicyclic aromatic molecule carrying nitrogen – are elusive. Here, we afford compelling evidence on the gas-phase formation of methylene amidogen (H 2 CN•) and cyanomethyl (H 2 CCN•) radicals via molecular beam studies and electronic structure calcula-tions. The modeling of the chemistries of Taurus Molecular Cloud (TMC-1) and Titan’s atmosphere contemplates a complex chain of reactions synthesizing pyridine, pyridinyl, and (iso)quinoline from H 2 CN• and H 2 CCN• at levels of up to 75%. Furthermore, this study affords unique entry points to precursors of DNA and RNA nitrogen bases in hydrocarbon-rich extraterrestrial environments thus changing the way we think about the origin of prebiotic molecules in our Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗

NASA NTRS1988

Turbulence modeling in hypersonic inlets

A study is conducted to analyze the performance of different turbulence models when applied to flow through a Mach 7.4 hypersonic inlet. The analysis, which is two-dimensional, is done by comparing computational results from a Parabolized Navier-Stokes code and a full Navier-Stokes code, with experimental data. The McDonald-Camarata (MC) and Baldwin-Lomax (BL) models were the two zero-equation models used in the study. The Turbulent Kinetic Energy (TKE) model was chosen as a representative higher order model. The MC model, when run with transition of flow, provides a solution which compares excellently with the data. Transition has a first order effect on the overall solution provided by the code. The BL model predicts separation of flow in the inlet, which contradicts experimental findings. The TKE model does not perform any better than the MC and BL models, despite the fact that it is a higher order turbulence model. The BL and TKE models predict transition in the inlet at a location which is much earlier than observed in the experiment. This may be attributed to the empirical constants used to determine the point of transition.

Ng, W. F.↗

NASA NTRS1977

Factors affecting ion kinetic temperature, number density, and containment time in the NASA Lewis bumpy-torus plasma

The degree of toroidal symmetry of the plasma, the number of midplane electrode rings, the configuration of electrode rings, and the location of the diagnostic instruments with respect to the electrode rings used to generate the plasma are discussed. Impurities were deliberately introduced into the plasma, and the effects of the impurity fraction on ion kinetic temperature and electron number density were observed. It is concluded that, if necessary precautions are taken, the plasma communicates extremely well along the magnetic field lines and displays a high degree of symmetry from sector to sector for a wide range of electrode ring configurations and operating conditions. Finally, some characteristic data taken under nonoptimized conditions are presented, which include the highest electron number density and the longest particle containment time (1.9 msec) observed. Also, evidence from a paired comparison test is presented which shows that the electric field acting along the minor radius of the toroidal plasma improves the plasma density and the calculated containment time more than an order of magnitude if the electric field points inward, relative to the values observed when it points (and pushes ions) radially outward.

Roth, J. R.↗

NASA NTRS2013

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species With Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for nearequilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion of the heating and is then compared to the total heating measured in flight.

Liechty, Derek S.↗

DOE OSTI2020

Coupling between interfacial charge and mechanical deformation at high temperatures in ceramics. Final report

Basic research, under this project, has led to new fundamental discoveries in the effect of electric field (and temperature) on diffusion kinetics, electroluminescence and electronic conductivity in ceramic materials. For example, yttria stabilized zirconia sinters in just seconds well below 1000°C, while, at the same time, exhibiting a highly non-linear rise in conductivity, and intense electroluminescence. The phenomenon has been shown to occur in nearly all ceramics whether they are ionic or electronic conductors, or insulators. The rise in conductivity leads to Joule heating but in-situ measurements with a platinum standard at APS and Brookhaven (NSLS-II) have shown that the temperatures are much below what would have been needed to sinter in a few seconds. We have explored the hypothesis that huge concentrations of Frenkel pairs of vacancies and interstitials are produced. While these point defects accelerate sintering (by the interstitials being transported to the pores and the vacancies to the grain boundaries), they may also produce electron-hole pairs which can impart high electrical conductivity or recombine to emit photons. Many measurements, including in-situ measurements of lattice expansion are consistent with this hypothesis. However, the energy barrier for generating Frenkels in zirconia for example is much larger than can be expected at temperature and electrical fields in these experiments. Most recent work suggests that short wavelength phonons play a critical role in the genesis of the “flash” phenomenon. For example, the Debye temperature has been shown to be the lower bound for initiating flash. Molecular dynamics simulations have shown that large concentrations of Frenkel pairs can form under the proliferation of “optical” phonons at the edge of the Brillouin zone, that is with a wavelength at or below the lattice parameter. This project has led to the development of new experimental techniques as well as a software hardware interface for process and control and data analysis in real time, with a time scale of less than one second. This approach spells a new paradigm in materials science research. It potential impact on high rate manufacturing is self evident.

36 MATERIALS SCIENCE↗

NASA NTRS2014

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species with Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for near-equilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion.

Liechty, Derek S.↗