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 577 records · Page 32

Validation of a Two-Phase CFD Model for Autogenous Pressurization and Expulsion

This paper presents a two-phase computational fluid dynamics (CFD) model for simulating autogenous pressurization and expulsion in a cryogenic liquid hydrogen storage tank, utilizing a volume-of-fluid (VOF) approach combined with the kinetics-based Schrage equation to accurately capture the phase boundary and interfacial heat, mass, and momentum transfer between vapor and liquid phases. The model is validated against experimental data from NASA's K-site tank facility, specifically experimental case 225, which involved autogenous pressurization and controlled expulsion of liquid hydrogen. Various turbulence models are evaluated to assess their influence on model accuracy. The CFD simulations successfully replicate key thermodynamic behaviors observed during the experiments, including pressure evolution, temperature profiles, and phase-change dynamics at the vapor-liquid interface. The predicted tank pressures and temperatures agree well with experimental data, whereas the pressurant mass predictions are within 16% of observed values. This study emphasizes the importance of selecting appropriate turbulence models to accurately simulate the complex flow and heat transfer phenomena during tank autogenous pressurization and expulsion. By improving the accuracy and reliability of CFD models for these processes, this research contributes to developing efficient cryogenic propellant management strategies for future space missions.

Computational Fluid Dynamics↗

Lewis number and Damkohler number effects in vort ex-flame interactions

A combined experimental-numerical study of the interaction of a two-dimensional vortex pair with a plane premixed laminar flame has been carried out. The selected geometry is two-dimensional in order to allow for the use of quantitative line-of-sight measurement techniques and for comparisons with two-dimensional direct numerical simulations. The experiment is used to identify possible effects of the Lewis number and of radiative heat losses, and direct numerical simulations reproducing the experimental conditions are used to investigate the role of the Lewis number, of heat losses, and of multi-step kinetics. The heat losses considered in this study are essentially radiative losses from the burnt gases and do not include conductive losses to the walls. This report presents the final results of this study of vortex-flame interactions. The main result is the creation of experimental data sets of vortex-flame interactions, which demonstrate the importance of the Lewis and Damkohler numbers to the flame and which allow for comparison with direct numerical simulations. In turn, the comparison between experimental and numerical results shows that the Lewis number effects and complex chemistry play significant roles, which can be reproduced using a two-step reaction mechanism.

Samaniego, Jean-Michael↗

A comparison of observed and simulated properties of sudden stratospheric warmings

Review of observational data and dynamical numerical simulations of stratospheric warmings. Classes of warmings, major and minor (major if poleward movement of planetary-scale thermal systems entails reversal of polar circulation at 10 mb or below), trajectories of warm cells, vertical and horizontal scale of warm-air systems, the time-scale of warming, initial zonal flow conditions prior to a warming, circulation reversals, and details of the energy budget before and after a warming are discussed. The 1963 and 1973 types of warmings are contrasted: the strong baroclinic conversion of eddy potential to eddy kinetic energy was not repeated in the latter, but both events were preceded by very large fluxes from the troposphere. Numerical model simulations by various authors are compared and evaluated.

Quiroz, R. S.↗

Project ARGO: Gas phase formation in simulated microgravity

The ARGO study investigated the reduced incidence of joint pain decompression sickness (DCS) encountered in microgravity as compared with an expected incidence of joint pain DCS experienced by test subjects in Earth-based laboratories (unit gravity) with similar protocols. Individuals who are decompressed from saturated conditions usually acquire joint pain DCS in the lower extremities. Our hypothesis is that the incidence of joint pain DCS can be limited by a significant reduction in the tissue gas micronuclei formed by stress-assisted nucleation. Reductions in dynamic and kinetic stresses in vivo are linked to hypokinetic and adynamic conditions of individuals in zero g. We employed the Doppler ultrasound bubble detection technique in simulated microgravity studies to determine quantitatively the degree of gas phase formation in the upper and lower extremities of test subjects during decompression. We found no evidence of right-to-left shunting through pulmonary vasculature. The volume of gas bubble following decompression was examined and compared with the number following saline contrast injection. From this, we predict a reduced incidence of DCS on orbit, although the incidence of predicted mild DCS still remains larger than that encountered on orbit.

