Search NASA⌕ Search

SEARCH · Search NASA

Results for “Enthalpy method”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1959

A General Method for Automatic Computation of Equilibrium Compositions and Theoretical Rocket Performance of Propellants

A general computer program for chemical equilibrium and rocket performance calculations was written for the IBM 650 computer with 2000 words of drum storage, 60 words of high-speed core storage, indexing registers, and floating point attachments. The program is capable of carrying out combustion and isentropic expansion calculations on a chemical system that may include as many as 10 different chemical elements, 30 reaction products, and 25 pressure ratios. In addition to the equilibrium composition, temperature, and pressure, the program calculates specific impulse, specific impulse in vacuum, characteristic velocity, thrust coefficient, area ratio, molecular weight, Mach number, specific heat, isentropic exponent, enthalpy, entropy, and several thermodynamic first derivatives.

Gordon, Sanford↗
NASA NTRS2008

Multi-wavelength Observations and Modelling of a Solar Flare

Aims: The aim of this work is to investigate the dynamic behavior of a C-class solar flare through the evolution of temperature, emission measure, energy loss and velocity. In particular, the variation of these properties with time are studied using multi-wavelength observations in combination with a recently developed 0-D hydrodynamic model. Methods: The temperature and emission measure evolution were studied using several instruments covering a wide range of temperatures - the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI, >5 MK), GOES-12 (5- 30 MK), the Transition Region and Coronal Explorer (TRACE 171 A, 1 MK) and the Coronal Diagnostic Spectrometer (CDS, 0.03-8 MK). The temperature and emission measure were analysed through the systematic cooling of flare plasma through the response functions of these instruments. These parameters were then investigated using the Enthalpy Based Thermal Evolution of Loops model (EBTEL). The Doppler shifts at both flare footpoints were analysed using five emission lines seen by CDS. Results: The flare began with clear evidence for pre-flare heating. Upflows of approx.90 km/s and low level emission, both observed in Fe XIX before the main impulsive phase were explained by pre-flare gentle chromospheric evaporation. During the main impulsive phase, the flare plasma was heated to a temperature of >13 MK in approximately 10 minutes. Explosive chromospheric evaporation was observed, driving upflows of approx.80 km/s in Fe XIX and simultaneous downflows of approx.20 km/s in He I and O v. At the peak of the Rare, conduction modelled by EBTEL was found to be the dominant loss mechanism, working efficiently to both lower the temperatures and drive gentle chromospheric evaporation. As the temperature fell below approx.8 MK, radiation became the dominant loss mechanism. During the final stages of the decay phase, downflowing plasma was observed at the footpoints in He I, O v and Mg x at velocities of up to approx.40 km/s, suggesting loop draining occurred. Conclusions. This is the first extensive study of the evolution of flare plasma using both spectroscopic and broad-band instruments in conjunction with a comprehensive hydrodynamic model. The flare began with pre-flare heating and then evolved following the predictions of the standard flare model. Detailed analysis of the plasma heating mechanisms was carried out and the heating function most consistent with observations was found to be Gaussian in shape. The simulations suggested that both direct heating and heating by a non-thermal beam played significant roles in this event.

Raftery, Claire L.↗
NASA NTRS1992

High heat flux measurements and experimental calibrations/characterizations

Recent progress in techniques employed in the measurement of very high heat-transfer rates in reentry-type facilities at the Arnold Engineering Development Center (AEDC) is described. These advances include thermal analyses applied to transducer concepts used to make these measurements; improved heat-flux sensor fabrication methods, equipment, and procedures for determining the experimental time response of individual sensors; performance of absolute heat-flux calibrations at levels above 2,000 Btu/cu ft-sec (2.27 kW/cu cm); and innovative methods of performing in-situ run-to-run characterizations of heat-flux probes installed in the test facility. Graphical illustrations of the results of extensive thermal analyses of the null-point calorimeter and coaxial surface thermocouple concepts with application to measurements in aerothermal test environments are presented. Results of time response experiments and absolute calibrations of null-point calorimeters and coaxial thermocouples performed in the laboratory at intermediate to high heat-flux levels are shown. Typical AEDC high-enthalpy arc heater heat-flux data recently obtained with a Calspan-fabricated null-point probe model are included.

