Results for “Kinetic simulations” — Page 11

NASA NTRS1999

Generation of Electron Suprathermal Tails in the Upper Solar Atmosphere: Implications for Coronal Heating

We present a mechanism for the generation of non-Maxwellian electron distribution function in the upper regions of the solar atmosphere in the presence of collisional damping. It is suggested that finite amplitude, low frequency, obliquely propagating electromagnetic waves can carry a substantial electric field component parallel to the mean magnetic field that can be significantly larger than the Dreicer electric field. This long wavelength electric fluctuation is capable of generating high frequency electron plasma oscillations and low frequency ion acoustic-like waves. The analysis has been performed using 1-1/2D Vlasov and PIC numerical simulations in which both electrons and ions are treated kinetically and self consistently. The simulation results indicate that high frequency electron plasma oscillations and low frequency ion acoustic-like waves are generated. The high frequency electron plasma oscillation drives electron plasma turbulence, which subsequently is damped out by the background electrons. The turbulence damping results in electron acceleration and the generation of non-Maxwellian suprathermal tails on time scales short compared to collisional damping. Bulk heating also occurs if the fluctuating parallel electric field is strong enough. This study suggests that finite amplitude, low frequency, obliquely propagating, electromagnetic waves can play a significant role in the acceleration and heating of the solar corona electrons and in the coupling of medium and small-scale phenomena.

Vinas, Adolfo F.↗

NASA NTRS2011

The Distribution of Chromosomal Aberrations in Human Cells Predicted by a Generalized Time-Dependent Model of Radiation-Induced Formation of Aberrations

New experimental data show how chromosomal aberrations for low- and high-LET radiation are dependent on DSB repair deficiencies in wild-type, AT and NBS cells. We simulated the development of chromosomal aberrations in these cells lines in a stochastic track-structure-dependent model, in which different cells have different kinetics of DSB repair. We updated a previously formulated model of chromosomal aberrations, which was based on a stochastic Monte Carlo approach, to consider the time-dependence of DSB rejoining. The previous version of the model had an assumption that all DSBs would rejoin, and therefore we called it a time-independent model. The chromosomal-aberrations model takes into account the DNA and track structure for low- and high-LET radiations, and provides an explanation and prediction of the statistics of rare and more complex aberrations. We compared the program-simulated kinetics of DSB rejoining to the experimentally-derived bimodal exponential curves of the DSB kinetics. We scored the formation of translocations, dicentrics, acentric and centric rings, deletions, and inversions. The fraction of DSBs participating in aberrations was studied in relation to the rejoining time. Comparisons of simulated dose dependence for simple aberrations to the experimental dose-dependence for HF19, AT and NBS cells will be made.

Ponomarev, Artem L.↗

NASA NTRS1976

Introductory study of the chemical behavior of jet emissions in photochemical smog

Jet aircraft emissions data from the literature were used as initial conditions for a series of computer simulations of photochemical smog formation in static air. The chemical kinetics mechanism used in these simulations was an updated version which contains certain parameters designed to account for hydrocarbon reactivity. These parameters were varied to simulate the reaction rate constants and average carbon numbers associated with the jet emissions. The roles of surface effects, variable light sources, NO/NO2 ratio, continuous emissions, and untested mechanistic parameters were also assessed. The results of these calculations indicate that the present jet emissions are capable of producing oxidant by themselves. The hydrocarbon/nitrous oxides ratio of present jet aircraft emissions is much higher than that of automobiles. These two ratios appear to bracket the hydrocarbon/nitrous oxides ratio that maximizes ozone production. Hence an enhanced effect is seen in the simulation when jet exhaust emissions are mixed with automobile emissions.

Whitten, G. Z.↗

NASA NTRS1971

Photoactivation of the manganese catalyst of O2 evolution. II A two-quantum mechanism.

A kinetic analysis is presented of the data obtained in previously described kinetic experiments relating to the photoactivation of the Mn complex required for photosynthetic O2 evolution. The results obtained from the computer simulation of the derived kinetic model compare favorably with the experimentally derived values. They also suggest that photoactivation can be rationalized in terms of a two-quantum process.

Radmer, R.↗

NASA NTRS1990

Relaxation processes in a turbulent compressible magnetofluid

The compressible extensions of time asymptotic relaxation states of incompressible two-dimensional magnetohydrodynamic turbulence are studied. A polytropic equation of state is used with viscous and resistive dissipation. The incompressible case is known to allow three distinct time asymptotic types of behavior: magnetic energy dominated relaxation, kinetic energy dominated relaxation, and cross helicity dominated relaxation. At low Mach numbers the incompressible scenario is reproducible from the compressible simulations, and compressibility plays only a secondary role. At moderate, but still subsonic, Mach numbers the distinct incompressible processes are still recognizable, but strong compressibility features dominate the high-wave-number regime of several simulations. In particular, the magnetic and kinetic energy dominated simulations display regions of strong acoustic turbulence near the dissipation scale.

