Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonequilibrium Systems”

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

The effect of fabrication on corrosion in aluminum 2195: Environmental and microstructural considerations

Aluminum alloys containing lithium are particularly attractive to the aerospace structural designer. Lithium's density is only 0.53 g/cc, thus an addition of one weight percent lithium not only increases yield strength, but decreases the density by almost three percent while increasing the modulus by over six percent. The fact that lithium improves these physical properties simultaneously has led to intense study and development of the alloy system. Heretofore, problems in large scale alloy production have retarded commercial development. During the last fifteen years, advances in production technology have rekindled interest in Al-Li alloys, and aluminum suppliers have developed many candidate aerospace materials. However, if these alloys are to be employed successfully, a more complete understanding of their nonequilibrium metallurgy is required. Peel and Starke have each pointed out that an understanding of the weldability of these alloys is a critical step in their implementation. This study addresses the critical lack of information on the environmental compatibility of welded Al 2195 components. Corrosion data for these systems is incomplete, particularly for welded materials exposed to sea water or sea water condensate.

Walsh, Daniel W.↗

Boundary Layer Stability and Laminar-Turbulent Transition Analysis with Thermochemical Nonequilibrium Applied to Martian Atmospheric Entry

As Martian atmospheric entry vehicles increase in size to accommodate larger payloads, transitional ow may need to be taken into account in the design of the heat shield in order to reduce heat shield mass. The mass of the Thermal Protection System (TPS) comprises a significant portion of the vehicle mass, and a reduction of this mass would result in fuel savings. The current techniques used to design entry shields generally assume fully turbulent flow when the vehicle is large enough to expect transitional flow, and while this worst-case scenario provides a greater factor of safety it may also result in overdesigned TPS and unnecessarily high vehicle mass. Greater accuracy in the prediction of transition would also reduce uncertainty in the thermal and aerodynamic loads. Stability analysis, using e(sup⁡ N) -based methods including Linear Stability Theory (LST) and the Parabolized Stability Equations (PSE), offers a physics-based method of transition prediction that has been thoroughly studied and applied in perfect gas flows, and to a more limited extent in reacting and nonequilibrium flows. These methods predict the amplification of a known disturbance frequency and allow identification of the most unstable frequency. Transition is predicted to occur at a critical amplification or N Factor, frequently determined through experiment and empirical correlations. The LAngley Stability and TRansition Analysis Code (LASTRAC), with modifications for thermochemically reacting flows and arbitrary gas mixtures, will be presented with LST results on a simulation of a high enthalpy CO2 gas wind tunnel test relevant to Martian atmospheric entry. The results indicate transition caused by modified Tollmien-Schlichting waves on the leeward side, which are predicted to be more stable and cause transition slightly downstream when thermochemical nonequilibrium is included in the stability analysis for the same mean flow solution.

Transition↗

Chemical nonequilibrium Navier-Stokes solutions for hypersonic flow over an ablating graphite nosetip

The general boundary conditions including mass and energy balances of chemically equilibrated or nonequilibrated gas adjacent to ablating surfaces have been derived. A computer procedure based on these conditions was developed and interfaced with the Navier-Stokes solver for predictions of the flow field, surface temperature, and surface ablation rates over re-entry space vehicles with ablating Thermal Protection Systems (TPS). The Navier-Stokes solver with general surface thermochemistry boundary conditions can predict more realistic solutions and provide useful information for the design of TPS. A test case with a proposed hypersonic test vehicle configuration and associated free stream conditions was developed. Solutions with various surface boundary conditions were obtained, and the effect of nonequilibrium gas as well as surface chemistry on surface heating and ablation rate were examined. The solutions of the GASP code with complete ablating surface conditions were compared with those of the ASC code. The direction of future work is also discussed.

Chen, Y. K.↗

A comparison of two algorithms for simulating collisionless systems

Two completely different simulation algorithms are compared by applying them to the same stellar dynamical problems: one is a self-consistent field (SCF) method for solving Poisson's equation and the other is a phase-space method for integrating the collisionless Boltzmann equation. We consider simulations of spherical stellar systems which are initially far from equilibrium and relax to their final states by gravitational collapse. The initial conditions consist of either uniform-density spheres or nonequilibrium models having Plummer density profiles, in which velocity dispersions are assigned according to given virial ratios. If a few tens of radial expansion terms with hundreds of thousands of particles are used in the SCF code, excellent agreement is found between the results it generates and those obtained with the phase-space solver, provided that a sufficiently large number of grid cells are employed with the latter. These findings imply that for simulating collisionless systmes over many dynamical times, the SCF approach based on sampling phase space is competitive with the approach treating phase space as a continuous fluid. The results of our tests make it possible to estimate the number of particles and basis functions required in situations like those modeled. Limitations of the SCF method and the choice of an optimal set of basis functions are also discussed.

