Search NASA⌕ Search

SEARCH · Search NASA

Results for “Collision rate coefficient”

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 163 records · Page 9

Formation of the 0.511 MeV line in solar flares

The slowing down and annihilation of positrons and the formation of positronium in a solar flare plasma are investigated to determine how the width of the 0.511-MeV line and its strength relative to the three-photon continuum from positronium decay depend on the temperature and density of the medium in which the positron comes to rest. The calculations are limited to the cases of annihilation in a completely ionized plasma, in a partially ionized plasma with an electron/neutron density ratio of 1.0 or 0.1, and in an atomic gas with a very small ion density. Thermally averaged rate coefficients are obtained for the free annihilation of positrons and for positronium formation through radiative recombination in the fully ionized plasma. Positronium formation rates and the resultant energy distributions of the positronium atoms in the partially ionized medium are determined by numerically solving the Fokker-Planck equation in a medium where the ambient free electrons have a Maxwell-Boltzmann distribution of finite temperature but the density of the medium is sufficiently low that positronium atoms decay without further collisions following their formation. A Monte Carlo calculation is performed for the positron energy loss, positronium formation through charge exchange, and positronium breakup in the weakly ionized medium. The energy distributions of decaying positronium atoms and the relative number of triplet to singlet positronium decays are evaluated for ion concentrations and densities characteristic of the solar photosphere.

Crannell, C. J.↗

Microscopic Simulation and Macroscopic Modeling for Thermal and Chemical Non-Equilibrium

This paper deals with the accurate microscopic simulation and macroscopic modeling of extreme non-equilibrium phenomena, such as encountered during hypersonic entry into a planetary atmosphere. The state-to-state microscopic equations involving internal excitation, de-excitation, dissociation, and recombination of nitrogen molecules due to collisions with nitrogen atoms are solved time-accurately. Strategies to increase the numerical efficiency are discussed. The problem is then modeled using a few macroscopic variables. The model is based on reconstructions of the state distribution function using the maximum entropy principle. The internal energy space is subdivided into multiple groups in order to better describe the non-equilibrium gases. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients. The modeling is completely physics-based, and its accuracy depends only on the assumed expression of the state distribution function and the number of groups used. The model makes no assumption at the microscopic level, and all possible collisional and radiative processes are allowed. The model is applicable to both atoms and molecules and their ions. Several limiting cases are presented to show that the model recovers the classical twotemperature models if all states are in one group and the model reduces to the microscopic equations if each group contains only one state. Numerical examples and model validations are carried out for both the uniform and linear distributions. Results show that the original over nine thousand microscopic equations can be reduced to 2 macroscopic equations using 1 to 5 groups with excellent agreement. The computer time is decreased from 18 hours to less than 1 second.

Modeling↗

The Effect of the Spin-Forbidden Co((sup 1) Sigma plus) plus O((sup 3) P) Yields CO2 (1 Sigma (sub G) plus) Recombination Reaction on Afterbody Heating of Mars Entry Vehicles

Vibrationally excited CO2, formed by two-body recombination from CO((sup 1) sigma plus) and O((sup 3) P) in the wake behind spacecraft entering the Martian atmosphere reaction, is potentially responsible for the higher than anticipated radiative heating of the backshell, compared to pre-flight predictions. This process involves a spin-forbidden transition of the transient triplet CO2 molecule to the longer-lived singlet. To accurately predict the singlet-triplet transition probability and estimate the thermal rate coefficient of the recombination reaction, ab initio methods were used to compute the first singlet and three lowest triplet CO2 potential energy surfaces and the spin-orbit coupling matrix elements between these states. Analytical fits to these four potential energy surfaces were generated for surface hopping trajectory calculations, using Tully's fewest switches surface hopping algorithm. Preliminary results for the trajectory calculations are presented. The calculated probability of a CO((sup 1) sigma plus) and O((sup 3) P) collision leading to singlet CO2 formation is on the order of 10 (sup -4). The predicted flowfield conditions for various Mars entry scenarios predict temperatures in the range of 1000 degrees Kelvin - 4000 degrees Kelvin and pressures in the range of 300-2500 pascals at the shoulder and in the wake, which is consistent with a heavy-particle collision frequency of 10 (sup 6) to 10 (sup 7) per second. Owing to this low collision frequency, it is likely that CO((sup 1) sigma plus) molecules formed by this mechanism will mostly be frozen in a highly nonequilibrium rovibrational energy state until they relax by photoemission.

CO2 recombination↗

Electron-nitrogen molecular collisions in high-temperature nonequilibrium air