Powell, Michael R.↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗

Topological Evolution of a Fast Magnetic Breakout CME in 3-Dimensions

W present the extension of the magnetic breakout model for CME initiation to a fully 3-dimensional, spherical geometry. Given the increased complexity of the dynamic magnetic field interactions in 3-dimensions, we first present a summary of the well known axisymmetric breakout scenario in terms of the topological evolution associated with the various phases of the eruptive process. In this context, we discuss the completely analogous topological evolution during the magnetic breakout CME initiation process in the simplest 3-dimensional multipolar system. We show that an extended bipolar active region embedded in an oppositely directed background dipole field has all the necessary topological features required for magnetic breakout, i.e. a fan separatrix surface between the two distinct flux systems, a pair of spine fieldlines, and a true 3-dimensional coronal null point at their intersection. We then present the results of a numerical MHD simulation of this 3-dimensional system where boundary shearing flows introduce free magnetic energy, eventually leading to a fast magnetic breakout CME. The eruptive flare reconnection facilitates the rapid conversion of this stored free magnetic energy into kinetic energy and the associated acceleration causes the erupting field and plasma structure to reach an asymptotic eruption velocity of greater than or approx. equal to 1100 km/s over an approx.15 minute time period. The simulation results are discussed using the topological insight developed to interpret the various phases of the eruption and the complex, dynamic, and interacting magnetic field structures.

Lynch, B. J.↗

Signatures of current loop coalescence in solar flares

The nonlinear coalescence instability of current carrying solar loops can explain many of the characteristics of the solar flares such as their impulsive nature, heating and high energy particle acceleration, amplitude oscillations of electromagnetic emission as well as the characteristics of 2-D microwave images obtained during a solar flare. The physical characteristics of the explosive coalescence of currents are presented in detail through computer simulation and theory. Canonical characteristics of the explosive coalescence are: (1) a large amount of impulsive increase of kinetic energies of electrons and ions; (2) simultaneous heating and acceleration of electrons and ions in high and low energy spectra; (3) ensuing quasi-periodic amplitude oscillations in fields and particle quantities; and (4) the double peak (or triple peak) structure in these profiles, participate in the coalescence process, yielding varieties of phenomena.

Sakai, J.↗

The performance of multilayer insulation in a rapidly depressurizing environment

The performance of multilayer insulation in a rapidly depressurizing environment is determined by the variation of heat transfer with internal pressure and the pressure history of the interstitial gas. Measurements of thermal performance were made on three multilayer insulation configurations at room temperature and steady state pressures from 10(exp -6) to 10(exp -1) torr. The heat transfer due to gas conduction alone correlated with the kinetic theory of gases in the molecular flow regime. Pressure histories were measured in a unique apparatus which simulated the depressurization rate of a boost vehicle. The pressure histories on both sides of two of the configurations were measured to bound the actual interstitial pressure. The results of the two types of measurements agreed with earlier work and were combined to make performance predictions using an actual ascent pressure history.

Efromson, Ronald A.↗

Mirror and ion cyclotron anisotropy instabilities in the magnetosheath

The kinetic properties of the ion cyclotron anisotropy and the mirror instabilities are investigated using one-dimensional hybrid simulations. It is found that, for moderate values of the ion beta and the proton temperature anisotropy, the two instabilities produce similar levels of turbulence. For high values of beta or temperature anisotropy, the ion cyclotron instability produces higher fluctuation levels of turbulence than does the mirror instability.

Mckean, M. E.↗

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

The AEDC aerospace chamber 7V: An advanced test capability for infrared surveillance and seeker sensors

