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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Core Physics and Kinetics Calculations for the Fissioning Plasma Core Reactor

Highly efficient, compact nuclear reactors would provide high specific impulse spacecraft propulsion. This analysis and numerical simulation effort has focused on the technical feasibility issues related to the nuclear design characteristics of a novel reactor design. The Fissioning Plasma Core Reactor (FPCR) is a shockwave-driven gaseous-core nuclear reactor, which uses Magneto Hydrodynamic effects to generate electric power to be used for propulsion. The nuclear design of the system depends on two major calculations: core physics calculations and kinetics calculations. Presently, core physics calculations have concentrated on the use of the MCNP4C code. However, initial results from other codes such as COMBINE/VENTURE and SCALE4a. are also shown. Several significant modifications were made to the ISR-developed QCALC1 kinetics analysis code. These modifications include testing the state of the core materials, an improvement to the calculation of the material properties of the core, the addition of an adiabatic core temperature model and improvement of the first order reactivity correction model. The accuracy of these modifications has been verified, and the accuracy of the point-core kinetics model used by the QCALC1 code has also been validated. Previously calculated kinetics results for the FPCR were described in the ISR report, "QCALC1: A code for FPCR Kinetics Model Feasibility Analysis" dated June 1, 2002.

Butler, C.↗

NASA NTRS1986

Oxidation and embrittlement of Ti-6Al-2Sn-4Zr-2Mo alloy

A thermogravimetric analysis monitoring study has been conducted for Ti-6Al-2Sn-4Zr-2Mo in the 593-760 C range, with a view to its total oxidation kinetics in laboratory air. In addition, tensile elongation was obtained for selected foil specimens after their exposure to simulated Space Shuttle reentry conditions. The total weight gain with time exhibited two parabolic stages separated by a transient region which is the result of a simultaneous increase in the parabolic growth rate constant for the oxide and in oxygen solubility at the oxide-metal interface. A model for the total oxidation kinetics accounting for the two individual components (oxide growth and solid solution formation) is proposed; the activation energy for oxygen diffusion in the alpha-solid solution is shown to approximate the activation energy for the degradation of tensile elongation of the alloy in its foil gage condition.

Shenoy, R. N.↗

NASA NTRS2023

Hybrid-Vlasov Simulation of Soft X-Ray Emissions at the Earth’S Dayside Magnetospheric Boundaries

Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause, the magnetosheath and the polar cusps by remote sensing techniques. The Solar wind–Magnetosphere–Ionosphere Link Explorer (SMILE) and Lunar Environment heliospheric X-ray Imager (LEXI) missions aim to obtain soft X-ray images of near-Earth space thanks to their Soft X-ray Imager (SXI) instruments. While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission, the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma. To investigate the possible signatures of ion-kinetic-scale processes in soft X-ray images, we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model. The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field. We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations, providing meridional and equatorial views. We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirror mode wave structures in the magnetosheath and flux transfer events (FTEs) at the magnetopause. Our results suggest that, although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers, mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images. For instance, a local increase by 30% in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12% enhancement in the line-of-sight- and time integrated soft X-ray emissivity originating from this region. Likewise, a proton density increase by 14% in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%. These are likely conservative estimates, given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind. These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale, transient phenomena occurring on the dayside.

magnetosphere↗

NASA NTRS2022

Hybrid Kinetic Model of the Interaction Between the Dense Plasma Clouds and Magnetospheric Plasma on Large Time and Spatial Scales, and Comparison With MMS Observations

We present a new simulation results of the cloud dynamics in the ambient magnetospheric plasma on the large time and spatial scales. It was assumed that these impulsive structures observed by the MMS spacecraft originally were created because of the reconnection at the magnetopause. Our new 3-D hybrid kinetic modeling on the large time and spatial scales captures several of these processes: an excitation of the electromagnetic waves (whistler and shear-Alfvén waves) and plasma instabilities (mirror and flute); a formation of shock waves, and collapsing diamagnetic cavity; particle acceleration. A strong overshoot in plasma density profile was observed in the modeling and MMS observation at the interface between the cloud and magnetospheric plasma. The cloud expansion into ambient magnetospheric plasma causes the flute waves connected with excitation of the Rayleigh-Taylor instability observed at the overshoot in plasma density profile across the external magnetic field. The modeling demonstrates a formation of the whistler waves at the initial stage which propagate in the external magnetic field direction. At the later stage, a formation of shear-Alfvén waves was observed

