Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic modeling”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗
NASA NTRS1995

Three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawnside and duskside of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H(+) fluxes at tens of keV, which are always overestimated. A newly invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, J(sub o)(1 + Ay(sup n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗
NASA NTRS1999

A Gas-Kinetic Scheme for Multimaterial Flows and Its Application in Chemical Reaction

This paper concerns the extension of the multicomponent gas-kinetic BGK-type scheme to multidimensional chemical reactive flow calculations. In the kinetic model, each component satisfies its individual gas-kinetic BGK equation and the equilibrium states of both components are coupled in space and time due to the momentum and energy exchange in the course of particle collisions. At the same time, according to the chemical reaction rule one component can be changed into another component with the release of energy, where the reactant and product could have different gamma. Many numerical test cases are included in this paper, which show the robustness and accuracy of kinetic approach in the description of multicomponent reactive flows.

Lian, Yongsheng↗
NASA NTRS1987

Grain dynamics in zero gravity

The dynamics of granular materials has proved difficult to model, primarily because of the complications arising from inelastic losses, friction, packing, and the effect of many grains being in contact simultaneously. One interesting limit for which it was recently possible to construct a theory is that where the grain-grain interactions are dominated by binary collisions. The kinetic model of granular systems if similar to the kinetic theory of gases, except that collisional energy losses are always present in the former and must be treated explicitly. Few granular materials on Earth are describable by this limiting model, since gravity tends to collapse the grains into a high-density state where Coulombic friction effects are dominant. The planned Space Station offers an unusual opportunity to test the kinetic grain model and to explore its predictions. Without gravity, the regime of low interparticle velocities, where an elastic description of the collision is still valid, is investigated. This will allow direct interpretation by dynamical computer simulations as well as by kinetic theory.

Werner, B. T.↗
NASA NTRS1990

Fracture-tough, corrosion-resistant bearing steels

The fundamental principles allowing design of stainless bearing steels with enhanced toughness and stress corrosion resistance has involved both investigation of basic phenomena in model alloys and evaluation of a prototype bearing steel based on a conceptual design exercise. Progress in model studies has included a scanning Auger microprobe (SAM) study of the kinetics of interfacial segregation of embrittling impurities which compete with the kinetics of alloy carbide precipitation in secondary hardening steels. These results can define minimum allowable carbide precipitation rates and/or maximum allowable free impurity contents in these ultrahigh strength steels. Characterization of the prototype bearing steel designed to combine precipitated austenite transformation toughening with secondary hardening shows good agreement between predicted and observed solution treatment response including the nature of the high temperature carbides. An approximate equilibrium constraint applied in the preliminary design calculations to maintain a high martensitic temperature proved inadequate, and the solution treated alloy remained fully austenitic down to liquid nitrogen temperature rather than transforming above 200 C. The alloy can be martensitically transformed by cryogenic deformation, and material so processed will be studied further to test predicted carbide and austenite precipitation behavior. A mechanistically-based martensitic kinetic model was developed and parameters are being evaluated from available kinetic data to allow precise control of martensitic temperatures of high alloy steels in future designs. Preliminary calculations incorporating the prototype stability results suggest that the transformation-toughened secondary-hardening martensitic-stainless design concept is still viable, but may require lowering Cr content to 9 wt. pct. and adding 0.5 to 1.0 wt. pct. Al. An alternative design approach based on strain-induced martensitic transformation during cryogenic forming, thus removing the high martensitic constraint, may permit alloy compositions offering higher fracture roughness.

Olson, Gregory B.↗
NASA NTRS1994

Direct numerical simulations of turbulent non-premixed methane-air flames modeled with reduced kinetics

Turbulent non-premixed stoichiometric methane-air flames modeled with reduced kinetics have been studied using the direct numerical simulation approach. The simulations include realistic chemical kinetics, and the molecular transport is modeled with constant Lewis numbers for individual species. The effect of turbulence on the internal flame structure and extinction characteristics of methane-air flames is evaluated. Consistent with earlier DNS with simple one-step chemistry, the flame is wrinkled and in some regions extinguished by the turbulence, while the turbulence is weakened in the vicinity of the flame due to a combination of dilatation and an increase in kinematic viscosity. Unlike previous results, reignition is observed in the present simulations. Lewis number effects are important in determining the local stoichiometry of the flame. The results presented in this work are preliminary but demonstrate the feasibility of incorporating reduced kinetics for the oxidation of methane with direct numerical simulations of homogeneous turbulence to evaluate the limitations of various levels of reduction in the kinetics and to address the formation of thermal and prompt NO(x).

Card, J. M.↗
NASA NTRS1977

Bhatnagar-Gross-Krook model solution of back-scattering of outgas flow from spherical spacecraft

The BGK kinetic model for a binary mixture of gases was used to predict the backscattering return flux of molecules outgassing from a spherical spacecraft. For simplicity, the outgassing molecules were considered to be a single species, with the ambient atmospheric molecules being the second species in the binary mixture. Two basic collision mechanisms were considered: 1) collisions involving only outgas molecules (self-scattering); and 2) collisions between outgas molecules and molecules in the ambient atmosphere (ambient-scattering). Results were obtained in the form of simple equations relating outgas return flux to parameters based on outgas characteristics, spacecraft size, and the characteristics of the ambient flow.