Ab initio calculation of vibrational excitation and deexcitation cross sections of N2 by low energy resonant electron impact have been carried out. The calculation includes initial target state in v = 0-12 and final target state in v = 0-17. The calculated cross sections are found to agree with available experimental data to within reported experimental error. Effect of target rotation has been investigated. Comparison of vibrational excitation cross sections at J = 0, 50, and 150 shows that the major effect of target rotation is to lower the resonance energy. Vibrational excitation and deexcitation rate coefficients are tabulated for electron temperatures of 1100-55,000 K. The present tabulation supercedes an earlier compilation by the Naval Research Laboratory that was derived from experimental data with a factor-of-two error. The tabulated rate coefficients are suited for flowfield modeling of aeroassisted orbital transfer vehicles.

Huo, W. M.↗

Collision cross sections and diffusion parameters for H and D in atomic oxygen

Modeling the behavior of H and D in planetary exospheres requires detailed knowledge of the differential scattering cross sections for all of the important neutral-neutral and ion-neutral collision processes affecting these species over their entire ranges of interaction energies. In the upper atmospheres of Earth, Venus, and other planets as well, the interactions of H and D with atomic oxygen determine the rates of diffusion of escaping hydrogen isotopes through the thermosphere, the velocity distributions of exospheric atoms that encounter the upper thermosphere, the lifetimes of exospheric orbiters with periapsides near the exobase, and the transfer of momentum in collisions with hot O. The nature of H-O and D-O collisions and the derivation of a data base consisting of phase shifts and the differential, total, and momentum transfer cross sections for these interactions in the energy range 0.001 - 10 eV are discussed. Coefficients of mutual diffusion and thermal diffusion factors are calculated for temperatures of planetary interest.

Hodges, R. R., Jr.↗

Cloud based tool for analysis of chemical kinetic mechanisms