A. S. Lipatov↗

NASA NTRS1993

Performance assessment of low pressure nuclear thermal propulsion

A low pressure nuclear thermal propulsion (LPNTP) system, which takes advantage of hydrogen dissociation/recombination, was proposed as a means of increasing engine specific impulse (Isp). The effect of hydrogen dissociation/recombination on LPNTP Isp is examined. A two-dimensional computer model was used to show that the optimum chamber pressure is approximately 100 psia (at a chamber temperature of 3,000 K), with an Isp approximately 15 s higher than at 1,000 psia. At high chamber temperatures and low chamber pressures, the increase in Isp is due to both lower average molecular weights caused by dissociation and added kinetic energy from monatomic hydrogen recombination. Monatomic hydrogen recombination increases the Isp more then hydrogen dissociation. Variations in the mole fraction of monatomic hydrogen are similar to variations in static pressure along the axial nozzle position. Most recombination occurs close to the nozzle throat. Practical variations in nozzle geometry have minimal impact on recombination. Other models which can simulate a wider range of nozzle designs should be used in the future. The uncertainty of the hydrogen kinetic reaction rates at high temperatures (approximately 3,000 K) affects the accuracy of the analysis and should be verified with simple bench tests.

Gerrish, H. P., Jr.↗

NASA NTRS1990

Magnetic reconnection in collisionless plasmas - Prescribed fields

The structure of the dissipation region during magnetic reconnection in collisionless plasma is investigated by examining a prescribed two-dimensional magnetic x line configuration with an imposed inductive electric field E(y). The calculations represent an extension of recent MHD simulations of steady state reconnection (Biskamp, 1986; Lee and Fu, 1986) to the collisionless kinetic regime. It is shown that the structure of the x line reconnection configuration depends on only two parameters: a normalized inductive field and a parameter R which represents the opening angle of the magnetic x lines.

Burkhart, G. R.↗

NASA NTRS1993

An investigation of small scales of turbulence in a boundary layer at high Reynolds numbers

The assumption that turbulence at large wave-numbers is isotropic and has universal spectral characteristics which are independent of the flow geometry, at least for high Reynolds numbers, has been a cornerstone of closure theories as well as of the most promising recent development in the effort to predict turbulent flows, viz. large eddy simulations. This hypothesis was first advanced by Kolmogorov based on the supposition that turbulent kinetic energy cascades down the scales (up the wave-numbers) of turbulence and that, if the number of these cascade steps is sufficiently large (i.e. the wave-number range is large), then the effects of anisotropies at the large scales are lost in the energy transfer process. Experimental attempts were repeatedly made to verify this fundamental assumption. However, Van Atta has recently suggested that an examination of the scalar and velocity gradient fields is necessary to definitively verify this hypothesis or prove it to be unfounded. Of course, this must be carried out in a flow with a sufficiently high Reynolds number to provide the necessary separation of scales in order unambiguously to provide the possibility of local isotropy at large wave-numbers. An opportunity to use our 12-sensor hot-wire probe to address this issue directly was made available at the 80'x120' wind tunnel at the NASA Ames Research Center, which is normally used for full-scale aircraft tests. An initial report on this high Reynolds number experiment and progress toward its evaluation is presented.

Wallace, James M.↗

NASA NTRS1994

A Monte Carlo study of photochemical model uncertainties affecting the calculation of subsonic perturbations in the upper troposphere

Monte Carlo methods are frequently applied to the evaluation of uncertainties in models with multiple inputs that themselves have associated imprecisions. In the case of photochemical models used to evaluate changes in O3 or OH, inputs analyzed include rate coefficients measured in the laboratory, and chemical and physical constituents measured in the atmosphere. The Monte Carlo method was used with the 1-dimensional GSFC tropospheric photochemical model to examine uncertainty propagation to calculation of NO(x) and the major oxidants (O3, OH) in the upper troposphere. In all cases chemical kinetics inputs are varied and the NO(x) perturbation of a subsonic fleet is simulated. A series of model runs is used to explore sensitivities of model-computed parameters to other parameters: heterogeneous processes, the uncertainty in upper tropospheric H2O vapor measurements, aircraft emissions at different latitudes.

Thompson, Anne M.↗

