Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermochemical nonequilibrium”

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 127 records · Page 7

Semi-Analytic Reconstruction of Flux in Finite Volume Formulations

Semi-analytic reconstruction uses the analytic solution to a second-order, steady, ordinary differential equation (ODE) to simultaneously evaluate the convective and diffusive flux at all interfaces of a finite volume formulation. The second-order ODE is itself a linearized approximation to the governing first- and second- order partial differential equation conservation laws. Thus, semi-analytic reconstruction defines a family of formulations for finite volume interface fluxes using analytic solutions to approximating equations. Limiters are not applied in a conventional sense; rather, diffusivity is adjusted in the vicinity of changes in sign of eigenvalues in order to achieve a sufficiently small cell Reynolds number in the analytic formulation across critical points. Several approaches for application of semi-analytic reconstruction for the solution of one-dimensional scalar equations are introduced. Results are compared with exact analytic solutions to Burger s Equation as well as a conventional, upwind discretization using Roe s method. One approach, the end-point wave speed (EPWS) approximation, is further developed for more complex applications. One-dimensional vector equations are tested on a quasi one-dimensional nozzle application. The EPWS algorithm has a more compact difference stencil than Roe s algorithm but reconstruction time is approximately a factor of four larger than for Roe. Though both are second-order accurate schemes, Roe s method approaches a grid converged solution with fewer grid points. Reconstruction of flux in the context of multi-dimensional, vector conservation laws including effects of thermochemical nonequilibrium in the Navier-Stokes equations is developed.

Gnoffo, Peter A.↗

A Study of Ablation-Flowfield Coupling Relevant to the Orion Heatshield

The coupled interaction between an ablating surface and the surrounding aerothermal environment is studied. An equilibrium ablation model is coupled to the LAURA flowfield solver, which allows the char ablation rate (m-dot(sub c)) to be computed as part of the flowfield solution. The wall temperature (T(sub w)) and pyrolysis ablation rate (m-dot(sub g)) may be specified by the user, obtained from the steady-state ablation approximation, or computed from a a material response code. A 32 species thermochemical nonequilibrium flowfield model is applied, which permits the treatment of C, H, O, N, and Si containing species. Coupled ablation cases relevant to the Orion heatshield are studied. These consist of diffusion-limited oxidation cases with Avcoat as the ablation material. The m-dot(sub c) values predicted from the developed coupled ablation analysis were compared with those obtained from a typical uncoupled ablation analysis. The coupled results were found to be as much as 50% greater than the uncoupled values. This is shown to be a result of the cumulative effect of the two fundamental approximations inherent in the uncoupled analysis.

Johnston, Christopher O.↗

Fully-Implicit Navier-Stokes (FIN-S)

FIN-S is a SUPG finite element code for flow problems under active development at NASA Lyndon B. Johnson Space Center and within PECOS: a) The code is built on top of the libMesh parallel, adaptive finite element library. b) The initial implementation of the code targeted supersonic/hypersonic laminar calorically perfect gas flows & conjugate heat transfer. c) Initial extension to thermochemical nonequilibrium about 9 months ago. d) The technologies in FIN-S have been enhanced through a strongly collaborative research effort with Sandia National Labs.

Kirk, Benjamin S.↗

Hypersonic Chemically Reacting Boundary-Layer Stability using LASTRAC

Hypersonic boundary layer transition is critical to the design of all hypersonic vehicles due to its effect on the heat transfer into the vehicle surface and potential drag enhancement or reduction during reentry. Boundary layer transition and boundary layer stability analysis under hypersonic conditions has been studied for decades, yet there is ample room for improved accuracy and further investigations into the relevant phenomena. In this work, we present a recent implementation of chemical equilibrium, finite-rate chemistry, and thermochemical nonequilibrium capabilities into LASTRAC, an existing well-established boundary-layer stability analysis code. Verification against existing numerical results in the literature are presented. LASTRAC was previously able to address calorically perfect flows. By using solutions of the Parabolized Stability Equations (PSE) with chemical and thermal nonequilibrium, we are able to investigate the effects of chemical and thermal nonequilibrium on a variety of phenomena including stationary crossflow instability on a swept wing and 2nd mode instabilities over a wedge.

