Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic simulations”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Prebiotic Processing induced by Comet and Meteor Impact

In their study of organic synthesis from impact shocks using the laser-induced-plasma (LIP) technique, McKay and Borucki(l) found that organic synthesis preferentially occurred in a reducing gas mixture rich in methane, and not in a mixture rich in carbon dioxide. This result means chemical models based on the thermodynamical equilibrium approach do not apply to shock chemistry. In this study, we employ the technique of reacting flow, i.e., chemical kinetics in a fluid flow, to simulate the chemistry occurring in LIP and in the wake region from comet or meteor impact. Three different air compositions have been used: (1) 1/3 CO2 and 2/3 H2, (2) pure CH4, and (3) 1/4 CH4, 1/4 CO2, and 1/2 H2O. The stoichiometric ratio of gas mixtures (1) and (3) are kept the same. For (1) we obtain equal mole fractions of CO and H2O as the major products and for (2) C2H2 is the major product. In both cases our results are in agreement with Ref. (1). For (3) we find an interesting case where the nature of chemicals produced to be critically dependent on the flowfield temperature. At the higher temperature part of the wake region, CO and H2O are the dominant products, whereas in the cooler region C2H2 is the dominant product. Further studies of these reactions, as well as for the gas mixture including N2, are being pursued.

Dateo, Christopher E.↗
NASA NTRS2012

Paramyxovirus Infection Mimics In Vivo Cellular Dynamics in Three-Demensional Human Bronchio-Epithelial Tissue-Like Assemblies