Ghosh, S.↗

NASA NTRS2012

Numerical Simulation of Rocket Exhaust Interaction with Lunar Soil

This technology development originated from the need to assess the debris threat resulting from soil material erosion induced by landing spacecraft rocket plume impingement on extraterrestrial planetary surfaces. The impact of soil debris was observed to be highly detrimental during NASA s Apollo lunar missions and will pose a threat for any future landings on the Moon, Mars, and other exploration targets. The innovation developed under this program provides a simulation tool that combines modeling of the diverse disciplines of rocket plume impingement gas dynamics, granular soil material liberation, and soil debris particle kinetics into one unified simulation system. The Unified Flow Solver (UFS) developed by CFDRC enabled the efficient, seamless simulation of mixed continuum and rarefied rocket plume flow utilizing a novel direct numerical simulation technique of the Boltzmann gas dynamics equation. The characteristics of the soil granular material response and modeling of the erosion and liberation processes were enabled through novel first principle-based granular mechanics models developed by the University of Florida specifically for the highly irregularly shaped and cohesive lunar regolith material. These tools were integrated into a unique simulation system that accounts for all relevant physics aspects: (1) Modeling of spacecraft rocket plume impingement flow under lunar vacuum environment resulting in a mixed continuum and rarefied flow; (2) Modeling of lunar soil characteristics to capture soil-specific effects of particle size and shape composition, soil layer cohesion and granular flow physics; and (3) Accurate tracking of soil-borne debris particles beginning with aerodynamically driven motion inside the plume to purely ballistic motion in lunar far field conditions. In the earlier project phase of this innovation, the capabilities of the UFS for mixed continuum and rarefied flow situations were validated and demonstrated for lunar lander rocket plume flow impingement under lunar vacuum conditions. Applications and improvements to the granular flow simulation tools contributed by the University of Florida were tested against Earth environment experimental results. Requirements for developing, validating, and demonstrating this solution environment were clearly identified, and an effective second phase execution plan was devised. In this phase, the physics models were refined and fully integrated into a production-oriented simulation tool set. Three-dimensional simulations of Apollo Lunar Excursion Module (LEM) and Altair landers (including full-scale lander geometry) established the practical applicability of the UFS simulation approach and its advanced performance level for large-scale realistic problems.

Liever, Peter↗

NASA NTRS2012

Motion Imagery and Robotics Application (MIRA): Standards-Based Robotics

This technology development originated from the need to assess the debris threat resulting from soil material erosion induced by landing spacecraft rocket plume impingement on extraterrestrial planetary surfaces. The impact of soil debris was observed to be highly detrimental during NASA s Apollo lunar missions and will pose a threat for any future landings on the Moon, Mars, and other exploration targets. The innovation developed under this program provides a simulation tool that combines modeling of the diverse disciplines of rocket plume impingement gas dynamics, granular soil material liberation, and soil debris particle kinetics into one unified simulation system. The Unified Flow Solver (UFS) developed by CFDRC enabled the efficient, seamless simulation of mixed continuum and rarefied rocket plume flow utilizing a novel direct numerical simulation technique of the Boltzmann gas dynamics equation. The characteristics of the soil granular material response and modeling of the erosion and liberation processes were enabled through novel first principle-based granular mechanics models developed by the University of Florida specifically for the highly irregularly shaped and cohesive lunar regolith material. These tools were integrated into a unique simulation system that accounts for all relevant physics aspects: (1) Modeling of spacecraft rocket plume impingement flow under lunar vacuum environment resulting in a mixed continuum and rarefied flow; (2) Modeling of lunar soil characteristics to capture soil-specific effects of particle size and shape composition, soil layer cohesion and granular flow physics; and (3) Accurate tracking of soil-borne debris particles beginning with aerodynamically driven motion inside the plume to purely ballistic motion in lunar far field conditions.

Martinez, Lindolfo↗

NASA NTRS1995

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗

NASA NTRS2003

Temporal Evolution of Nanostructures in a Model Nickel-Base Superalloy: Experiments and Simulations

