Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heat equation”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Redesigned Human Metabolic Simulator

A design has been formulated for a proposed improved version of an apparatus that simulates atmospheric effects of human respiration by introducing controlled amounts of carbon dioxide, water vapor, and heat into the air. Denoted a human metabolic simulator (HMS), the apparatus is used for testing life-support equipment when human test subjects are not available. The prior version of the HMS, to be replaced, was designed to simulate the respiratory effects of as many as four persons. It exploits the catalytic combustion of methyl acetate, for which the respiratory quotient (the molar ratio of carbon dioxide produced to oxygen consumed) is very close to the human respiratory quotient of about 0.86. The design of the improved HMS provides for simulation of the respiratory effects of as many as eight persons at various levels of activity. The design would also increase safety by eliminating the use of combustion. The improved HMS (see figure) would include a computer that would exert overall control. The computer would calculate the required amounts of oxygen removal, carbon dioxide addition, water addition, and heat addition by use of empirical equations for metabolic profiles of respiration and heat. A blower would circulate air between the HMS and a chamber containing a life-support system to be tested. With the help of feedback from a mass flowmeter, the blower speed would be adjusted to regulate the rate of flow according to the number of persons to be simulated and to a temperature-regulation requirement (the air temperature would indirectly depend on the rate of flow, among other parameters). Oxygen would be removed from the circulating air by means of a commercially available molecular sieve configured as an oxygen concentrator. Oxygen, argon, and trace amounts of nitrogen would pass through a bed in the molecular sieve while carbon dioxide, the majority of nitrogen, and other trace gases would be trapped by the bed and subsequently returned to the chamber. If, as recommended, the oxygen concentrator were of a rotating twelve-bed design, then variations in the product stream could be made very small. Carbon dioxide would be added directly to the circulating air by simple injection from a supply tank. The rate of injection would be maintained at the required rate by use of a mass flowmeter/controller. In the same way, nitrogen would be added to make up for the small amount of nitrogen lost through the oxygen concentrator. Water vapor would be added to the circulating air by heating the corresponding required flow of water to steam in a heat exchanger. More heat, required to complete the simulation of the thermal effect of respiration, would be added through another heat exchanger. Heat would be supplied to both heat exchangers via a hot-oil loop.

Duffield, Bruce↗
DOE OSTI2024

A Reactor Scale-Up Methodology from Lab-Scale to Pilot-Scale Operations: Numerical Modeling of THFA Dehydration to DHP in Packed-Bed Reactors

This manuscript discusses developing a model-based scale-up methodology for a successful technology transfer of gas-phase catalytic reactors from a lab-scale to a pilot-scale operation. The manuscript demonstrates the methodology for gas-phase dehydration of tetrahydrofurfuryl alcohol (THFA) to dihydropyran (DHP) process over commercial Al 2 O 3 catalysts. A two-dimensional reactor model was developed using COMSOL Multiphysics 6.1 software. The model solves heat and mass transport equations in bed-scale and particle scales simultaneously. This powerful feature enables accurate prediction of the heat and mass transfer limitations in pilot-scale reactors, if any exists. Further, the model uses isothermal lab-scale experimental data to derive and validate the reaction chemistry, flow fields and boundary conditions. The model was then scaled-up to project conversion, selectivity, yield and formation rate of DHP in a pilot-scale reactor. The results highlight the complex nature of chemistry, heat, and mass transfer effects in lab-scale and pilot-scale reactors. The model results inform the possible operational limitations of the pilot-scale reactor and design strategies to improve process efficiency. Although the scale-up approach is explained through the THFA dehydration process, the methodology is applicable to any catalytic packed-bed reactor models for a successful process scale-up.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1962

Numerical Analysis of the Transient Response of Advanced Thermal Protection Systems for Atmospheric Entry

Equations for the transfer of heat through thermal protection shields are derived in finite difference form. These equations are applicable to charring ablators, impregnated ceramics, subliming ablators, heat sinks, and insulating materials and have been programed for solution on a high-speed digital computer. In the program, thermal properties can be functions of temperature. Provision is made for analysis of heat shields subjected to simultaneous convective and radiative heat inputs. Some typical results are presented. Limited comparisons with experimental results are made.

THERMAL PROTECTION↗
DOE OSTI2021

The interpretation of temperature and salinity variables in numerical ocean model output and the calculation of heat fluxes and heat content

