Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERNAL ENERGY”

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 235 records · Page 13

Numerical Modelling of the Deep Impact Mission Experiment

NASA s Deep Impact Mission (launched January 2005) will provide, for the first time ever, insights into the interior of a comet (Tempel 1) by shooting a approx.370 kg projectile onto the surface of a comets nucleus. Although it is usually assumed that comets consist of a very porous mixture of water ice and rock, little is known about the internal structure and in particular the constitutive material properties of a comet. It is therefore difficult to predict the dimensions of the excavated crater. Estimates of the crater size are based on laboratory experiments of impacts into various target compositions of different densities and porosities using appropriate scaling laws; they range between 10 s of meters up to ~250 m in diameter [1]. The size of the crater depends mainly on the physical process(es) that govern formation: Smaller sizes are expected if (1) strength, rather than gravity, limits crater growth; and, perhaps even more crucially, if (2) internal energy losses by pore-space collapse reduce the coupling efficiency (compaction craters). To investigate the effect of pore space collapse and strength of the target we conducted a suite of numerical experiments and implemented a novel approach for modeling porosity and the compaction of pores in hydrocode calculations.

Wuennemann, K.↗

Advanced Electrical Materials and Components Development: An Update

The primary means to develop advanced electrical components is to develop new and improved materials for magnetic components (transformers, inductors, etc.), capacitors, and semiconductor switches and diodes. This paper will give an update of the Advanced Power Electronics and Components Technology being developed by the NASA Glenn Research Center for use in future Power Management and Distribution subsystems used in space power systems for spacecraft and lunar and planetary surface power. The initial description and status of this technology program was presented two years ago at the First International Energy Conversion Engineering Conference held at Portsmouth, Virginia, August 2003. The present paper will give a brief background of the previous work reported and a summary of research performed the past several years on soft magnetic materials characterization, dielectric materials and capacitor developments, high quality silicon carbide atomically smooth substrates, and SiC static and dynamic device characterization under elevated temperature conditions. The rationale for and the benefits of developing advanced electrical materials and components for the PMAD subsystem and also for the total power system will also be briefly discussed.

Schwarze, Gene E.↗

Mestastable State Population in Laser Induced Plasmas

Laser induced plasma has been used as a source of neutrals and ions in the study of astrophysical plasmas. The purity of state of this source is essential in the determination of collision parameters such as the charge transfer rate coefficients between ions and neutrals. We will show that the temperature of the laser induced plasma is a rapidly decreasing function of time. The temperature is initially high but cools off rapidly through collisions with the expanding plasma electrons as the plasma recombines and streams into the vacuum. This rapid expansion of the plasma, similar to a supersonic jet, drastically lowers the internal energy of the neutrals and ions.

Kwong, V. H. S.↗

Microfabrication of a Segmented-Involute-Foil Regenerator, Testing in a Sunpower Stirling Convertor and Supporting Modeling and Analysis

Under Phase II of a NASA Research Award contract, a prototype nickel segmented-involute-foil regenerator was microfabricated via LiGA and tested in the NASA/Sunpower oscillating-flow test rig. The resulting figure-of-merit was about twice that of the approx.90% porosity random-fiber material currently used in the small 50-100 W Stirling engines recently manufactured for NASA. That work was reported at the 2007 International Energy Conversion Engineering Conference in St. Louis, was also published as a NASA report, NASA/TM-2007-2149731, and has been more completely described in a recent NASA Contractor Report, NASA/CR-2007-2150062. Under a scaled-back version of the original Phase III plan, a new nickel segmentedinvolute- foil regenerator was microfabricated and has been tested in a Sunpower Frequency-Test-Bed (FTB) Stirling convertor. Testing in the FTB convertor produced about the same efficiency as testing with the original random-fiber regenerator. But the high thermal conductivity of the prototype nickel regenerator was responsible for a significant performance degradation. An efficiency improvement (by a 1.04 factor, according to computer predictions) could have been achieved if the regenerator been made from a low-conductivity material. Also the FTB convertor was not reoptimized to take full advantage of the microfabricated regenerator's low flow resistance; thus the efficiency would likely have been even higher had the FTB been completely reoptimized. This report discusses the regenerator microfabrication process, testing of the regenerator in the Stirling FTB convertor, and the supporting analysis. Results of the pre-test computational fluid dynamics (CFD) modeling of the effects of the regenerator-test-configuration diffusers (located at each end of the regenerator) is included. The report also includes recommendations for accomplishing further development of involute-foil regenerators from a higher-temperature material than nickel.