Kidd, Carl T.↗
NASA NTRS2002

Computational Thermochemistry of Jet Fuels and Rocket Propellants

The design of new high-energy density molecules as candidates for jet and rocket fuels is an important goal of modern chemical thermodynamics. The NASA Glenn Research Center is home to a database of thermodynamic data for over 2000 compounds related to this goal, in the form of least-squares fits of heat capacities, enthalpies, and entropies as functions of temperature over the range of 300 - 6000 K. The chemical equilibrium with applications (CEA) program written and maintained by researchers at NASA Glenn over the last fifty years, makes use of this database for modeling the performance of potential rocket propellants. During its long history, the NASA Glenn database has been developed based on experimental results and data published in the scientific literature such as the standard JANAF tables. The recent development of efficient computational techniques based on quantum chemical methods provides an alternative source of information for expansion of such databases. For example, it is now possible to model dissociation or combustion reactions of small molecules to high accuracy using techniques such as coupled cluster theory or density functional theory. Unfortunately, the current applicability of reliable computational models is limited to relatively small molecules containing only around a dozen (non-hydrogen) atoms. We propose to extend the applicability of coupled cluster theory- often referred to as the 'gold standard' of quantum chemical methods- to molecules containing 30-50 non-hydrogen atoms. The centerpiece of this work is the concept of local correlation, in which the description of the electron interactions- known as electron correlation effects- are reduced to only their most important localized components. Such an advance has the potential to greatly expand the current reach of computational thermochemistry and thus to have a significant impact on the theoretical study of jet and rocket propellants.

Crawford, T. Daniel↗
NASA NTRSDate not supplied

Evaluating Novel Approaches for In Situ Analysis in the Search for Extraterrestrial Life

Universal approaches to detect extraterrestrial life are required to determine whether life exists outside of earth. A significant limitation of current detection strategies is the absence of a reliable, universal approach. They do not provide mechanisms to detect life that has evolved distinctly from life on Earth. We hypothesize that the most basic, essential features of life that can be detected are that they decrease their internal entropy at the expense of free energy obtained from their surroundings and that they produce heat associated with metabolic reactions. We propose that by using a thermodynamic analysis of the energy associated with these fundamental components of life, one can readily quantify the degree of free energy difference between life forms or their remnants and their abiotic surroundings. Analysis of the free energy entropy term (and thus structural complexity) and enthalpy term (as a manifestation of the energy release caused by metabolic reactions) may provide definitive detection of the presence of life. As a part of our analysis we review data available in the literature about the thermodynamics of microbial growth to assess the extent microbial metabolism and growth may contribute to negative entropy of maintenance (metabolism) and structural complexity. Additionally, we evaluate the different mechanisms by which structural complexity can be measured using calorimetry. An experimental framework is considered using concepts from previous work that could be used to further validate the calorimetric methods for life detection by potentially differentiating between life forms or their remnants from an abiotic surrounding. Results obtained demonstrate the applications of this approach and indicate that modification of a microcalorimetric technique has potential for incorporation into a space exploration platform as a life-detection payload.

Life Detection, Calorimetry↗
NASA NTRS2020

NASA Ames Thermophysics Ground Test Facilities Supporting Future Planetary Atmospheric Entry

A review of the current facility capabilities for testing Thermal Protection Systems and quantifying their entry environments at NASA Ames Research Center is presented based on the expected targets of interest to the Planetary Science and Astrobiology communities. While the operational capabilities of these facilities are generally considered sufficient for supporting future missions to targets of interest, expanded ground test capabilities such as larger sample sizes, flight-relevant gas mixtures, dusty environments, and flight-relevant shear/pressure combinations would reduce future entry vehicle design uncertainties and applied margins. These reduced uncertainties may translate into reduced entry vehicle masses, increased robustness, and decreased operational risks during entry phases for science missions. Expanded ground test capabilities would also offer the ability to study material failure modes in environments even more representative of flight than are currently achievable. A list of desired future test capabilities is presented along with suggestions of possible methods of achieving each. The main recommendation of this paper is the undertaking of a detailed study of the benefits and associated costs of each of these expanded capabilities to determine the best future path.

