Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonequilibrium dynamics”

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 37 records · Page 2

Comparison of methodologies for describing relaxation in nonequilibrium gaseous systems

The heat transfer process in hypervelocity vehicles is dominated by nonequilibrium gas dynamics. One model used in computational fluid dynamics (CFD) codes to predict hypervelocity heat transfer is the 'two-temperature' model. An analysis has been made to test the validity of the two-temperature model for predicting another nonequilibrium phenomenon, sound absorption and deviation of signal speed in a high temperature gas. It is found that the two temperature model's prediction capabilities degenerate with increasing temperature. These results are felt to have significance concerning the two-temperature's ability to predict heat transfer in hypervelocity flows.

Schreiber, Willard C.↗

New outlook on control of crystalline and chemical perfection during growth of silicon

Significant progress has been made in our understanding of the Czochralski crystal growth process with the realization that the incorporation of oxygen into silicon is directly related to the internal gettering capability of wafers during device fabrication. It was also recognized that the electronic properties exhibited by silicon during various stages of device fabrication were significantly affected by the thermal history of the silicon during the post-growth cool-down period. Turbulent melt convection, induced by unavoidable destabilizing thermal gradients, was found to interfere with homogeneous dopant (and oxygen) incorporation and to influence markedly the dynamics of nonequilibrium point defects in the solidified silicon matrix during the cool-down period. In view of the unavoidability of destabilizing thermal gradients in conventional crystal growth configurations, melt stabilization through the application of magnetic fields is generally considered a viable approach. Control of heat input to the melt through heat pipe systems, as used during growth of germanium, cannot be applied to industrial growth of silicon. Recent studies have shown that heat exchange systems located coaxially about a growing crystal can be used to stabilize and control not only heat transfer in the grown crystal, but also in the melt adjacent to the solidification interface.

Witt, A. F.↗

A study of numerical methods for hyperbolic conservation laws with stiff source terms

The proper modeling of nonequilibrium gas dynamics is required in certain regimes of hypersonic flow. For inviscid flow this gives a system of conservation laws coupled with source terms representing the chemistry. Often a wide range of time scales is present in the problem, leading to numerical difficulties as in stiff systems of ordinary differential equations. Stability can be achieved by using implicit methods, but other numerical difficulties are observed. The behavior of typical numerical methods on a simple advection equation with a parameter-dependent source term was studied. Two approaches to incorporate the source term were utilized: MacCormack type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. Various comparisons over a wide range of parameter values were made. In the stiff case where the solution contains discontinuities, incorrect numerical propagation speeds are observed with all of the methods considered. This phenomenon is studied and explained.

Leveque, R. J.↗

Theoretical investigation of non-equilibrium chemistry and optical radiation in hypersonic flow fields

Future space vehicles returning from distant missions or high earth orbits may enter the upper regions of the atmosphere and use aerodynamic drag to reduce their velocity before they skip out of the atmosphere and enter low earth orbit. The Aeroassist Flight Experiment (AFE) is designed to explore the special problems encountered in such entries. A computer code was developed to calculate the radiative transport along line-or-sight in the general 3-D flow field about an arbitrary entry vehicle, if the temperatures and species concentrations along the line-of-sight are known. The radiative heating calculation at the stagnation point of the AFE vehicle along the entry trajectory was performed, including a detailed line-by-line accounting of the radiative transport in the vacuum ultraviolet (below 200 nm) by the atomic N and O lines. A method was developed for making measurements of the haze particles in the Titan atmosphere above 200 km altitude. Several other tasks of a continuing nature, to improve the technical ability to calculate the nonequilibrium gas dynamic flow field and radiative heating of entry vehicles, were completed or advanced.

Whiting, Ellis E.↗

A fully coupled flow simulation around spacecraft in low earth orbit

The primary objective of this investigation is to provide a full flow simulation of a spacecraft in low earth orbit (LEO). Due to the nature of the environment, the simulation includes the highly coupled effects of neutral particle flow, free stream plasma flow, nonequilibrium gas dynamics effects, spacecraft charging and electromagnetic field effects. Emphasis is placed on the near wake phenomenon and will be verified in space by the Wake Shield Facility (WSF) and developed for application to Space Station conditions as well as for other spacecraft. The WSF is a metallic disk-type structure that will provide a controlled space platform for highly accurate measurements. Preliminary results are presented for a full flow around a metallic disk.

Justiz, C. R.↗

Scale-Similar Models for Large-Eddy Simulations

Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Sarghini, F.↗

On the Dynamics of Computing a Chemically Relaxed Nonequilibrium Flow

In order to gain insights into the strong dependence of numerical solutions on initial data for finite time steps, a set of nonlinear test problems rich enough to capture the behavior of difference schemes were recently identified and the numerical basins of attraction for these problems were computed using commonly used time discretizations in CFD. Our study revealed a wealth of surprisingly nonlinear behavior of numerical schemes that were not observed before, in particular for the implicit time discretizations that are commonly used in CFD. The goal of this work is to apply these tools to study a practical model from non-equilibrium flowfield relaxation. This type of problem arises in chemically nonequilibrium hypersonic flows such as in a shock tube experiment or an expanding nozzle. Here we consider a reacting mixture of (N2, N) for an inviscid one-dimensional steady model. Preliminary numerical results indicate that, aside from the possibility of spurious numerical solutions being introduced by the time discretizations, limitations on the model for physical or accurate solutions may also play a part in the dynamics observed.

Sweby, P. K.↗

Overview of X-38 Hypersonic Wind Tunnel Data and Comparison with Numerical Results

