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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Magnetotail Structure and its Internal Particle Dynamics During Northward IMF

This study uses Global magnetohydrodynamic (MHD) simulations driven by solar wind data along with Geotail observations of the magnetotail to investigate the magnetotail's response to changes in the interplanetary magnetic field (IMF); observed events used in the study occurred on March 29, 1993 and February 9, 1995. For events from February 9, 1995, we also use the time-dependent MHD magnetic and electric fields and the large-scale kinetic (LSK) technique to examine changes in the Geotail ion velocity distributions. Our MHD simulation shows that on March 29, 1993, during a long period of steady northward IMF, the tail was strongly squeezed and twisted around the Sun-Earth axis in response to variations in the IMF B(sub y) component. The mixed (magnetotail and magnetosheath) plasma observed by Geotail results from the spacecraft's close proximity to the magnetopause and its frequent crossings of this boundary. In our second example (February 9, 1995) the IMF was also steady and northward, and in addition had a significant B(sub y) component. Again the magnetotail was twisted, but not as strongly as on March 29, 1993. The Geotail spacecraft, located approximately 30 R(sub E) downtail, observed highly structured ion distribution functions. Using the time-dependent LSK technique, we investigate the ion sources and acceleration mechanisms affecting the Geotail distribution functions during this interval. At 1325 UT most ions are found to enter the magnetosphere on the dusk side earthward of Geotail with a secondary source on the dawn side in the low latitude boundary layer (LLBL). A small percentage come from the ionosphere. By 1347 UT the majority of the ions come from the dawn side LLBL. The distribution functions measured during the later time interval are much warmer, mainly because particles reaching the spacecraft from the dawn side are affected by nonadiabatic scattering and acceleration in the neutral sheet.

Ashour-Abdalla, M.↗
NASA NTRS2011

Comparison of DSMC and CFD Solutions of Fire II Including Radiative Heating

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. These flows may also contain significant radiative heating. To prepare for these missions, NASA is developing the capability to simulate rarefied, ionized flows and to then calculate the resulting radiative heating to the vehicle's surface. In this study, the DSMC codes DAC and DS2V are used to obtain charge-neutral ionization solutions. NASA s direct simulation Monte Carlo code DAC is currently being updated to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced Quantum-Kinetic chemistry model, and to include electronic energy levels as an additional internal energy mode. The Fire II flight test is used in this study to assess these new capabilities. The 1634 second data point was chosen for comparisons to be made in order to include comparisons to computational fluid dynamics solutions. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid. It is shown that there can be quite a bit of variability in the vibrational temperature inferred from DSMC solutions and that, from how radiative heating is computed, the electronic temperature is much better suited for radiative calculations. To include the radiative portion of heating, the flow-field solutions are post-processed by the non-equilibrium radiation code HARA. Acceptable agreement between CFD and DSMC flow field solutions is demonstrated and the progress of the updates to DAC, along with an appropriate radiative heating solution, are discussed. In addition, future plans to generate more high fidelity radiative heat transfer solutions are discussed.

Liechty, Derek S.↗
NASA NTRS1996

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗
NASA NTRS1986

Comparative energetics diagnosis of the observed and simulated blockings in the Northern Hemisphere

The blocking phenomena, in particular winter blocking, are stuided using energetic diagnostics. Sequences of blockings during the FGGE winter and summer in observed and simulated circulations, and energy transformations associated with blockings are examined. It is observed that the simulation of the Northern Hemisphere circulation correlates better with the observed circulation during summer than during winter. It is noted that the observed winter blockings are maintained by the supply of kinetic energy to the ultralong waves from the shorter cyclone scales through nonlinear wave-wave interaction; in the simulated winter blocking, the energy converted at the cyclone scale cascades down to the shorter wave range; and the summer blockings are supported by the baroclinic conversion of long waves in both the observed and simulated circulations.

Kung, Ernest C.↗
NASA NTRS2017

Electron Drift Resonance in the MHD-Coupled Comprehensive Inner Magnetosphere-Ionosphere Model