arc jet↗
NASA NTRS2000

Convection in Slab and Spheroidal Geometries

Three-dimensional numerical simulations of compressible turbulent thermally driven convection, in both slab and spheroidal geometries, are reviewed and analyzed in terms of velocity spectra and mixing-length theory. The same ideal gas model is used in both geometries, and resulting flows are compared. The piecewise-parabolic method (PPM), with either thermal conductivity or photospheric boundary conditions, is used to solve the fluid equations of motion. Fluid motions in both geometries exhibit a Kolmogorov-like k(sup -5/3) range in their velocity spectra. The longest wavelength modes are energetically dominant in both geometries, typically leading to one convection cell dominating the flow. In spheroidal geometry, a dipolar flow dominates the largest scale convective motions. Downflows are intensely turbulent and up drafts are relatively laminar in both geometries. In slab geometry, correlations between temperature and velocity fluctuations, which lead to the enthalpy flux, are fairly independent of depth. In spheroidal geometry this same correlation increases linearly with radius over the inner 70 percent by radius, in which the local pressure scale heights are a sizable fraction of the radius. The effects from the impenetrable boundary conditions in the slab geometry models are confused with the effects from non-local convection. In spheroidal geometry nonlocal effects, due to coherent plumes, are seen as far as several pressure scale heights from the lower boundary and are clearly distinguishable from boundary effects.

Porter, David H.↗
NASA NTRS2009

Simulation of RCC Crack Growth Due to Carbon Oxidation in High-Temperature Gas Environments

The carbon wall oxidation technique coupled with a CFD technique was employed to study the flow in the expanding crack channel caused by the oxidation of the channel carbon walls. The recessing 3D surface morphing procedure was developed and tested in comparison with the arcjet experimental results. The multi-block structured adaptive meshing was used to model the computational domain changes due to the wall recession. Wall regression rates for a reinforced carbon-carbon (RCC) samples, that were tested in a high enthalpy arcjet environment, were computationally obtained and used to assess the channel expansion. The test geometry and flow conditions render the flow regime as the transitional to continuum, therefore Navier-Stokes gas dynamic approach with the temperature jump and velocity slip correction to the boundary conditions was used. The modeled mechanism for wall material loss was atomic oxygen reaction with bare carbon. The predicted channel growth was found to agree with arcjet observations. Local gas flow field results were found to affect the oxidation rate in a manner that cannot be predicted by previous mass loss correlations. The method holds promise for future modeling of materials gas-dynamic interactions for hypersonic flight.

Titov, E. V.↗
NASA NTRS2001

MHD Energy Bypass Scramjet Engine

Revolutionary rather than evolutionary changes in propulsion systems are most likely to decrease cost of space transportation and to provide a global range capability. Hypersonic air-breathing propulsion is a revolutionary propulsion system. The performance of scramjet engines can be improved by the AJAX energy management concept. A magneto-hydro-dynamics (MHD) generator controls the flow and extracts flow energy in the engine inlet and a MHD accelerator downstream of the combustor accelerates the nozzle flow. A progress report toward developing the MHD technology is presented herein. Recent theoretical efforts are reviewed and ongoing experimental efforts are discussed. The latter efforts also include an ongoing collaboration between NASA, the US Air Force Research Laboratory, US industry, and Russian scientific organizations. Two of the critical technologies, the ionization of the air and the MHD accelerator, are briefly discussed. Examples of limiting the combustor entrance Mach number to a low supersonic value with a MHD energy bypass scheme are presented, demonstrating an improvement in scramjet performance. The results for a simplified design of an aerospace plane show that the specific impulse of the MHD-bypass system is better than the non-MHD system and typical rocket over a narrow region of flight speeds and design parameters. Equilibrium ionization and non-equilibrium ionization are discussed. The thermodynamic condition of air at the entrance of the engine inlet determines the method of ionization. The required external power for non-equilibrium ionization is computed. There have been many experiments in which electrical power generation has successfully been achieved by magneto-hydrodynamic (MHD) means. However, relatively few experiments have been made to date for the reverse case of achieving gas acceleration by the MHD means. An experiment in a shock tunnel is described in which MHD acceleration is investigated experimentally. MHD has several potential aerospace applications. The first is to improve the performance of hypersonic air-breathing engines for space launch and cruise vehicles. The second is to improve the performance of a high enthalpy wind tunnel. The third is to control a hypersonic vehicle. With such applications in mind, theoretical and experiments are being conducted at the NASA Ames Research Center to develop the MHD technology.