An advanced sensor test capability is now operational at the Air Force Arnold Engineering Development Center (AEDC) for calibration and performance characterization of infrared sensors. This facility, known as the 7V, is part of a broad range of test capabilities under development at AEDC to provide complete ground test support to the sensor community for large-aperture surveillance sensors and kinetic kill interceptors. The 7V is a state-of-the-art cryo/vacuum facility providing calibration and mission simulation against space backgrounds. Key features of the facility include high-fidelity scene simulation with precision track accuracy and in-situ target monitoring, diffraction limited optical system, NIST traceable broadband and spectral radiometric calibration, outstanding jitter control, environmental systems for 20 K, high-vacuum, low-background simulation, and an advanced data acquisition system.

Simpson, W. R.↗

Kinetic physics of the mirror instability

Kinetic mechanisms for the growth and saturation of the mirror instability are described using one-dimensional hybrid simulations. Two parameter regimes are considered. In the first regime, a relatively small ion anisotropy excites a slowly growing instability that produces small-amplitude waves; most ions respond to the waves as an adiabatic fluid. In the second regime a large anisotropy excites a rapidly growing instability that generates large-amplitude waves; the response of many ions in this case is nonadiabatic. The difference in ion response is due to the relative importance of two ion populations, resonant and nonresonant. Resonant ions, those ions with low velocities parallel to the background magnetic field, contribute to the growth of the instability as a result of their gyrointeractions with the noncoplanar component of the waves electric field and respond to the mirror waves nonadiabatically. Nonreasonant ions, those with large parallel velocities, respond as an adiabatic fluid. In both regimes, ion anisotropy is reduced by means of the magnetic mirror force; in the second regime, the anisotropy reduction reduces the free energy available for wave growth and leads to the saturation of the mirror instability.

Mckean, M. E.↗

Modeling the decay of energy containing eddies: A source of solar wind heating

To understand the solar wind heating and acceleration mechanisms one needs to understand the decay of energy containing eddies. With this goal in mind, attempts have been made to extend the fluid dynamic phenomenology of large scale quasi-equilibrium to the case of magnetohydrodynamics. Matthaeus et al. have proposed a model for the inhomogeneous transport and decay of five mean variables, namely, two mean square Elsasser variables z(exp 2) (sub +/-) their correlation lengths, and the difference between the kinetic and magnetic energies. We test the validity of this model in the simplified case of homogeneous turbulence simulated in a periodic box. We propose a class of models and show that they may fit the simulation satisfactorily. Analytic solutions of this class of model reveal their inherent properties and demonstrate the difficulties associated with finite cross helicity. It is noted that adjustments are required to make the simplest models, which are based upon isotropic turbulence, scale properly with respect to the strength of the mean magnetic field. This can be interpreted as due to anisotropic turbulence, which can be modelled by simple parameterization in the phenomenology.

Hossain, M.↗

Assessment of the Subgrid-Scale Models at Low and High Reynolds Numbers