Relativistic electrons in the outer radiation belt are highly dynamic and respond to interplanetary solar wind structures interacting with the Earth's magnetic field. A known mechanism dictating electron dynamics is the drift-resonant interaction with ultralow frequency (ULF) waves. The present work simulates the ring current and radiation belt electron populations in the bounce-averaged, kinetic Comprehensive Inner Magnetosphere-Ionosphere model coupled with the Block Adaptive Tree Solar Wind Roe-type Upwind Scheme global magnetospheric magnetohydrodynamic (MHD) code using an idealized ULF wave solar wind density driver. ULF waves generated with 10 min periods (at 1.67 mHz frequencies) in the MHD model are characterized and the corresponding energization of electrons and radial transport of electron phase space density is presented. The drift-resonant electron energy is determined in the simulation and is consistent with the electron resonance conditions in dipolar magnetic fields. The present results will be an important component of understanding inner magnetospheric dynamics and how these inner magnetospheric populations interact with ULF waves resulting from interplanetary solar wind structures.

Komar, C. M.↗
NASA NTRS1993

Role of secondary instability theory and parabolized stability equations in transition modeling

In modeling the laminar-turbulent transition region, the designer depends largely on benchmark data from experiments and/or direct numerical simulations that are usually extremely expensive. An understanding of the evolution of the Reynolds stresses, turbulent kinetic energy, and quantifies in the transport equations like the dissipation and production is essential in the modeling process. The secondary instability theory and the parabolized stability equations method are used to calculate these quantities, which are then compared with corresponding quantities calculated from available direct numerical simulation data for the incompressible boundary-layer flow of laminar-turbulent transition conditions. The potential of the secondary instability theory and the parabolized stability equations approach in predicting these quantities is discussed; results indicate that inexpensive data that are useful for transition modeling in the early stages of the transition region can be provided by these tools.

El-Hady, Nabil M.↗
NASA NTRS2010

Oxidation of SiC Fiber-Reinforced SiC Matrix Composites with a BN Interphase

SiC-fiber reinforced SiC matrix composites with a BN interphase were oxidized in reduced oxygen partial pressures of oxygen to simulate the environment for hypersonic vehicle leading edge applications. The constituent fibers as well as composite coupons were oxidized in oxygen partial pressures ranging from 1000 ppm O2 to 5% O2 balance argon. Exposure temperatures ranged from 816 C to 1353 C (1500 F to 2450 F). The oxidation kinetics of the coated fibers were monitored by thermogravimetric analysis (TGA). An initial rapid transient weight gain was observed followed by parabolic kinetics. Possible mechanisms for the transient oxidation are discussed. One edge of the composite coupon seal coat was ground off to simulate damage to the composite which allowed oxygen ingress to the interior of the composite. Oxidation kinetics of the coupons were characterized by scanning electron microscopy since the weight changes were minimal. It was found that sealing of the coupon edge by silica formation occurred. Differences in the amount and morphology of the sealing silica as a function of time, temperature and oxygen partial pressure are discussed. Implications for use of these materials for hypersonic vehicle leading edge materials are summarized.

Opila, Elizabeth↗
NASA NTRS1980

A numerical study of a marine subtropical stratus cloud layer and its stability

A model for numerical simulation of stratus cloud layers is constructed by combining a second-order closure, turbulent transfer model with a thermal radiative transfer model. The turbulent transfer model allows water vapor saturation. The combined turbulence-radiation model is applied to both a horizontally uniform one-dimensional case and a horizontally nonuniform two-dimensional case. In the latter, the dynamics of mesoscale circulations are also incorporated. Results of the two-dimensional simulation show that the layer cloud instability occurs where the sea surface temperature is high and the large-scale subsidence is weak. The simulated instability is analyzed in view of an instability criterion, the eddy kinetic energy budget, and evaporative cooling near the cloud top.

Moeng, C.-H.↗
NASA NTRS1988

The evolution of cross helicity in driven/dissipative two-dimensional magnetohydrodynamics

The paper presents a series of incompressible two-dimensional simulations of driven/dissipative MHD turbulence where the amount of correlation between the kinetic and magnetic forcing is regulated, thereby controlling the amount of cross helicity injection. It is shown that correlated forcing provides a strong source of magnetofluid cross helicity. The clear appearance of the 'minority species effect,' which is the most striking and systematic effect during this driven dynamic alignment process, is discussed.

Ghosh, S.↗
NASA NTRS1988

On the prediction of the behavior of the velocity, temperature and humidity fields above a thermally stratified air-water interface

A turbulent kinetic energy-dissipation (k-epsilon) model is employed to numerically simulate the behavior of the mean velocity, temperature and humidity fields, along with such second-order statistics as Reynolds stresses, turbulent heat and mass transfer correlation, and bulk transfer coefficients in the boundary layer above wind-generated laboratory waves. The k and epsilon equations have been properly modified to account for thermal stratification, the presence of free surface, and favorable pressure gradient effects. Predictions are compared with available data.