The size and complexity of chemical kinetic models used to simulate combustion of hydrocarbon fuels continue to increase with our improved understanding of the underlying physical mechanisms, along with the improved ability to utilize such models as the result of increased computational power and efficiency. As mechanisms grow beyond thousands of species and tens of thousands of reactions, it becomes increasingly difficult to manually check for errors and physical inconsistencies. In this work, we present several automated methods to check for such issues in the specification of chemical reaction models. First, we demonstrate how discontinuities in thermodynamic data can cause simulation difficulties manifested as long wall-clock times and convergence failures. To correct this type of problem, we describe an automated method for refitting thermodynamic parameters. We also outline several methods to check the timescales of reaction rate coefficients to ensure physical consistency. All the methods are made available through a web-based tool (https://combustiontools.llnl.gov) intended to aid mechanism developers and users to improve the accuracy and performance of fuel models used by the combustion community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantification of Uncertainty and Risk Sensitivity for Safety of Emerging Operations

The growing need to develop and deploy small unmanned aerial vehicles (sUAVs) for various applications in the airspace necessitates reliable tools to accurately predict the flight trajectories of the sUAVs. The knowledge of the predicted trajectories help decision makers anticipate potential conflict, assess the risk, and take appropriate risk mitigation actions. In addition, uncertainties in vehicle models, weather, and controller action further highlights the need for reliable prediction tools. In this project, the application of mixed sparse grid-based quadrature and generalized polynomial chaos(gPC) expansion method for uncertainty quantification and collision assessment in air traffic consisting of fixed-wing small unmanned aerial vehicles (sUAV) was studied. From the results obtained, it can be concluded that this provides a reliable framework to carry out quantitative conflict assessment in an unmanned air traffic, which when employed, can improve the functionalities of the unmanned traffic management system. It was observed that the results from the gPC expansion framework developed in the project can be utilized to conduct rapid probabilistic collision assessment for near real-time unmanned traffic management in the airspace. From the vehicle models, position updates, and wind-field data, a priori gPC based 3-σcon-fidence ellipses can provide estimates of potential conflict at some future instants. The computational costs scaled linearly when the uncertain inputs were fewer. Further, the largest allowable distribution of para-metric uncertainties that leads to the smallest risk of collision in traffic of small unmanned aerial vehicles could be calculated. The time of closest approach between two sUAVs can be established paving way for development of proactive mitigation strategies. The separation between the sUAVs was found to be most significantly affected by uncertainties in the maximum available thrusts, zero-lift drag coefficients, and wing planform areas of the sUAVs. The study of uncertain wind-fields indicated that a heterogeneous traffic mix resulted in an increased probability of conflict. Increased measurement update rate reduced the uncertain-ties in the trajectories of the vehicles, further reducing the probability of conflict but rapid updates of all vehicles in the airspace poses a stringent communication limitation. The gPC framework also provided the means to analyze vehicle impact (crash region) due to loss of control resulting from actuator failure in sUAS traffic, essentially to predict impact and crash zones for representative vehicles. The predicted regions when compared with non-participant density, provides a means to develop an early mitigation strategy, should the sUAV detect an imminent actuator failure.

Rajnish Bhusal↗

Evolution of inclined galactic gas disks. I - A cloud-fluid approach

In this paper, numerical models are developed to follow the time-dependent behavior of inclined gas disks in nonspherical galaxies. Differential equations are developed which describe the time rate of change of the orbits of mass elements composing a thin, smoothly warping and twisting gas disk. These equations include multiple moments due to the nonsphericity of the galactic mass distribution, Coriolis effects arising from determination of orbits in a reference frame which tumbles with the surface figure of the galaxy, and viscous forces arising from cloud-cloud collisions. A cloud-fluid approach is used in determining the viscous forces. Physical arguments are used to determine the maximum permissible local value of the coefficient of kinematic viscosity. For small initial inclinations, analytic expressions are obtained for the inclination as a function of time and for the settling time, and the features of the solution are discussed.

Steiman-Cameron, Thomas Y.↗

Analysis of the upper atmosphere CO2(nu-2) vibrational temperatures retrieved from ATMOS/Spacelab 3 observations

An analysis of the upper atmospheric (80-116 km) CO2(nu-2) vibrational temperatures retrieved from atmospheric trace molecule spectroscopy (ATMOS) experiment Spacelab 3 spectra by using a nonlocal thermodynamic equilibrium (non-LTE) radiative transfer model is presented. Thermal collisions with atmospheric atomic oxygen keep this vibrational state very close to LTE up to around 100 km. Above this height, the different deviations from LTE retrieved from ATMOS/Spacelab 3 spectra for the Northern and the Southern Hemispheres are explained in terms of this collisional process and in terms of the different kinetic temperature profiles measured at those locations. From these simultaneous observations of the kinetic and CO2(nu-2) vibrational temperatures, a deactivation rate of CO2(nu-2) by O(3P) has been derived which leads to a rate coefficient value between 3 and 6 x 10 exp -12 cu cm/s and favors an independent or negative temperature dependence rate constant for the atmospheric temperature range. Cooling rates induced by the CO2 15-micron fundamental band in the upper mesosphere and lower thermosphere were derived from the simultaneous kinetic temperature, CO2(nu-2) vibrational temperature, and CO2 concentration, as measured by ATMOS/Spacelab 3, and found to be a factor of between 5 and 10 times larger than those generally accepted until very recently.

Lopez-Puertas, M.↗

Vibrational nonequilibrium effects on diatomic dissociation rates

The collision-induced dissociation rate of diatomic molecules from a ladder of rotational and anharmonic vibrational states is developed, and the correction for vibrational nonequilibrium is considered. The result is similar to an analytic correction derived by Hammerling et al. (1959) for harmonic oscillators. An empirical correction algorithm suggested by Park (1987, 1990) gives similar results when vibrational temperature is comparable to kinetic temperature but underestimates the dissociation rate when vibrational temperature is small compared with the kinetic temperature. This algorithm uses an effective temperature in the experimentally determined Arrhenius expression for the rate coefficient, which is a weighted average of the vibrational and kinetic temperature, whereas theory indicates that kinetic temperature should appear only in the exponential term of the Arrhenius expression. Nevertheless, an effective temperature can always be found that will numerically duplicate the proper rate coefficient at any given condition, but a constant weighting factor cannot be expected to provide this. However, the algorithm can he adjusted to give reasonable results over a range of conditions if the geometric weighting factor is taken to be a simple linear function of the ratio of vibrational to kinetic temperature in the gas.

Hansen, C. F.↗

Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions

A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.

plasma↗

Event generator for jet tomography in electron-ion collisions

We develop the first event generator, the electron-heavy-ion-jet-interaction-generator (eHIJING), for the jet tomography study of electron-ion collisions. In this generator, energetic jet partons produced from the initial hard scattering undergo multiple collisions with the nuclear target. The collision rate is proportional to the transverse-momentum-dependent (TMD) gluon density in the nucleus, which is given by a simple model inspired by the physics of gluon saturation. Medium-modified QCD splitting functions within the higher-twist (HT) and generalized higher-twist (GHT) frameworks are utilized to simulate parton showering in the nuclear medium that takes into account the non-Abelian Landau-Pomeranchuck-Midgal interference effect. Employing eHIJING, we revisit hadron production in semi-inclusive deep inelastic scattering (SIDIS) as measured by EMC, HERMES, and recent CLAS experiments. eHIJING with both GT and GHT frameworks gives reasonably good descriptions of these experimental data. Predictions for experiments at the future electron-ion colliders are also provided. It is demonstrated that future measurements of the transverse momentum broadening of single hadron spectra can be used to map out the two-dimensional kinematic (Q 2 , x B ) dependence of the jet transport coefficient $\hat{q}$ in cold nuclear matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Estimated Probability of a Cervical Spine Injury During an ISS Mission

Introduction: The Integrated Medical Model (IMM) utilizes historical data, cohort data, and external simulations as input factors to provide estimates of crew health, resource utilization and mission outcomes. The Cervical Spine Injury Module (CSIM) is an external simulation designed to provide the IMM with parameter estimates for 1) a probability distribution function (PDF) of the incidence rate, 2) the mean incidence rate, and 3) the standard deviation associated with the mean resulting from injury/trauma of the neck. Methods: An injury mechanism based on an idealized low-velocity blunt impact to the superior posterior thorax of an ISS crewmember was used as the simulated mission environment. As a result of this impact, the cervical spine is inertially loaded from the mass of the head producing an extension-flexion motion deforming the soft tissues of the neck. A multibody biomechanical model was developed to estimate the kinematic and dynamic response of the head-neck system from a prescribed acceleration profile. Logistic regression was performed on a dataset containing AIS1 soft tissue neck injuries from rear-end automobile collisions with published Neck Injury Criterion values producing an injury transfer function (ITF). An injury event scenario (IES) was constructed such that crew 1 is moving through a primary or standard translation path transferring large volume equipment impacting stationary crew 2. The incidence rate for this IES was estimated from in-flight data and used to calculate the probability of occurrence. The uncertainty in the model input factors were estimated from representative datasets and expressed in terms of probability distributions. A Monte Carlo Method utilizing simple random sampling was employed to propagate both aleatory and epistemic uncertain factors. Scatterplots and partial correlation coefficients (PCC) were generated to determine input factor sensitivity. CSIM was developed in the SimMechanics/Simulink environment with a Monte Carlo wrapper (MATLAB) used to integrate the components of the module. Results: The probability of generating an AIS1 soft tissue neck injury from the extension/flexion motion induced by a low-velocity blunt impact to the superior posterior thorax was fitted with a lognormal PDF with mean 0.26409, standard deviation 0.11353, standard error of mean 0.00114, and 95% confidence interval [0.26186, 0.26631]. Combining the probability of an AIS1 injury with the probability of IES occurrence was fitted with a Johnson SI PDF with mean 0.02772, standard deviation 0.02012, standard error of mean 0.00020, and 95% confidence interval [0.02733, 0.02812]. The input factor sensitivity analysis in descending order was IES incidence rate, ITF regression coefficient 1, impactor initial velocity, ITF regression coefficient 2, and all others (equipment mass, crew 1 body mass, crew 2 body mass) insignificant. Verification and Validation (V&V): The IMM V&V, based upon NASA STD 7009, was implemented which included an assessment of the data sets used to build CSIM. The documentation maintained includes source code comments and a technical report. The software code and documentation is under Subversion configuration management. Kinematic validation was performed by comparing the biomechanical model output to established corridors.

