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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Thermal Behaviour of Unusual Local-Scale Surface Features on Vesta

On Vesta, the region of the infrared spectrum beyond approximately 3.5 micrometers is dominated by the thermal emission of the asteroid's surface, which can be used to determine surface temperature by means of temperature-retrieval algorithms. The thermal behavior of areas of unusual albedo seen at the local scale can be related to physical properties that can provide information about the origin of those materials. Dawn's Visible and Infrared Mapping Spectrometer (VIR) hyperspectral cubes are used to retrieve surface temperatures, with high accuracy as long as temperatures are greater than 180 K. Data acquired in the Survey phase (23 July through 29 August 2011) show several unusual surface features: 1) high-albedo (bright) and low-albedo (dark) material deposits, 2) spectrally distinct ejecta, 3) regions suggesting finer-grained materials. Some of the unusual dark and bright features were re-observed by VIR in the subsequent High-Altitude Mapping Orbit (HAMO) and Low-Altitude Mapping Orbit (LAMO) phases at increased pixel resolution. To calculate surface temperatures, we applied a Bayesian approach to nonlinear inversion based on the Kirchhoff law and the Planck function. These results were cross-checked through application of alternative methods. Here we present temperature maps of several local-scale features that were observed by Dawn under different illumination conditions and different local solar times. Some bright terrains have an overall albedo in the visible as much as 40% brighter than surrounding areas. Data from the IR channel of VIR show that bright regions generally correspond to regions with lower thermal emission, i.e. lower temperature, while dark regions correspond to areas with higher thermal emission, i.e. higher temperature. This behavior confirms that many of the dark appearances in the VIS mainly reflect albedo variations. In particular, it is shown that during maximum daily insolation, dark features in the equatorial region may rise to temperatures greater than 270 K. However, individual features may show different thermal behaviours, as a result of differences in composition and/or structure (e.g, average grain size of the surface regolith, porosity, etc.). To complement the temperature and near-infrared emissivity derived from the infrared spectra, a separate work is devoted to calculate thermal inertia and other thermal properties using theoretical models which solve the heat equation for airless bodies, and model the distribution of temperatures due to surface roughness variations.

Tosi, F.↗
NASA NTRS2017

Pressurization of a Flightweight, Liquid Hydrogen Tank: Evaporation and Condensation at a Liquid Vapor Interface

Evaporation and condensation at a liquidvapor interface is important for long-term, in-space cryogenic propellant storage. Yet the current understanding of interfacial physics does not predict behavior or evaporation condensation rates. The proposed paper will present a physical model, based on the 1-D Heat equation and Schrages equation which demonstrates thin thermal layers at the fluidvapor interface.

Cryogenic hydrogen↗
NASA NTRS2017

Pressurization of a Flightweight, Liquid Hydrogen Tank: Evaporation and Condensation at a Liquid Vapor Interface

Evaporation and condensation at a liquid-vapor interface is important for long-term, in-space cryogenic propellant storage. Yet the current understanding of inter-facial physics does not consistently predict behavior of evaporation or condensation rates. The proposed paper will present a physical model, based on the 1-D Heat equation and Schrage's equation, which demonstrates thin thermal layers at the fluid vapor interface.

cryogenic hydrogen↗
NASA NTRS2017

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗
NASA NTRS1992

A near-wall two-equation model for turbulent heat fluxes

A near-wall two-equation model for turbulent heat fluxes is derived from the temperature variance and its dissipation-rate equations and the assumption of gradient transport. Only incompressible flows with non-buoyant heat transfer are considered. The near-wall asymptotics of each term in the exact equations are examined and used to derive near-wall correction functions that render the modeled equations consistent with these behavior. Thus modeled, the equations are used to calculate fully-developed pipe and channel flows with heat transfer. It is found that the proposed two-equation model yields asymptotically correct near-wall behavior for the normal heat flux, the temperature variance and its near-wall budget and correct limiting wall values for these properties compared to direct simulation data and measurements obtained under different wall boundary conditions.

Sommer, T. P.↗
NASA NTRS1988

Regressed relations for forced convection heat transfer in a direct injection stratified charge rotary engine

Currently, the heat transfer equation used in the rotary combustion engine (RCE) simulation model is taken from piston engine studies. These relations have been empirically developed by the experimental input coming from piston engines whose geometry differs considerably from that of the RCE. The objective of this work was to derive equations to estimate heat transfer coefficients in the combustion chamber of an RCE. This was accomplished by making detailed temperature and pressure measurements in a direct injection stratified charge (DISC) RCE under a range of conditions. For each specific measurement point, the local gas velocity was assumed equal to the local rotor tip speed. Local physical properties of the fluids were then calculated. Two types of correlation equations were derived and are described in this paper. The first correlation expresses the Nusselt number as a function of the Prandtl number, Reynolds number, and characteristic temperature ratio; the second correlation expresses the forced convection heat transfer coefficient as a function of fluid temperature, pressure and velocity.

Lee, Chi M.↗
NASA NTRS1988

Regressed relations for forced convection heat transfer in a direct injection stratified charge rotary engine

Currently, the heat transfer equation used in the rotary combustion engine (RCE) simulation model is taken from piston engine studies. These relations have been empirically developed by the experimental input coming from piston engines whose geometry differs considerably from that of the RCE. The objective of this work was to derive equations to estimate heat transfer coefficients in the combustion chamber of an RCE. This was accomplished by making detailed temperature and pressure measurements in a direct injection stratified charge (DISC) RCE under a range of conditions. For each specific measurement point, the local gas velocity was assumed equal to the local rotor tip speed. Local physical properties of the fluids were then calculated. Two types of correlation equations were derived and are described in this paper. The first correlation expresses the Nusselt number as a function of the Prandtl number, Reynolds number, and characteristic temperature ratio; the second correlation expresses the forced convection heat transfer coefficient as a function of fluid temperature, pressure and velocity.