Papadimitrakis, Yiannis A.↗
NASA NTRS1990

Growth morphology with anisotropic surface kinetics

The morphological evolution of crystals growing from an incongruent vapor phase is studied using a Monte Carlo model, and the full range of growth morphologies is recovered. The diffusion in the bulk nutrient and the anisotropy in the interface kinetics are morphologically destabilizing and stabilizing, respectively. For a given set of simulation parameters and lattice symmetries there is a critical size, which scales linearly with the mean free path in the vapor, beyond which a crystal cannot retain its stable, macroscopically faceted growth shape. Surface diffusion stabilizes faceted growth on the shorter scale of the mean surface diffusion length. In simulations with a uniform drift superimposed on the random walk nutrient transport, crystal faces oriented toward the drift show enhanced morphological stability compared to the purely diffusive situation. Rotational drifts with periodic reversal of direction are morphologically stabilizing for all crystal facets.

Xiao, Rong-Fu↗
NASA NTRS1997

Physics of Boundaries and their Interactions in Space Plasmas

This report describes the work done by SciberNet, Inc. during the month of October. We are working on the further refinement of the model used in our large-scale hybrid simulations of the magnetopause. Specifically, we are experimenting with several ways of modeling the effects of cold magnetospheric ions into our simulations. In addition, we are preparing two presentations for the upcoming Fall AGU highlighting the results of these simulations. We have also made progress in our development of a new kinetic linear code which we are using to study the linear properties of the Kelvin-Helmholtz instability at the magnetopause. We have extended the code from the electrostatic limit to the fully electromagnetic regime and are currently in the process of debugging and testing the code. Finally, we have made several test runs with our 2-D hybrid code for the magnetopause. The inflow-outflow boundary conditions are working properly. However, there are issues related to the setup and evolution of the original equilibrium that we are still trying to resolve. Finally, we are preparing several presentations for the upcoming Fall AGU.

Omidi, Nojan↗
NASA NTRS2004

Microgravity Investigation of Crew Reactions in 0-G (MICRO-G)

