Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Rarefied Flows”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modeling Heat Pipe Startup And Noncondensable Gases In Sockeye

A one-dimensional gas mixture flow model was developed and implemented in the heat pipe code Sockeye to model the effects of non-condensable gases. Additionally a startup model based on the dusty gas model was implemented to model the transition from rarefied gas dynamics to continuum flow, which occurs during the frozen startup of high-temperature heat pipes. Multiple startup models and non-condensable gas models were tested against experimental data for sodium heat pipes, showing excellent agreement. Additionally, the newly developed gas mixture model for modeling non-condensable gas is further tested with a theoretical case study with arbitrary heating configurations. Finally, several recommendations and conclusions are made from the studies in this work to guide future heat pipe modeling efforts.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Direct simulation of stagnation streamline flow for hypersonic reentry

The direct simulation Monte Carlo (DSMC) method is applied to simulate one-dimensional flow along the stagnation streamline. The freestream conditions considered are those encountered by the nose region of the Space Shuttle Orbiter during the hypersonic reentry. The range of altitudes (75 to 92 km) considered in the present calculations covers continuum to the less rarefied portion of the transition flow regime. The calculations account for thermal as well as the chemical nonequilibrium effects. The attention is focused on the flow structure along the stagnation streamline for different specified shock locations at a given altitude. The effect of shock location on the stagnation-point gas composition and heat transfer is analyzed for specific altitude; consequently, these results are appropriate for hemispherical bodies with different nose radii. Finally, comparison is made between the present DSMC results and the viscous shock-layer calculations, which show good agreement at the lowest altitude.

Dogra, Virendra K.↗
NASA NTRS2008

A Modification and Analysis of Lagrangian Trajectory Modeling and Granular Dynamics of Lunar Dust Particles

A previously developed mathematical model is amended to more accurately incorporate the effects of lift and drag on single dust particles in order to predict their behavior in the wake of high velocity gas flow. The model utilizes output from a CFD or DSMC simulation of exhaust from a rocket nozzle hot gas jet. An extension of the Saffman equation for lift based on the research of McLaughlin (1991) and Mei (1992) is used, while an equation for the Magnus force modeled after the work of Oesterle (1994) and Tsuji et al (1985) is applied. A relationship for drag utilizing a particle shape factor (phi = 0.8) is taken from the work of Haider and Levenspiel (1989) for application to non-spherical particle dynamics. The drag equation is further adjusted to account for rarefaction and compressibility effects in rarefied and high Mach number flows according to the work of Davies (1945) and Loth (2007) respectively. Simulations using a more accurate model with the correction factor (Epsilon = 0.8 in a 20% particle concentration gas flow) given by Richardson and Zaki (1954) and Rowe (1961) show that particles have lower ejection angles than those that were previously calculated. This is more prevalent in smaller particles, which are shown through velocity and trajectory comparison to be more influenced by the flow of the surrounding gas. It is shown that particles are more affected by minor changes to drag forces than larger adjustments to lift forces, demanding a closer analysis of the shape and behavior of lunar dust particles and the composition of the surrounding gas flow.

Long, Jason M.↗
NASA NTRSDate not supplied

Applications of PLIF within Supersonic/Hypersonic Flows relevant to Space Technology and Exploration

We summarize two recent applications of planar laser-induced fluorescence using seeded nitric oxide (NO-PLIF) within two different large-scale facilities for space technology and exploration. The first measurement campaign described is the NO-PLIF measurements performed within the 31-in Mach 10 air tunnel at NASA Langley Research Center for flow visualization, velocity, and temperature measurements of the wake produced by a hypersonic blunt-body. PLIF within this wake-flow was performed to support the developments of HIAD (hypersonic inflatable aerodynamic decelerator) technology, including the recent LOFTID flight mission. The second campaign is the NO-PLIF measurements performed within a 20-ft vacuum chamber at the NASA Marshall Space Flight Center for flow visualization of a Mach 5 plume impinging on a flat plate under rarefied and continuum conditions. PLIF within this supersonic flow was performed to support future lunar and Martian lander missions by investigating the plume impingement behavior. While both of these measurement campaigns were performed using 10 Hz PLIF laser systems, current efforts include utilization of tunable high-repetition-rate lasers for NO-PLIF measurements up to kHz-MHz rates.

Planar laser-induced fluorescence↗
NASA NTRS1993

Simulated rarefied entry of the Galileo probe into the atmosphere of Jupiter