Ibrahim, Mounir B.↗

Cleaning Carbon Nanotubes by Use of Mild Oxygen Plasmas

Experiments have shown that it is feasible to use oxygen radicals (specifically, monatomic oxygen) from mild oxygen plasmas to remove organic contaminants and chemical fabrication residues from the surfaces of carbon nanotubes (CNTs) and metal/CNT interfaces. A capability for such cleaning is essential to the manufacture of reproducible CNT-based electronic devices. The use of oxygen radicals to clean surfaces of other materials is fairly well established. However, previously, cleaning of CNTs and of graphite by use of oxygen plasmas had not been attempted because both of these forms of carbon were known to be vulnerable to destruction by oxygen plasmas. The key to success of the present technique is, apparently, to ensure that the plasma is mild . that is to say, that the kinetic and internal energies of the oxygen radicals in the plasma are as low as possible. The plasma oxygen-radical source used in the experiments was a commercial one marketed for use in removing hydrocarbons and other organic contaminants from vacuum systems and from electron microscopes and other objects placed inside vacuum systems. In use, the source is installed in a vacuum system and air is leaked into the system at such a rate as to maintain a background pressure of .0.56 torr (.75 Pa). In the source, oxygen from the air is decomposed into monatomic oxygen by radio-frequency excitation of a resonance of the O2 molecule (N2 is not affected). Hence, what is produced is a mild (non-energetic) oxygen plasma. The oxygen radicals are transported along with the air molecules in the flow created by the vacuum pump. In the experiments, exposure to the oxygen plasma in this system was shown to remove organic contaminants and chemical fabrication residues from several specimens. Many high-magnification scanning electron microscope (SEM) images of CNTs were taken before and after exposure to the oxygen plasma. As in the example shown in the figure, none of these images showed evidence of degradation of CNT structures.

Petkov, Mihail↗

Thermodynamic Properties of Nitrogen Including Liquid and Vapor Phases from 63K to 2000K with Pressures to 10,000 Bar

Tables of thermodynamic properties of nitrogen are presented for the liquid and vapor phases for temperatures from the freezing line to 2000K and pressures to 10,000 bar. The tables include values of density, internal energy, enthalpy, entropy, isochoric heat capacity, isobaric heat capacity velocity of sound, the isotherm derivative, and the isochor derivative. The thermodynamic property tables are based on an equation of state, P=P (p,T), which accurately represents liquid and gaseous nitrogen for the range of pressures and temperatures covered by the tables. Comparisons of property values calculated from the equation of state with measured values for P-p-T, heat capacity, enthalpy, latent heat, and velocity of sound are included to illustrate the agreement between the experimental data and the tables of properties presented here. The coefficients of the equation of state were determined by a weighted least squares fit to selected P-p-T data and, simultaneously, to isochoric heat capacity data determined by corresponding states analysis from oxygen data, and to data which define the phase equilibrium criteria for the saturated liquid and the saturated vapor. The vapor pressure equation, melting curve equation, and an equation to represent the ideal gas heat capacity are also presented. Estimates of the accuracy of the equation of state, the vapor pressure equation, and the ideal gas heat capacity equation are given. The equation of state, derivatives of the equation, and the integral functions for calculating derived thermodynamic properties are included.

Thermodynamic properties↗

Analysis of Photospheric Convection Cells with SDO/HMI

Supergranulation is a component of solar convection that assists in the outward transportation of internal energy. Supergranule cells are approximately 35 Mm across, have lifetimes on the order of a day and have divergent horizontal velocities of around 300 m/s, a factor of 10 higher than their central radial components. While they have been observed using Doppler methods for around half a century, their existence is also observed in other datasets such as magnetograms and Ca II K images. These datasets clearly show the influence of supergranulation on solar magnetism and how the local field is organized by the flows of supergranule cells. The Heliospheric and Magnetic Imager (HMI) aboard SDO is making fresh observations of convection phenomena at a higher cadence and a higher resolution that should make granular features visible. Granulation and supergranulation characteristics can now be compared within the same datasets, which may lead to further understanding of any mutual influences. The temporal and spatial enhancements of HMI will also reduce the noise level within studies of convection so that more detailed studies of their characteristics may be made. We present analyses of SDO/HMI Dopplergrams that provide new estimates of convection cell sizes, lifetimes, and velocity flows, as well as the rotation rates of the convection patterns across the solar disk. We make comparisons with previous data produced by MDI, as well as from data simulations.