Kline, H. L.↗

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↗

A Multi-Architecture Approach for Implicit Computational Fluid Dynamics on Unstructured Grids

High-performance computing (HPC) architectures are trending toward manycore paradigms such as graphics processing units (GPUs). Approximately half of the top 100 publicly disclosed supercomputers in the world utilize GPU accelerators for performance. This is in contrast to a decade ago, where there were only a few such machines in the top 100. It is not currently possible to compile and run legacy central processing unit (CPU) software efficiently on GPUs without significant refactoring. Though a number of frameworks offering performance portability exist, none offer a standardized specification that is supported by all major hardware vendors. Additionally, experiences show that obtaining a high percentage of peak performance often requires architecture-specific code. This work details a pragmatic multi-architecture computational fluid dynamics library focused on aerospace problems across the speed range from low subsonic to hypersonic flows involving thermochemical nonequilibrium. A thin abstraction layer above NVIDIA CUDA C++ is utilized, which enables primarily single-source software currently capable of running efficiently on multicore CPUs, NVIDIA GPUs, AMD GPUs, and Intel GPUs. Results on various problems of interest across the speed range are presented and performance is compared between various architectures.

GPU↗

A Multi-Architecture Approach for Implicit Computational Fluid Dynamics on Unstructured Grids

High-performance computing (HPC) architectures are trending toward manycore paradigms such as graphics processing units (GPUs). Approximately half of the top 100 publicly disclosed supercomputers in the world utilize GPU accelerators for performance. This is in contrast to a decade ago, where there were only a few such machines in the top 100. It is not currently possible to compile and run legacy central processing unit (CPU) software efficiently on GPUs without significant refactoring. Though a number of frameworks offering performance portability exist, none offer a standardized specification that is supported by all major hardware vendors. Additionally, experiences show that obtaining a high percentage of peak performance often requires architecture-specific code. This work details a pragmatic multi-architecture computational fluid dynamics library focused on aerospace problems across the speed range from low subsonic to hypersonic flows involving thermochemical nonequilibrium. A thin abstraction layer above NVIDIA CUDA C++ is utilized, which enables primarily single-source software currently capable of running efficiently on multicore CPUs, NVIDIA GPUs, AMD GPUs, and Intel GPUs. Results on various problems of interest across the speed range are presented and performance is compared between various architectures.

GPU↗

Aeroheating Environments of Aerocapture Systems for Uranus Orbiters

Aeroheating environments for an aerocapture system enabling flagship-class science at Uranus are presented. Applicability of a low lift-to-drag entry vehicle aeroshell with flight heritage for Martian entries is considered in the context of an end-to-end Uranus orbiter and probe mission design using aerocapture for orbit insertion. A feasible trajectory space for aerocapture orbit insertion based on various launch opportunities and interplanetary trajectories yield atmosphere-relative entry velocities between 22 km/s and 31 km/s with maximum freestream densities on the order of 10−5 kg/m3. A preliminary study of the aerothermal environments for the entire trajectory space is discussed followed by a more detailed study of the nominal design trajectories. Design conditions are based on laminar, thermochemical nonequilibrium Navier-Stokes flowfield simulations with coupled radiation transport for full lift-up (deep) and lift-down (shallow) design trajectories. Convective heating is found to be the dominant heating mode, with radiation accounting for only 1 % to 5 % of the total heat flux for the design trajectories. The peak design conditions are found to be 422 W/cm2 for heat flux, 10.5 kPa for pressure, 149 Pa for shear stress, and 88.1 kJ/cm2 for the total heat load.

Uranus↗

Aeroheating Environments of Aerocapture Systems for Uranus Orbiters