Respiratory syncytial virus and parainfluenza virus cause severe respiratory disease, especially in infants, children and the elderly. An in vitro model that accurately mimics infection of the human respiratory epithelium (HRE) would facilitate vaccine development greatly. Monolayer cultures traditionally used to study these viruses do not accurately and precisely differentiate the replication efficiencies of wild type and attenuated viruses. Therefore, we engineered novel three-dimensional (3D) tissue-like assemblies (TLAs) of human broncho-epithelial (HBE) cells to produce a more physiologically relevant in vitro model of the HRE. TLAs resemble HRE structurally and by expression of differentiated epithelial cell markers. Most significantly, wild type viruses exhibited a clear growth advantage over attenuated strains in TLAs unlike monolayer cultures. In addition, the TLAs responded to virus infection by secreting pro-inflammatory mediators similar to the respiratory epithelia of infected children. These characteristics make the TLA model a valuable platform technology to develop and evaluate live, attenuated respiratory virus vaccine candidates for human use. Respiratory virus diseases, the most frequent and least preventable of all infectious diseases, range in severity from the common cold to severe bronchiolitis and pneumonia . Two paramyxoviruses, respiratory syncytial virus (RSV) and parainfluenza virus type 3 (PIV3), are responsible for a majority of the most severe respiratory diseases of infants and young children. RSV causes 70% of all bronchiolitis cases and is a major cause of morbidity and mortality worldwide, especially in infants. PIV3 causes 10-15% of bronchiolitis and pneumonia during infancy, second only to RSV, and 40% of croup in infants To date, licensed vaccines are not available to prevent these respiratory diseases. At present, traditional monkey kidney (Vero and LLC-MK2) and human (HEp-2) tissue culture cells and small animal models (mouse, cotton rat, guinea pig, ferret, and hamster) fail to accurately imitate viral replication and human disease states (8). Lacking an authentic model has impeded the development and evaluation of live, attenuated vaccine candidates. Development of a physiologically relevant in vitro tissue culture model that reproduces characteristics of the HRE, the primary target of RSV and PIV3, would aid in predicting clinical attenuation and safety of vaccine candidates. Successful tissue engineering of a 3D human intestinal model using novel NASA technology inspired the development of a tri-culture 3D model for the HRE. Sequential layering of primary mesenchymal cells (comprised of normal human fibroblasts and endothelial cells) followed by BEAS-2B epithelial cells derived from human bronchi and tracheae were recapitulated on Cultisphere and/or cytodex3 microcarriers in cylindrical vessels that rotate horizontally creating an organized epithelial structure. Horizontal rotation randomizes the gravity vector modeling aspects of microgravity. Mesenchymal and epithelial cells grown under these conditions reproduce the structural organization, multi-cellular complexity, and differentiation state of the HRE. The opportunity to study respiratory viruses in a nasal epithelium model is invaluable because the most promising respiratory virus vaccine candidates are live attenuated viruses for intranasal administration. Here we characterize the interactions of respiratory viruses and epithelial cells grown under modeled microgravity in comparison to gravity-ladened monolayers. 3D HBE TLAs and traditional monolayers (2D) are infected at 35 C, the upper temperature of the upper HRE, to simulate in vivo infection conditions. Growth kinetics of wild type (wt) RSV and PIV3 viruses were compared in 2D and 3D cells to that of strains attenuated in humans or rhesus macaques. This novel 3D HBE model also offers an opportunity to study whether the epithelial cell function, especially in host defenses recapitulated by mimicking the structural organization of the HRE. In vivo, airway epithelial cells play a significant and dynamic role in host defense by blocking paracellular permeability and modulating airway function through cellular interactions or tight junctions. As regulators of the innate immune response, epithelial cells constitutively express cytokines, chemokines, and colony stimulating factors including RANTES, IL-8, IL-6, GM-CSF, and G-CSF for proactive host defense. In response to viral infection, epithelial cells induce potent immuno-modulatory and pro-inflammatory cytokines that recruit phagocytic and inflammatory cells to clear the virus and enhance protection. Although disease pathogenesis is classically attributed to the cytopathic effects of the pathogen, severe disease states associated with RSV and PIV3 are attributed to the inflammatory response, especially in infants. RSV is a potent inducer of cytokines and pro-inflammatory mediators in epithelial cells in vivo. A differentiated human epithelial model independent of the complete functional immune system will help elucidate the role of epithelial cells in respiratory disease. We reported here, virus and host cell interactions in 3D HBE TLAs are similar to that in vivo. Because the epithelial cell organization of the TLAs impacts not only the expression of airway epithelial characteristics, but also cellular communication, the TLAs represent a more physiologically relevant model of the HRE than BEAS-2B or other non-tumour monolayer models of respiratory disease. As a result, wild type respiratory viruses have a clear growth advantage over attenuated viruses in TLAs unlike traditional monolayers. In addition, the TLAs respond to wild type virus infection by secreting pro-inflammatory mediators characteristic of infected HRE. TLAs expressing microbial defense mechanisms provide an excellent model to study the interactions of respiratory pathogens with their host and to identify the innate immunity mediators. Therefore, 3D HBE TLAs offer advantages for the study of respiratory viruses and the development of viral vaccine candidates.

Deatly, Anne M.↗
NASA NTRS2018

Collisionless Relaxation of the Ion Ring Distribution in Space Plasma

Energetic processes often produce transversely-heated angular distributions of the magnetized core (lowest energy) plasma. This characteristic is found in solar wind ion pickup, resulting from cometary or interstellar gas ionization, in Earths' ionosphere, and with hot ions formed around the Space Transportation System during gas releases. We investigate the thermalization of O+ ion pickup using the 2.5D hybrid simulation method (with fluid electrons and kinetic ions) of the ion pickup (ring) distributions, formed in the auroral ionosphere, with a range of ring velocities and thermal to magnetic pressure ratios. We find that in the unstable collisonless regime the anisotropy of the non-thermal distribution produces the ion-cyclotron instability, and the nonlinear relaxation is accompanied by wave-particle scattering that results in an emitted power of EMIC waves. We conclude that ionospheric pickup thermalization is slow due to the small ring speed compared to the thermal and Alfven speeds, while in the solar wind and other space plasmas regions with larger ion-ring velocity the collisionless relaxation and thermalization is rapid in terms of O+ ion gyro-period.

