Engineering PapersSearch

SEARCH · Engineering Papers

Results for “mixture model”

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.

Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching

Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.

line chilldown

Experimental Effervescence and Freezing Point Depression Measurements of Nitrogen in Liquid Methane-Ethane Mixtures

NASA is designing an unmanned submarine to explore the depths of the hydrocarbon-rich seas on Saturn’s moon Titan. Data from Cassini indicates that the Titan north polar environment sustains stable seas of variable concentrations of ethane, methane, and nitrogen, with a surface temperature near 93 K. The submarine must operate autonomously, study atmosphere/sea exchange, interact with the seabed, hover at the surface or any depth within the sea, and be capable of tolerating variable hydrocarbon compositions. Currently, the main thermal design concern is the effect of effervescence on submarine operation, which affects the ballast system, science instruments, and propellers. Twelve effervescence measurements on various liquid methane-ethane compositions with dissolved gaseous nitrogen are thus presented from 1.5 bar to 4.5 bar at temperatures from 92 K to 96 K to simulate the conditions of the seas. After conducting effervescence measurements, two freezing point depression measurements were conducted. The freezing liquid line was depressed more than 15 K below the triple point temperatures of pure ethane (90.4 K) and pure methane (90.7 K). Experimental effervescence measurements will be used to compare directly with effervescence modeling to determine if changes are required in the design of the thermal management system as well as the propellers.

Freezing Point Depression

Metallurgical and Thermal Properties of 57Ni-40Ti-3Hf (wt%)

Selected properties of 57Ni-40Ti-3Hf (wt%) were measured to further develop this material for use in aerospace rolling-element bearings. The average grain size and microindentation hardness of the material at room temperature were 90 μm and 549 HV, respectively. From room temperature to 500 °C, electrical resistivity increased from 8.9×10–7 to 11.3×10–7 Ω⋅m, thermal expansion coefficients varied from 1.1×10–5 to 1.2×10–5 mm/mm⋅°C, specific heat increased from 0.44 to 0.53 W⋅sec/g⋅K, thermal diffusivity rose from approximately 0.003 to 0.005 cm2/sec and thermal conductivity increased from approximately 0.09 to 0.19 W/cm⋅K. The behavior of the material can be explained by well-known metallurgical theory. These results are expected to provide a basis for design and modeling of future aerospace components.

rule of mixtures

RTE: A Computer Code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Mohammad H N Naraghi

A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations

Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.

Radiation

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG

Prediction of prompt NO(x) in hydrocarbon air flames

The gas turbine industry is directing particular attention to very low NOx combustors, whether for aircraft or land based CCGT systems. These low NOx combustors frequently use liquid fuels or natural gas burning under very lean premixed conditions with air or under rich-lean conditions, although only the first case is studied here. In land based systems, diluted steam or nitrogen are sometimes injected into the combustion chambers to reduce flame temperature. The NOx emissions from such systems are the product of three chemical mechanisms which are interrelated: the hydrocarbon prompt NO, the thermal NO (extended Zeldovich mechanism), and the nitrous oxide route to NO. Formation of NO2 from NO also occurs, as well as emission of carbon monoxide and unburnt hydrocarbons. When the fuel-oxidant proportion decreases towards leaner conditions, flame temperatures are lowered, resulting in the total NO being reduced and the thermal-NO contribution greatly diminished to the benefit of the remaining two mechanisms of NO formation. While knowledge of the elementary reactions and their chemical kinetics concerning methane and simple hydrocarbons combustion has existed for a number of years, its use for computer modeling is limited to simple flow dynamics configurations. Nevertheless, understanding of such combusting flows under a wide range of experimental conditions allows for analogies or speculations with more complex actual systems. Such understanding can be achieved by means of one dimensional laminar premixed flame modeling, with a full chemical mechanism which incorporates the three routes of NO formation. Complementary to this understanding is the modeling of the actual combustion system using a full description of the fluid dynamics coupled with a reduced chemical scheme, which is then compared against the first model. The objective of this investigation is to evaluate the relative importance of the three mechanisms of NO formation in lean premixed methane-air combustion with increasing pressure using the one dimensional plug flow package PREMIX, and to test the validity of a three dimensional model with a global chemical mechanism against the one dimensional model in the atmospheric pressure case. Methane is chosen because it is the only mechanism which is reasonably well known and is a good guide to the behavior of other hydrocarbons. The mixture ratio chosen is richer than that in lean gas turbines, but the combustion of this mixture with the low preheat gives realistic gas turbine final flame temperatures. Conditions of NO2 formation are also analyzed in the one dimensional model and results extrapolated to the case of gas turbines.