Flow properties and aerodynamics are computed with a direct simulation Monte Carlo (DSMC) method for rarefied entry of the Galileo Probe into the atmosphere of Jupiter. Accurate predictions of vehicle drag coefficients are needed in order to assess atmospheric properties from the onboard Atmosphere Structure Experiment where highly-sensitive accelerometers will measure the drag force to within 10-6 barr during the initial entry phase at high altitudes. The corresponding flow rarefraction extends from the free molecule limit to the near continuum transition regime (Re less than 1000). Simulation results indicate that C(sub D) varies from 2.1 at the free molecule limit down to 1.6 at Re(infinity) = 1,000. Temperatures, densities, and internal energies throughout the flow field were also computed at each altitude ranging from 735 km to 353 km above the 1 barr level in the Jovian atmosphere. Surface heating and temperatures of the probe were computed directly in the DSMC code by assuming radiative equilibrium. Material response was re-asssessed accurately during entry by accounting for conductivity, heat capacity, and pyrolysis which led to surface material mass efflux several times that of the freestream mass influx. The simulation also accounted for the quantum nature of the rotational energy mode of the dominant atmospheric species H2 through partial internal excitation in the freestream gas.

Hass, Brian L.↗
NASA NTRS2002

Simulations of Ground and Space-Based Oxygen Atom Experiments

Fast, pulsed atomic oxygen sources are a key tool in ground-based investigations of spacecraft contamination and surface erosion effects. These technically challenging ground-based studies provide a before and after picture of materials under low-earth-orbit (LEO) conditions. It would be of great interest to track in real time the pulsed flux from the source to the surface sample target and beyond in order to characterize the population of atoms and molecules that actually impact the surface and those that make it downstream to any coincident detectors. We have performed simulations in order to provide such detailed descriptions of these ground-based measurements and to provide an assessment of their correspondence to the actual LEO environment. Where possible we also make comparisons to measured fluxes and erosion yields. To perform the calculations we use a detailed description of a measurement beam and surface geometry based on the W, pulsed apparatus at Montana State University. In this system, a short pulse (on the order of 10 microseconds) of an O/O2 beam impacts a flat sample about 40 cm downstream and slightly displaced &om the beam s central axis. Past this target, at the end of the beam axis is a quadrupole mass spectrometer that measures the relative in situ flux of 0102 to give an overall normalized erosion yield. In our simulations we use the Direct Simulation Monte Carlo (DSMC) method, and track individual atoms within the atomic oxygen pulse. DSMC techniques are typically used to model rarefied (few collision) gas-flows which occur at altitudes above approximately 110 kilometers. These techniques are well suited for the conditions here, and multi-collision effects that can only be treated by this or a similar technique are included. This simulation includes collisions with the surface and among gas atoms that have scattered from the surface. The simulation also includes descriptions of the velocity spread and spatial profiles of the O/O2 beam obtained from separate measurements. These computations use basic engineering models for the gas-gas and gas-surface scattering and focus on the influence of multi-collision effects. These simulations characterize many important quantities of interest including the actual flux of atoms that reach the surface, the energy distribution of this flux, as well as the direction of the velocity of the flux that strikes the surface. These quantities are important in characterizing the conditions which give rise to measured surface erosion. The calculations also yield time- snapshots of the pulse as it impacts and flows around the surface. These snapshots reveal the local environment of gas near the surface for the duration of the pulse. We are also able to compute the flux of molecules that travel downstream and reach the spectrometer, and we characterize their velocity distribution. The number of atoms that reach the spectrometer can in fact be influenced by the presence of the surface due to gas-gas collisions from atoms scattered h m the surface, and it will generally be less than that with the surface absent. This amounts to an overall normalization factor in computing erosion yields. We discuss these quantities and their relationship to the gas-surf$ce interaction parameters. We have also performed similar calculations corresponding to conditions (number densities, temperatures, and velocities) of low-earth orbit. The steady-state nature and lower overall flux of the actual space environment give rise to differences in the nature of the gas-impacts on the surface from those of the ground-based measurements using a pulsed source.

Minton, T. K.↗
NASA NTRS1991

Hypersonic merged layer blunt body flows with wakes

An attempt is made here to understand the basic physics of the flowfield with wake on a blunt body of revolution under hypersonic rarefied conditions. A merged layer model of flow is envisioned. Full steady-state Navier-Stokes equations in spherical polar coordinate system are computed from the surface with slip and temperature jump conditions to the free stream by the Accelerated Successive Replacement method of numerical integration. Analysis is developed for bodies of arbitrary shape, but actual computations have been carried out for a sphere and sphere-cone body. Particular attention is paid to set the limit of the onset of separation, wake closure, shear-layer impingement, formation and dissipation of the shocks in the flowfield. Validity of the results is established by comparing the present results for sphere with the corresponding results of the SOFIA code in the common region of their validity and with the experimental data.

Jain, Amolak C.↗
NASA NTRS1989

Rarefied-continuum gas dynamics transition for SUMS project