Ofman, Leon↗
NASA NTRSDate not supplied

Examination of Mars2020 Shock-Layer Conditions via Infrared Laser Absorption Spectroscopy of CO2 and CO

Amid-infrared laser absorption diagnostic was deployed to study a simulated Mars2020 shock layer in the Electric Arc Shock Tube (EAST) facility at NASA’s Ames Research Center. Rapid RF-diplexing techniques enabled quantitative temperature and number density measurements of CO2 and CO with μs-resolution over an incident shock velocity range of 1.39 – 3.75 km/s. Two interband cascade lasers were utilized at 4.17 and 4.19 μm to resolve eight CO2 asymmetric stretch fundamental band (𝜈3) transitions from two different vibrational levels: 0000 (ground state) and 0110 (first excited bending mode). The probed rotational levels span across J” = 58 to J” = 140. Results are compared to DPLR simulations of the shock layer using kinetic mechanisms of Johnston et al. and Cruden et al. At shock velocities below 3.1 km/s, the agreement between the measurements and the Johnston mechanism is typically within 5% for temperature and within 10% for number density. At shock velocities above 3.1 km/s, the CO2 measurement becomes sensitive to a thin boundary layer and corrections of this effect are presented. On test cases with enough energy to dissociate CO2, a quantum cascade laser scanned the P(2, 20), P(0, 31), and P(3, 14) transitions of the CO fundamental band at 4.98 μm. CO formation rate is measured to be close to the Johnston kinetic mechanism at low velocities, and then trending towards the Cruden kinetic mechanism at high velocities. On a few low velocity test cases, rovibrational relaxation of the Martian atmosphere is probed with μs resolution.

Laser Absorption↗
NASA NTRS1994

Coalescence of two current loops with a kink instability simulated by a three-dimensional electromagnetic particle code

We have studied the dynamics of a coalescence of current loops using three-dimensional electromagnetic (EM) particle simulation code. Our focus is the investigation of such kinetic processes as energy trasnfer, heating particles, and electromagnetic emissions associated with a current loop coalescence which cannot be studied by MHD simulations. First, the two loops undergo a pinching oscillation due to a pressure imbalance between the inside and outside of the current loop. During the pinching oscillation, a kinetic kink instability is excited and electrons in the loops are heated perpendicularly to an ambient magnetic field. Next, the two current loops collide and coalesce, while at the same time a helical structure grows further. Subsequently, the perturbed current, which is due to these helically bunched electrons, can drive a whistler instability. It should be noted in this case that the whistler wave is excited by the kinetic kink instability and not a beam instability. After the coalescence of two helical loops, tilting motions can be observed in the direction of left-hand rotation, and the helical structure will relax resulting in strong plasma heating mostly in the direction perpendicular to the ambient magnetic field. It is also shown that high-frequency electromagnetic waves can be emitted from the region where the two loops coalesce and propagate strongly in the direction of the electron drift velocity. These processes may be important in understanding heating mechansims for coronal loops as well as radio wave emission mechanisms from active regions of solar plasmas.

Nishikawa, K.-I.↗
NASA NTRS1993

Computation of large-scale statistics in decaying isotropic turbulence

We have performed large-eddy simulations of decaying isotropic turbulence to test the prediction of self-similar decay of the energy spectrum and to compute the decay exponents of the kinetic energy. In general, good agreement between the simulation results and the assumption of self-similarity were obtained. However, the statistics of the simulations were insufficient to compute the value of gamma which corrects the decay exponent when the spectrum follows a k(exp 4) wave number behavior near k = 0. To obtain good statistics, it was found necessary to average over a large ensemble of turbulent flows.

Chasnov, Jeffrey R.↗
NASA NTRS1992

Reduction of iron-bearing lunar minerals for the production of oxygen