The temporal evolution of the nanostructure of a model Ni-base superalloy (Ni-5.2 at.% Al-14.2 at.% Cr) is studied experimentally employing three-dimensional atom-probe (3DAP) microscopy in conjunction with kinetic Monte Carlo (KMC) simulations at 600 C. It is demonstrated that not only can the mean compositions of individual gamma' (Ni3Al with the Li2 structure) precipitates be measured but the Ni, Al, and Cr concentration profiles within the precipitates can also be determined for precipitates with a mean radius ( ) as small as 0.85 nm. The three asymptotic tinie dependencies of the Lifshitz-Slyzov-Wagner (LSW) theory of coarsening (Ostwald ripening) are measured and found to deviate from its theoretical predictions: possible explanations for these discrepancies are discussed. At 0.25 hr. there is 3DAP microscope evidence for the presence of precipitates of another nickel-rich phase. ="Ni3Cr" (Ni3Cr(1-x)Al(x)), which exhibits short-range order (SRO) and that is metastable with respect to Ni3Al. This metastable phase is also found by KMC simulations and has the composition Ni3Cr(1-x)Al(x), which is Ni-2.91 at.% Al-21.98 at.% Cr at 16 hours. Our results demonstrate that the decomposition of gamma the primary gamma (FCC) phase results in the concurrent formation of an ordered phase and a disordered phase by 0.25 hours.

Sudbrack, Chantal K.↗

NASA NTRS2022

Spectral Power-law Formation by Sequential Particle Acceleration in Multiple Flare Magnetic Islands

We present a first-principles model of pitch-angle and energy distribution function evolution as particles are sequentially accelerated by multiple flare magnetic islands. Data from magnetohydrodynamic (MHD) simulations of an eruptive flare/coronal mass ejection provide ambient conditions for the evolving particle distributions. Magnetic islands, which are created by sporadic reconnection at the self-consistently formed flare current sheet, contract and accelerate the particles. The particle distributions are evolved using rules derived in our previous work. In this investigation, we assume that a prescribed fraction of particles sequentially "hops" to another accelerator and receives an additional boost in energy and anisotropy. This sequential process generates particle number spectra that obey an approximate power law at mid-range energies and presents low- and high-energy breaks. We analyze these spectral regions as functions of the model parameters. We also present a fully analytic method for forming and interpreting such spectra, independent of the sequential acceleration model. The method requires only a few constrained physical parameters, such as the percentage of particles transferred between accelerators, the energy gain in each accelerator, and the number of accelerators visited. Our investigation seeks to bridge the gap between MHD and kinetic regimes by combining global simulations and analytic kinetic theory. The model reproduces and explains key characteristics of observed flare hard X-ray spectra as well as the underlying properties of the accelerated particles. Our analytic model provides tools to interpret high-energy observations for missions and telescopes, such as RHESSI, FOXSI, NuSTAR, Solar Orbiter, EOVSA, and future high-energy missions.

magnetic reconnection↗

NASA NTRS2006

Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model

The global dynamics of the ionized and neutral gases in the environment of Io plays an important role in the interaction of Jupiter s corotating magnetospheric plasma with Io. Stationary simulations of this problem have already been done using the magnetohydrodynamics (MHD) and the electrodynamics approaches. One of the major results of recent simplified two-fluid model simulations [Saur, J., Neubauer, F.M., Strobel, D.F., Summers, M.E., 2002. J. Geophys. Res. 107 (SMP5), 1-18] was the production of the structure of the double-peak in the magnetic field signature of the Io flyby. These could not be explained before by standard MHD models. In this paper, we present a hybrid simulation for Io with kinetic ions and fluid electrons. This method employs a fluid description for electrons and neutrals, whereas for ions a particle approach is used. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. Our model may provide a much more accurate description for the ion dynamics than previous approaches and allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. The first results of such a simulation of the dynamics of ions in Io s environment are discussed in this paper. Comparison with the Galileo IO flyby results shows that this approach provides an accurate physical basis for the interaction and can therefore naturally reproduce all the observed salient features.

Lipatov, Alexander S.↗

NASA NTRS1975

Predicted exhaust emissions from a methanol and jet fueled gas turbine combustor

A computer model of a gas turbine combustor has been used to predict the kinetic combustion and pollutant formation processes for methanol and simulated jet fuel. Use of the kinetic reaction mechanisms has also allowed a study of ignition delay and flammability limit of these two fuels. The NOX emissions for methanol were predicted to be from 69 to 92% lower than those for jet fuel at the same equivalence ratio which is in agreement with experimentally observed results. The high heat of vaporization of methanol lowers both the combustor inlet mixture temperatures and the final combustion temperatures. The lower combustion temperatures lead to low NOX emissions while the lower inlet mixture temperatures increase methanol's ignition delay. This increase in ignition delay dictates the lean flammability limit of methanol to be 0.8, while jet fuel is shown to combust at 0.4.

Adelman, H. G.↗

NASA NTRS2018

A High-Order Kinetic Energy Conserving Scheme for Compressible Large-Eddy Simulation

No abstract available

compressible↗

NASA NTRS2012

Multi-Scale Modeling of Liquid Phase Sintering Affected by Gravity: Preliminary Analysis