Aeroheating environments for an aerocapture system enabling flagship-class science at Uranus are presented. Applicability of a low lift-to-drag entry vehicle aeroshell with flight heritage for Martian entries is considered in the context of an end-to-end Uranus orbiter and probe mission design using aerocapture for orbit insertion. A feasible trajectory space for aerocapture orbit insertion based on various launch opportunities and interplanetary trajectories yield atmosphere-relative entry velocities between 22 km/s and 31 km/s with maximum freestream densities on the order of 10 −5 kg/m 3 . A preliminary study of the aerothermal environments for the entire trajectory space is discussed followed by a more detailed study of the nominal design trajectories. Design conditions are based on laminar, thermochemical nonequilibrium Navier-Stokes flowfield simulations with coupled radiation transport for full lift-up (deep) and lift-down (shallow) design trajectories. Convective heating is found to be the dominant heating mode, with radiation accounting for only 1 % to 5 % of the total heat flux for the design trajectories. The peak design conditions are found to be 422 W/cm 2 for heat flux, 10.5 kPa for pressure, 149 Pa for shear stress, and 88.1 kJ/cm 2 for the total heat load.

Aerocapture↗

Commentary on the 1985 NASA/Vanderbilt Symposium on Future Hypervelocity Flight Requirements

The discussion, started in Semiannual Status Report Number 1, on aerothermal problems of hypervelocity flight and experiments that may lead to significant improvements in analytical/computational predictive methods, continues. The commentary is based on presentations made by speakers at a symposium on this subject held in December 1985. Symposium participants focused on the serious deficiencies that exist in knowledge of real-gas, nonequilibrium thermochemical-kinetic processes, catalytic processes, surface and shock slip, gas/surface interaction, boundary layer transition, and vortical leeside flows under hypervelocity conditions. Programs of laboratory research and computations leading toward in-flight experiments were recommended. Feasibility of appropriate measurement techniques for the flight environment was assessed and problems for study in that area identified. A synopsis of the oral presentations is given.

Potter, J. L.↗

A parametric study of dissociation and ionization models at 12 km/sec

Thermochemical nonequilibrium-solution dependence on available models for the chemical reaction rates is examined. Solutions from the Kang and Dunn (1973) reaction-rate set, the Park rate set of 1987, and the Park rate set of 1991 are compared. The blunt-nosed, axisymmetric geometry considered is a 60-deg sphere cone with nose radius of 1.07 m and cicular aft skirt. The nonequilibrium test case is 12 km/sec entry into the earth's atmosphere at 80 km altitude. The model variations are implemented into the Langley aerothermodynamics upwind relaxation algorithm code. While variations in the reaction rates have no effect on the surface pressure distribution and little effect on the convective heating, the effect on degree of ionization and radiative heating can be a factor of three.

Mitcheltree, R. A.↗

Effects of Nonequilibrium at Edge of Boundary Layer on Convective Heat Transfer to a Blunt Body

This investigation is a continuation of a previous study on nonequilibrium convective heat transfer to a blunt body. In the previous study, for relatively high Reynolds number flows, it was found that: nonequilibrium convective heat transfer to a blunt body is not strongly dependent on freestream parameters, provided that the thermochemical equilibrium is reached at the edge of boundary layer; and successful testing of convective heat transfer in an arc-jet environment is possible by duplicating the surface pressure and total enthalpy. The nonequilibrium convective heat transfer computations are validated against the results of Fay and Riddell/Goulard theory. Present work investigates low Reynolds number conditions which are typical in an actual arc-jet flow environment. One expects that there will be departures from the Fay and Riddell/Goulard result since certain assumptions of the classical theory are not satisfied. These departures are of interest because the Fay and Riddell/Goulard formulas are extensively used in arc-jet testing (e.g., to determine the enthalpy of the flow and the catalytic efficiency of heat shield materials). For practical sizes of test materials, density of the test flow (and Reynolds number) in an arc-jet is such that thermochemical equilibrium may not be reached at the edge of boundary layer. For blunt body flows of nitrogen and air, computations will be presented to show the effects of thermochemical nonequilibrium at the boundary layer edge on nonequilibrium heat transfer.

Goekcen, Tahir↗

Modeling of Laser Vaporization and Plume Chemistry in a Boron Nitride Nanotube Production Rig