The kinetics and mechanism of the reduction of simulants of the iron-bearing lunar minerals olivine ((Fe,Mg)2SiO4), pyroxene ((Fe,Mg,Ca)SiO3), and ilmenite (FeTiO3) are investigated, extending previous work with ilmenite. Fayalite is reduced by H2 at 1070 K to 1480 K. A layer of mixed silica glass and iron forms around an unreacted core. Reaction kinetics are influenced by permeation of hydrogen through this layer and a reaction step involving dissociated hydrogen. Reaction mechanisms are independent of Mg content. Augite, hypersthene, and hedenbergite are reduced in H2 at the same temperatures. The products are iron metal and lower iron silicates mixed throughout the mineral. Activation energy rises with calcium content. Ilmenite and fayalite are reduced with carbon deposited on partially reduced minerals via the CO disproportionation reaction. Reduction with carbon is rapid, showing the carbothermal reduction of lunar minerals is possible.

Massieon, Charles↗
NASA NTRS2010

Spacecraft Solar Particle Event (SPE) Shielding: Shielding Effectiveness as a Function of SPE model as Determined with the FLUKA Radiation Transport Code

Analysis of both satellite and surface neutron monitor data demonstrate that the widely utilized Exponential model of solar particle event (SPE) proton kinetic energy spectra can seriously underestimate SPE proton flux, especially at the highest kinetic energies. The more recently developed Band model produces better agreement with neutron monitor data ground level events (GLEs) and is believed to be considerably more accurate at high kinetic energies. Here, we report the results of modeling and simulation studies in which the radiation transport code FLUKA (FLUktuierende KAskade) is used to determine the changes in total ionizing dose (TID) and single-event environments (SEE) behind aluminum, polyethylene, carbon, and titanium shielding masses when the assumed form (i. e., Band or Exponential) of the solar particle event (SPE) kinetic energy spectra is changed. FLUKA simulations have fully three dimensions with an isotropic particle flux incident on a concentric spherical shell shielding mass and detector structure. The effects are reported for both energetic primary protons penetrating the shield mass and secondary particle showers caused by energetic primary protons colliding with shielding mass nuclei. Our results, in agreement with previous studies, show that use of the Exponential form of the event

Koontz, Steve↗
NASA NTRS1995

A Large-Particle Monte Carlo Code for Simulating Non-Linear High-Energy Processes Near Compact Objects

High-energy radiation processes in compact cosmic objects are often expected to have a strongly non-linear behavior. Such behavior is shown, for example, by electron-positron pair cascades and the time evolution of relativistic proton distributions in dense radiation fields. Three independent techniques have been developed to simulate these non-linear problems: the kinetic equation approach; the phase-space density (PSD) Monte Carlo method; and the large-particle (LP) Monte Carlo method. In this paper, we present the latest version of the LP method and compare it with the other methods. The efficiency of the method in treating geometrically complex problems is illustrated by showing results of simulations of 1D, 2D and 3D systems. The method is shown to be powerful enough to treat non-spherical geometries, including such effects as bulk motion of the background plasma, reflection of radiation from cold matter, and anisotropic distributions of radiating particles. It can therefore be applied to simulate high-energy processes in such astrophysical systems as accretion discs with coronae, relativistic jets, pulsar magnetospheres and gamma-ray bursts.

Stern, Boris E.↗
NASA NTRS1989

Excitation of kinetic Alfven waves by resonant mode conversion and longitudinal heating of magnetized plasmas

The excitation of the kinetic Alfven wave by resonant mode conversion and longitudinal heating of the plasma by the kinetic Alfven wave were demonstrated on the basis of a macroscale particle simulation. The longitudinal electron current was shown to be cancelled by the ions. The kinetic Alfven wave produced an ordered motion of the plasma particles in the wave propagation direction. The electrons were pushed forward along the ambient magnetic field by absorbing the kinetic Alfven wave through the Landau resonance.

Tanaka, Motohiko↗
NASA NTRS2009