NASA NTRS1999

Compaction as the Origin of the Unusual Craters on the Asteroid Mathilde

Asteroid Mathilde has been pummeled by at least five giant impacts (Figure 1). Previous experience with cratering suggests Mathilde's giant craters should each be surrounded by kilometer-deep blankets of ejecta, i.e. material excavated during the impact events. Curiously, there appears to be very little ejecta around Mathilde's craters; they show no evidence of filling by ejecta from adjacent large craters. A previous explanation for the missing ejecta, based on computer simulations, is that Mathilde's unusually high porosity (50 +/- 20%) confines the deposited impact kinetic energy to a localized volume, and produces excavation velocities so high (greater than approximately 20m/s) that nearly all ejecta escape Mathilde's gravitational field. Here we report on laboratory experiments in a highly porous material that give a different explanation. The crater is formed primarily by compaction, not excavation. The small amount of material that is lofted has velocities and ranges so small that nearly all of it is re-deposited within the crater bowl, thereby sparing neighboring craters from ejecta in-filling. This peculiar style of cratering implies that highly porous asteroids are minor contributors of meteorites, because essentially no ejecta escape these asteroids.

Housen, Kevin R.↗

NASA NTRS2004

Two Formal Gas Models For Multi-Agent Sweeping and Obstacle Avoidance

The task addressed here is a dynamic search through a bounded region, while avoiding multiple large obstacles, such as buildings. In the case of limited sensors and communication, maintaining spatial coverage - especially after passing the obstacles - is a challenging problem. Here, we investigate two physics-based approaches to solving this task with multiple simulated mobile robots, one based on artificial forces and the other based on the kinetic theory of gases. The desired behavior is achieved with both methods, and a comparison is made between them. Because both approaches are physics-based, formal assurances about the multi-robot behavior are straightforward, and are included in the paper.

Kerr, Wesley↗

NASA NTRS2010

Testing of an Integrated Reactor Core Simulator and Power Conversion System with Simulated Reactivity Feedback

A Direct Drive Gas-Cooled (DDG) reactor core simulator has been coupled to a Brayton Power Conversion Unit (BPCU) for integrated system testing at NASA Glenn Research Center (GRC) in Cleveland, Ohio. This is a closed-cycle system that incorporates an electrically heated reactor core module, turboalternator, recuperator, and gas cooler. Nuclear fuel elements in the gas-cooled reactor design are replaced with electric resistance heaters to simulate the heat from nuclear fuel in the corresponding fast spectrum nuclear reactor. The thermodynamic transient behavior of the integrated system was the focus of this test series. In order to better mimic the integrated response of the nuclear-fueled system, a simulated reactivity feedback control loop was implemented. Core power was controlled by a point kinetics model in which the reactivity feedback was based on core temperature measurements; the neutron generation time and the temperature feedback coefficient are provided as model inputs. These dynamic system response tests demonstrate the overall capability of a non-nuclear test facility in assessing system integration issues and characterizing integrated system response times and response characteristics.

Bragg-Sitton, Shannon M.↗

NASA NTRS2009

Testing of an Integrated Reactor Core Simulator and Power Conversion System with Simulated Reactivity Feedback

A Direct Drive Gas-Cooled (DDG) reactor core simulator has been coupled to a Brayton Power Conversion Unit (BPCU) for integrated system testing at NASA Glenn Research Center (GRC) in Cleveland, OH. This is a closed-cycle system that incorporates an electrically heated reactor core module, turbo alternator, recuperator, and gas cooler. Nuclear fuel elements in the gas-cooled reactor design are replaced with electric resistance heaters to simulate the heat from nuclear fuel in the corresponding fast spectrum nuclear reactor. The thermodynamic transient behavior of the integrated system was the focus of this test series. In order to better mimic the integrated response of the nuclear-fueled system, a simulated reactivity feedback control loop was implemented. Core power was controlled by a point kinetics model in which the reactivity feedback was based on core temperature measurements; the neutron generation time and the temperature feedback coefficient are provided as model inputs. These dynamic system response tests demonstrate the overall capability of a non-nuclear test facility in assessing system integration issues and characterizing integrated system response times and response characteristics.

Bragg-Sitton, Shannon M.↗

NASA NTRS2003

Dynamical Study of Prebiotic Processing by Comet Impact