Abstract. The international Thermodynamic Equation of Seawater 2010 (TEOS-10) defined the enthalpy and entropy of seawater, thus enabling the global ocean heat content to be calculated as the volume integral of the product of in situ density, ρ, and potential enthalpy, h0 (with reference sea pressure of 0 dbar). In terms of Conservative Temperature, Θ, ocean heat content is the volume integral of ρcp0Θ, where cp0 is a constant “isobaric heat capacity”. However, many ocean models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) as well as all models that contributed to earlier phases, such as CMIP5, CMIP3, CMIP2, and CMIP1, used EOS-80 (Equation of State – 1980) rather than the updated TEOS-10, so the question arises of how the salinity and temperature variables in these models should be physically interpreted, with a particular focus on comparison to TEOS-10-compliant observations. In this article we address how heat content, surface heat fluxes, and the meridional heat transport are best calculated using output from these models and how these quantities should be compared with those calculated from corresponding observations. We conclude that even though a model uses the EOS-80, which expects potential temperature as its input temperature, the most appropriate interpretation of the model's temperature variable is actually Conservative Temperature. This perhaps unexpected interpretation is needed to ensure that the air–sea heat flux that leaves and arrives in atmosphere and sea ice models is the same as that which arrives in and leaves the ocean model. We also show that the salinity variable carried by present TEOS-10-based models is Preformed Salinity, while the salinity variable of EOS-80-based models is also proportional to Preformed Salinity. These interpretations of the salinity and temperature variables in ocean models are an update on the comprehensive Griffies et al. (2016) paper that discusses the interpretation of many aspects of coupled Earth system models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

Equatorial response to northeasterly cold surges as inferred from satellite cloud imagery

Daily convective activity in equatorial regions during Dec.-Jan., 1974-1975 was estimated from NOAA satellite data and composited in relation to the phase of cold surge episodes which then occured. Eastward and westward propagating cloud clusters displayed phase speeds of about 5 and 10 m/sec, respectively, and were generated near southwestern Borneo after a cold surge. Kelvin wave responses were indicated by confinement of the eastward-moving components to the equatorial region. Further, an analysis of symmetric and asymmetric characteristics showed that components symmetric with respect to the equator displayed a larger response amplitude, which was taken to mean that strongly trapped waves were generated from heat sources near the equator.

Lau, K.-M.↗
NASA NTRSDate not supplied

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗
NASA NTRS1973

An equation of state for oxygen and nitrogen

Recent measurements of the thermodynamic properties of oxygen and nitrogen have been used to develop a single equation of state for both fluids; with appropriate coefficients, the equation describes the pressure-density-temperature surface for each fluid within the estimated experimental uncertainty of the data. The formulation uses least squares curve fitting of the equation of state to the available experimental data (including pressure-density-temperature data, saturation data, and measured isochoric heat capacities). New vapor pressure equations and equations for the ideal gas heat capacities of oxygen and nitrogen are included for use in calculation of derived thermodynamic properties.

Jacobsen, R. T.↗
NASA NTRS1983

Character and stability of axisymmetric thermal convection in spheres and spherical shells

The influence of shell size and mode of heating on the behavior and stability of axisymmetric, infinite Prandtl number convection in a spherical geometry is studied. Heating from within and below features convection onset governed by a self-adjoint system of equations and boundary conditions. For heating only from within or from below, linearized equations and boundary conditions are non-self-adjoint. Identification of the parameter which initiates the departure from self-adjointness, together with the properties of the self-adjoint solution, provide a basis for calculating the heat transfer characteristics of the non-self-adjoint situations. The investigations are an effort to develop a model for heat transfer in planetary interiors. Further development of the technique by modifying the Galerkin method by the introduction of diagonal mode truncation is suggested to permit the consideration of higher values of the Rayleigh numbers, i.e., those more commensurate with terrestrial planet mantles.

Zebib, A.↗
NASA NTRS1980

Surface temperatures in sliding systems - A finite element analysis

Finite element equations are developed for studying surface temperatures resulting from frictional heating in sliding systems. The equations include the effect of velocity of moving components, an effect which is found to be quite significant, even at low sliding velocities. A program was written using the equations and it was applied to the study of surface temperatures in two different sliding systems: dry or boundary lubricated sleeve bearings and a labyrinth gas path seal configuration. Very good agreement was achieved between analytical predictions using the program and experimental temperature measurements. The program was used to study the influence of various material parameters on surface temperatures in the two sliding systems.

