Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal equilibrium”

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 73 records · Page 4

Navier-Stokes simulation of real gas flows in nozzles

Air flow in a hypersonic nozzle causes real gas effects due to reaction among the species constituting air. Such reactions may be in chemical equilibrium or in chemical nonequilibrium. Here using the CSCM upwind scheme for the compressible Navier-Stokes equations, the real gas flowfield in an arcjet nozzle is computed for both the equilibrium case and the nonequilibrium case. A hypersonic nozzle flow arising from a pebble bed heated plenum is also computed for the equilibrium situation. Between the equilibrium cases, the chemistry is treated by two different schemes and comments are made as to computational complexity. For the nonequilibrium case, a full set of seventeen reactions and full implicit coupling of five species with gasdynamics is employed to compute the flowfield. For all cases considered here the gas is assumed to be a calorically imperfect mixture of ideal gases in thermal equilibrium.

Nagaraj, N.↗

Electromagnetic fluctuations for anisotropic media and the generalized Kirchhoff's law

In this paper the polarimetric emission parameters for anisotropic media are derived using the generalized Kirchhoff's law for media with a uniform temperature and the fluctuation-dissipation theory for media with a temperature profile. Both finite-size objects and half-space media are considered. When the object has a uniform temperature across its body, the Kirchhoff's law, based on the condition of energy conservation in thermal equilibrium is generalized to obtain the emission parameters of an anisotropic medium, which can be interpreted as the absorptivity or the absorption cross section of the complementary object with a permittivity that is the transpose of the original object. When the medium has a nonuniform temperature distribution, the fluctuation-dissipation theory is applied for deriving the covariances between vector components of the thermal currents and, consequently, the covariances of the polarizations of electric fields radiated by the thermal currents. To verify the formulas derived from the fluctuation-dissipation theory, we let the temperature of the object be a constant and show that the results reduce to those obtained from the generalized Kirchhoff's law.

Yueh, Simon H.↗

Miniature Arcs for Synthesis of Carbon Nanotubes in Microgravity

Although many methods are available for producing single-walled carbon nanotubes (SWNTs), the conventional carbon arc process remains the most popular due to its simplicity and large production rate. In the carbon arc, SWNTs are catalytically synthesized by rapidly evaporating a graphite anode impregnated with NiN metal catalyst from which the nanotubes grow in an inert atmosphere. However, high temperatures inside the carbon arc generate strong buoyancy driven convection, and it is hypothesized that the non-uniform environment created by this flow has a large effect on the growth and morphology of the SWNTs. To study the effect of buoyancy on the arc process, a miniature carbon arc apparatus was developed to synthesize SWNTs in a microgravity environment substantially free from these strong convective flows. The reactor was operated for either 2.2 or 5 seconds during free-fall in the drop towers at the NASA Glenn Research Center. Two apparatus designs differing mainly in their production rate and power capacity were investigated. The first consisted of a miniaturized carbon arc employing a 1 mm diameter graphite anode and powered by a 0.54 F capacitor bank charged to 65 V. The second, larger apparatus employed a 4 mm diameter anode and was powered by a portable battery pack capable of providing in excess of 300 amps at 30 volts to the arc for the duration of a 5 second drop. Initial results indicated that transient heating is a very large effect in the short-duration drop tower carbon arcs, and thermal equilibrium of the arc plasma, buffer gas, and apparatus was not attained during the short microgravity periods. In addition, removal of the buoyant convection by the microgravity now allowed clear observation of large jets of evaporated carbon vapor streaming from the anode and mixing with the inert buffer gas. The initial mixing of these jets with the cold buffer gas combined with the thermal transient made it difficult to establish a uniform high temperature environment around the arc in the 2.1 to 5 second microgravity time interval, and even with a very high-powered arc, the arc region was cooler than in continuously operated arcs. Despite these difficulties, the miniature arc produced SWNTs in microgravity. However, given the large thermal transient to overcome, no dramatic difference in sample yield or composition was noted between normal gravity and q2-,andL%econd long microgravity runs.

Alford, J. M.↗