Comets and meteoroids that bombarded the Earth, other planets and moons are considered possible deliverers of prebiotic materials manufactured in space. Simultaneously, chemical processing can be initiated by the large kinetic energy imparted to the planetary atmosphere during comet/meteoroid entry. The efficiency of organic synthesis and the diversity of products by impact shock are known to decrease as the reducing power of the atmosphere decreases. It is generally accepted that rich organic products are produced in a methane atmosphere whereas a carbon dioxide atmosphere is reported to yield a dearth of products In order to understand the details of impact chemistry and how it depends on the atmospheric composition, we carried out simulations of the chemistry initiated by comet/meteoroid impact upon a planetary atmosphere using different atmospheric compositions. The simulations were done by solving the set of coupled equations for mass, momentum, and energy conservations, chemical kinetics, and transport, that describe a high-energy impact shock, subsequent expansion and cooling of the hot shocked gas by mixing with the ambient gas, and the eventual steady state composition.

Huo, Winifred M.↗

NASA NTRS1985

Simulation and modeling of homogeneous, compressed turbulence

Low Reynolds number homogeneous turbulence undergoing low Mach number isotropic and one-dimensional compression was simulated by numerically solving the Navier-Stokes equations. The numerical simulations were performed on a CYBER 205 computer using a 64 x 64 x 64 mesh. A spectral method was used for spatial differencing and the second-order Runge-Kutta method for time advancement. A variety of statistical information was extracted from the computed flow fields. These include three-dimensional energy and dissipation spectra, two-point velocity correlations, one-dimensional energy spectra, turbulent kinetic energy and its dissipation rate, integral length scales, Taylor microscales, and Kolmogorov length scale. Results from the simulated flow fields were used to test one-point closure, two-equation models. A new one-point-closure, three-equation turbulence model which accounts for the effect of compression is proposed. The new model accurately calculates four types of flows (isotropic decay, isotropic compression, one-dimensional compression, and axisymmetric expansion flows) for a wide range of strain rates.

Wu, C. T.↗

NASA NTRSDate not supplied

Comparison between Hayabusa 2 Spectral Measurements and Simulations

Recent Hayabusa 2 spectral measurements are compared with state-of-the-art shock-layer radiation simulations resulting from the LAURA/HARA code suite. These simulations include coupled ablation, which accounts for the injection of ablation products into the flowfield, and coupled radiation, which accounts for radiative energy loss in the flowfield. To enable the cou- pled ablation simulations, a best-estimate model is developed for Hayabusa’s carbon-phenolic ablator based on available published information. The comparison between the simulations and measurements focuses on three atomic nitrogen lines, two atomic oxygen lines, one atomic hydrogen line, and the CN Violet band system. Parametric uncertainties are evaluated for these features using a recently developed capability in LAURA/HARA. These uncertainties capture the impact of flowfield and radiation modeling uncertainties on the observed radiation simulations, and include the impact of coupled radiation and ablation product emission. The resulting parametric uncertainties range from 50 to over 100% for the atomic lines, and 80% for CN Violet. Uncertainties in flowfield kinetics and non-Boltzmann modeling provide the dominant contribution to these uncertainties. For the atomic nitrogen and oxygen lines, com- parisons between the nominal simulations and measurements show agreement well within these simulation uncertainty bounds, with the measurements within 20% of the nominal simulation over most of the trajectory, and within 10% at peak emission. This excellent agreement is un- precedented for observed radiation measurements of atomic lines, where previous Stardust and Hayabusa 1 comparisons were significantly worse. The improved agreement for these Hayabusa 2 comparisons is both the result of improved measurement quality and enhanced flowfield and radiation modeling. For the ablation product emission from the 656 nm atomic hydrogen line, the agreement between measurements and simulations is within ±20% for all trajectory points except one outlier. This good comparison is surprising considering this species depends on the ablation rate and is shown to have strong wake contribution, which both contribute to the large ≈±70% parametric uncertainty. For the ablation product emission from CN Violet, the simulations over-predict the measurements by over a factor of two early in the trajectory, with the comparison improving later in the trajectory. This disagreement for CN Violet is consistent with the ≈±80% parametric uncertainty bounds evaluated for the simulations. To utilize the excellent agreement between the atomic line simulations and measurements in the assessment of the radiative heating margin for a future flight vehicle, the flowfield property binning (FPB) approach is applied. Using the Mars Sample Return (MSR) Earth Entry System (EES) as an example, the combination of the excellent Hayabusa comparisons and the FPB analysis justify nearly a 10% reduction in the radiative heating margin from an assumed baseline margin of 30%. This direct quantitative link between observed radiation measurements and the radia- tive heating margin of future flight vehicles demonstrates the value of these observed radiation measurements.