A multi-scale simulation concept taking into account impact of gravity on liquid phase sintering is described. The gravity influence can be included at both the micro- and macro-scales. At the micro-scale, the diffusion mass-transport is directionally modified in the framework of kinetic Monte-Carlo simulations to include the impact of gravity. The micro-scale simulations can provide the values of the constitutive parameters for macroscopic sintering simulations. At the macro-scale, we are attempting to embed a continuum model of sintering into a finite-element framework that includes the gravity forces and substrate friction. If successful, the finite elements analysis will enable predictions relevant to space-based processing, including size and shape and property predictions. Model experiments are underway to support the models via extraction of viscosity moduli versus composition, particle size, heating rate, temperature and time.

Olevsky, Eugene↗

NASA NTRSDate not supplied

Solid-state Precipitation in Molecular Dynamics: KMC-MD Hybrid Simulations

A critical drawback of molecular dynamics (MD) simulations when modeling physical properties is the limitation to short timescales. Solid state precipitation is one behavior that is not significantly observable at timescales accessible to MD, and thus has been difficult to simulate at an MD scale. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation during MD simulation implemented in LAMMPS. This method incorporates the atomic neighbor hopping behavior utilized by KMC simulations into a MD simulation with a variable local lattice configuration. By accelerating nearest-neighbor atomic hopping with this technique, solid-state precipitation can be observed at an MD timescale. Utilizing theory based on persistent-embryo methods for liquid-solid precipitation in MD we demonstrate the calculation of critical nucleus size for the gamma’ precipitation out of gamma phase in the Ni-Al system.

Molecular Dynamics↗

NASA NTRSDate not supplied

Solid-state Precipitation in Molecular Dynamics: kMC-MD Hybrid Simulations

A critical drawback of molecular dynamics (MD) simulations when modeling physical properties is the limitation to short timescales. Solid state precipitation is one behavior that is not significantly observable at timescales accessible to MD, and thus has been difficult to simulate at an MD scale. Taking inspiration from Kinetic Monte-Carlo (KMC) methods for simulating diffusion on a rigid lattice, we have developed a hybrid KMC-MD approach to enable simulation of solid-state precipitation during MD simulation implemented in LAMMPS. This method incorporates the atomic neighbor hopping behavior utilized by KMC simulations into a MD simulation with a variable local lattice configuration. By accelerating nearest-neighbor atomic hopping with this technique, solid-state precipitation can be observed at an MD timescale. Utilizing theory based on persistent-embryo methods for liquid-solid precipitation in MD we demonstrate the calculation of critical nucleus size for the gamma’ precipitation out of gamma phase in the Ni-Al system.

Molecular Dynamics↗

NASA NTRS1992

Measurement and analysis of a small nozzle plume in vacuum

Measurements of Pitot pressure and flow angle were made in the plume of a nozzle flowing nitrogen and exhausting to a vacuum. The measurements were compared to results from a numerical simulation of the flow that was based on kinetic theory and used the direct-simulation Monte Carlo (DSMC) method. Numerical results were compared with measurements made in the plume at various axial and radial stations. Total pressure measurements were made with Pitot tubes sized for specific regions of the plume. Flow angle measurements were made with a conical probe. The measurement area for flow angle extended to 160 mm (5 exit diameters) downstream of the nozzle exit plane and radially to 60 mm (1.9 exit diameters) from the plume axis. The total pressure measurements extended 480 mm (16 exit diameters) downstream and radially to 60 mm. Comparisons of computed results from the DSMC method with measurements of flow angle displayed improved agreement with increasing distance from the exit plane. Pitot pressures computed from the DSMC method were in reasonably good agreement with experimental results over the entire measurement area.

Penko, Paul F.↗

NASA NTRS1993

Measurement and analysis of a small nozzle plume in vacuum

Pitot pressures and flow angles are measured in the plume of a nozzle flowing nitrogen and exhausting to a vacuum. Total pressures are measured with Pitot tubes sized for specific regions of the plume and flow angles measured with a conical probe. The measurement area for total pressure extends 480 mm (16 exit diameters) downstream of the nozzle exit plane and radially to 60 mm (1.9 exit diameters) off the plume axis. The measurement area for flow angle extends to 160 mm (5 exit diameters) downstream and radially to 60 mm. The measurements are compared to results from a numerical simulation of the flow that is based on kinetic theory and uses the direct-simulation Monte Carlo (DSMC) method. Comparisons of computed results from the DSMC method with measurements of flow angle display good agreement in the far-field of the plume and improve with increasing distance from the exit plane. Pitot pressures computed from the DSMC method are in reasonably good agreement with experimental results over the entire measurement area.

Penko, P. F.↗