Brooker, John E.↗

Charge Transfer Calculations and Database for Astrophysics

A variety of theoretical approaches, having different but overlapping energy ranges of applicability, are applied to investigate charge transfer processes for collisions of atomic ions with atoms and molecules. The methods include quantal molecular-orbital close-coupling, classical trajectory Monte Carlo, and continuum distorted wave methods. Recent collision systems studied include S(+4) + H, S(+4) + He, N(+7) + He, H2O, CO, and CO2, O(+q)(q = 1 - 8) + H, H2, and S(+q)(q = 1 - 16) + H2. The database effort is concentrating on astrophysically important reactions of atomic ions X(+q)(X=H-Zn, q=1-4, and selected higher charges) with H, He, various metal atoms, H2, and other selected molecular targets. Existing data, much of it produced by us, has been compiled and critically evaluated. Data for many reactions missing in the literature are estimated using the multichannel Landau-Zener approximation. Fits to cross sections and rate coefficients using standard functions are provided as well as tabulations of the raw data. The database is available on the World Wide Web at cfadc.phy.ornl.gov/astro/ps/data.

Wang, J. G.↗

Nonadiabatically Driven Quantum Interference Effects in the Ultracold K + KRb → Rb + K 2 Chemical Reaction

The K + KRb → Rb + K 2 chemical reaction is the first ultracold atom–diatom chemical reaction for which experimental results have been reported for temperatures below 1 μK more than a decade ago. The reaction occurs through coupling with an excited electronic state that is accessible even in the ultracold limit. A previous quantum dynamics study, excluding nonadiabatic effects, has reported a rate coefficient that is about 35% below the experimental value. Here, we report the first nonadiabatic quantum dynamics study of this reaction and obtain rate coefficients in better agreement with experiments. Our results show that short-range dynamics mediated by coupling with the excited electronic state introduces quantum interference effects that influence both the state-to-state rate coefficients and the overall reaction rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solar observations and atomic data for the 3s2 1s0-3s3p 3P1 transition in S v