Mehta, Unmeel B.↗
NASA NTRSDate not supplied

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗
NASA NTRS1994

Shock-tunnel combustor testing for hypersonic vehicles

Proposed configurations for the next generation of transatmospheric vehicles will rely on air breathing propulsion systems during all or part of their mission. At flight Mach numbers greater than about 7 these engines will operate in the supersonic combustion ramjet mode (scramjet). Ground testing of these engine concepts above Mach 8 requires high pressure, high enthalpy facilities such as shock tunnels and expansion tubes. These impulse, or short duration facilities have test times on the order of a millisecond, requiring high speed instrumentation and data systems. One such facility ideally suited for scramjet testing is the NASA-Ames 16-Inch shock tunnel, which over the last two years has completed a series of tests for the NASP (National Aero-Space Plane) program at simulated flight Mach numbers ranging from 12-16. The focus of the experimental programs consisted of a series of classified tests involving a near-full scale hydrogen fueled scramjet combustor model in the semi-free jet method of engine testing whereby the compressed forebody flow ahead of the cowl inlet is reproduced (see appendix A). The AIMHYE-1 (Ames Integrated Modular Hypersonic Engine) test entry for the NASP program was completed in April 1993, while AIMHYE-2 was completed in May 1994. The test entries were regarded as successful, resulting in some of the first data of its kind on the performance of a near full scale scramjet engine at Mach 12-16. The data was distributed to NASP team members for use in design system verification and development. Due to the classified nature of the hardware and data, the data reports resulting from this work are classified and have been published as part of the NASP literature. However, an unclassified AIAA paper resulted from the work and has been included as appendix A. It contains an overview of the test program and a description of some of the important issues.

Loomis, Mark P.↗
NASA NTRS2017

Evaluation of the Uncertainty in JP-7 Kinetics Models Applied to Scramjets

One of the challenges of designing and flying a scramjet-powered vehicle is the difficulty of preflight testing. Ground tests at realistic flight conditions introduce several sources of uncertainty to the flow that must be addressed. For example, the scales of the available facilities limit the size of vehicles that can be tested and so performance metrics for larger flight vehicles must be extrapolated from ground tests at smaller scales. To create the correct flow enthalpy for higher Mach number flows, most tunnels use a heater that introduces vitiates into the flow. At these conditions, the effects of the vitiates on the combustion process is of particular interest to the engine designer, where the ground test results must be extrapolated to flight conditions. In this paper, the uncertainty of the cracked JP-7 chemical kinetics used in the modeling of a hydrocarbon-fueled scramjet was investigated. The factors that were identified as contributing to uncertainty in the combustion process were the level of flow vitiation, the uncertainty of the kinetic model coefficients and the variation of flow properties between ground testing and flight. The method employed was to run simulations of small, unit problems and identify which variables were the principal sources of uncertainty for the mixture temperature. Then using this resulting subset of all the variables, the effects of the uncertainty caused by the chemical kinetics on a representative scramjet flow-path for both vitiated (ground) and nonvitiated (flight) flows were investigated. The simulations showed that only a few of the kinetic rate equations contribute to the uncertainty in the unit problem results, and when applied to the representative scramjet flowpath, the resulting temperature variability was on the order of 100 K. Both the vitiated and clean air results showed very similar levels of uncertainty, and the difference between the mean properties were generally within the range of uncertainty predicted.