Orion Landing Simulation Eight Soil Model Comparison

LS-DYNA finite element simulations of a rigid Orion Crew Module (CM) were used to investigate the CM impact behavior on eight different soil models. Ten different landing conditions, characterized by the combination of CM vertical and horizontal velocity, hang angle, and roll angle were simulated on the eight different soils. The CM center of gravity accelerations, pitch angle, kinetic energy, and soil contact forces were the outputs of interest. The simulation results are presented, with comparisons of the CM behavior on the different soils. The soils analyzed in this study can be roughly categorized as soft, medium, or hard, according to the CM accelerations that occur when landing on them. The soft group is comprised of the Carson Sink Wet soil and the Kennedy Space Center (KSC) Low Density Dry Sand. The medium group includes Carson Sink Dry, the KSC High Density In-Situ Moisture Sand and High Density Flooded Sand, and Cuddeback B. The hard soils are Cuddeback A and the Gantry Unwashed Sand. The softer soils were found to produce lower peak accelerations, have more stable pitch behavior, and to be less sensitive to the landing conditions. This investigation found that the Cuddeback A soil produced the highest peak accelerations and worst stability conditions, and that the best landing performance was achieved on the KSC Low Density Dry Sand.

Mark, Stephen D.↗
NASA NTRS1984

Studies of new media radiation induced laser

Various lasants were investigated especially, 2-iodohepafluoropropane (i-C3F7I) for the direct solar pumped lasers. Optical pumping of iodine laser was achieved using a small flashlamp. Using i-C3F7I as a laser gain medium, threshold inversion density, small signal gain, and laser performance at the elevated temperature were measured. The experimental results and analysis are presented. The iodine laser kinetics of the C3F7I and IBr system were numerically simulated. The concept of a direct solar-pumped laser amplifier using (i-C3F7I) as the laser material was evaluated and several kinetic coefficients for i-C3F7I laser system were reexamined. The results are discussed.

Han, K. S.↗
NASA NTRSDate not supplied

State-Specific Kinetic Modeling for Predictions of Radiative Heating in H 2 /He Entry Flows

A 23-species State-to-State (StS) thermochemical model for H 2 /He mixtures is developed and used to simulate 1-D shocks at conditions relevant for ice and gas giant entry flows. To implement this StS model, a literature review of the state-specific excitation and ionization rates of atomic H is first performed. While electron impact rates from various sources are found to be in relatively good agreement, large discrepancies are found in the limited data available on heavy particle impact rates. To investigate the impact of these discrepancies on the predicted shock and radiance profiles, simulations are performed for three different versions of the kinetic model, each with a different set of heavy particle impact rates from the literature. Of these three models, two are found to over-predict and one is found to under-predict the radiance values seen in corresponding experimental data from the NASA Ames Electric Arc Shock Tube (EAST) facility. A sensitivity analysis to the electron and heavy particle impact rates is then performed, from which an improved kinetic model with modified rates is proposed. This modified kinetic model is found to reproduce the total integrated radiance seen in the EAST experiments accurately. This suggests that the associated radiative heating that would be encountered by an entry probe in these environments would be well predicted by this kinetic model. However, some discrepancies are still found i the relative contributions to the total radiance from different spectral features, suggesting that additional improvements to the kinetic model still need to be made.

Outer Planet↗
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↗
NASA NTRS1993

Direct numerical simulation of isotropic turbulence interacting with a weak shock wave

Direct numerical simulations are used to investigate the interaction of isotropic quasi-incompressible turbulence with a weak shock wave. A linear analysis of the interaction is conducted for comparison with the simulations. Both the simulations and the analysis show that turbulence is enhanced during the interaction. Turbulent kinetic energy and transverse vorticity components are amplified, and turbulent lengthscales are decreased. It is suggested that the amplification mechanism is primarily linear. Simulations also showed a rapid evolution of turbulent kinetic energy just downstream of the shock, a behavior not reproduced by the linear analysis. Analysis of the budget of the turbulent kinetic energy transport equation shows that this behavior can be attributed to the pressure transport term. Multiple compression peaks were found along the mean streamlines at locations where the local shock thickness had increased significantly.