Robertson, S. J.↗
NASA NTRS2002

Fundamentals of Plasma Sails Propulsion Concept

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.

Khazanov, G. V.↗
NASA NTRS2024

Fundamental Experimental Tests and Modeling of LOX/CH4 Engines at High Pressures

In recent years, the space travel industry has grown exponentially, resulting in the need for a low-cost, efficient rocket engine fuel. As such, there has been renewed interest in utilizing liquid natural gas (LNG) as it is less likely to soot than kerosene-based fuels and has widespread availability and low cost of use as compared to the traditional RP-X and even liquid methane (LCH4) fuels due to the reduced need for refinement. However, current literature does not fully cover varying blends of natural gas and the effects that impurities have on natural gas/methane fuels, with most studies being confined to low-pressure applications. Therefore, new experimental ignition delay time measurements at rocket engine relevant pressures are essential to ensure current chemical kinetic models capture the ignition behavior at these elevated conditions and to determine the optimal purity levels for a reliable, cost-efficient, aerospace-grade LCH4/LNG fuel. The current work explores various blends of natural gas/methane fuel with various impurities, including higher hydrocarbons as well as nitrogen-carriers. A shock tube study was carried out to study the ignition delay times and carbon monoxide time histories of these natural gas blends utilizing chemiluminescence and laser absorption spectroscopy at 20, 50 and 100 atm over temperatures ranging from 1400-1700K and at an equivalence ratio of 1. The experimental data was then compared to two different chemical kinetic mechanisms: the industry-standard GRI 3.0 and the in-house UCF 2022. The data is used to improve chemical kinetic mechanisms and modeling of LNG fuels for rocket engine applications.

Methane↗
NASA NTRS2002

Origins of Solar Systems Program

The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments.

O'D. Alexander, Conel↗
NASA NTRS2003

Final Summary of Research Report to the National Aeronautics and Space Administration Origins of Solar Systems Program

The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments. Cosmic spherules present an opportunity to test our evaporation models under flash heating conditions that would be difficult to simulate experimentally. However, there is surprisingly little known about the isotopic compositions of silicate cosmic spherules, and a number of questions need to be addressed. Is the range of compositions they exhibit due to evaporation? If they are, are the relative volatilities consistent with the models/experiments and are the isotopic fractionations consistent with Rayleigh conditions? For instance, do the alkalis and S evaporate prior to significant melting so that conditions did not meet the Rayleigh criteria of rapid diffusion? If so, their isotopic fractionation might be considerably suppressed. Could this mechanism of K loss apply to chondrule formation? The Fe isotopic fractionation during evaporation of silicates has not been measured, so cosmic spherules might provide a clue to whether FeO diffusion is fast enough to maintain Rayleigh conditions during evaporation. And so on.

O'D. Alexander, Conel↗
NASA NTRS2013

NASA Models of Space Radiation Induced Cancer, Circulatory Disease, and Central Nervous System Effects

The risks of late effects from galactic cosmic rays (GCR) and solar particle events (SPE) are potentially a limitation to long-term space travel. The late effects of highest concern have significant lethality including cancer, effects to the central nervous system (CNS), and circulatory diseases (CD). For cancer and CD the use of age and gender specific models with uncertainty assessments based on human epidemiology data for low LET radiation combined with relative biological effectiveness factors (RBEs) and dose- and dose-rate reduction effectiveness factors (DDREF) to extrapolate these results to space radiation exposures is considered the current "state-of-the-art". The revised NASA Space Risk Model (NSRM-2014) is based on recent radio-epidemiology data for cancer and CD, however a key feature of the NSRM-2014 is the formulation of particle fluence and track structure based radiation quality factors for solid cancer and leukemia risk estimates, which are distinct from the ICRP quality factors, and shown to lead to smaller uncertainties in risk estimates. Many persons exposed to radiation on earth as well as astronauts are life-time never-smokers, which is estimated to significantly modify radiation cancer and CD risk estimates. A key feature of the NASA radiation protection model is the classification of radiation workers by smoking history in setting dose limits. Possible qualitative differences between GCR and low LET radiation increase uncertainties and are not included in previous risk estimates. Two important qualitative differences are emerging from research studies. The first is the increased lethality of tumors observed in animal models compared to low LET radiation or background tumors. The second are Non- Targeted Effects (NTE), which include bystander effects and genomic instability, which has been observed in cell and animal models of cancer risks. NTE's could lead to significant changes in RBE and DDREF estimates for GCR particles, and the potential effectiveness of radiation mitigator's. The NSRM- 2014 approaches to model radiation quality dependent lethality and NTE's will be described. CNS effects include both early changes that may occur during long space missions and late effects such as Alzheimer's disease (AD). AD effects 50% of the population above age 80-yr, is a degenerative disease that worsens with time after initial onset leading to death, and has no known cure. AD is difficult to detect at early stages and the small number of low LET epidemiology studies undertaken have not identified an association with low dose radiation. However experimental studies in mice suggest GCR may lead to early onset AD. We discuss modeling approaches to consider mechanisms whereby radiation would lead to earlier onset of occurrence of AD. Biomarkers of AD include amyloid beta (A(Beta)) plaques, and neurofibrillary tangles (NFT) made up of aggregates of the hyperphosphorylated form of the micro-tubule associated, tau protein. Related markers include synaptic degeneration, dentritic spine loss, and neuronal cell loss through apoptosis. Radiation may affect these processes by causing oxidative stress, aberrant signaling following DNA damage, and chronic neuroinflammation. Cell types to be considered in multi-scale models are neurons, astrocytes, and microglia. We developed biochemical and cell kinetics models of DNA damage signaling related to glycogen synthase kinase-3(Beta) (GSK3(Beta)) and neuroinflammation, and considered multi-scale modeling approaches to develop computer simulations of cell interactions and their relationships to A(Beta) plaques and NFTs. Comparison of model results to experimental data for the age specific development of A(Beta) plaques in transgenic mice will be discussed.

