Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonequilibrium physics”

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

Point-implicit relaxation strategies for viscous, hypersonic flows

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's symmetric total variation diminishing scheme. Newton relaxation of the fully coupled equation set is employed on a cell-to-cell basis. Under-relaxation of the inviscid and over-relaxation of the viscous contributions to the residual are implemented. Computational work is easily partitioned among many processors in an asynchronous, dynamic mode for convergence acceleration. An overview of the physical models employed herein for thermochemical nonequilibrium is included. Several test cases and comparisons with experimental data are presented involving hypersonic flow over blunt bodies which illustrate the qualitative and quantitative capabilities of this approach.

Gnoffo, P. A.↗

Hypersonic flow - Physical issues and computational simulation

Computational results are presented on hypersonic flows, with emphasis on thermal and chemical nonequilibrium effects and the role played by grids. The continuum equations and their boundary conditions are presented, and the direct simulation Monte Carlo method is described. The importance of using the correct computational and physical models in describing hypersonic nonequilibrium flows is demonstrated.

Hassan, H. A.↗

Nonequilibrium, large-amplitude MHD fluctuations in the solar wind

Compressible MHD simulations in one dimension with three-dimensional vectors are used to investigate a number of processes relevant to problems in interplanetary physics. The simulations indicate that a large-amplitude nonequilibrium (e.g., linearly polarized) Alfvenic wave, which always starts with small relative fluctuations in the magnitude B of the magnetic field, typically evolves to flatten the magnetic profile in most regions. Under a wide variety of conditions B and the density rho become anticorrelated on average. If the mean magnetic field is allowed to decrease in time, the point where the transverse magnetic fluctuation amplitude delta B(sub T) is greater than the mean field B(sub 0) is not special, and large values of delta B(sub T)/B(sub 0) do not cause the compressive thermal energy to increase remarkably or the wave energy to dissipate at an unusually high rate. Nor does the 'backscatter' of the waves that occurs when the sound speed is less than the Alfven speed result, in itself, in substantial energy dissipation, but rather primarily in a phase change between the magnetic and velocity fields. For isolated wave packets the backscatter does not occur for any of the parameters examined; an initial radiation of acoustic waves away from the packet establishes a stable traveling structure. Thus these simulations, although greatly idealized compared to reality, suggest a picture in which the interplanetary fluctuations should have small deltaB and increasingly quasi-pressure balanced compressive fluctuations, as observed, and in which the dissipation and 'saturation' at delta B(sub T)/B(sub 0) approximately = 1 required by some theories of wave acceleration of the solar wind do not occur. The simulations also provide simple ways to understand the processes of nonlinear steepening and backscattering of Alfven waves and demonstrate the existence of previously unreported types of quasi-steady MHD states.

Roberts, D. Aaron↗

Probe and radar electron temperatures in an isotropic nonequilibrium plasma.

Electron temperatures measured by electrostatic probes and radar backscatter are distinct physical quantities, the temperature from each technique determined from a different moment of the electron-distribution function. Numerical inequality of temperatures results from a non-Maxwellian electron-distribution function or, equivalently, from a nonequilibrium electron plasma. Probe and backscatter electron temperatures are studied for low- and high-energy (isotropic) distortions of the distribution function. The nonequilibrium plasma generally produces higher probe than backscatter temperatures; however, the temperature difference is small for distortions due to realistic photoelectron populations. If the ionosphere is in a highly nonequilibrium state, probe and backscatter temperatures would differ from the temperature characterizing the average electron kinetic energy, and a single temperature applicable to a variety of physical processes would no longer exist.

Hoegy, W. R.↗

Updated Stagnation-Point Aeroheating Correlations for Mars Entry

The objective of this work was to develop new stagnation-point aeroheating engineering correlations for Mars entry vehicles. New convective and radiative heating relations were formulated over a wide range of entry conditions. These relations were based on computational fluid dynamics simulations informed by recent experimental testing and modeling enhancements. For convective heating, a relation similar to the Sutton-Graves relation was developed that is accurate to within ±25% across the domain of relevant Mars entry conditions. For radiative heating, two different functional forms were developed: exponential and polynomial relations. Each of these relation, which represents two different levels of complexity for implementation, were able to predict the radiative heating within about ±25%. The correlations were tested by applying them to the Mars Pathfinder entry trajectory to demonstrate their applicability. These new correlations are significant improvement over existing relations in terms of the accuracy, domain of applicability, and the captured physics, including CO2 inferred radiation and thermochemical nonequilibrium.

Thomas K West, IV↗