Valerie Dupont

Al/Cl2 molten salt battery

Molten salt battery has been developed with theoretical energy density of 5.2 j/kg (650 W-h/lb). Battery, which operates at 150 C, can be used in primary mode or as rechargeable battery. Battery has aluminum anode and chlorine cathode. Electrolyte is mixture of AlCl3, NaCl, and some alkali metal halide such as KCl.

Giner, J.

The Future of Oaks in the Santa Monica Mountains: A Case Study in Using Remote Sensing Data for Species Distributions Models

The Woolsey Fire began on November 8, 2018, and lasted for almost two weeks, during which it burned almost 100,000 acres of valuable landscape and habitat, including a vast area of woodland. The persistence of key woodland species provides aesthetic, monetary, and ecological value to the landscape through carbon sequestration, air temperature moderation, and erosion mitigation, among other ecosystem services. This study investigated the impact of the Woolsey Fire on native woodland species distributions and identified areas suitable for restoration within the Santa Monica Mountains National Recreation Area. The team partnered with the Resource Conservation District of the Santa Monica Mountains; National Park Service, Santa Monica Mountains National Recreation Area; California Department of Parks and Recreation, Los Angeles County Division; County of Los Angeles Fire Department, Prevention Services Bureau, Forestry Division; County of Los Angeles Department of Regional Planning; and the University of Montana. The Earth observations used include data from Landsat 8 Operational Land Imager, NASA ER-2 Jet Airborne Visible InfraRed Imaging Spectrometer, Shuttle Radar Topography Mission, and RapidEye. The team produced maps of burn severity from the Woolsey Fire, its impact on plant species distributions, and habitat suitability projections for 2050 and 2099 to assist partners in prioritizing areas for restoration. A plant community classification was successfully created using Multiple Endmember Spectral Mixture Analysis (MESMA). Overall accuracy was assessed at 90.54% by comparing the classification to validation pixels derived from ground truth information provided by our partners.

Roger Ly

Calculations of the Moon's Thermal History at Different Concentrations of Radioactive Elements, Taking Into Account Differentiation on Melting