Williams, Peter E.↗

Updated Chemical Kinetics and Sensitivity Analysis Code

An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗

Modeling, Measurements, and Fundamental Database Development for Nonequilibrium Hypersonic Aerothermodynamics

The design of entry vehicles requires predictions of aerothermal environment during the hypersonic phase of their flight trajectories. These predictions are made using computational fluid dynamics (CFD) codes that often rely on physics and chemistry models of nonequilibrium processes. The primary processes of interest are gas phase chemistry, internal energy relaxation, electronic excitation, nonequilibrium emission and absorption of radiation, and gas-surface interaction leading to surface recession and catalytic recombination. NASAs Hypersonics Project is advancing the state-of-the-art in modeling of nonequilibrium phenomena by making detailed spectroscopic measurements in shock tube and arcjets, using ab-initio quantum mechanical techniques develop fundamental chemistry and spectroscopic databases, making fundamental measurements of finite-rate gas surface interactions, implementing of detailed mechanisms in the state-of-the-art CFD codes, The development of new models is based on validation with relevant experiments. We will present the latest developments and a roadmap for the technical areas mentioned above

Bose, Deepak↗

Unresolved Problems by Shock Capturing: Taming the Overheating Problem

The overheating problem, first observed by von Neumann [1] and later studied extensively by Noh [2] using both Eulerian and Lagrangian formulations, remains to be one of the unsolved problems by shock capturing. It is historically well known to occur when a flow is under compression, such as when a shock wave hits and reflects from a wall or when two streams collides with each other. The overheating phenomenon is also found numerically in a smooth flow undergoing rarefaction created by two streams receding from each other. This is in contrary to one s intuition expecting a decrease in internal energy. The excessive amount in the temperature increase does not reduce by refining the mesh size or increasing the order of accuracy. This study finds that the overheating in the receding flow correlates with the entropy generation. By requiring entropy preservation, the overheating is eliminated and the solution is grid convergent. The shock-capturing scheme, as being practiced today, gives rise to the entropy generation, which in turn causes the overheating. This assertion stands up to the convergence test.

Liou, Meng-Sing↗

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↗

Microscopic Simulation and Macroscopic Modeling for Thermal and Chemical Non-Equilibrium

This paper deals with the accurate microscopic simulation and macroscopic modeling of extreme non-equilibrium phenomena, such as encountered during hypersonic entry into a planetary atmosphere. The state-to-state microscopic equations involving internal excitation, de-excitation, dissociation, and recombination of nitrogen molecules due to collisions with nitrogen atoms are solved time-accurately. Strategies to increase the numerical efficiency are discussed. The problem is then modeled using a few macroscopic variables. The model is based on reconstructions of the state distribution function using the maximum entropy principle. The internal energy space is subdivided into multiple groups in order to better describe the non-equilibrium gases. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients. The modeling is completely physics-based, and its accuracy depends only on the assumed expression of the state distribution function and the number of groups used. The model makes no assumption at the microscopic level, and all possible collisional and radiative processes are allowed. The model is applicable to both atoms and molecules and their ions. Several limiting cases are presented to show that the model recovers the classical twotemperature models if all states are in one group and the model reduces to the microscopic equations if each group contains only one state. Numerical examples and model validations are carried out for both the uniform and linear distributions. Results show that the original over nine thousand microscopic equations can be reduced to 2 macroscopic equations using 1 to 5 groups with excellent agreement. The computer time is decreased from 18 hours to less than 1 second.

Modeling↗

Calculation of State Specific Rate Coefficients for Non-Equilibrium Hypersonics Applications: from H(Psi) = E(Psi) to k(T) = A *exp(-E(sub a)/RT)