The neutral atomic phases of the interstellar medium

We calculate the thermal equilibrium gas temperature of the diffuse interstellar medium. Our method incorporates a new photoelectric heating rate from small grains and polycyclic aromatic hydrocarbons (PAHs) that accounts for a size distribution of particles extending from 100 to 3 A radius. We also include a detailed treatment of the ionization rates and heating due to the soft X-ray background and due to cosmic rays. Phase diagrams (thermal pressure P versus hydrogen density n) are presented for gas that is illuminated by local interstellar far-ultraviolet (FUV) and X-ray radiation fields. A stable two-phase medium is produced with thermal pressure in the range P/k approximately = to 10(exp 3-4) K/cc. We demonstrate that photoelectric heating from PAHs dominates in the warm neutral phase (WNM) and cold neutral phase (CNM). If the C II (158 micrometers cooling per hydrogen nucleus in the solar neighborhood represents an average value for the Galaxy, we predict L(sub CII) approximately = to 7 x 10(exp 7) solar luminosities from the CNM in the Galaxy, comparable to that observed by the Cosmic Background Explorer (COBE). We discuss the dependence of the results on absorbing column density, gas phase abundances, dust abundances and metallicity, FUV field, and the X-ray radiation field. These results will be useful in modeling the multiphase structure of high-velocity clouds in the halo, the interstellar matter (ISM) at other galactocentric radii, and the ISM in external galaxies and galactic nuclei.

Wolfire, M. G.↗

Squeezed states and graviton-entropy production in the early universe

Squeezed states are a very useful framework for the quantum treatment of tensor perturbations (i.e. gravitons production) in the early universe. In particular, the non equilibrium entropy growth in a cosmological process of pair production is completely determined by the associated squeezing parameter and is insensitive to the number of particles in the initial state. The total produced entropy may represent a significant fraction of the entropy stored today in the cosmic blackbody radiation, provided pair production originates from a change in the background metric at a curvature scale of the Planck order. Within the formalism of squeezed thermal states it is also possible to discuss the stimulated emission of gravitons from an initial thermal bath, under the action of the cosmic gravitational background field. We find that at low energy the graviton production is enhanced, if compared with spontaneous creation from the vacuum; as a consequence, the inflation scale must be lowered, in order not to exceed the observed CMB quadrupole anisotropy. This effect is important, in particular, for models based on a symmetry-breaking transition which require, as initial condition, a state of thermal equilibrium at temperatures higher than the inflation scale and in which inflation has a minimal duration.

Giovannini, Massimo↗

Containerless Ripple Turbulence

One of the longest standing unsolved problems in physics relates to the behavior of fluids that are driven far from equilibrium such as occurs when they become turbulent due to fast flow through a grid or tidal motions. In turbulent flows the distribution of vortex energy as a function of the inverse length scale [or wavenumber 'k'] of motion is proportional to 1/k(sup 5/3) which is the celebrated law of Kolmogorov. Although this law gives a good description of the average motion, fluctuations around the average are huge. This stands in contrast with thermally activated motion where large fluctuations around thermal equilibrium are highly unfavorable. The problem of turbulence is the problem of understanding why large fluctuations are so prevalent which is also called the problem of 'intermittency'. Turbulence is a remarkable problem in that its solution sits simultaneously at the forefront of physics, mathematics, engineering and computer science. A recent conference [March 2002] on 'Statistical Hydrodynamics' organized by the Los Alamos Laboratory Center for Nonlinear Studies brought together researchers in all of these fields. Although turbulence is generally thought to be described by the Navier-Stokes Equations of fluid mechanics the solution as well as its existence has eluded researchers for over 100 years. In fact proof of the existence of such a solution qualifies for a 1 M$ millennium prize. As part of our NASA funded research we have proposed building a bridge between vortex turbulence and wave turbulence. The latter occurs when high amplitude waves of various wavelengths are allowed to mutually interact in a fluid. In particular we have proposed measuring the interaction of ripples [capillary waves] that run around on the surface of a fluid sphere suspended in a microgravity environment. The problem of ripple turbulence poses similar mathematical challenges to the problem of vortex turbulence. The waves can have a high amplitude and a strong nonlinear interaction. Furthermore, the steady state distribution of energy again follows a Kolmogorov scaling law; in this case the ripple energy is distributed according to 1/k (sup 7/4). Again, in parallel with vortex turbulence ripple turbulence exhibits intermittency. The problem of ripple turbulence presents an experimental opportunity to generate data in a controlled, benchmarked system. In particular the surface of a sphere is an ideal environment to study ripple turbulence. Waves run around the sphere and interact with each other, and the effect of walls is eliminated. In microgravity this state can be realized for over 2 decades of frequency. Wave turbulence is a physically relevant problem in its own right. It has been studied on the surface of liquid hydrogen and its application to Alfven waves in space is a source of debate. Of course, application of wave turbulence perspectives to ocean waves has been a major success. The experiment which we plan to run in microgravity is conceptually straightforward. Ripples are excited on the surface of a spherical drop of fluid and then their amplitude is recorded with appropriate photography. A key challenge is posed by the need to stably position a 10cm diameter sphere of water in microgravity. Two methods are being developed. Orbitec is using controlled puffs of air from at least 6 independent directions to provided the positioning force. This approach has actually succeeded to position and stabilize a 4cm sphere during a KC 135 segment. Guigne International is using the radiation pressure of high frequency sound. These transducers have been organized into a device in the shape of a dodecahedron. This apparatus 'SPACE DRUMS' has already been approved for use for combustion synthesis experiments on the International Space Station. A key opportunity presented by the ripple turbulence data is its use in driving the development of codes to simulate its properties.

Putterman, Seth↗

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗

The aeronomy of vibrationally excited ozone

Theoretical calculations show that above 80 km in the earth's atmosphere the production of vibrationally excited ozone by chemical processes leads to number densities which are usually larger than those expected for local thermodynamic equilibrium. Quenching of highly excited molecules produced in O+O2+M, O3+M provided a significant source of the lower lying states above the mesopause while the 9.6 microns emission of O3 (0,0,1) was a major sink. Analysis of available laboratory results implied that reactions involving excited ozone play a significant role in the global ozone balance despite the relatively small abundance of the molecule. However, this effect is implicit in many of the rate coefficients currently used in stratospheric calculations. In the upper mesosphere and lower thermosphere, where the excited state populations differ from those for thermal equilibrium, published reaction rate data are not necessarily applicable to aeronomic calculations.

Frederick, J. E.↗

Measurement of Vibrational Non-Equilibrium in a Supersonic Freestream Using Dual-Pump CARS

Measurements have been conducted at the University of Virginia Supersonic Combustion Facility of the flow in a constant area duct downstream of a Mach 2 nozzle, where the airflow has first been heated to approximately 1200 K. Dual-pump CARS was used to acquire rotational and vibrational temperatures of N2 and O2 at two planes in the duct at different downstream distances from the nozzle exit. Wall static pressures in the nozzle are also reported. With a flow of clean air, the vibrational temperature of N2 freezes at close to the heater stagnation temperature, while the O2 vibrational temperature is about 1000 K. The results are well predicted by computational fluid mechanics models employing separate "lumped" vibrational and translational/rotational temperatures. Experimental results are also reported for a few percent steam addition to the air and the effect of the steam is to bring the flow to thermal equilibrium.

Cutler, Andrew D.↗

Nonequilibrium Supersonic Freestream Studied Using Coherent Anti-Stokes Raman Spectroscopy

Measurements were conducted at the University of Virginia Supersonic Combustion Facility of the flow in a constant-area duct downstream of a Mach 2 nozzle. The airflow was heated to approximately 1200 K in the facility heater upstream of the nozzle. Dual-pump coherent anti-Stokes Raman spectroscopy was used to measure the rotational and vibrational temperatures of N2 and O2 at two planes in the duct. The expectation was that the vibrational temperature would be in equilibrium, because most scramjet facilities are vitiated air facilities and are in vibrational equilibrium. However, with a flow of clean air, the vibrational temperature of N2 along a streamline remains approximately constant between the measurement plane and the facility heater, the vibrational temperature of O2 in the duct is about 1000 K, and the rotational temperature is consistent with the isentropic flow. The measurements of N2 vibrational temperature enabled cross-stream nonuniformities in the temperature exiting the facility heater to be documented. The measurements are in agreement with computational fluid dynamics models employing separate lumped vibrational and translational/rotational temperatures. Measurements and computations are also reported for a few percent steam addition to the air. The effect of the steam is to bring the flow to thermal equilibrium, also in agreement with the computational fluid dynamics.

Cutler, Andrew D.↗

Convective Melting of Particles in Flow Under Microgravity Conditions

Study of melting of dispersed or packed solid particles in a fluid under gravity and microgravity conditions provides benchmark information for many engineering applications such as material processing, environmental assessment and protection and space fire protection. During such processes, packed or dispersed solid particles are interacting with fluid flow at above-melting temperatures. By unmasking the buoyancy effects in coupled flow, phase change and heat transport phenomena, a better understanding of melting rate in non-thermal equilibrium, convective conditions can be studied. A series of flight experiments were conducted onboard the NASA KC-135 Microgravity Research Airplane. The Particle Melting in Plow (PMF) module was designed to allow flow through the initially packed ice particles at controlled temperature and velocity. To achieve this, a close-loop flow system was designed. Video images were taken to record the visualization of the melting process, from which a time variation of packed particle thickness distribution at different times can be obtained by the image analysis method. The fluid temperature distribution within the melting zone is measured by thermocouples. An infrared camera was mounted from the top of the test section to record the ice-water thermal images at a given location. The results from thermal images yield local temperature variation between melting solid and liquid and local Stephan number. Typical results for a number of cases are presented. The mathematical model, describing mass, energy, and momentum balance equations for the liquid and solid phases, is presented. It is found that melting rate is influenced mainly by the ratio of Reynolds number (based on the initial particle diameter) to the Froud number, and Stephan number. At the absence of gravity, Froud number approaches zero, Reynolds number and Stephan number become dominant factors governing the melting rate. The numerically determined results are compared with the experimental ones. It is found that the discrepancy between the predicted and measured melting rate is largely due to the inaccuracy in the constitutive equations for effective thermal physical properties, such as effective thermal conductivity and diffusivity, and transport properties, such as particle interaction coefficient, and local heat transfer coefficient of particles.

Stewart, C.↗

Athermal Energy Loss from X-rays Deposited in Thin Superconducting Films on Solid Substrates

When energy is deposited in a thin-film cryogenic detector, such as from the absorption of an X-ray, an important feature that determines the energy resolution is the amount of athermal energy that can be lost to the heat bath prior to the elementary excitation systems coming into thermal equilibrium. This form of energy loss will be position-dependent and therefore can limit the detector energy resolution. An understanding of the physical processes that occur when elementary excitations are generated in metal films on dielectric substrates is important for the design and optimization of a number of different types of low temperature detector. We have measured the total energy loss in one relatively simple geometry that allows us to study these processes and compare measurements with calculation based upon a model for the various di.erent processes. We have modeled the athermal phonon energy loss in this device by finding an evolving phonon distribution function that solves the system of kinetic equations for the interacting system of electrons and phonons. Using measurements of device parameters such as the Debye energy and the thermal di.usivity we have calculated the expected energy loss from this detector geometry, and also the position-dependent variation of this loss. We have also calculated the predicted impact on measured spectral line-shapes, and shown that they agree well with measurements. In addition, we have tested this model by using it to predict the performance of a number of other types of detector with di.erent geometries, where good agreement is also found.

Kozorezov, Alexander G.↗

Behavior of Langmuir Probes in Non-Equilibrium Plasmas

Langmuir probes are diagnostic tools used to determine electron temperature, number density, and plasma potential. Irving Langmuir first used an electrostatic probe in the 1920s to find these characteristics in ionized gases. Single, double, and triple Langmuir probes are commonly used in plasma diagnostics because of their relative simplicity. In the single probe, a swept voltage is applied between the probe tip and circuit common to acquire a waveform showing the collected current as a function of applied voltage. A double Langmuir probe consists of two tips, both inserted into the plasma, with a voltage applied between them. As this voltage is swept, a current-voltage characteristic is measured. In a triple probe three probe tips are electrically coupled to each other with constant non-swept voltages applied between each of the tips. The voltages are selected to represent three points on the single Langmuir probe I-V curve. Elimination of the voltage sweep makes it possible to measure time-varying plasma properties in transient plasmas. Triple Langmuir probe measurements have been widely employed for various types of plasmas, including pulsed and time-varying plasmas such as those seen in pulsed plasma thrusters (PPTs), dense plasma focus devices, plasma flows, and fusion experiments. The typical Langmuir probe analysis for determining electron temperature and number density of the plasma (for a single, double, or triple Langmuir probe) includes an assumption that the plasma is in thermal equilibrium. While the this assumption may be justified for some applications, it is unlikely that it is fully justifiable for pulsed and time-varying plasmas or for the entire time a plasma device is in use. In the present work, we model the responses of Langmuir probes as they are inserted into a range of simple equilibrium and non-equilibrium plasmas. We return to basic governing equations of probe current collection and compute the current to the probes for a distribution function consisting of two Maxwellian distributions with different temperatures (the two-temperature Maxwellian). A variation of this method is also employed, where one of the Maxwellians is offset from zero (in velocity space) to add a suprathermal beam of electrons to the tail of the main Maxwellian distribution (the bump-on-the-tail distribution function). For a range of parameters in these non-Maxwellian distributions, we compute the current collection to the probes. Comparing the distribution function that was assumed a priori with the plasma density and temperature one would infer when applying standard probe theory to analyze the collected currents serves to illustrate the effect a non- Maxwellian plasma would have on results interpreted using the equilibrium probe current collection theory, allowing us to state the magnitudes of these deviations as a function of the assumed distribution function properties.

Polzin, Kurt A.↗

An "Inefficient Fin" Non-Dimensional Parameter to Measure Gas Temperatures Efficiently

A gas containment vessel that is not in thermal equilibrium with the bulk gas can affect its temperature measurement. The physical nature of many gas dynamics experiments often makes the accurate measurement of temperature a challenge. The environment itself typically requires that the thermocouple be sheathed, both to protect the wires and hot junction of the instrument from their environment, and to provide a smooth, rigid surface for pressure sealing of the enclosure. However, that enclosure may also be much colder than the gas to be sensed, or vice-versa. Either way, the effect of such gradients is to potentially skew the temperature measurements themselves, since heat may then be conducted by the instrument. Thermocouple designers traditionally address this problem by insulating the sheath from the thermocouple leads and hot junction as much as possible. The thermocouple leads are typically packed in a ceramic powder inside the sheath, protecting them somewhat from temperature gradients along the sheath, but there is no effective mechanism to shield the sheath from the enclosure body itself. Standard practice dictates that thermocouples be used in installations that do not present large thermal gradients along the probe. If this conduction dominates heat transfer near the tip of the probe, then temperature measurements may be expected to be skewed. While the same problem may be experienced in the measurement of temperature at various points within a solid in a gradient, it tends to be aggravated in the measurements of gas temperature, since heat transfer dependent on convection is often less efficient than conduction along the thermocouple. The proposed solution is an inefficient fin thermocouple probe. Conventional wisdom suggests that in many experiments where gas flows through an enclosure (e.g., flow in pipe, manifold, nozzle, etc.), the thermocouple be introduced flush to the surface, so as not to interfere with the flow. In practice, however, many such experiments take place where the flow is already turbulent, so that a protruding thermocouple probe has a negligible effect on the flow characteristics. The key question then becomes just how far into the flow should a thermocouple protrude in order to properly sense the gas temperature at that point. Modeling the thermocouple as an "inefficient fin" directly addresses this question. The appropriate assumptions in this case are: one-dimensional conduction along the fin; steady-state, constant, and homogeneous thermal conductivity; negligible radiation; and a uniform, constant heat transfer coefficient over the probe surface. It is noted that the nature of the ceramic-filled probe makes the key assumption of homogeneous thermal conductivity that much more conservative.

Lemieux, Patrick↗

Navier-Stokes and viscous shock-layer solutions for radiating hypersonic flows

Results are presented from the Navier-Stokes and viscous shock-layer (VSL) calculations with nonequilibrium and equilibrium chemistry, respectively. These calculations contain coupling to the Aerotherm radiation code RAD. A simplified form of the electron energy equation is used to obtain an electron temperature in the Navier-Stokes calculations. The radiation in the flowfield is calculated using this temperature. The Navier-Stokes code is used at high altitude only, whereas the VSL code is employed for the entire entry period to make estimates of the radiative and convective heating to the Fire II vehicle. Results from the Navier-Stokes code have also been compared with the predictions of Lee and Kawamura, who used gray-gas radiation model and thin-layer Navier-Stokes equations. Quite good agreement is obtained between the measured and computed values of radiative and convective heating from the VSL code in th medium-to-low altitude flight regime of the Fire II vehicle. At high altitudes, the Navier-Stokes calculations considerably overpredict the Fire II flight data for radiative intensity. This is attributed to the deficiencies in the Aerotherm radiation model when used for low-density flight conditions. This model contains the thermal equilibrium assumption and precludes accounting for the collision-limiting phenomenon at high altitudes. Present Navier-Stokes calculations highlight the effect of these assumptions on radiative heating calculations for such conditions.

Gupta, Roop N.↗

Predictive Modeling of Carbon Ablators Using Micro and Macro-Scale Modeling

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗

Predictive Modeling of Carbon Ablators

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗

Abundant Cool Magnesium-Rich Pyroxene Crystals in Comet Hale-Bopp

Modeling of the observed dust emission from Comet Hale-Bopp over a large range of heliocentric distances (2.8 AU - 0.93 AU -1.7 AU) led to the discovery of Mg-rich pyroxene crystals in the coma These pyroxene crystals are apparent in the 10 micron spectrum only when the comet is close to perihelion (r(sub h) = 1.2 AU) because they are cooler than the other silicate minerals. The pyroxene crystals are cooler than the other species because of their high Mg-content. They do not absorb as efficiently as the other silicate minerals. Given the same high Mg content of Mg/(Mg + Fe) = 0.9, radiative equilibrium computations show that pyroxene crystals are expected to be 150 K cooler than olivine crystals. The pyroxene crystals are also about 10x more abundant than the other silicate mineral species. Their high Mg content and relatively large abundance are in agreement with the preponderance of pyroxene interplanetary dust particles (IDPs) and the recent reanalysis of the PUMA-I flyby of Comet Halley. Before Hale-Bopp, only olivine crystals were detected spectroscopically in comets, probably because the pyroxene crystals are less optically active, hence significantly cooler and harder to detect in contrast to the warmer silicate species. Determining the relative abundances of silicate minerals depends on their Mg contents. If the pyroxene crystals in Comet Hale-Bopp are solar nebula condensates, then they probably had to form during the early FU Orionis epoch when the inner disk was hot enough and be transported out to the region of formation of icy planetesimals without being reheated. Reheating events appear to reincorporate Fe back into the crystals or form Fe-rich rims, which are not consistent with the high-Mg-content crystals. The condensation of Mg-rich pyroxene crystals is expected from solar nebula thermal equilibrium computations. However, their subsequent transport to the outer solar nebula unaltered has yet to be theoretically demonstrated. The discovery of Mg-rich crystals in Comet Hale-Bopp and in AGB stars opens the possibility that these crystals are relic interstellar grains. One-third of IDPs have been shown to have significant deuterium enrichments, thus indicating that they contain presolar material. By spectroscopic analogy to IDPs the Mg-rich pyroxene crystals in Comet Hale-Bopp may be presolar grains. If so, then the comet contains largely ISM silicates. ISM grains may have been the dominant source of dust in the outer early solar nebula.

Wooden, D. H.↗