Lee, Chi M.↗
NASA NTRS2004

Sensitivity Equation Derivation for Transient Heat Transfer Problems

The focus of the paper is on the derivation of sensitivity equations for transient heat transfer problems modeled by different discretization processes. Two examples will be used in this study to facilitate the discussion. The first example is a coupled, transient heat transfer problem that simulates the press molding process in fabrication of composite laminates. These state equations are discretized into standard h-version finite elements and solved by a multiple step, predictor-corrector scheme. The sensitivity analysis results based upon the direct and adjoint variable approaches will be presented. The second example is a nonlinear transient heat transfer problem solved by a p-version time-discontinuous Galerkin's Method. The resulting matrix equation of the state equation is simply in the form of Ax = b, representing a single step, time marching scheme. A direct differentiation approach will be used to compute the thermal sensitivities of a sample 2D problem.

Hou, Gene↗
DOE OSTI2020

Multiphase representation of hotspot formation and growth

Multiphase reacting flow approach has been adopted for modeling of hotspot formation and growth on sub-zonal scale. Mechanical response of the high explosive (HE) is modeled by the Visco-Plastic SCRAM model which includes plastic work, crack growth, and frictional heating along the cracks. In this work a sub-zonal regions with high temperature, hotspots, are modeled as a new phase. Formation and growth of the hotspot phase from the cold phase is modeled as a phase change. Each phase has its own temperature and pressure. Chemical reactions also progress in both phases. The growth of hotspots are determined by the solution of the heat conduction equation with appropriate heat sources due to friction and plastic heating. Sample calculations using the BABBO code show successful detonation calculations with hotspot formation and growth.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1973

An equation of state for oxygen and nitrogen

Recent measurements of thermodynamic properties of oxygen and nitrogen have provided data necessary for development of a single equation of state for both fluids. Data are available in summary report and two-part detailed study on thermodynamic properties of oxygen and nitrogen. Same data are used to develop vapor-pressure equation and heat-capacity equation.

Jacobsen, R. T.↗
NASA NTRS1981

Numerical treatment of nonlinear latent heat boundary conditions at moving interfaces in genuine two dimensional solidification problems

The proposed method for the treatment of two dimensional solidification problems is based on quasilinearization of the transformed heat conduction equation and latent heat condition at the interface and an iterative sequence in which these are solved simultaneously. Modern algorithms for solving such sparse systems mean that most of the storage advantage of other methods are reduced and the speed of solution can be improved.

Beckett, P. M.↗
NASA NTRS1993

Computational solution of the defect stream-function equation for nonequilibrium turbulent boundary layers

The derivation of the accurate, second-order, almost linear, approximate equation governing the defect stream function for nonequilibrium compressible turbulent boundary layers is reviewed. The similarity of this equation to the heat conduction equation is exploited in the development of an unconditionally stable, tridiagonal computational method which is second-order accurate in the marching direction and fourth-order accurate in the surface-normal direction. Results compare well with experimental data. Nonlinear effects are shown to be small. This two-dimensional method is simple and has been implemented on a programmable calculator.

Barnwell, Richard W.↗
NASA NTRS2001

Thermal Conductivity of Ceramic Thermal Barrier and Environmental Barrier Coating Materials

Thermal barrier and environmental barrier coatings (TBC's and EBC's) have been developed to protect metallic and Si-based ceramic components in gas turbine engines from high temperature attack. Zirconia-yttria based oxides and (Ba,Sr)Al2Si2O8(BSAS)/mullite based silicates have been used as the coating materials. In this study, thermal conductivity values of zirconia-yttria- and BSAS/mullite-based coating materials were determined at high temperatures using a steady-state laser heat flux technique. During the laser conductivity test, the specimen surface was heated by delivering uniformly distributed heat flux from a high power laser. One-dimensional steady-state heating was achieved by using thin disk specimen configuration (25.4 mm diam and 2 to 4 mm thickness) and the appropriate backside air-cooling. The temperature gradient across the specimen thickness was carefully measured by two surface and backside pyrometers. The thermal conductivity values were thus determined as a function of temperature based on the 1-D heat transfer equation. The radiation heat loss and laser absorption corrections of the materials were considered in the conductivity measurements. The effects of specimen porosity and sintering on measured conductivity values were also evaluated.

Zhu, Dong-Ming↗
NASA NTRS2006

Heat Transfer In High-Temperature Multilayer Insulation

The combined radiation/conduction heat transfer in high-temperature multilayer insulations for typical reentry of reusable launch vehicles from low Earth orbit was investigated experimentally and numerically. The high-temperature multilayer insulation investigated consisted of gold-coated reflective foils separated by alumina fibrous insulation spacers. The steady-state heat transfer through four multilayer insulation configurations was investigated experimentally over the temperature range of 300-1300 K and environmental pressure range of 1.33 10(exp -5)-101.32 kPa. It was shown that including the reflective foils reduced the effective thermal conductivity compared to fibrous insulation sample at 1.5 times the density of the multilayer sample. A finite volume numerical model was developed to solve the governing combined radiation/conduction heat transfer equations. The radiation heat transfer in the fibrous insulation spacers was modeled using the modified two-flux approximation assuming anisotropic scattering and gray medium. The numerical model was validated by comparison with steady-state experimental data. The root mean square deviation between the predicted and measured effective thermal conductivity of the samples was 9.5%.

Daryabeigi, Kamran↗