Development of High-Fidelity Physics-Based Models to describe hypersonic flight through the atmospheres of Earth and Mars is underway at NASA Ames Research Center. The goal is to construct chemistry models of the collisional and radiative processes that occur in the bow shock and boundary layers of spacecraft during atmospheric entry that are free of empiricism. In this talk I will discuss our philosophy and describe some of our progress. Topics to be covered include thermochemistry, internal energy relaxation, collisional dissociation and radiative emission and absorption. For this work we start by solving the Schrodinger equation to obtain accurate interaction potentials and radiative properties. Then we invoke classical mechanics to compute state-specific heavy particle collision cross sections and reaction rate coefficients. Finally, phenomenological rate coefficients and relaxation times are determined from master equation solutions.

aerothermodynamics↗

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↗

Object-Oriented/Data-Oriented Design of a Direct Simulation Monte Carlo Algorithm

Over the past decade, there has been much progress towards improved phenomenological modeling and algorithmic updates for the direct simulation Monte Carlo (DSMC) method, which provides a probabilistic physical simulation of gas Rows. These improvements have largely been based on the work of the originator of the DSMC method, Graeme Bird. Of primary importance are improved chemistry, internal energy, and physics modeling and a reduction in time to solution. These allow for an expanded range of possible solutions In altitude and velocity space. NASA's current production code, the DSMC Analysis Code (DAC), is well-established and based on Bird's 1994 algorithms written in Fortran 77 and has proven difficult to upgrade. A new DSMC code is being developed in the C++ programming language using object-oriented and data-oriented design paradigms to facilitate the inclusion of the recent improvements and future development activities. The development efforts on the new code, the Multiphysics Algorithm with Particles (MAP), are described, and performance comparisons are made with DAC.

Liechty, Derek S.↗

Extension of the Viscous Collision Limiting Direct Simulation Monte Carlo Technique to Multiple Species

There are many flows fields that span a wide range of length scales where regions of both rarefied and continuum flow exist and neither direct simulation Monte Carlo (DSMC) nor computational fluid dynamics (CFD) provide the appropriate solution everywhere. Recently, a new viscous collision limited (VCL) DSMC technique was proposed to incorporate effects of physical diffusion into collision limiter calculations to make the low Knudsen number regime normally limited to CFD more tractable for an all-particle technique. This original work had been derived for a single species gas. The current work extends the VCL-DSMC technique to gases with multiple species. Similar derivations were performed to equate numerical and physical transport coefficients. However, a more rigorous treatment of determining the mixture viscosity is applied. In the original work, consideration was given to internal energy non-equilibrium, and this is also extended in the current work to chemical non-equilibrium.

Liechty, Derek S.↗

Development of the US3D Code for Advanced Compressible and Reacting Flow Simulations

Aerothermodynamics and hypersonic flows involve complex multi-disciplinary physics, including finite-rate gas-phase kinetics, finite-rate internal energy relaxation, gas-surface interactions with finite-rate oxidation and sublimation, transition to turbulence, large-scale unsteadiness, shock-boundary layer interactions, fluid-structure interactions, and thermal protection system ablation and thermal response. Many of the flows have a large range of length and time scales, requiring large computational grids, implicit time integration, and large solution run times. The University of Minnesota NASA US3D code was designed for the simulation of these complex, highly-coupled flows. It has many of the features of the well-established DPLR code, but uses unstructured grids and has many advanced numerical capabilities and physical models for multi-physics problems. The main capabilities of the code are described, the physical modeling approaches are discussed, the different types of numerical flux functions and time integration approaches are outlined, and the parallelization strategy is overviewed. Comparisons between US3D and the NASA DPLR code are presented, and several advanced simulations are presented to illustrate some of novel features of the code.

CFD↗

Multilayer Insulation Ascent Venting Model

The thermal and venting transient experienced by tank-applied multilayer insulation (MLI) in the Earth-to-orbit environment is very dynamic and not well characterized. This new predictive code is a first principles-based engineering model which tracks the time history of the mass and temperature (internal energy) of the gas in each MLI layer. A continuum-based model is used for early portions of the trajectory while a kinetic theory-based model is used for the later portions of the trajectory, and the models are blended based on a reference mean free path. This new capability should improve understanding of the Earth-to-orbit transient and enable better insulation system designs for in-space cryogenic propellant systems.

Tramel, R. W.↗