Accurate SGS models must be capable of correctly representing the energy transfer between GS and SGS. Recent direct assessment of the energy transfer carried out using direct numerical simulation (DNS) data for wall-bounded flows revealed that the energy exchange is not unidirectional. Although GS kinetic energy is transferred to the SGS (forward scatter (F-scatter) on average, SGS energy is also transferred to the GS. The latter energy exchange (backward scatter (B-scatter) is very significant, i.e., the local energy exchange can be backward nearly as often as forward and the local rate of B-scatter is considerably higher than the net rate of energy dissipation.

Horiuti, K.↗

Turbulent Premixed Flame Propagation in Microgravity

A facility in which turbulent Couette flow could be generated in a microgravity environment was designed and built. To fit into the NASA Lewis drop tower the device had to be very compact. This means that edge effects and flow re-circulation were expected to affect the flow. The flow was thoroughly investigated using LDV and was found to be largely two dimensional away from the edges with constant turbulence intensities in the core. Slight flow asymmetries are introduced by the non symmetric re-circulation of the fluid outside the test region. Belt flutter problems were remedied by adding a pair of guide plates to the belt. In general, the flow field was found to be quite similar to previously investigated Couette flows. However, turbulence levels and associated shear stresses were higher. This is probably due to the confined re-circulation zone reintroducing turbulence into the test section. An estimate of the length scales in the flow showed that the measurements were able to resolve nearly all the length scales of interest. Using a new LES method for subgrid combustion it has been demonstrated that the new procedure is computational feasible even on workstation type environment. It is found that this model is capable of capturing the propagation of the premixed names by resolving the flame in the LES grid within 2-3 grid points. In contrast, conventional LES results in numerical smearing of the flame and completely inaccurate estimate of the turbulent propagation speed. Preliminary study suggests that there is observable effect of buoyancy in the 1g environment suggesting the need for microgravity experiments of the upcoming experimental combustion studies. With the cold flow properties characterized, an identical hot flow facility is under construction. It is assumed that the turbulence properties ahead of the flame in this new device will closely match the results obtained here. This is required since the hot facility will not enable LDV measurements. The reacting flow facility is also being constructed with planned drop tower experiments in mind. Therefore, issues related to safety and structural integrity are being take into account. Further development of the numerical model will also be carried out to include finite-rate kinetics for representative premixed cases. More detail analysis of the flame structure and propagation nature will be investigated. Simulations will also be compared to the flame properties observed in the experiments.

Menon, S.↗

Self-Pressurization and Spray Cooling Simulations of the Multipurpose Hydrogen Test Bed (MHTB) Ground-Based Experiment

This paper presents a CFD (computational fluid dynamics) model for simulating the self-pressurization of a large scale liquid hydrogen storage tank. In this model, the kinetics-based Schrage equation is used to account for the evaporative and condensing interfacial mass flows. Laminar and turbulent approaches to modeling natural convection in the tank and heat and mass transfer at the interface are compared. The flow, temperature, and interfacial mass fluxes predicted by these two approaches during tank self-pressurization are compared against each other. The ullage pressure and vapor temperature evolutions are also compared against experimental data obtained from the MHTB (Multipuprpose Hydrogen Test Bed) self-pressurization experiment. A CFD model for cooling cryogenic storage tanks by spraying cold liquid in the ullage is also presented. The Euler- Lagrange approach is utilized for tracking the spray droplets and for modeling interaction between the droplets and the continuous phase (ullage). The spray model is coupled with the VOF (volume of fluid) model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. Droplet ullage heat and mass transfer are modeled. The flow, temperature, and interfacial mass flux predicted by the model are presented. The ullage pressure is compared with experimental data obtained from the MHTB spray bar mixing experiment. The results of the models with only droplet/ullage heat transfer and with heat and mass transfer between the droplets and ullage are compared.

cryogenic fluid storage↗

Atmosphere Modeling and Performance Sensitivity for the Mars Sample Return Earth Entry System

The Capture, Containment, and Return System (CCRS) mission is a key element of the joint NASA-European Space Agency (ESA) planned Mars Sample Return (MSR) Campaign. The CCRS assembled Earth Entry System (EES) will enter the mission’s final segment in its Approach, Entry, Descent, and Landing (AEDL) Phase. The EES AEDL aims to deliver a highly reliable, safe, and accurate return while maintaining strict containment assurance targets established by the campaign. As currently designed, the EES would be the first fully passive sample return capsule with no parachute or onboard control system, prompting a significant effort in Earth atmosphere characterization and modeling. Earth’s atmosphere, specifically winds, have a strong influence on EES flight mechanics and landing footprint during its free fall landing. This paper describes Earth atmosphere modeling, atmosphere characterization, and performance sensitivities incorporated into the teams’s efforts to ensure AEDL success. By utilizing high resolution balloon radiosondes, analyzing wind structural and distributional compositions, and investigating flight mechanics sensitivities, the AEDL atmosphere team has been able to better understand and simulate local wind conditions at the Utah Test and Training Range (UTTR). Utilizing tools such as horizontal turbulent kinetic energy, vertical wind shear, or integrated wind, the team have been able to reveal valuable information about wind profiles in a deeper context than previously conducted, ultimately improving understanding of AEDL flight mechanics sensitivity and EES design.

Kaustubh Ray↗