High resolution solar observations of the S v intersystem line at 1199.18 A are available from Skylab. This line is potentially useful as a density diagnostic for high density plasmas expected in solar flares. S v lines are also prominent in solar spectra at wavelengths below 1000 A. Collision strengths and radiative decay rates are calculated for levels of the configurations 3 s2, 3s3p, 3p2, and 3s3d. Level populations for the five lowest energy levels have been calculated as a function of electron density. These calculations are carried out assuming the temperature at which S v is most abundant in solar plasmas. The calculated population of the 3s3p 3P1 level reaches a pseudo-Boltzmann equilibrium at a density which is about five times higher than is derived from solar spectra and previous density determinations. It is believed that the discrepancy is due to a significant contribution from resonances to the 3s2-3s3p 3P excitation rate coefficient. It is suggested that the contribution of the resonances to the S v excitation rate coefficient can be estimated from a comparison of the experimental results and the theoretical calculations.

Feldman, U.↗

Search for same-charge top-quark pair production in pp collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for the production of top-quark pairs with the same electric charge ($tt$ or $\overline{tt}$) is presented. The analysis uses proton-proton collision data at $\sqrt{s}$ = 13 TeV, recorded by the ATLAS detector at the Large Hadron Collider, corresponding to an integrated luminosity of 140 fb −1 . Events with two same-charge leptons and at least two b-tagged jets are selected. Neural networks are employed to define two selections sensitive to additional couplings beyond the Standard Model that would enhance the production rate of same-sign top-quark pairs. No significant signal is observed, leading to an upper limit on the total production cross-section of same-sign top-quark pairs of 1.6 fb at 95% confidence level. Corresponding limits on the three Wilson coefficients associated with the $\mathcal{O}^{(1)}_{tu}$, $\mathcal{O}^{(1)}_{Qu}$, and $\mathcal{O}^{(8)}_{Qu}$ operators in the Standard Model Effective Field Theory framework are derived.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Minijet quenching in a concurrent jet + hydro evolution and the nonequilibrium quark-gluon plasma

Minijets, created by perturbative hard QCD collisions at moderate energies, can represent a significant portion of the total multiplicity of a heavy-ion collision event. Since their transverse momenta are initially larger than the typical saturation scale describing the bulk of the equilibrating quark-gluon plasma (QGP), they ought to be described through the physics of parton energy loss. Indeed, their typical stopping distances are larger than the usual hydrodynamization time, so they do not in general hydrodynamize at the same pace than the bulk of the collision. Therefore, in general minijets cannot be described solely by a unique preequilibrium stage that bridges the initial, overoccupied glasma state, with the hydrodynamical evolution. In this work we make use of a new concurrent minijet+hydrodynamic framework in which the properties of the hydrodynamically evolving QGP are modified due to the injection of energy and momentum from the minijets. We study the system for different choices of the minimum transverse momentum associated with minijet production. In order to achieve a realistic description of charged particle multiplicity, the amount of entropy associated to the low- x initial state needs to be reduced. Moreover, the fact that the injected momentum from the randomly oriented minijets is not correlated with the spatial gradients of the system reduces overall flow, and the value of the QGP transport coefficients needs to be reduced accordingly in order to describe the measured flow coefficients in experiments. They are, in effect, an important new source of fluctuations, resulting in a spikier, notably modified hydrodynamical evolution when compared to the scenario in which the presence of minijets is ignored. We avow that their abundance makes it necessary to include their physics in holistic descriptions of heavy-ion collisions. We discuss the impact of the minijets on a number of observables, such as p T spectra and p T -differential flow v n for a wide range of centrality classes. In contrast to elliptic, triangular or quadrangular flow, here we find that directed flow, v 1 , has the strongest potential to discriminate between different minijet production rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