Cucinotta, Francis A.↗
NASA NTRS1995

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗
NASA NTRS1984

Theoretcial studies of solar-pumped lasers

A method of pumping a COhZ laser by a hot cavity was demonstrated. The cavity, heated by solar radiation, should increase the efficiency of solar pumped lasers used for energy conversion. Kinetic modeling is used to examine the behavior of such a COhZ laser. The kinetic equations are solved numerically vs. time and, in addition, steady state solutions are obtained analytically. The effect of gas heating filling the lower laser level is included. The output power and laser efficiency are obtained as functions of black body temperature and gas ratios (COhZ-He-Ar) and pressures. The values are compared with experimental results.

Harries, W. L.↗
NASA NTRS1989

Conservation equations and physical models for hypersonic air flows in thermal and chemical nonequilibrium

The conservation equations for simulating hypersonic flows in thermal and chemical nonequilibrium and details of the associated physical models are presented. These details include the curve fits used for defining thermodynamic properties of the 11 species air model, curve fits for collision cross sections, expressions for transport properties, the chemical kinetics models, and the vibrational and electronic energy relaxation models. The expressions are formulated in the context of either a two or three temperature model. Greater emphasis is placed on the two temperature model in which it is assumed that the translational and rotational energy models are in equilibrium at the translational temperature, T, and the vibrational, electronic, and electron translational energy modes are in equilibrium at the vibrational temperature, T sub v. The eigenvalues and eigenvectors associated with the Jacobian of the flux vector are also presented in order to accommodate the upwind based numerical solutions of the complete equation set.

Peter A. Gnoffo↗
NASA NTRS2000

Modeling of HiPco Process for Carbon Nanotube Production

High-pressure carbon monoxide (HiPco) reactor, developed at Rice University, is used to produce single-walled carbon nanotubes (SWNT) from gas-phase reactions of iron carbonyl and nickel carbonyl in carbon monoxide at high pressures (10 - 100 atm). Computational modeling is used to better understand the HiPco process. In the present model, decomposition of the precursor, metal cluster formation and growth, and carbon nanotube growth are addressed. Decomposition of precursor molecules is necessary to initiate metal cluster formation. The metal clusters serve as catalysts for carbon nanotube growth. Diameter of metal clusters and number of atoms in these clusters are some of the essential information for predicting carbon nanotube formation and growth, which is then modeled by Boudouard reaction (2CO ---> C(s) + CO2) with metal catalysts. The growth kinetic model is integrated with a two-dimensional axisymmetric reactor flow model to predict reactor performance.

Gokcen, T.↗
NASA NTRS2002

Paralinear Oxidation of Silicon Nitride in a Water Vapor/Oxygen Environment

Three silicon nitride materials were exposed to dry oxygen flowing at 0.44 cm/s at temperatures between 1200 and 1400 C. Reaction kinetics were measured with a continuously recording microbalance. Parabolic kinetics were observed. When the same materials were exposed to a 50% H2O - 50% O2 gas mixture flowing at 4.4 cm/s, all three types exhibited paralinear kinetics. The material is oxidized by water vapor to form solid silica. The protective silica is in turn volatilized by water vapor to form primarily gaseous Si(OH)4. Nonlinear least squares analysis and a paralinear kinetic model were used to determine both parabolic and linear rate constants from the kinetic data. Volatilization of the protective silica scale can result in accelerated consumption of Si3N4. Recession rates under conditions more representative of actual combustors are compared to the furnace data.

Fox, Dennis S.↗