Hozumi, Shunsuke↗

Overview of the Dragonfly Entry Aerosciences Measurements (DrEAM) Suite

NASA Ames and Langley arepartnering with DLR to propose acomprehensive instrumentation suite known as the Dragonfly Entry Aerosciences Meas-urements (DrEAM). DrEAM being the first competed mission to fly EDL instrumentation as part of NASA’s Engineering Science Investigation (ESI).DrEAM will provide key aerothermodynamicdata and performance analysisfor Dragonfly’s forebody and backshell therma lprotection system (TPS),and also includes a DLR-provided Data Acquisition System(DAS).Titan’s atmosphere predominantly consists of nitro-gen (~98% by mole) with small amounts of methane(~2% by mole)and other trace gases. CN is a strong ra-diator and is found in nonequilibriumconcentrationsfor Titan entry. The accurate modeling of nonequilibrium CN radiation has proven to be a difficult task. Prompted by the Huygensmission, many experimental campaigns and analyses were performed to better understand the aerothermal environments experience by the probe dur-ing Titan entry[1].However, the Huygens probe carried no heatshieldinstrumentation. Therefore,the DrEAM instrumentation suite will significantly advance the state-of-the-art not only by documenting theenviron-ment and performance of Dragonfly’sentry system but also by making keymeasurements in Titan’s atmos-phere for thefirst time, thus providing new benchmark dataapplicable to entry science more generally. Current Measurement Goals: Aerothermal envi-ronments and TPS responsewill be measured using sen-sors similarto the Mars Entry, Descent, and Landing In-strumentation2 (MEDLI2) Integrated Sensor Plug(MISP) and the COMbined Aerothermaland Radiome-ter Sensor (COMARS)suite[2], with the latter supplied by DLR.For MEDLI2, MISP usedembedded thermo-couples (TCs) todirectly measure in-depth temperature of theTPS at several locations,which can also be used to infer surface environmentsvia inverse analysis. For DrEAM, the MISPstyle plugswill be known as Drag-onfly Sensors for Aero-Thermal Reconstruction (Drag-STR)plugs. On Schiaparelli, the COMARSsuite in-cludedthree total surface-mounted heatflux sensors, three pressure sensors, and one radiometer. For DrEAM, the COMARS package will be known as COmbined Sensor System for Titan Atmosphere (COSSTA). Atmospheric density measurements and capsuleaerodynamic data will be obtained throughthe onboard Inertial Measurement Unit (IMU),supple-mented by pressure transducers similar tothose used by the MEDLI Mars Entry AtmosphericData System (MEADS).The DrEAM pressure sensors will be known as Dragonfly Atmospheric Flight Transducers(DrAFT) Both DragSTR and DrAFThave flight heritagefrom MISP and MEADS on the MSLand Mars 2020mis-sions, and the COMARS suite successfullyflew on the ESA Schiaparelli EDMlander. A preliminary layout of the sensors is shown in Fig. 1. Because Dragonflyuses the same aeroshell provider (i.e., LockheedMartin) and materials for the TPS,with what are expected to be sim-ilar thicknesses as MSL and Mars 2020on both the heat shield and backshell, theDrEAMinstrumentation will look to utilize the same techniques and processes as de-veloped byMEDLI and MEDLI2for vehicle integra-tion. This commonality alsoenables DrEAMto lever-age the extensiveground test qualifications performed forMEDLI and MEDLI2 and claim substantial heritagefor this system.

A Brandis↗

Conservation equations and physical models for hypersonic air flows over the aeroassist flight experiment vehicle

The code development and application program for the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA), with emphasis directed toward support of the Aeroassist Flight Experiment (AFE) in the near term and Aeroassisted Space Transfer Vehicle (ASTV) design in the long term is reviewed. LAURA is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3-D, viscous, hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over an AFE wind tunnel model. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium over the full scale AFE configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗

Containerless processing of Nb-Ge alloys in a long drop tube

The thirty-two meter drop tube at the Marshall Space Flight Center was used to study the effect of zero gravity containerless processing on the structure and properties of materials. The concept involves the suppression of heterogeneous nucleation of solid in liquid and, therefore, solidification accompanied by large degrees of undercooling. Under these conditions metastable phases can be formed or, at the very least, unique nonequilibrium microstructures (containing equilibrium phases) with unique properties can be produced. The drop tube solidification was applied to niobium base alloys with emphasis on the Nb-Ge binary system in an effort to produce metastable phases with high superconducting transition temperatures in bulk specimens. In the past, only lower Ge alloys (Nb-13 a/o, Nb-18 a/o, and Nb-22 a/o) could be undercooled. Higher Ge alloys (e.g., Nb-25 a/o Ge and Nb-27 a/o Ge) can now be undercooled on a routine basis.

Bayuzick, R. J.↗

The lower main sequence and the nature of secondary stars in ultracompact binaries

The possible nature of the secondary stars in ultracompact binary stellar systems (with orbital periods of less than about 1 hr) are systematically investigated. Using a simplified stellar evolution code, which assumes isentropic stellar models, nearly 3000 separate models for hydrogen-burning main-sequence stars with masses less than 0.3 solar mass are generated. The effects of the (homogeneous) chemical composition on the properties of such stars and in, particular, on the minimum main-sequence mass are explored in detail. It is found that this minimum mass is a sensitive function of the hydrogen content and can be as small as 0.035 solar mass. The properties of fully degenerate low-mass stars, as well as low-mass stars that are not in thermal equilibrium are also studied. In particular, it is shown that the thermal time scale of such nonequilibrium stars may exceed 10 billion yr. The results also confirm the existence of a second branch of the main sequence and shed new light on the thermal instability of this branch. The available observational information on the three known ultracompact binary systems (4U 1626-67, G61-29, 4U 1916-05) is summarized, and the results of the stellar model calculations are combined with the empirical results to place constraints on the properties of the secondary stars.

Rappaport, S.↗

Upwind-biased, point-implicit relaxation strategies for hypersonic flowfield simulations on supercomputers

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, 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 Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The relaxation strategy is well suited for computers employing either vector or parallel architectures, and the relation between computer architecture and algorithm is emphasized. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large. block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over a blunt body. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km.) over the Aeroassist Flight Experiment (AFE) configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗

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.↗

Prospects for Nuclear Electric Propulsion Using Closed-Cycle Magnetohydrodynamic Energy Conversion

Nuclear electric propulsion (NEP) has long been recognized as a major enabling technology for scientific and human exploration of the solar system, and it may conceivably form the basis of a cost-effective space transportation system suitable for space commerce. The chief technical obstacles to realizing this vision are the development of efficient, high-power (megawatt-class) electric thrusters and the development of low specific mass (less than 1 kg/kWe) power plants. Furthermore, comprehensive system analyses of multimegawatt class NEP systems are needed in order to critically assess mission capability and cost attributes. This Technical Publication addresses some of these concerns through a systematic examination of multimegawatt space power installations in which a gas-cooled nuclear reactor is used to drive a magnetohydrodynamic (MHD) generator in a closed-loop Brayton cycle. The primary motivation for considering MHD energy conversion is the ability to transfer energy out of a gas that is simply too hot for contact with any solid material. This has several intrinsic advantages including the ability to achieve high thermal efficiency and power density and the ability to reject heat at elevated temperatures. These attributes lead to a reduction in system specific mass below that obtainable with turbine-based systems, which have definite solid temperature limits for reliable operation. Here, the results of a thermodynamic cycle analysis are placed in context with a preliminary system analysis in order to converge on a design space that optimizes performance while remaining clearly within established bounds of engineering feasibility. MHD technology issues are discussed including the conceptual design of a nonequilibrium disk generator and opportunities for exploiting neutron-induced ionization mechanisms as a means of increasing electrical conductivity and enhancing performance and reliability. The results are then used to make a cursory examination of piloted Mars missions during the 2018 opportunity.

Litchford, R. J.↗

A hybrid multigrid technique for computing steady-state solutions to supersonic flows