This program is to develop an analytic method for reducing SUMS data for the determination of the undisturbed atmosphere conditions ahead of the shuttle along its descending trajectory. It is divided into an internal flow problem, an external flow problem and their matching conditions. Since the existing method of Direct Simulation Monte Carlo (DSMC) failed completely for the internal flow problem, the emphasis is on the internal flow of a highly non-equilibrium, rarefied air through a short tube of a diameter much less than the gaseous mean free path. A two fluid model analysis of this internal flow problem has been developed and studied with typical results illustrated. A computer program for such an analysis and a technical paper published in Lecture Notes in Physics No. 323 (1989) are included as Appendices 3 and 4. A proposal for in situ determination of the surface accommodation coefficients sigma sub t and sigma e is included in Appendix 5 because of their importance in quantitative data reduction. A two fluid formulation for the external flow problem is included as Appendix 6 and a review article for AIAA on Hypersonic propulsion, much dependent on ambient atmospheric density, is also included as Appendix 7.

Cheng, Sin-I↗
NASA NTRS1991

Thermomolecular effect on pressure measurements with orifices in transitional flow

Results are presented of a study carried out in order to verify a version of the method (presented at the 7th International Rarefied Gas Dynamics Symposium) for calculating the influence of thermomolecular flow on the pressures measured by means of orifices in walls exposed to rarefied gases. In these experiments, the 'orifice effect' is examined under conditions that are significantly different from those for which the semiempirical (SE) method of Kinslow and Potter (1971) was shown to be successful. Comparisons between the predicted results with flight data and the results obtained in the laboratory, combined with limited DSMC calculations, indicate that the SE method remains a useful tool which is relatively easy to apply and which gives good results when an appropriate energy accommodation coefficient is used.

Potter, J. L.↗
DOE OSTI2022

Modeling rarefied gas chemistry with QuiPS, a novel quasi-particle method

The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Overall, variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.

42 ENGINEERING↗
NASA NTRSDate not supplied

Towards Plume Impingement Modeling in Space Environments

After 30 years human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to act as a home base for Lunar exploration. A large space station must be assembled while in orbit, using a “piecemeal” approach. In this context, it means that the different modules will arrive at different times and attach to what is already in service. The ISS provides an excellent example of this approach, and the proposed Lunar Gateway will undergo a similar assembly process. This assembly is achieved via “docking” maneuvers between modules, which are made possible by sequential firings of the onboard reaction control system (RCS) thrusters. They work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking approach. The issue is that the gas from these thrusters' forms flow structures described as “plumes” and can impinge onto the outer surfaces of the space station, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion. All mechanisms that can damage the space station and must be avoided. Both permanent and visiting modules will have these RCS thruster exhaust impingement problems. Accurately and efficiently modeling these plume is an involved multi-physics calculation but also an important tool when designing the control algorithms of approaching modules. This poster presents progress towards this simulation on two fronts. First is the verification of OpenFOAM for rarefied plume impingement calculations by direct comparison to published DAC cases. Second is the estimation of plume impingement strikes over the time scale of an entire docking event. This is done by using a simple plume source flow model and a prescribed motion visualizer—coded in Python. Together these tools help push NASA's capabilities for simulating these plume impingement effects.

Rarefied Flows↗
NASA NTRS1966

Omegatron studies of a skimmed beam system.

Continuum source molecular beam facility, examining effect of source-skimmer distance variation on flux and spatial distribution of beam in collimation chamber downstream

MOLECULAR BEAM↗
NASA NTRS1990

Bipropellant engine plume analyses for the Mars Observer scientific spacecraft

This paper discusses the methodology utilized to characterize exhaust plume expansion of the 490 N bipropellant engine and provides an analysis of exhaust plune impingement effects on the Mars Observer oxidizer tank structure. Nozzle and plume continuum flowfields are computed with a method of characteristics solution while transitional and rarefied regions of the exhaust gas flow were defined by a direct simulation Monte Carlo method. Results of this study provide sufficient data for an assessment of spacecraft thermal impacts due to plume heating effects and thrust losses from exhaust plume impingement.

Rattenni, L., Jr.↗
NASA NTRS1991

An engineering assessment of gas/surface interactions in free-molecular aerodynamics

Data available on the gas/surface interaction (GSI) phenomenon, particularly in the rarefied to near-free-molecular flow regime are collected, categorized, and analyzed with a purpose of developing a GSI model that could be used as a guide to spacecraft designers. The study shows that there are not enough useful data for building a high-confidence GSI model. However, sufficient results are obtained to suggest that continued reliance on the diffuse GSI model is inappropriate, particularly in the near-free-molecular regime.

Knox, E. C.↗
NASA NTRSDate not supplied

ANALYSIS OF THE MSL/MEDLI ENTRY DATA WITH COUPLED CFD AND MATERIAL RESPONSE.

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heat-shield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material. PICA is a lightweight carbon fiber/polymeric resin material that offers out-standing performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP). The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on Open-FOAM (PATO) using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA, is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Mars Science Laboratory↗
NASA NTRSDate not supplied

Analysis of MSL/MEDLI Entry Data with Coupled CFD and Material Response

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on OpenFOAM (PATO) [4,5,6] using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA [7], is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Thermal Protection Systems↗