Flow in a pressurized, vapor condensation (PVC) boron nitride nanotube (BNNT) production rig is modeled. A laser provides a thermal energy source to the tip of a boron ber bundle in a high pressure nitrogen chamber causing a plume of boron-rich gas to rise. The buoyancy driven flow is modeled as a mixture of thermally perfect gases (B, B2, N, N2, BN) in either thermochemical equilibrium or chemical nonequilibrium assuming steady-state melt and vaporization from a 1 mm radius spot at the axis of an axisymmetric chamber. The simulation is intended to define the macroscopic thermochemical environment from which boron-rich species, including nanotubes, condense out of the plume. Simulations indicate a high temperature environment (T > 4400K) for elevated pressures within 1 mm of the surface sufficient to dissociate molecular nitrogen and form BN at the base of the plume. Modifications to Program LAURA, a finite-volume based solver for hypersonic flows including coupled radiation and ablation, are described to enable this simulation. Simulations indicate that high pressure synthesis conditions enable formation of BN vapor in the plume that may serve to enhance formation of exceptionally long nanotubes in the PVC process.

Gnoffo, Peter A.↗

Computation of nonequilibrium viscous flows in arc-jet wind tunnel nozzles

The axisymmetric Navier-Stokes equations are numerically solved for air flows in converging-diverging conical nozzles. A nonequilibrium formulation with a two-temperature thermochemical model has been employed for vibrationally excited, partially dissociated and ionized air. Evaluations of the formulation is made through comparisons of the computed flow quantities with the experimental data obtained in expanding arc-jet nozzle flows. Two sets of experimental data obtained in conical nozzles with geometric exit-to-throat area ratios of 625 and 8100 cover reservoir pressures ranging from 5 to 95 atm and enthalpies from 5 to 10 MJ/kg. The computed results are in good agreement with local mass flux, Pitot pressure, and static wall pressure measurements at the nozzle exit. For low pressure reservoir conditions, comparisons of computations with the experimental data indicate that the wall catalycity is very important in predicting species mole fraction at the nozzle exit since gas phase chemistry is frozen. The current computations overpredict the atomic oxygen concentrations which in turn results in underprediction of the vibrational temperature.

Gokcen, Tahir↗

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for high speed lunar return experiments conducted in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The nonequilibrium radiative transport equations are solved via NASA's in-house radiation code NEQAIR using flow field input from US3D flow solver. The post-shock flow properties for the 10 km/s Earth entry conditions are computed using the stagnation line of a blunt-body and a full facility CFD (Computational Fluid Dynamics) simulation of the EAST shock tube. The shocked gas in the blunt-body flow achieves a thermochemical equilibrium away from the shock front whereas EAST flow exhibits a nonequilibrium behavior due to strong viscous dissipation of the shock by boundary layer. The full-tube flow calculations capture the influence of the boundary layer on the shocked gas state and provide a realistic fluid dynamic input for the radiative predictions. The integrated radiance behind the shock is calculated in NEQAIR for wavelength regimes from Vacuum-UltraViolet (VUV) to InfraRed (IR), which are pertinent to the emission characteristics of high enthalpy shock waves in air. These radiance profiles are validated against corresponding EAST shots. The full-tube simulations successfully predict a sharp radiance peak at the shock front which gets smeared in the test data due to the spatial resolution in the measurements. The full facility based radiance behind the shock shows a slightly better match with the test data in the VUV and Red spectral regions, as compared to that from a blunt-body based predictions. The UV radiance is very similar for both geometries and under-predicts the test behavior. The IR test data matches better with the blunt-body based predictions where the full-tube simulations show a significant over-prediction.

Radiation↗

Computational procedure for evaluation of space shuttle TPS requirements

The rudiments of a computational procedure are presented for predicting heat transfer rates to the windward pitchplane of a shuttle vehicle under conditions of nonequilibrium flow. The procedure involves computations in three regions of the flow field, namely, the stagnation region, inviscid flow region, and boundary layer region. The thermochemical state in the inviscid flow was calculated with the aid of correlations developed from shadowgraphs of shuttle configurations. The flow in the boundary layer and stagnation regions was computed with a nonequilibrium viscous layer code. For boundary layers, the outer edge thermochemical state was determined by matching the mass flow in the boundary layer with the mass flow represented by each inviscid streamline. Typical solutions are presented to show the effects of surface pressure, catalycity, and homogeneous chemical reactions on flow field variables and surface heat transfer rates.

Tong, H.↗