There is a need for a human factors, technology-based bioastronautics research effort to develop an integrated system that reduces risk and provides scientific knowledge of astronaut-induced loads and motions during long-duration missions on the International Space Station (ISS), which will lead to appropriate countermeasures. The primary objectives of the Microgravity Investigation of Crew Reactions in 0-G (MICRO-GI research effort are to quantify astronaut adaptation and movement as well as to model motor strategies for differing gravity environments. The overall goal of this research program is to improve astronaut performance and efficiency through the use of rigorous quantitative dynamic analysis, simulation and experimentation. The MICRO-G research effort provides a modular, kinetic and kinematic capability for the ISS. The collection and evaluation of kinematics (whole-body motion) and dynamics (reacting forces and torques) of astronauts within the ISS will allow for quantification of human motion and performance in weightlessness, gathering fundamental human factors information for design, scientific investigation in the field of dynamics and motor control, technological assessment of microgravity disturbances, and the design of miniaturized, real-time space systems. The proposed research effort builds on a strong foundation of successful microgravity experiments, namely, the EDLS (Enhanced Dynamics Load Sensors) flown aboard the Russian Mir space station (19961998) and the DLS (Dynamic Load Sensors) flown on Space Shuttle Mission STS-62. In addition, previously funded NASA ground-based research into sensor technology development and development of algorithms to produce three-dimensional (3-0) kinematics from video images have come to fruition and these efforts culminate in the proposed collaborative MICRO-G flight experiment. The required technology and hardware capitalize on previous sensor design, fabrication, and testing and can be flight qualified for a fraction of the cost of an initial spaceflight experiment. Four dynamic load sensors/restraints are envisioned for measurement of astronaut forces and torques. Two standard ISS video cameras record typical astronaut operations and prescribed IVA motions for 3-D kinematics. Forces and kinematics are combined for dynamic analysis of astronaut motion, exploiting the results of the detailed dynamic modeling effort for the quantitative verification of astronaut IVA performance, induced-loads, and adaptive control strategies for crewmember whole-body motion in microgravity. This comprehensive effort, provides an enhanced human factors approach based on physics-based modeling to identify adaptive performance during long-duration spaceflight, which is critically important for astronaut training as well as providing a spaceflight database to drive countermeasure design.

Newman, Dava↗
NASA NTRS1991

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

This research is involved with the implementations of advanced computational schemes based on large eddy simulations (LES) and direct numerical simulations (DNS) to study the phenomenon of mixing and its coupling with chemical reactions in compressible turbulent flows. In the efforts related to LES, a research program was initiated to extend the present capabilities of this method for the treatment of chemically reacting flows, whereas in the DNS efforts, focus was on detailed investigations of the effects of compressibility, heat release, and nonequilibrium kinetics modeling in high speed reacting flows. The efforts to date were primarily focussed on simulations of simple flows, namely, homogeneous compressible flows and temporally developing hign speed mixing layers. A summary of the accomplishments is provided.

Givi, Peyman↗
NASA NTRS2003

Integrated Power and Attitude Control for a Spacecraft with Flywheels and Control Moment Gyroscopes

A law is designed for simultaneous control of the orientation of an Earth-pointing spacecraft, the energy stored by counter-rotating flywheels, and the angular momentum of the flywheels and control moment gyroscopes used together as all integrated set of actuators for attitude control. General. nonlinear equations of motion are presented in vector-dyadic form, and used to obtain approximate expressions which are then linearized in preparation for design of control laws that include feedback of flywheel kinetic energy error as it means of compensating for damping exerted by rotor bearings. Two flywheel 'steering laws' are developed such that torque commanded by all attitude control law is achieved while energy is stored or discharged at the required rate. Using the International Space Station as an example, numerical simulations are performed to demonstrate control about a torque equilibrium attitude and illustrate the benefits of kinetic energy error feedback.

Roithmayr, Carlos M.↗
NASA NTRS2003

Dynamics and Control of Attitude, Power, and Momentum for a Spacecraft Using Flywheels and Control Moment Gyroscopes

Several laws are designed for simultaneous control of the orientation of an Earth-pointing spacecraft, the energy stored by counter-rotating flywheels, and the angular momentum of the flywheels and control moment gyroscopes used together as an integrated set of actuators for attitude control. General, nonlinear equations of motion are presented in vector-dyadic form, and used to obtain approximate expressions which are then linearized in preparation for design of control laws that include feedback of flywheel kinetic energy error as a means of compensating for damping exerted by rotor bearings. Two flywheel steering laws are developed such that torque commanded by an attitude control law is achieved while energy is stored or discharged at the required rate. Using the International Space Station as an example, numerical simulations are performed to demonstrate control about a torque equilibrium attitude, and illustrate the benefits of kinetic energy error feedback. Control laws for attitude hold are also developed, and used to show the amount of propellant that can be saved when flywheels assist the CMGs. Nonlinear control laws for large-angle slew maneuvers perform well, but excessive momentum is required to reorient a vehicle like the International Space Station.

Roithmayr, Carlos M.↗
NASA NTRSDate not supplied

Surrogate Fuel Chemistry from ReaxFF Molecular Dynamics

Reactive molecular dynamics with ReaxFF simulations are used to study the chemistry of four component jet fuel mixtures of n-dodecane, iso-octane, n-propylbenzene, and mesitylene. ReaxFF does not require knowledge of the chemistry a priori providing insight into reactions otherwise not seen in established chemical kinetic mechanisms. The results show that a difference in the local chemical environment can change the fuel decomposition pathways across the surrogate mixtures where the same fuel species are used. These simulations highlight the need to consider cross-chemical interactions in chemical kinetic mechanism development and the importance of resonance stabilized polyaromatic hydrocarbon reactions for soot growth and formation. A method for predicting the threshold sooting index from the carbon species formed at the end of the simulations is established. The results of the simulations are consistent with measurements of the threshold sooting index with the predicted values within 20% of the measured data

Surrogate Fuels↗
NASA NTRS1994

Kinetics of phase transformation in glass forming systems

The objectives of this research were to (1) develop computer models for realistic simulations of nucleation and crystal growth in glasses, which would also have the flexibility to accomodate the different variables related to sample characteristics and experimental conditions, and (2) design and perform nucleation and crystallization experiments using calorimetric measurements, such as differential scanning calorimetry (DSC) and differential thermal analysis (DTA) to verify these models. The variables related to sample characteristics mentioned in (1) above include size of the glass particles, nucleating agents, and the relative concentration of the surface and internal nuclei. A change in any of these variables changes the mode of the transformation (crystallization) kinetics. A variation in experimental conditions includes isothermal and nonisothermal DSC/DTA measurements. This research would lead to develop improved, more realistic methods for analysis of the DSC/DTA peak profiles to determine the kinetic parameters for nucleation and crystal growth as well as to assess the relative merits and demerits of the thermoanalytical models presently used to study the phase transformation in glasses.

Ray, Chandra S.↗