Vibration-dissociation coupling in nonequilibrium flows

The purpose of the study is to evaluate all important physical phenomena contributing to the vibration-dissociation problem by developing a numerical technique for the simulation of the vibrational relaxation and dissociation/recombination of N(2) in an isothermal box. A set of coupled master equations representing each vibrational level of N(2) is solved, and both V-V and V-T exchanges between molecules and atoms are included, along with the effects of dissociation and recombustion from/to each level. It is observed that the inverse range parameter has a strong impact on the vibrational exchange rates and relaxation. In a heating case, the lower vibrational levels relax through a series of equilibrium distributions, while in a cooling case, a population inversion is formed in the upper levels. Other cases of heating and cooling will be simulated, and together with this work, used to develop a matrix of nonequilibrium relaxation runs.

Landrum, D. Brian↗

The development of an explicit thermochemical nonequilibrium algorithm and its application to compute three dimensional AFE flowfields

This study presents a three-dimensional explicit, finite-difference, shock-capturing numerical algorithm applied to viscous hypersonic flows in thermochemical nonequilibrium. The algorithm employs a two-temperature physical model. Equations governing the finite-rate chemical reactions are fully-coupled to the gas dynamic equations using a novel coupling technique. The new coupling method maintains stability in the explicit, finite-rate formulation while allowing relatively large global time steps. The code uses flux-vector accuracy. Comparisons with experimental data and other numerical computations verify the accuracy of the present method. The code is used to compute the three-dimensional flowfield over the Aeroassist Flight Experiment (AFE) vehicle at one of its trajectory points.

Palmer, Grant↗

Thermochemical nonequilibrium in atomic hydrogen at elevated temperatures

A numerical study of the nonequilibrium flow of atomic hydrogen in a cascade arc was performed to obtain insight into the physics of the hydrogen cascade arc. A rigorous mathematical model of the flow problem was formulated, incorporating the important nonequilibrium transport phenomena and atomic processes which occur in atomic hydrogen. Realistic boundary conditions, including consideration of the wall electrostatic sheath phenomenon, were included in the model. The governing equations of the asymptotic region of the cascade arc were obtained by writing conservation of mass and energy equations for the electron subgas, an energy conservation equation for heavy particles and an equation of state. Finite-difference operators for variable grid spacing were applied to the governing equations and the resulting system of strongly coupled, stiff equations were solved numerically by the Newton-Raphson method.

Scott, R. K.↗

A shock-layer theory based on thirteen-moment equations and DSMC calculations of rarefied hypersonic flows

Grad's thirteen-moment equations are applied to the flow behind a bow shock under the formalism of a thin shock layer. Comparison of this version of the theory with Direct Simulation Monte Carlo calculations of flows about a flat plate at finite attack angle has lent support to the approach as a useful extension of the continuum model for studying translational nonequilibrium in the shock layer. This paper reassesses the physical basis and limitations of the development with additional calculations and comparisons. The streamline correlation principle, which allows transformation of the 13-moment based system to one based on the Navier-Stokes equations, is extended to a three-dimensional formulation. The development yields a strip theory for planar lifting surfaces at finite incidences. Examples reveal that the lift-to-drag ratio is little influenced by planform geometry and varies with altitudes according to a 'bridging function' determined by correlated two-dimensional calculations.

Cheng, H. K.↗

DSMC simulation of ionized rarefied flows

Recently a renewed interest has been exhibited in understanding contamination and space environmental effect in the LEO range of 200 to 600 km as well as in low density microelectronics manufacturing technologies. Realistic simulations must model the physics of the highly coupled effects of neutral and charged particle flows, thermodynamic nonequilibrium, surface charging, and electromagnetic field effects. The computational requirements are enormous for this level of modeling and are almost impossible on current single processor supercomputers. An effort was initiated to develop this capability on a massively parallel MIMD supercomputers. The Wake Shield Facility experiment will be used for test case calculations. Preliminary results indicate that parallel supercomputers may be the only computers capable of simulating the required level of physics and spatial resolution. Calculations which contain from 1 to 8 million simulation particles are presented.

Bartel, Timothy J.↗

Sub-T/g/ annealing studies of advanced epoxy-matrix graphite-fiber-reinforced composites