Norris, A. T.↗
NASA NTRS2000

DSMC Simulation of Separated Flows About Flared Bodies at Hypersonic Conditions

This paper describes the results of a numerical study of interacting hypersonic flows at conditions that can be produced in ground-based test facilities. The computations are made with the direct simulation Monte Carlo (DSMC) method of Bird. The focus is on Mach 10 flows about flared axisymmetric configurations, both hollow cylinder flares and double cones. The flow conditions are those for which experiments have been or will be performed in the ONERA R5Ch low-density wind tunnel and the Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock (LENS) tunnel. The range of flow conditions, model configurations, and model sizes provides a significant range of shock/shock and shock/boundary layer interactions at low Reynolds number conditions. Results presented will highlight the sensitivity of the calculations to grid resolution, contrast the differences in flow structure for hypersonic cold flows and those of more energetic but still low enthalpy flows, and compare the present results with experimental measurements for surface heating, pressure, and extent of separation.

Moss, James N.↗
NASA NTRSDate not supplied

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗
NASA NTRS1993

Containerless measurements on liquids at high temperatures

The application of containerless techniques for measurements of the thermophysical properties of high temperature liquids is reviewed. Recent results obtained in the materials research laboratories at Intersonics are also presented. Work to measure high temperature liquid properties is motivated by both the need for reliable property data for modeling of industrial processes involving molten materials and generation of data form basic modeling of materials behavior. The motivation for this work and examples of variations in thermophysical property values from the literature are presented. The variations may be attributed to changes in the specimen properties caused by chemical changes in the specimen and/or to measurement errors. The two methods used to achieve containerless conditions were aeroacoustic levitation and electromagnetic levitation. Their qualities are presented. The accompanying slides show the layout of levitation equipment and present examples of levitated metallic and ceramic specimens. Containerless techniques provide a high degree of control over specimen chemistry, nucleation and allow precise control of liquid composition to be achieved. Effects of minor additions can thus be measured in a systematic way. Operation in reduced gravity enables enhanced control of liquid motion which can allow measurement of liquid transport properties. Examples of nucleation control, the thermodynamics of oxide contamination removal, and control of the chromium content of liquid aluminum oxide by high temperature containerless processes are presented. The feasibility of measuring temperature, emissivity, liquidus temperature, enthalpy, surface tension, density, viscosity, and thermal diffusivity are discussed in the final section of the paper.

Weber, Richard↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗
NASA NTRS2006

BLIMPK/Streamline Surface Catalytic Heating Predictions on the Space Shuttle Orbiter

This paper describes the results of an analysis of localized catalytic heating effects to the U.S. Space Shuttle Orbiter Thermal Protection System (TPS). The analysis applies to the High-temperature Reusable Surface Insulation (HRSI) on the lower fuselage and wing acreage, as well as the critical Reinforced Carbon-Carbon on the nose cap, chin panel and the wing leading edge. The object of the analysis was to use a modified two-layer approach to predict the catalytic heating effects on the Orbiter windward HRSI tile acreage, nose cap, and wing leading edge assuming localized highly catalytic or fully catalytic surfaces. The method incorporated the Boundary Layer Integral Matrix Procedure Kinetic (BLIMPK) code with streamline inputs from viscous Navier-Stokes solutions to produce heating rates for localized fully catalytic and highly catalytic surfaces as well as for nominal partially catalytic surfaces (either Reinforced Carbon-Carbon or Reaction Cured Glass) with temperature-dependent recombination coefficients. The highly catalytic heating results showed very good correlation with Orbiter Experiments STS-2, -3, and -5 centerline and STS-5 wing flight data for the HRSI tiles. Recommended catalytic heating factors were generated for use in future Shuttle missions in the event of quick-time analysis of damaged or repaired TPS areas during atmospheric reentry. The catalytic factors are presented along the streamlines as well as a function of stagnation enthalpy so they can be used for arbitrary trajectories.

Marichalar, Jeremiah J.↗