A NASA team of engineers has been organized to design a crew return vehicle for returning International Space Station crew members from orbit. The hypersonic characteristics of this X-23/X-2&4 derived crew return vehicle (designated X-38) are being evaluated in various wind tunnels in support of this effort. Aerodynamic data has been acquired in three NASA hypersonic facilities at Mach 20, and Mach 6. Computational Fluid Dynamics tools have been applied at the appropriate wind tunnel conditions to make comparisons with portions of this data. Experimental data from the Mach 6 Air and CF4 facilities illustrate a net positive pitching moment increment due to density ratio, as well as increased elevon effectiveness. Chemical nonequilibrium computational fluid dynamics solutions at flight conditions reinforce this conclusion.

Campbell, Charles H.↗

Stagnation point nonequilibrium radiative heating and the influence of energy exchange models

A nonequilibrium radiative heating prediction method has been used to evaluate several energy exchange models used in nonequilibrium computational fluid dynamics methods. The radiative heating measurements from the FIRE II flight experiment supply an experimental benchmark against which different formulations for these exchange models can be judged. The models which predict the lowest radiative heating are found to give the best agreement with the flight data. Examination of the spectral distribution of radiation indicates that despite close agreement of the total radiation, many of the models examined predict excessive molecular radiation. It is suggested that a study of the nonequilibrium chemical kinetics may lead to a correction for this problem.

Hartung, Lin C.↗

Stagnation Point Nonequilibrium Radiative Heating and the Influence of Energy Exchange Models

A nonequilibrium radiative heating prediction method has been used to evaluate several energy exchange models used in nonequilibrium computational fluid dynamics methods. The radiative heating measurements from the FIRE II flight experiment supply an experimental benchmark against which different formulations for these exchange models can be judged. The models which predict the lowest radiative heating are found to give the best agreement with the flight data. Examination of the spectral distribution of radiation indicates that despite close agreement of the total radiation, many of the models examined predict excessive molecular radiation. It is suggested that a study of the nonequilibrium chemical kinetics may lead to a correction for this problem.

Hartung, Lin C.↗

Gravitational potential energy of the earth: A spherical harmonic approach

A spherical harmonic equation for the gravitational potential energy of the earth is derived for an arbitrary density distribution by conceptually bringing in mass-elements from infinity and building up the earth shell upon spherical shell. The zeroth degree term in the spherical harmonic equation agrees with the usual expression for the energy of a radial density distribution. The second degree terms give a maximum nonhydrostatic energy in the mantle and crust of -2.77 x 10 to the twenty-ninth power ergs, an order of magnitude. If the earth is assumed to be a homogeneous viscous oblate spheroid relaxing to an equilibrium shape, then a lower limit to the mantle viscosity of 1.3 x 10 to the twentieth power poises is found by assuming the total geothermal flux is due to viscous dissipation. If the nonequilibrium figure is dynamically maintained by the earth acting as a heat engine at one per cent efficiency, then the viscosity is ten to the twenty second power poises, a number preferred by some as the viscosity of the mantle.

Rubincam, D. P.↗

Gravitational potential energy of the earth - A spherical harmonic approach

A spherical harmonic equation for the gravitational potential energy of the earth is derived for an arbitrary density distribution by conceptually bringing in mass-elements from infinity and building up the earth shell upon spherical shell. The zeroth degree term in the spherical harmonic expansion agrees with the usual expression for the energy of a radial density distribution. The second degree terms give a maximum nonhydrostatic energy in the crust and mantle of -2.77 x 10 to the 29th ergs, an order of magnitude below McKenzie's (1966) estimate. McKenzie's result stems from mathematical error. Our figure is almost identical with Kaula's (1963) estimate of the minimum shear strain energy in the mantle, a not unexpected result on the basis of the virial theorem. If the earth is assumed to be a homogeneous viscous oblate spheroid relaxing to an equilibrium shape, then a lower limit to the mantle viscosity of 1.3 x 10 to the 20th P is found by assuming that the total geothermal flux is due to viscous dissipation of energy. This number is almost six orders of magnitude below MacDonald's (1966) estimate of the viscosity and removes his objection to convection. If the nonequilibrium figure is dynamically maintained by the earth acting as a heat engine at 1% efficiency, then the viscosity is 10 to the 22nd P, a number preferred by Cathles (1975) and Peltier and Andrew (1976) as the viscosity of the mantle.

Rubincam, D. P.↗

Aerodynamics of Stardust Sample Return Capsule

Successful return of interstellar dust and cometary material by the Stardust Sample Return Capsule requires an accurate description of the Earth entry vehicle's aerodynamics. This description must span the hypersonic-rarefied, hypersonic-continuum, supersonic, transonic, and subsonic flow regimes. Data from numerous sources are compiled to accomplish this objective. These include Direct Simulation Monte Carlo analyses, thermochemical nonequilibrium computational fluid dynamics, transonic computational fluid dynamics, existing wind tunnel data, and new wind tunnel data. Four observations are highlighted: 1) a static instability is revealed in the free-molecular and early transitional-flow regime due to aft location of the vehicle s center-of-gravity, 2) the aerodynamics across the hypersonic regime are compared with the Newtonian flow approximation and a correlation between the accuracy of the Newtonian flow assumption and the sonic line position is noted, 3) the primary effect of shape change due to ablation is shown to be a reduction in drag, and 4) a subsonic dynamic instability is revealed which will necessitate either a change in the vehicle s center-of-gravity location or the use of a stabilizing drogue parachute.

Mitcheltree, R. A.↗