Annealing at temperatures below the glass transition temperature (sub-T/g/ annealing) in structural epoxy network glasses is found to affect long-term mechanical properties of symmetrically (+ or - 45 deg) reinforced 16-lamina graphite/epoxy composites. Post-cured specimens of Thornel 300 graphite/Narmco 5208 epoxy are sub-T(g) annealed at 140 C for about 10, 100, 1000, 10,000, and 100,000 min with a prior quenching above T(g). The ultimate tensile strength, strain-to-break, and static toughness of the composite materials are found to decrease as a function of sub-T(g) annealing time, while no weight loss is observed. Time-dependent changes in properties are explained by a physical aging concept based on free-volume decreases in the nonequilibrium glassy state of network epoxies, and results imply that changes in composite properties may occur with the service time.

Kong, E. S. W.↗

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models.The wall geometry and flow conditions of the DNS are representative of the experimental data for a Mach 4.9 turbulent boundary layer that was tested on a two-dimensional planar convexwall model in the high-speed blow-down wind tunnel located at the National Aerothermochemistry Laboratory at Texas A&M University. The favorable pressure gradient is induced bythe streamwise surface curvature. The DNS results will be compared with experimental measurements and RANS computations for the same flow conditions and wall geometry. The DNS datasets will also be analyzed to provide an improved physical understanding of hypersonic boundary-layer turbulence subject to mechanical nonequilibrium, including the validity of Morkovin’s hypothesis and the various available compressibility transformations, as well as to obtain turbulence statistics that are relevant to RANS modeling, including the distribution ofeddy viscosity, the budgets of turbulent kinetic energy and Reynolds stresses, and the velocity-temperature correlations. The RANS predictions will be obtained by using both the commonly used one- and two-equation models that assume a constant turbulent Prandtl number and an algebraic, variable Prandtl number model for the turbulent energy flux.

turbulent boundary layers↗

Nonequilibrium material effects on the behavior of polymeric composite matrices and their related composites

The effects of physical aging on the material properties of some linear and network macromolecular glasses are discussed. The free volume concept is used to describe this behavior. The effect of physical aging on properties of some uniaxial graphite/fiber epoxy resin composites is investigated using stress relaxation in both tensile and flexural modes. The matrix polymers used were resins both of which are based on a 4,4-methylenedianiline derivative of epichlorohydrin with diamino diphenyl sulfone (DDS) as the curing agent. The matrix resin, as used in the practical application in composites, not fully cured and the glass transition of the network was dependent on the curing schedule. The physical aging of the bulk crosslinked epoxy was found to depend on the annealing temperature, and the T sub g of the resin. The physical aging of the composite, monitored by the stress relaxation method, was found to be dependent on the testing direction.

Wilkes, G. L.↗

An analysis of artificial viscosity effects on reacting flows using a spectral multi-domain technique

Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.

Macaraeg, M. G.↗

Analysis of artificial viscosity effects on reacting flows using a spectral multidomain technique

Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.

Macaraeg, Michele G.↗

Direct simulation of rarefied hypersonic flows

As the capability of the space transportation vehicles (STV's) expand to meet the requirements for future space exploration and utilization, the effects of rarefied hypersonic flows will play a more significant role in defining the aerodynamic and aerothermodynamic performance of STV's. This is particularly true of the low lift/drag aeroassisted STV's where aerobraking occurs at relatively high altitudes and high velocity. Because of the limitations of the continuum description as expressed by the Navier-Stokes equations and the difficulties of solving the Boltzmann equation, the particle of molecular approach has been developed over the last three decades for modeling rarefied gas effects. The direct simulation Monte Carlo (DSMC) method of Bird is the most used method today for simulating rarefied flows. The DSMC method provides a direct physical simulation as opposed to a numerical solution of a set of model equations. This is accomplished by developing phenomenological models of the relevant physical events. The DSMC method accounts for translational, thermal, chemical, and radiative nonequilibrium effects. The general features of the DSMC method, the numerical requirements for obtaining meaningful results, the modeling used to simulate high temperature gas effects, and applications of the method to calculate the flow about an aeroassist flight experiment vehicle (AFE) are reviewed. The AFE simulates a geosynchronous return while entering the Earth's upper atmosphere at approximately 10 km/s. Results obtained using a general 3-D code are presented for the more rarefied portion of the atmospheric encounter (altitudes of 200 to 100 km) emphasizing surface, flowfield, and aerodynamic characteristics of the AFE. Finally, results obtained using axisymmetric and 1-D versions of the code are presented for lower altitude conditions.

Moss, James N.↗

Non equilibrium flow in an arc-jet

A method for simulating the flow of a nitrogen plasma in thermodynamic nonequilibrium is described. The method relies on the modeling of the relevant physical and chemical phenomena and on a fully coupled implicit finite volume upwind numerical approximation of the governing equations. The results agree with experimental data.

Rostand, P.↗