Recently, Li and Sanders have introduced a class of finite difference schemes to approximate generally discontinuous solutions to hyperbolic systems of conservation laws. These equations have the form together with relevant boundary conditions. When modelling hypersonic spacecraft reentry, the differential equations above are frequently given by the compressible Euler equations coupled with a nonequilibrium chemistry model. For these applications, steady state solutions are often sought. Many tens (to hundreds) of super computer hours can be devoted to a single three space dimensional simulation. The primary difficulty is the inability to rapidly and reliably capture the steady state. In these notes, we demonstrate that a particular variant from the schemes presented can be combined with a particular multigrid approach to capture steady state solutions to the compressible Euler equations in one space dimension. We show that the rate of convergence to steady state coming from this multigrid implementation is vastly superior to the traditional approach of artificial time relaxation. Moreover, we demonstrate virtual grid independence. That is, the rate of convergence does not depend on the degree of spatial grid refinement.

Sanders, Richard↗

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗

951 Gaspra - A pre-Galileo estimate of its surface evolution

An effort is made to quantitatively model the collisional environment and surface evolution of asteroid 951 Gaspra prior to the flyby of the Galileo spacecraft later this year. Gaspra's average intrinsic collisional probability is calculated to be 5.4 x 10 to the -18th/sq km per year with a root mean square impact velocity of 5.0 km/sec. The projectile population and impact strength for Gaspra are highly uncertain and a range of values is considered. Based on these models, Gaspra is unlikely to be a primordial object surviving from the beginning of the solar system. It is most likely a second or multiple generation fragment forming about 100 million yr ago. For a low impact flux, the crater size distribution directly reveals the size distribution of the impacting population. For a nonequilibrium surface, model-dependent isochrons for estimating the age of the surface are given. Lower limits are set by the model for Gaspra's age (about 1 billion yr) and flux for the case of a high impact flux producing an equilibrium surface.

Namiki, Noriyuki↗

NASA Ames's electric arc-driven shock tube facility and research on nonequilibrium phenomena in low density hypersonic flows

Basic requirements for a ground test facility simulating low density hypersonic flows are discussed. Such facilities should be able to produce shock velocities in the range of 10-17 km/sec in an initial pressure of 0.010 to 0.050 Torr. The facility should be equipped with diagnostics systems to be able to measure the emitted radiation, characteristic temperatures and populations in various energy levels. In the light of these requirements, NASA Ames's electric arc-driven low density shock tube facility is described and available experimental diagnostics systems and computational tools are discussed.

Sharma, Surendra P.↗

Investigation of thermal protection systems effects on viscid and inviscid flow fields for manned entry systems

Procedures and methods for predicting aerothermodynamic heating to delta orbiter shuttle vehicles were reviewed. A number of approximate methods were found to be adequate for large scale parameter studies, but are considered inadequate for final design calculations. It is recommended that final design calculations be based on a computer code which accounts for nonequilibrium chemistry, streamline spreading, entropy swallowing, and turbulence. It is further recommended that this code be developed with the intent that it can be directly coupled with an exact inviscid flow field calculation when the latter becomes available. A nonsimilar, equilibrium chemistry computer code (BLIMP) was used to evaluate the effects of entropy swallowing, turbulence, and various three dimensional approximations. These solutions were compared with available wind tunnel data. It was found study that, for wind tunnel conditions, the effect of entropy swallowing and three dimensionality are small for laminar boundary layers but entropy swallowing causes a significant increase in turbulent heat transfer. However, it is noted that even small effects (say, 10-20%) may be important for the shuttle reusability concept.

Bartlett, E. P.↗

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗

Transient Ablation of Teflon Hemispheres

For high-speed entry of space vehicles into atmospheric environments, ablation is a practical method for alleviating severe aerodynamic heating. Several studies have been undertaken on steady or quasi-steady ablation. However, ablation is a very complicated phenomenon in which a nonequilibrium chemical process is associated with an aerodynamic process that involves changes in body shape with time. Therefore, it seems realistic to consider that ablation is an unsteady phenomenon. In the design of an ablative heat-shield system, since the ultimate purpose of the heat shield is to keep the internal temperature of the space vehicle at a safe level during entry, the transient heat conduction characteristics of the ablator may be critical in the selection of the material and its thickness. This note presents an experimental study of transient ablation of Teflon, with particular emphasis on the change in body shape, the instantaneous internal temperature distribution, and the effect of thermal expansion on ablation rate.

Arai, Norio↗