ON THE RELAXATION OF GASES TOWARD CONTINUUM FLOW
Gas dynamics - relaxation of gases towards continuum flow
SEARCH · Search NASA
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.
Gas dynamics - relaxation of gases towards continuum flow
We proposed to create a single computational code incorporating methods that can model both rarefied and continuum flow to enable the efficient simulation of flow about space craft and high altitude hypersonic aerospace vehicles. The code was to use a single grid structure that permits a smooth transition between the continuum and rarefied portions of the flow. Developing an appropriate computational boundary between the two regions represented a major challenge. The primary approach chosen involves coupling a four-speed Lattice Boltzmann model for the continuum flow with the DSMC method in the rarefied regime. We also explored the possibility of using a standard finite difference Navier Stokes solver for the continuum flow. With the resulting code we will ultimately investigate three-dimensional plume impingement effects, a subject of critical importance to NASA and related to the work of Drs. Forrest Lumpkin, Steve Fitzgerald and Jay Le Beau at Johnson Space Center. Below is a brief background on the project and a summary of the results as of the end of the grant.
Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.
Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.
An efficient particle simulation technique, developed for use on vector architecture based supercomputers for studying hypersonic rarefied gas flows is employed to simulate the complex wake generated by Mach six flow over a 10 deg half-angle wedge for freestream Reynolds numbers of 1780 and 3560. Data obtained are compared against higher Reynolds number experimental results. Simulations utilized as many as 10 to the 5th computational cells and 10 to the 7th simulated particles having power-law interaction potentials. A code performance of 1.8-2.4 microsec of Cray-2 CPU time to process a single particle per timestep is achieved. Diffuse adiabatic and isothermal wedge surface models are used in this investigation. Although the wedge geometry is two-dimensional, the simulation incorporates a width-wise direction, resulting in a three-dimensional computation.
Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo (DSMC) method, and the continuum approach solves the appropriate equations of fluid flow. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code that is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuum of a previous study that considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry model employed.
Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo method (DSMC), and the continuum approach solves the Euler equations. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code which is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuation of a previous study which considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry models employed.
Rarefied gas dynamics - flow regime between continuum flow and free molecular flow
Heat transfer phenomena of rarefied gas flows is discussed based on a literature survey of analytical and experimental rarefied gas dynamics. Subsonic flows are emphasized for the purposes of meteorological thermometry in the high atmosphere. The heat transfer coefficients for three basic geometries are given in the regimes of free molecular flow, transition flow, slip flow, and continuum flow. Different types of heat phenomena, and the analysis of theoretical and experimental data are presented. The uncertainties calculated from the interpolation rule compared with the available experimental data are discussed. The recovery factor for each geometry in subsonic rarefied flows is also given.
A new thermochemical nonequilibrium formulation for hypersonic transitional flows of air is presented. Air is assumed to have five chemical species (N2, O2, NO, N, O) and three temperatures corresponding to the translational, rotational, and vibrational modes of energy. In the present study, the no-slip boundary conditions are replaced by slip boundary conditions to extend the range of the Navier-Stokes equations to high-speed low-density flows.
Shock tube determination of drag coefficient of small spherical particles accelerating in laminar, non-reacting, incompressible continuum flow
Observations of the interplanetary medium are considered together with the structure of the solar atmosphere and the kinetic properties of the interplanetary plasma. Methods of gas kinetics and continuum flows are examined, taking into account microscopic and macroscopic representations, the kinetic flow equations, the formal solution of the Vlasov equation, and the continuum flow equations. The collective particle behavior of that interplanetary gas is discussed along with the hydrodynamic coronal expansion, a wave-pump problem, the free expansion phenomenon, and a generalized free expansion problem.
Results of an investigation to determine the full scale drag coefficient in the high speed, low density regime of the Viking lander capsule 1 entry vehicle are presented. The principal flight data used in the study were from onboard pressure, mass spectrometer, and accelerometer instrumentation. The hypersonic continuum flow drag coefficient was unambiguously obtained from pressure and accelerometer data; the free molecule flow drag coefficient was indirectly estimated from accelerometer and mass spectrometer data; the slip flow drag coefficient variation was obtained from an appropriate scaling of existing experimental sphere data. Comparison of the flight derived drag hypersonic continuum flow regime except for Reynolds numbers from 1000 to 100,000, for which an unaccountable difference between flight and ground test data of about 8% existed. The flight derived drag coefficients in the free molecule flow regime were considerably larger than those previously calculated with classical theory. The general character of the previously determined temperature profile was not changed appreciably by the results of this investigation; however, a slightly more symmetrical temperature variation at the highest altitudes was obtained.
Aerodynamic heating of spacecraft during reentry, investigating characteristics for low density rarefied flow and high density continuum flow regimes
The turbulence structure of weakly ionized, continuum plasma is analyzed with emphasis on the turbulent fluctuation and the transport properties of the ionized species. The mean gradients of the velocity and ionized-species concentrations are considered to be uniform. It is assumed that the turbulence Reynolds number is small and that the turbulence field is homogeneous in order to close the dynamic equations of turbulence. It was found that the ratio of the Debye shielding length to the local integral scale of turbulence must be sufficiently small in order for the detailed quasi-neutrality to prevail. The condition of the detailed quasi-neutrality is necessary, in addition to the equal mean concentrations of electrons and ions, in order for the turbulence properties involving the electrons to be the same as those involving the ions.
The flow field in a cometary coma varies from free molecular flow to transition flow and finally to continuum flow. General discussions of these flow regimes with two-phase effects are given. Some numerical results of similarity solution of a strong shock propagation in a two-phase flow are presented.
Impact pressure probe response characteristics in free-molecular and continuum flows with transition Knudsen numbers
The reentry flow field and thermal environment around the straight wing shuttle orbiter vehicle were determined. Both rarefied and continuum flow fields and associated heating rates on various configurations representative of the orbiter at high angle of attack were calculated. Rarefied flow fields and heating rates were computed by the Monte Carlo direct simulation technique for altitudes above 82.3 km. Continuum inviscid flow fields were calculated by 2-D unsteady and 3-D steady finite difference/artificial viscosity methods and also by a 2-D shock layer analysis technique. Viscous flow fields and heating rates in the continuum regime were computed by a boundary layer integral matrix method for laminar flow and by an aerodynamic surface heating technique for turbulent flow. Shapes considered in the study included flat plates (representing the underside of the orbiter fuselage or the wing MAC), orbiter fuselage cross sections, orbiter wing airfoils, and 3-D orbiter configurations, all at high angle of attack (40 - 60 deg). The theoretical results showed good agreement with measured pressure and heat transfer data.