Shock-Layer Radiative Heating↗

NASA NTRS2018

Tuning Process for the Modified Magnussen Combustion Model

In the application of CFD to turbulent reacting flows, one of the main limitations to predictive accuracy is the chemistry model. Using a full or skeletal kinetics model may provide good predictive ability, however, at considerable computational cost. Adding the ability to account for the interaction between turbulence and chemistry improves the overall fidelity of a simulation but adds to this cost. An alternative is the use of simple models, such as the Magnussen model, which has negligible computational overhead, but lacks general predictive ability except for cases that can be tuned to the flow being solved. In this paper, a technique will be described that allows the tuning of the Magnussen model for an arbitrary fuel and flow geometry without the need to have experimental data for a particular case. The tuning is based on comparing the results of the Magnussen model and full finite-rate chemistry when applied to perfectly and partially stirred reactor sim- ulations. In addition, a modification to the Magnussen model is proposed that allows the upper kinetic limit for the reaction rate to be set, giving better physical agreement with full kinetic mechanisms. In order to improve the agreement with flame temperatures, the thermal properties of the product species is adjusted to better match the mixture proper- ties of the full mechanism. The combustion model is then applied to the simulation of a representative scramjet flowpath, and the results compared to experimental data and other kinetic models. This procedure allows a simple reacting model to be used in a predictive manner, and affords significant savings in computational costs for CFD simulations.

Norris, A. T.↗

NASA NTRS2022

After DART: Using the first full-scale test of a kinetic impactor to inform a future planetary defense mission

NASA’s Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ~10% and characterization of the Didymos binary system. But DART is a single experiment; how could these results be used in a future planetary defense necessity involving a different asteroid? We examine what aspects of Dimorphos’s response to kinetic impact will be constrained by DART results; how these constraints will help refine knowledge of the physical properties of asteroidal materials and predictive power of impact simulations; what information about a potential Earth impactor could be acquired before a deflection effort; and how design of a deflection mission should be informed by this understanding. We generalize the momentum enhancement factor β, showing that a particular direction-specific β will be directly determined by the DART results, and that a related direction-specific β is a figure of merit for a kinetic impact mission. The DART β determination constrains the ejecta momentum vector, which, with hydrodynamic simulations, constrains the physical properties of Dimorphos’s near-surface. In a hypothetical planetary defense exigency, extrapolating these constraints to a newly discovered asteroid will require Earth-based observations and benefit from in-situ reconnaissance. We show representative predictions for momentum transfer based on different levels of reconnaissance and discuss strategic targeting to optimize the deflection and reduce the risk of a counterproductive deflection in the wrong direction.

Asteroids↗

NASA NTRS2011

Two-Dimensional Analysis of Conical Pulsed Inductive Plasma Thruster Performance

A model of the maximum achievable exhaust velocity of a conical theta pinch pulsed inductive thruster is presented. A semi-empirical formula relating coil inductance to both axial and radial current sheet location is developed and incorporated into a circuit model coupled to a momentum equation to evaluate the effect of coil geometry on the axial directed kinetic energy of the exhaust. Inductance measurements as a function of the axial and radial displacement of simulated current sheets from four coils of different geometries are t to a two-dimensional expression to allow the calculation of the Lorentz force at any relevant averaged current sheet location. This relation for two-dimensional inductance, along with an estimate of the maximum possible change in gas-dynamic pressure as the current sheet accelerates into downstream propellant, enables the expansion of a one-dimensional circuit model to two dimensions. The results of this two-dimensional model indicate that radial current sheet motion acts to rapidly decouple the current sheet from the driving coil, leading to losses in axial kinetic energy 10-50 times larger than estimations of the maximum available energy in the compressed propellant. The decreased available energy in the compressed propellant as compared to that of other inductive plasma propulsion concepts suggests that a recovery in the directed axial kinetic energy of the exhaust is unlikely, and that radial compression of the current sheet leads to a loss in exhaust velocity for the operating conditions considered here.

Hallock, A. K.↗