Search NASA⌕ Search

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.

At least 55 records · Page 3

Reynolds number influences in aeronautics

Reynolds number, a measure of the ratio of inertia to viscous forces, is a fundamental similarity parameter for fluid flows and therefore, would be expected to have a major influence in aerodynamics and aeronautics. Reynolds number influences are generally large, but monatomic, for attached laminar (continuum) flow; however, laminar flows are easily separated, inducing even stronger, non-monatomic, Reynolds number sensitivities. Probably the strongest Reynolds number influences occur in connection with transitional flow behavior. Transition can take place over a tremendous Reynolds number range, from the order of 20 x 10(exp 3) for 2-D free shear layers up to the order of 100 x 10(exp 6) for hypersonic boundary layers. This variability in transition behavior is especially important for complex configurations where various vehicle and flow field elements can undergo transition at various Reynolds numbers, causing often surprising changes in aerodynamics characteristics over wide ranges in Reynolds number. This is further compounded by the vast parameterization associated with transition, in that any parameter which influences mean viscous flow development (e.g., pressure gradient, flow curvature, wall temperature, Mach number, sweep, roughness, flow chemistry, shock interactions, etc.), and incident disturbance fields (acoustics, vorticity, particulates, temperature spottiness, even electro static discharges) can alter transition locations to first order. The usual method of dealing with the transition problem is to trip the flow in the generally lower Reynolds number wind tunnel to simulate the flight turbulent behavior. However, this is not wholly satisfactory as it results in incorrectly scaled viscous region thicknesses and cannot be utilized at all for applications such as turbine blades and helicopter rotors, nacelles, leading edge and nose regions, and High Altitude Long Endurance and hypersonic airbreathers where the transitional flow is an innately critical portion of the problem.

Bushnell, Dennis M.↗

Heat Transfer from Cylinders in Transition From Slip Flow to Free-Molecule Flow