Lee, Sangsan↗
NASA NTRS2004

Chemical models for simulating single-walled nanotube production in arc vaporization and laser ablation processes

Chemical kinetic models for the nucleation and growth of clusters and single-walled carbon nanotube (SWNT) growth are developed for numerical simulations of the production of SWNTs. Two models that involve evaporation and condensation of carbon and metal catalysts, a full model involving all carbon clusters up to C80, and a reduced model are discussed. The full model is based on a fullerene model, but nickel and carbon/nickel cluster reactions are added to form SWNTs from soot and fullerenes. The full model has a large number of species--so large that to incorporate them into a flow field computation for simulating laser ablation and arc processes requires that they be simplified. The model is reduced by defining large clusters that represent many various sized clusters. Comparisons are given between these models for cases that may be applicable to arc and laser ablation production. Solutions to the system of chemical rate equations of these models for a ramped temperature profile show that production of various species, including SWNTs, agree to within about 50% for a fast ramp, and within 10% for a slower temperature decay time.

Nanotubes, Carbon/chemistry↗
NASA NTRS2010

Numerical Simulation of Vitiation Effects on a Hydrogen-Fueled Dual-Mode Scramjet

The Wind-US computational fluid dynamics (CFD) flow solver was used to simulate dual-mode direct-connect ramjet/scramjet engine flowpath tests conducted in the University of Virginia (UVa) Supersonic Combustion Facility (SCF). The objective was to develop a computational capability within Wind-US to aid current hypersonic research and provide insight to flow as well as chemistry details that are not resolved by instruments available. Computational results are compared with experimental data to validate the accuracy of the numerical modeling. These results include two fuel-off non-reacting and eight fuel-on reacting cases with different equivalence ratios, split between one set with a clean (non-vitiated) air supply and the other set with a vitiated air supply (12 percent H2O vapor). The Peters and Rogg hydrogen-air chemical kinetics model was selected for the scramjet simulations. A limited sensitivity study was done to investigate the choice of turbulence model and inviscid flux scheme and led to the selection of the k-epsilon model and Harten, Lax and van Leer (for contact waves) (HLLC) scheme for general use. Simulation results show reasonably good agreement with experimental data and the overall vitiation effects were captured.

Vyas, Manan A.↗
NASA NTRS2003

MHD Turbulence at Moderate Magnetic Reynolds Number

In the present article we will consider the decay of MHD turbulence under the influence of a strong external magnetic field at moderate magnetic Reynolds numbers. Typical values of R(sub m) that are considered here range from R(sub m) approx. 0.1 to R(sub m) approx. 20. As a comparison, the initial kinetic Reynolds number common to all our simulations is Re(sub L) = 199. This means that the range of Prandtl numbers explored is 5 x 10(exp -4) to 10(exp -1). Our motivation is mainly to exhibit how the transition from the QS approximation to FMHD occurs. At the lowest values of R(sub m) studied here, the QS approximation is shown to model the flow faithfully. However, for the higher values of R(sub m) considered, it is clearly inadequate but can be replaced by another approximation which will be referred to as the Quasi-Linear (QL) approximation. Another objective of the present study is to describe how variations in the magnetic Reynolds number (while maintaining all other parameters constant) affect the dynamics of the flow. This complements past studies where variations in either the strength of the external magnetic field or the kinetic Reynolds number were considered. This article is organized as follows. In section 2 we recall the definition of the quasi-static approximation. Section 3 is devoted to the description of the numerical experiments performed using the quasi-static approximation and full MHD. In section 4 we describe the quasi-linear approximation and test it numerically against full MHD. A concluding summary is given in section 5.

Knaepen, B.↗