An analytical model for the prediction of viscoseal performance in the regime from continuum to free-molecule flow
Predicting seal leakage rates using modified Hodgon model for rarefied gas viscoseal
Patterson, K. E.↗
SEARCH · Search NASA
Results for “CONTINUUM FLOW”
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.
An experiment for Shuttle aerodynamic force coefficient determination from inflight dynamical and atmospheric measurements
A two-phase experiment is proposed which utilizes the Shuttle Orbiter and its unique series of repeated entries into the earth's atmosphere as an airborne in situ aerodynamic testing laboratory. The objective of the experiment is to determine static aerodynamic force coefficients, first of the orbiter, and later of various entry configurations throughout the high speed flight regime, including the transition from free molecule to continuum fluid flow. The objective will be accomplished through analysis of inflight measurements from both shuttle-borne and shuttle-launched instrumented packages. Results are presented to demonstrate the feasibility of such an experiment.
Compton, H. R.↗
Zonally-decoupled DSMC solutions of hypersonic blunt body wake flows
Direct simulation Monte Carlo (DSMC) solutions are presented for the hypersonic flow behind a blunt body in which the wake region is solved in a zonally-decoupled manner. The forebody flow is solved separately using either a DSMC or a Navier-Stokes method, and the forebody exit plane solution is specified as the inflow condition to the decoupled DSMC solution of the wake region. Results are presented for a 70-deg, blunted cone at flow conditions that can be accommodated in existing low-density wind tunnels with the Knudsen number based on base diameter ranging from 0.03 to 0.001. The zonally-decoupled solutions show good agreement with fully-coupled DSMC solutions of the wake flow densities and velocities. The wake closure predicted by the zonally-decoupled solutions is in better agreement with fully-coupled results than that predicted by a fully-coupled Navier-Stokes method indicating the need to account for rarefaction in the wake for the cases considered. The combined use of Navier-Stokes for the forebody with a decoupled DSMC solution for the wake provides an efficient method for solving transitional blunt-body flows where the forebody flow is continuum and the wake is rarefied.
Wilmoth, Richard G.↗
HPTAM, a two-dimensional Heat Pipe Transient Analysis Model, including the startup from a frozen state
A two-dimensional Heat Pipe Transient Analysis Model, 'HPTAM,' was developed to simulate the transient operation of fully-thawed heat pipes and the startup of heat pipes from a frozen state. The model incorporates: (a) sublimation and resolidification of working fluid; (b) melting and freezing of the working fluid in the porous wick; (c) evaporation of thawed working fluid and condensation as a thin liquid film on a frozen substrate; (d) free-molecule, transition, and continuum vapor flow regimes, using the Dusty Gas Model; (e) liquid flow and heat transfer in the porous wick; and (f) thermal and hydrodynamic couplings of phases at their respective interfaces. HPTAM predicts the radius of curvature of the liquid meniscus at the liquid-vapor interface and the radial location of the working fluid level (liquid or solid) in the wick. It also includes the transverse momentum jump condition (capillary relationship of Pascal) at the liquid-vapor interface and geometrically relates the radius of curvature of the liquid meniscus to the volume fraction of vapor in the wick. The present model predicts the capillary limit and partial liquid recess (dryout) in the evaporator wick, and incorporates a liquid pooling submodel, which simulates accumulation of the excess liquid in the vapor core at the condenser end.
Tournier, Jean-Michel↗
Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods
Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.
Chung, Chang-Hong↗
Numerical investigations in the backflow region of a vacuum plume
The objective of this research is to numerically simulate the vacuum plume flow field in the backflow region of a low thrust nozzle exit. In space applications, the low thrust nozzles are used as a propulsion device to control the vehicle attitude, or to maneuver the vehicle flight trajectory. When the spacecraft is deployed in the orbit or cruising in a planetary mission, the vacuum plume is created behind the nozzle exit (so called backflow region), by the exhausting gas of the propulsion system or by venting internal gas to the extremely low density ambient. The low density vacuum plume flow regions cover the continuum, transitional and free molecular flow regimes, which were characterized by the Knudsen number K(sub n), K(sub n) = lambda(sub m)/L where lambda(sub m) is the mean free path of the gas molecules and L is the characteristic length of the flow field. The transitional regime is defined by 0.01 is less than or equal to K(sub n) is less than or equal to 10. The conventional Navier-Stokes equations are valid only in the flow region close to the nozzle exit since the validity of the Navier-Stokes equations fails asymptotically as the Knudsen number increases. The vacuum plume characteristics prediction is primarily a problem of transitional aerodynamics.
Liaw, Goang-Shin↗
DSMC Simulations of Disturbance Torque to ISS During Airlock Depressurization
The primary attitude control system on the International Space Station (ISS) is part of the United States On-orbit Segment (USOS) and uses Control Moment Gyroscopes (CMG). The secondary system is part of the Russian On orbit Segment (RSOS) and uses a combination of gyroscopes and thrusters. Historically, events with significant disturbances such as the airlock depressurizations associated with extra-vehicular activity (EVA) have been performed using the RSOS attitude control system. This avoids excessive propulsive "de-saturations" of the CMGs. However, transfer of attitude control is labor intensive and requires significant propellant. Predictions employing NASA's DSMC Analysis Code (DAC) of the disturbance torque to the ISS for depressurization of the Pirs airlock on the RSOS will be presented [1]. These predictions were performed to assess the feasibility of using USOS control during these events. The ISS Pirs airlock is vented using a device known as a "T-vent" as shown in the inset in figure 1. By orienting two equal streams of gas in opposite directions, this device is intended to have no propulsive effect. However, disturbance force and torque to the ISS do occur due to plume impingement. The disturbance torque resulting from the Pirs depressurization during EVAs is estimated by using a loosely coupled CFD/DSMC technique [2]. CFD is used to simulate the flow field in the nozzle and the near field plume. DSMC is used to simulate the remaining flow field using the CFD results to create an in flow boundary to the DSMC simulation. Due to the highly continuum nature of flow field near the T-vent, two loosely coupled DSMC domains are employed. An 88.2 cubic meter inner domain contains the Pirs airlock and the T-vent. Inner domain results are used to create an in flow boundary for an outer domain containing the remaining portions of the ISS. Several orientations of the ISS solar arrays and radiators have been investigated to find cases that result in minimal disturbance torque. Figure 1 shows surface pressure contours on the ISS and a plane of number density contours for a particular case.
Lumpkin, F. E., III↗
A Fluid Model for the Interaction of the Solar Wind and the Geomagnetic Field
Solar wind and geomagnetic field interaction in terms of continuum theory of fluid flow
SOLAR WIND↗
Numerical Simulation of Rocket Exhaust Interaction with Lunar Soil
This technology development originated from the need to assess the debris threat resulting from soil material erosion induced by landing spacecraft rocket plume impingement on extraterrestrial planetary surfaces. The impact of soil debris was observed to be highly detrimental during NASA s Apollo lunar missions and will pose a threat for any future landings on the Moon, Mars, and other exploration targets. The innovation developed under this program provides a simulation tool that combines modeling of the diverse disciplines of rocket plume impingement gas dynamics, granular soil material liberation, and soil debris particle kinetics into one unified simulation system. The Unified Flow Solver (UFS) developed by CFDRC enabled the efficient, seamless simulation of mixed continuum and rarefied rocket plume flow utilizing a novel direct numerical simulation technique of the Boltzmann gas dynamics equation. The characteristics of the soil granular material response and modeling of the erosion and liberation processes were enabled through novel first principle-based granular mechanics models developed by the University of Florida specifically for the highly irregularly shaped and cohesive lunar regolith material. These tools were integrated into a unique simulation system that accounts for all relevant physics aspects: (1) Modeling of spacecraft rocket plume impingement flow under lunar vacuum environment resulting in a mixed continuum and rarefied flow; (2) Modeling of lunar soil characteristics to capture soil-specific effects of particle size and shape composition, soil layer cohesion and granular flow physics; and (3) Accurate tracking of soil-borne debris particles beginning with aerodynamically driven motion inside the plume to purely ballistic motion in lunar far field conditions. In the earlier project phase of this innovation, the capabilities of the UFS for mixed continuum and rarefied flow situations were validated and demonstrated for lunar lander rocket plume flow impingement under lunar vacuum conditions. Applications and improvements to the granular flow simulation tools contributed by the University of Florida were tested against Earth environment experimental results. Requirements for developing, validating, and demonstrating this solution environment were clearly identified, and an effective second phase execution plan was devised. In this phase, the physics models were refined and fully integrated into a production-oriented simulation tool set. Three-dimensional simulations of Apollo Lunar Excursion Module (LEM) and Altair landers (including full-scale lander geometry) established the practical applicability of the UFS simulation approach and its advanced performance level for large-scale realistic problems.
Liever, Peter↗
Photometric Studies of Heat Flow at the Photosphere
Continuum photometry is carried out and the results of comparing these observations with models of photospheric heat flow are described. The main results are: (1) a possible detection of weak bright rings around some spot penumbrae; (2) no evidence is found for large scale photospheric brightness inhomogeneities exceeding 2-3 K which places tighter constraints on models of global scale convection; (3) supergranular scale continuum structures observed across the photosphere appear mainly due to random clumping of granules; (4) the one case observed of a sunspot emergence shows no thermal shadow exceeding 1.5 K rms one day prior to umbra appearance; (5) network and faculae are found to show a small excess brightness even at mu = 1, so detection of faculae at mu = 1 by differential photometry indicates a gentler temperature gradient near tau = 1 in the facular (relative to cell) atmosphere; (6) the limb darkening study shows no significant global variations to within 0.1% rms.
Foukal, P.↗
Hypersonic rarefied wake characterization
Results of a numerical study using the direct simulation Monte Carlo (DSMC) method are presented for hypersonic rarefied flow over an aeroassisted space transfer vehicle (ASTV). The emphasis of the study is the characterization of the near wake region which includes the ASTV payload. The study covered the transitional flow regime from near continuum to free molecular. Calculations show that the character of the near wake is significantly affected by the presence of the payload. Flow separation occurs when an afterbody is present throughout the transitional flow regime. In contrast, when no afterbody is present, no separation is observed until the flow approaches continuum.
Brewer, E. B.↗
Some Developments of the Equilibrium Particle Simulation Method for the Direct Simulation of Compressible Flows
The direct simulation Monte Carlo (DSMC) method is the established technique for the simulation of rarefied gas flows. In some flows of engineering interest, such as occur for aero-braking spacecraft in the upper atmosphere, DSMC can become prohibitively expensive in CPU time because some regions of the flow, particularly on the windward side of blunt bodies, become collision dominated. As an alternative to using a hybrid DSMC and continuum gas solver (Euler or Navier-Stokes solver) this work is aimed at making the particle simulation method efficient in the high density regions of the flow. A high density, infinite collision rate limit of DSMC, the Equilibrium Particle Simulation method (EPSM) was proposed some 15 years ago. EPSM is developed here for the flow of a gas consisting of many different species of molecules and is shown to be computationally efficient (compared to DSMC) for high collision rate flows. It thus offers great potential as part of a hybrid DSMC/EPSM code which could handle flows in the transition regime between rarefied gas flows and fully continuum flows. As a first step towards this goal a pure EPSM code is described. The next step of combining DSMC and EPSM is not attempted here but should be straightforward. EPSM and DSMC are applied to Taylor-Couette flow with Kn = 0.02 and 0.0133 and S(omega) = 3). Toroidal vortices develop for both methods but some differences are found, as might be expected for the given flow conditions. EPSM appears to be less sensitive to the sequence of random numbers used in the simulation than is DSMC and may also be more dissipative. The question of the origin and the magnitude of the dissipation in EPSM is addressed. It is suggested that this analysis is also relevant to DSMC when the usual accuracy requirements on the cell size and decoupling time step are relaxed in the interests of computational efficiency.
Macrossan, M. N.↗
Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm
In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Herein, several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a “dry diffusion” approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.
97 MATHEMATICS AND COMPUTING↗
Continuum theories for fluid-particle flows: Some aspects of lift forces and turbulence
A general framework is outlined for the modeling of fluid particle flows. The momentum exchange between the constituents embodies both lift and drag forces, constitutive equations for which can be made explicit with reference to known single particle analysis. Relevant results for lift are reviewed, and invariant representations are posed. The fluid and particle velocities and the particle volume fraction are then decomposed into mean and fluctuating parts to characterize turbulent motions, and the equations of motion are averaged. In addition to the Reynolds stresses, further correlations between concentration and velocity fluctuations appear. These can be identified with turbulent transport processes such as eddy diffusion of the particles. When the drag force is dominant, the classical convection dispersion model for turbulent transport of particles is recovered. When other interaction forces enter, particle segregation effects can arise. This is illustrated qualitatively by consideration of turbulent channel flow with lift effects included.
Mctigue, David F.↗
NASA's hypersonic fluid and thermal physics program (Aerothermodynamics)
This survey paper gives an overview of NASA's hypersonic fluid and thermal physics program (recently renamed aerothermodynamics). The purpose is to present the elements of, example results from, and rationale and projection for this program. The program is based on improving the fundamental understanding of aerodynamic and aerothermodynamic flow phenomena over hypersonic vehicles in the continuum, transitional, and rarefied flow regimes. Vehicle design capabilities, computational fluid dynamics, computational chemistry, turbulence modeling, aerothermal loads, orbiter flight data analysis, orbiter experiments, laser photodiagnostics, and facilities are discussed.
Graves, R. A.↗
Study of hypersonic flow past sharp cones
Calculations using the direct simulation Monte Carlo (DSMC) method of Bird for flow past sharp cones in the near continuum to free molecule flow regime are presented and compared with experiment. It is found that results are sensitive to the grid and the interaction potential. Moreover, the time counter method was found to be as accurate as other methods when the solution is grid independent. Finally, the results show that the effects of the wake on the forebody surface properties are minimal.
Taylor, Jeff C.↗
Comparison of continuum and particle simulations of expanding rarefied flows
Comparisons of Navier-Stokes solutions and particle simulations for a simple two-dimensional model problem at a succession of altitudes are performed in order to assess the importance of rarefaction effects on the base flow region. In addition, an attempt is made to include 'Burnett-type' extensions to the Navier-Stokes constitutive relations. The model geometry consists of a simple blunted wedge with a 0.425 meter nose radius, a 70 deg cone half angle, a 1.7 meter base length, and a rounded shoulder. The working gas is monatomic with a molecular weight and viscosity similar to air and was chosen to focus the study on the continuum and particle methodologies rather than the implementation of thermo-chemical modeling. Three cases are investigated, all at Mach 29, with densities corresponding to altitudes of 92 km, 99 km, and 105 km. At the lowest altitude, Navier-Stokes solutions agree well with particle simulations. At the higher altitudes, the Navier-Stokes equations become less accurate. In particular, the Navier-Stokes equations and particle method predict substantially different flow turning angle in the wake near the after body. Attempts to achieve steady continuum solutions including 'Burnett-type' terms failed. Further research is required to determine whether the boundary conditions, the equations themselves, or other unknown causes led to this failure.
Lumpkin, Forrest E., III↗
Nitric oxide formation by meteoroids in the upper atmosphere
The process of nitric oxide formation during atmospheric entry of meteoroids is analyzed theoretically. An ablating meteoroid is assumed to be a point source in a uniform flow with a continuum regime evolving in its wake. The amount of nitric oxide produced by high-temperature reactions of air in the continuum regime is calculated by numerical integration of chemical-rate equations. This is accomplished by assuming that flow properties are constant across the reacting region, its radius being determined from considerations of shock-wave formation and molecular diffusion. The results, when summed over the observed mass, velocity, and entry-angle distributions of meteoroids, provide annual global production rates of nitric oxide as a function of altitude. The peak production is found to occur between 90- and 100-km altitude, the total annual rate being around 40 million kg.
Menees, G. P.↗