Calculations of the thermal history of the moon were carried out by solving the thermal conductivity equation for the case in which the heat sources are the long-lived radioactive elements Th, U, and K 40 . The concentrations of these elements were adjusted to give four variations of the heat flow: 1.35 x 10 -8 and 0.91 x 10 -6 cal cm -2 s -1 (I and I**), 0.61 x 10 -6 (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x 10 -6 (variant III, the chondrite model of Urey and MacDonald). In the same calculations, we considered layering of the differentiated material with transport to the surface of the radioactive elements after the temperature of the layer rose to 200 K above the melting temperature, which is considered in five variants that differ in the amount of transported radioactive elements: 100 percent ( n = 1), 80 percent n = 0.8), 40 percent, 20 percent ( n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calculations show that the interior of the Moon was heated to melting during the first 0.7 to 2.3 x 10 9 years. The maximum fusion involved practically the entire Moon to a distance from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10 9 years ago (I, I**), or 2.5 to 3.0 x 10 9 years ago (II, III) and continued for 1 to 2 x 10 9 years. Today the Moon is cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of radioactive elements, a heat flow of 0.9 to 0.95 x 10 -6 cal cm -2 s -1 , which agrees with the radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data.

O I Ornatskaya

Studies of Flame Structure in Microgravity

The present research endeavor is concerned with gaining fundamental understanding of the configuration, structure, and dynamics of laminar premixed and diffusion flames under conditions of negligible effects of gravity. Of particular interest is the potential to establish and hence study the properties of spherically- and cylindrically-symmetric flames and their response to external forces not related to gravity. For example, in an earlier experimental study of the burner-stabilized cylindrical premixed flames, the possibility of flame stabilization through flow divergence was established, while the resulting one-dimensional, adiabatic, stretchless flame also allowed an accurate means of determining the laminar flame speeds of combustible mixtures. We have recently extended our studies of the flame structure in microgravity along the following directions: (1) Analysis of the dynamics of spherical premixed flames; (2) Analysis of the spreading of cylindrical diffusion flames; (3) Experimental observation of an interesting dual luminous zone structure of a steady-state, microbuoyancy, spherical diffusion flame of air burning in a hydrogen/methane mixture environment, and its subsequent quantification through computational simulation with detailed chemistry and transport; (4) Experimental quantification of the unsteady growth of a spherical diffusion flame; and (5) Computational simulation of stretched, diffusionally-imbalanced premixed flames near and beyond the conventional limits of flammability, and the substantiation of the concept of extended limits of flammability. Motivation and results of these investigations are individually discussed.

C K Law

Dissociation and Recombination Effects on the Performance of Pulse Detonation Engines

This paper summarizes major theoretical results for pulse detonation engine performance taking into account real gas chemistry, as well as significant performance differences resulting from the presence of ram and compression heating. An unsteady CFD analysis, as well as a thermodynamic cycle analysis, was conducted in order to determine the actual and the ideal performance for an air-breathing pulse detonation engine (PDE) using either a hydrogen-air or ethylene-air mixture over a flight Mach number range from 0 to 4. The results clearly elucidate the competitive regime of PDE application relative to ramjets and gas turbines.

Povinelli, Louis A.

Petrogenesis of Lunar Rocks: Rb-Sr Constraints and Lack of H2O

Rb and Sr isotopic data and other chemical data indicate major lunar differentiation at about 4.6 AE (AE = 10 9 years) and very limited subsequent differentiation. The constraints of limited differentiation after 4.6 AE and the apparent lack of H20 on the Moon, when applied to the derivation and petrogenesis of lunar samples, suggest the following: (1) soil samples, breccias, metaclastic rocks, and feldspathic basalts represent mixtures of repeatedly modified clastic material, which was ultimately derived from materials formed during the ~ 4.6 AE differentiation; and (2) mare basalts crystallized from melts which formed by partial melting, and which developed without equilibration between the melt and crystalline residuum.

Arden L Albee

Thermal Cycle and Capacitance Testing of Mezzo Industries Phase Change Material Heat Capacitor

A phase change material heat capacitor prototype, designed and built under a NASA Small Business Innovation Research grant by Mezzo Technologies, was tested through coolant thermal cycling within the expected flow rates and temperature range of the Orion spacecraft’s propylene glycol-water mixture coolant system. Testing was performed with the phase change material contained under sealed conditions after a degassing procedure and exposed to atmosphere without prior degassing. The useful heat storage capacitance was measured and compared to Orion spacecraft requirements and theoretical N-Pentadecane storage capacity.

Propylene Glycol Water

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

stagnation heat flux

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

inviscid theory

Boundary Layer Analysis of Shock Tube Flows

Shock tubes offer a controlled environment to reproduce kinetic and radiative phenomena characteristic of atmospheric entry flows under ground-test conditions. The boundary layer developing behind the incident shock wave determines the available test time, influences particle residence times important for similarity scaling, and can affect radiative energy transport. In this study, we couple a quasi-1D space marcher with the compressible boundary-layer equations to numerically compute the post-shock flow in a shock frame of reference for various test gas mixtures representative of different planetary atmospheres. The solvers are individually verified against CFD simulations and analytical correlations available in the literature. Coupled solutions are computed for a finite-rate chemistry in the boundary layer and a non-catalytic isothermal wall. Results yield refined estimates of the maximum separation distance, as well as insights into concentration profiles of relevant species within the boundary layer.

Andrea Fagnani

The Effects of Gravity on Combustion and Structure Formation During Combustion Synthesis in Gasless Systems

There have been relatively few publications examining the role of gravity during combustion synthesis (CS), mostly involving thermite systems. The main goal of this research was to study the influence of gravity on the combustion characteristics of heterogeneous gasless systems. In addition, some aspects of microstructure formation processes which occur during gasless CS were also studied. Four directions for experimental investigation have been explored: (1) the influence of gravity force on the characteristic features of heterogeneous combustion wave propagation (average velocity, instantaneous velocities, shape of combustion front); (2) the combustion of highly porous mixtures (with porosity greater than that for loose powders), which cannot be obtained in normal gravity; (3) the effect of gravity on sample expansion during combustion, in order to produce highly porous materials under microgravity conditions; and (4) the effect of gravity on the structure formation mechanism during the combustion synthesis of poreless composite materials.

Arvind Varma