Kennedy, F. E., Jr.↗
NASA NTRS2013

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗
NASA NTRS1976

Heating of the solar chromosphere and corona. I - Generalized inhomogeneous wave equation for magnetoacoustic motions

The generalized inhomogeneous wave equation that governs magnetoacoustic, vortical, and thermal motions in compressible fluids and that is applicable to the problem of heating of the solar chromosphere and corona is obtained. The effects of kinematic and bulk viscosity, heat conduction, Joule dissipation, and magnetic diffusivity are included. Under the usual assumptions, the generalized wave equation reduces to the well-known equations of Lighthill, Kulsrud, Phillips, and others. The major problems encountered in applying Lighthill's (1952) mechanism to sound generation in turbulent media are reviewed for both the subsonic and supersonic cases.

Anand, S. P. S.↗
NASA NTRS2003

Real-Gas Effects on Binary Mixing Layers

This paper presents a computational study of real-gas effects on the mean flow and temporal stability of heptane/nitrogen and oxygen/hydrogen mixing layers at supercritical pressures. These layers consist of two counterflowing free streams of different composition, temperature, and density. As in related prior studies reported in NASA Tech Briefs, the governing conservation equations were the Navier-Stokes equations of compressible flow plus equations for the conservation of total energy and of chemical- species masses. In these equations, the expressions for heat fluxes and chemical-species mass fluxes were derived from fluctuation-dissipation theory and incorporate Soret and Dufour effects. Similarity equations for the streamwise velocity, temperature, and mass fractions were derived as approximations to the governing equations. Similarity profiles showed important real-gas, non-ideal-mixture effects, particularly for temperature, in departing from the error-function profile, which is the similarity solution for incompressible flow. The temperature behavior was attributed to real-gas thermodynamics and variations in Schmidt and Prandtl numbers. Temporal linear inviscid stability analyses were performed using the similarity and error-function profiles as the mean flow. For the similarity profiles, the growth rates were found to be larger and the wavelengths of highest instability shorter, relative to those of the errorfunction profiles and to those obtained from incompressible-flow stability analysis. The range of unstable wavelengths was found to be larger for the similarity profiles than for the error-function profiles

Okong'o, Nora↗
NASA NTRS1998

Heat Transfer on a Film-Cooled Rotating Blade Using a Two Equation Turbulence Model

A three-dimensional Navier-Stokes code has been used to compare the heat transfer coefficient on a film-cooled, rotating turbine blade. The blade chosen is the ACE rotor with five rows containing 93 film cooling holes covering the entire span. This is the only film-cooled rotating blade over which experimental data is available for comparison. Over 2.278 million grid points are used to compute the flow over the blade including the tip clearance region, using Coakley's q-omega turbulence model. Results are also compared with those obtained by Garg and Abhari (1997) using the zero-equation Baldwin-Lomax (B-L) model. A reasonably good comparison with the experimental data is obtained on the suction surface for both the turbulence models. At the leading edge, the B-L model yields a better comparison than the q-omega model. On the pressure surface, however, the comparison between the experimental data and the prediction from either turbulence model is poor. A potential reason for the discrepancy on the pressure surface could be the presence of unsteady effects due to stator-rotor interaction in the experiments which are not modeled in the present computations. Prediction using the two-equation model is in general poorer than that using the zero-equation model, while the former requires at least 40% more computational resources.

Garg, Vijay K.↗
NASA NTRS2010

Analysis of Flame Deflector Spray Nozzles in Rocket Engine Test Stands

The development of a unified tightly coupled multi-phase computational framework is described for the analysis and design of cooling spray nozzle configurations on the flame deflector in rocket engine test stands. An Eulerian formulation is used to model the disperse phase and is coupled to the gas-phase equations through momentum and heat transfer as well as phase change. The phase change formulation is modeled according to a modified form of the Hertz-Knudsen equation. Various simple test cases are presented to verify the validity of the numerical framework. The ability of the methodology to accurately predict the temperature load on the flame deflector is demonstrated though application to an actual sub-scale test facility. The CFD simulation was able to reproduce the result of the test-firing, showing that the spray nozzle configuration provided insufficient amount of cooling.

Sachdev, Jai S.↗