Over 600 measured heat-transfer coefficients in the transition from slip to free-molecule flow have been correlated by using the Nusselt number Nu as a function of the Knudsen Kn and Reynolds Re (or Mach M) numbers. The experimental range for these heat-transfer data from transverse cylinders in air corresponds to the following dimensionless groups: M, 0.10 to 0.90; Re, 0.03 to 11.5; Kn, 0.10 to 5.0. The total air temperature T(sub t) was maintained constant at 80 F, but wire temperature was Varied from 150 to 580 F. At Kn=0.10, Nu extrapolates smoothly into slip-flow empirical curves that show Nu as a function of Re and M or Kn. The correlation gradually changes from the square root of Re(sub t) dependence characteristic of continuum flow to first-power Re dependence as Kn increases (decreasing Re). At the experimental limit Kn ft 5.0, the Nu data correlate with a mean fractional error of 413 percent by the prediction of free-molecule-flow theory. In comparing experimental results with theory, an accommodation coefficient of 0.57+/-0.07 was inferred from the heat-transfer data, which were obtained with etched tungsten wire in air. The wire recovery temperature T(sub e) was measured and compared with existing data and theory in terms of a ratio eta(equivalent to T(sub e)/T(sub t). The results can be divided into three groups by Kn criteria: For Kn less than 2.01, eta is independent of Kn, and eta decreases from 1.0 to 0.97 as M increases from 0 to 0.90; for 2.0 less than Kn less than 5.0, eta is a function of both Kn and M in this transition region to fully developed free-molecule flow; and for Kn greater than 5.0, eta predicted by free-molecule-flow theory is observed and increases from 1.0 to 1.08 as M increases from 0 to 0.90, again independent of Kn. Therefore, these T(sub e) data provide a guide to the boundary of fully developed free-molecule flow, which is.inferred from this research to exist for Kn greater than 5.0. This boundary criterion is substantiated by other published data on T(sub e) at supersonic speeds.

Cybulski, Ronald J.↗

Electrical conductivity and velocity of highly ionized plasma flows - Theory and experiment.

Use of an immersible, three-coil, magnetic-induction probe, previously tested in a low-density supersonic argon jet, to measure electrical conductivity and velocity profiles of a highly ionized high-density nitrogen jet in the continuum flow regime where effects due to probe bow shocks and boundary layers might not be negligible. Measured centerline values of electrical conductivity and velocity were compared with predictions based on a theoretical analysis previously developed to study the gas as it expanded adiabatically and inviscidly from an equilibrium sonic state to the nozzle exit. The resulting numerical exit plane values for electron density and electron temperature were then substituted into the Spitzer-Haerm conductivity formula to compute a theoretical conductivity which agreed within 40% of the measured conductivity, while the calculated and experimental velocity values differed by as much as 50%. The lack of agreement was attributed to the possible use of invalid assumptions and boundary conditions in the computer analysis or to the unknown effects of shocks on the probe data.

Vendell, E. W.↗

A Preliminary Theoretical Study of the Expansion Tube, a New Device for Producing High-Enthalpy Short-Duration Hypersonic Gas Flows

A preliminary theoretical study was conducted of a new technique for producing high-enthalpy gas flows. The device considered utilizes an unsteady expansion process for the purpose of total-enthalpy multiplication. Analyses were conducted for both perfect and real air in equilibrium, assuming idealized diaphragm bursts, centered expansion waves, and continuum flow. Results of the study showed the expansion tube capable of outperforming the conventional shock tunnel by a factor of approximately 2 in velocity for the same test-section ambient density and pressure. The degree of dissociation is low at all phases of the thermodynamic cycle in the expansion tube; thus the test-section air has a good possibility of being in equilibrium. Maximum pressures involved in the cycle for duplicating a typical reentry from, a lunar mission are low enough for existing pressure-vessel techniques. Both the known and anticipated advantages and disadvantages of this new concept are also discussed.

TUBING↗

Application of OpenFOAM to Plume Impingement in Space Environments

After 30 years of continuous human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. Travelling beyond low earth orbit requires NASA to learn how humans can live in Deep Space environments – beyond the protection of Earth’s magnetosphere and at distances from Earth that prevent a quick return in case of trouble. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to be put in orbit around the moon to act as a home base for Lunar exploration for NASA astronauts. The Gateway Lunar outpost will be built incrementally, via modules which will arrive at separate times and dock to the existing structure. The incremental addition of Gateway modules, and the docking of visiting vehicles, is achieved via a sequence of firings from the approaching body’s onboard reaction control system (RCS) thrusters to achieve the required approach trajectory. The typical hypergolic chemical RCS thrusters work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking process. The exhaust gas from the RCS thrusters form plumes that expand into the vacuum of space and can impinge onto the outer surfaces of the Lunar Gateway, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion - all mechanisms that can damage the Lunar Gateway and must be minimized. Both permanent and visiting modules will have this RCS thruster exhaust impingement problem. This research aims to establish existing OpenFOAM solvers as a methodology for improving simulation techniques of rocket exhaust plume impingement in space environments. The flow structure of a plume in a space environment is complex; a plume that originates from a hypergolic chemical RCS thruster and expands into a vacuum will experience several regimes of rarefication. This range includes the continuum flow in the rocket nozzle through the fully rarefied free molecular flow further from the nozzle. The flow physics is different at these two extremes, and as such, the simulation approach for plumes is generally divided into a traditional computational fluid dynamics (CFD) simulation in and near the nozzle which is coupled to a subsequent direct simulation Monte Carlo (DSMC) simulation. At this time, the scope of this research is developing, verifying, and validating a method using existing solvers in the OpenFOAM framework for performing coupled CFD/DSMC calculations to determine the extent of plume impingement loading on generic space structures. This presentation will detail code-to-code comparisons between the hyStrath dsmcFoam+ solver, developed using OpenFOAM and available as open-source, and NASA’s in-house DSMC Analysis Code (DAC). Comparisons to several open-source publication findings using DAC [3,4] are presented, and advantages of using an OpenFOAM based solver are also discussed. The presentation concludes with a discussion of future work, and a plan for coupling the dsmcFoam+ solver with CFD simulations of chemical rocket engines for unified coupled plume simulation.

DSMC↗

Modeling heat pipe startup and noncondensable gases in Sockeye

For this work, a one-dimensional gas mixture flow model was developed and implemented in the heat pipe code Sockeye to model the effects of noncondensable 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 and noncondensable gas models were tested against experimental data for sodium heat pipes, showing excellent agreement. Additionally, the newly developed gas mixture model for modeling noncondensable 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↗

Aerothermodynamics of Pyrolizing Surfaces in Hypersonic Rarefied Flows

Direct simulation Monte Carlo (DSMC) calculations of rarefied flows about entry bodies typically employ a fixed surface temperature or a radiative-equilibrium energy balance to compute that temperature. Such boundary conditions neglect any effects of heat capacitance and heat conduction in the spacecraft heat shield and, therefore, provide an upper bound for the surface temperature. Such calculations also neglect pyrolysis from the heat shield which can be significant for a high-energy incident flow at very low densities. Accurate prediction of both heating and aerodynamic forces requires including pyrolysis and surface heat transfer in the models for gas-surface interaction employed in DSMC methods. Although these physical models have long appeared in various continuum flow calculation codes, they have only recently appeared in DSMC codes which are required to simulate rarefied flows during entry at high altitudes. In the current implementation, routines from the widely distributed Charring Material Thermal Response and Ablation (CMA) program are coupled into a DSMC code to calculate the one-dimensional heat transfer into the carbon phenolic heat shield at each point on a vehicle surface. Temperature-dependent material properties, surface re-radiation, and in-depth pyrolysis were included in the calculation, but surface ablation was neglected. Sample calculations for entry of the Galileo probe into the atmosphere of Jupiter demonstrate that including pyrolysis in the model leads to significant differences in predicted aerodynamics. Granted, the drag coefficient does not depend strongly on the surface temperature which can itself be significantly below the radiative equilibrium value during entry. However, the surface mass flux due to pyrolysis of the material is significant once the probe drops to altitudes characterized by transition flow. This leads to a noticeable increase in drag and a decrease in heating compared to a body without pyrolysis.

Haas, Brian L.↗

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↗

Basic governing equations for the flight regimes of aeroassisted orbital transfer vehicles

The basic governing equations for the low-density, high-enthalpy flow regimes expected in the shock layers over the heat shields of the proposed aeroassisted orbital transfer vehicles are derived by combining and extending existing theories. The conservation equations are derived from gas kinetic principles for a four-component ionized gas consisting of neutral molecules, neutral atoms, singly ionized ions, and electrons, assuming a continuum flow. The differences among translational-rotational, vibrational, and electron temperatures are accounted for, as well as chemical nonequilibrium and electric-charge separation. Expressions for convective and viscous fluxes, transport properties, and the terms representing interactions among various energy modes are given explicitly. The expressions for the rate of electron-vibration energy transfer, which violates the Landau-Teller conditions, is derived by solving the system of master equations accounting for the multiple-level transitions.

Lee, J.-H.↗

Basic Governing Equations for the Flight Regimes of Aeroassisted Orbital Transfer Vehicles

The basic governing equations for the low-density, high-enthalpy flow regimes expected in the shock layers over the heat shields of the proposed aeroassisted orbital transfer vehicles are derived by combining and extending existing theories. The conservation equations are derived from gas kinetic principles for a four-component ionized gas consisting of neutral molecules, neutral atoms, singly ionized ions, and electrons, assuming a continuum flow. The differences among translational-rotational, vibrational, and electron temperatures are accounted for, as well as chemical nonequilibrium and electric-charge separation. Expressions for convective and viscous fluxes, transport properties, and the terms representing interactions among various energy modes are explicitly given. The expressions for the rate of electron-vibration energy transfer, which violates the Landau-Teller conditions, are derived by solving the system of master equations accounting for the multiple-level transitions.

Lee, Jong-Hun↗

Decoupled predictions of radiative heating in air using a particle simulation method

A particle simulation technique with decoupled radiation is used to estimate the radiative emission along the stagnation streamline and the radiative heating at the stagnation point of a blunt-nosed vehicle entering the earth's atmosphere at hypersonic speed. The direct simulation Monte Carlo (DSMC) method is used to compute the fluid mechanics of the weakly ionized flow. The radiative emission is computed using the NEQAIR computer code with the radiation decoupled from the flowfield solution. It is concluded that the new models make it possible to decrease the predicted total radiative heating at the stagnation point of the vehicle by a factor of 15. A comparison of the DSMC approach with a continuum flow model shows that the total predicted radiative heating estimates agree within a factor of 2.

Boyd, Iain D.↗

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗

Efficient Three-Dimensional Direct Simulation Monte Carlo for Complex Geometry Problems

The simulation of flowfields in the transition flow regime is notoriously difficult with high demands on computer resources (CPU time and storage) and user expertise/labor. This paper describes a new, efficient code which has been developed to simulate high Knudsen number flowfields in three dimensions about bodies of arbitrarily complex geometry. The algorithm has been tested over a wide range of conditions, from free molecular to near-continuum flow regimes, for slender and blunt bodies, for re-entry vehicles and spacecraft. A series of validation tests have been conducted using both wind-tunnel measurements and flight data.

Rault, Didier F. G.↗

Plume-Free Stream Interaction Heating Effects During Orion Crew Module Reentry

During reentry of the Orion Crew Module (CM), vehicle attitude control will be performed by firing reaction control system (RCS) thrusters. Simulation of RCS plumes and their interaction with the oncoming flow has been difficult for the analysis community due to the large scarf angles of the RCS thrusters and the unsteady nature of the Orion capsule backshell environments. The model for the aerothermal database has thus relied on wind tunnel test data to capture the heating effects of thruster plume interactions with the freestream. These data are only valid for the continuum flow regime of the reentry trajectory. A Direct Simulation Monte Carlo (DSMC) analysis was performed to study the vehicle heating effects that result from the RCS thruster plume interaction with the oncoming freestream flow at high altitudes during Orion CM reentry. The study was performed with the DSMC Analysis Code (DAC). The inflow boundary conditions for the jets were obtained from Data Parallel Line Relaxation (DPLR) computational fluid dynamics (CFD) solutions. Simulations were performed for the roll, yaw, pitch-up and pitch-down jets at altitudes of 105 km, 125 km and 160 km as well as vacuum conditions. For comparison purposes (see Figure 1), the freestream conditions were based on previous DAC simulations performed without active RCS to populate the aerodynamic database for the Orion CM. Other inputs to the analysis included a constant Orbital reentry velocity of 7.5 km/s and angle of attack of 160 degrees. The results of the study showed that the interaction effects decrease quickly with increasing altitude. Also, jets with highly scarfed nozzles cause more severe heating compared to the nozzles with lower scarf angles. The difficulty of performing these simulations was based on the maximum number density and the ratio of number densities between the freestream and the plume for each simulation. The lowest altitude solutions required a substantial amount of computational resources (up to 1800 processors) to simulate approximately 2 billion molecules for the refined (adapted) solutions.

Marichalar, J.↗

Inferences on the Emplacement Dynamics of Martian Impact Crater Ejecta: Constraints from Mola Topography

Lobate ejecta deposits surround many of the younger impact craters on Mars. Viking Orbiter images indicate the distal parts of the ejecta blankets of these lobate craters are characterized by ramparts. In the absence of detailed topographic data for characterizing the topology of these apparently fluidized ejecta deposits, physical models have relied upon their morphologic characteristics. The most widely accepted model for the formation of such rampart ejecta deposits on Mars invokes fluidization of the ejecta to produce one or more viscous flow lobes. The availability of high-precision topographic data from the Mars Orbiter Laser Altimeter [4,51 facilitates a more quantitative examination of the physical processes involved in the formation of rampart ejecta deposits on Mars. Here we investigate the emplacement constraints that can be developed from the dimensions, topography, and morphology of martian rampart craters. The primary assumptions we have adopted are: (1) the ejecta blanket is emplaced as a continuum flow over the martian surface, rather than an airfall deposit, and (2) that the observable dimensions of the deposits are indicative of flow dimensions during emplacement.

Garvin, J. B.↗

Fuel-injector/air-swirl characterization

The objectives of this program are to establish an experimental data base documenting the behavior of gas turbine engine fuel injector sprays as the spray interacts with the swirling gas flow existing in the combustor dome, and to conduct an assessment of the validity of current analytical techniques for predicting fuel spray behavior. Emphasis is placed on the acquisition of data using injector/swirler components which closely resemble components currently in use in advanced aircraft gas turbine engines, conducting tests under conditions that closely simulate or closely approximate those developed in actual combustors, and conducting a well-controlled experimental effort which will comprise using a combination of low-risk experiments and experiments requiring the use of state-of-the-art diagnostic instrumentation. Analysis of the data is to be conducted using an existing, TEACH-type code which employs a stochastic analysis of the motion of the dispersed phase in the turbulent continuum flow field.

Mcvey, J. B.↗

Ablation of silicate particles in high-speed continuum and transition flow with application to the collection of interplanetary dust particles

The ablative deceleration of spheres in the continuum and slip regimes is studied using spherical 7.1-micron-diam soda-lime glass particles launched from vacuum at about 4500 m/sec speed through a 13-micron-thick plastic film into a capture chamber containing Xe at 0.1 or 0.2 atm pressure and 295 K temperature. The results of SEM examinations of the collected ablated particles showed that the ratio of the ablated-particle radius (Rf) to the initial radius (R0) increased with gas pressure (from Rf/R0 about 0.67 at 0.1 atm, to about 0.88 at 0.2 atm). A model was developed to describe the ablation and deceleration of spheres in high-speed continuum and slip flow. The pressure dependence predicted by the model agreed with experimental results.

Rulison, Aaron J.↗

Wind tunnel mass spectrometer design study

The sampling probe to measure the free stream gas composition in a hypersonic wind tunnel is described. Using a free expansion of the continuum flow, a series of sampling inlet-pumping stages is evaluated. Mass analysis is carried out with a quadrupole mass spectrometer employing a crossed electron beam, molecular beam ion source, and an off-axis electron multiplier. The spatial distribution of the sampled molecular flux and the pumping requirements are evaluated for varying inlet conditions of gas density and flow velocity. The signal-to-noise ratio of the molecular beam formed by this series of coaxial inlets is also determined. The signal-to-noise ratio of the unscattered flux is sufficiently high that an unambiguous determination of the beam composition and of the free stream gas properties can be made. The magnitude of the detected signal is sufficiently high that measurement of trace constituents two orders of magnitude below the primary beam signal can be made.

Hueser, J. E.↗