Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Heat Transfer in High Temperature Multilayer Insulation

High temperature multilayer insulations have been investigated as an effective component of thermal-protection systems for atmospheric re-entry of reusable launch vehicles. Heat transfer in multilayer insulations consisting of thin, gold-coated, ceramic reflective foils and Saffil(TradeMark) fibrous insulation spacers was studied both numerically and experimentally. A finite volume numerical thermal model using combined conduction (gaseous and solid) and radiation in porous media was developed. A two-flux model with anisotropic scattering was used for radiation heat transfer in the fibrous insulation spacers between the reflective foils. The thermal model was validated by comparison with effective thermal conductivity measurements in an apparatus based on ASTM standard C201. Measurements were performed at environmental pressures in the range from 1x10(exp -4) to 760 torr over the temperature range from 300 to 1300 K. Four multilayer samples with nominal densities of 48 kg/cu m were tested. The first sample was 13.3 mm thick and had four evenly spaced reflective foils. The other three samples were 26.6 mm thick and utilized either one, two, or four reflective foils, located near the hot boundary with nominal foil spacing of 1.7 mm. The validated thermal model was then used to study relevant design parameters, such as reflective foil spacing and location in the stack-up and coating of one or both sides of foils.

Daryabeigi, Kamran↗
DOE OSTI2024

Analytical Modeling of Metal Foam Composite Phase Change Materials (PCM) in Thermal Energy Storage Using Asymptotic Analysis

The use of phase change materials (PCMs) for thermal energy storage can release or absorb a significant amount of latent heat during the freezing or melting process, offering a higher energy storage density. One of the main drawbacks of PCMs is their low thermal conductivity, resulting in poor thermal performance. Recent research has attempted to enhance heat transfer and increase the thermal conductivity of PCMs, including the use of metal foams. However, modeling the metal foam composite PCM using conventional methods is computationally expensive. This paper proposes an asymptotic solution for a Stefan-like problem subject to a convective boundary for outward solidification in a hollow cylinder, capable of predicting the freeze-melt cycle of the metal foam composite PCM. Specifically, three temporal regimes and four spatial layers are considered in the asymptotic analysis for each phase change process. The thermal conductivity is calculated by a theoretical three-dimensional tetrakaidecahedron model, while other thermophysical properties are obtained using the method of volume averaging. The results are verified with numerical data and validated against experimental data in the literature. The presented analytical modeling framework could have the potential to be applied to other types of composite PCMs with considerably lower computational costs compared with conventional methods.

analytical model↗
NASA NTRS1985

Methods for heat transfer and temperature field analysis of the insulated diesel phase 2 progress report

This report describes work done during Phase 2 of a 3 year program aimed at developing a comprehensive heat transfer and thermal analysis methodology for design analysis of insulated diesel engines. The overall program addresses all the key heat transfer issues: (1) spatially and time-resolved convective and radiative in-cylinder heat transfer, (2) steady-state conduction in the overall structure, and (3) cyclical and load/speed temperature transients in the engine structure. During Phase 2, radiation heat transfer model was developed, which accounts for soot formation and burn up. A methodology was developed for carrying out the multi-dimensional finite-element heat conduction calculations within the framework of thermodynamic cycle codes. Studies were carried out using the integrated methodology to address key issues in low heat rejection engines. A wide ranging design analysis matrix was covered, including a variety of insulation strategies, recovery devices and base engine configurations. A single cylinder Cummins engine was installed at Purdue University, and it was brought to a full operational status. The development of instrumentation was continued, concentrating on radiation heat flux detector, total heat flux probe, and accurate pressure-crank angle data acquisition.

Morel, T.↗
NASA NTRS1975

Expansion of chromospheric matter in the gradual phase of solar flares

Interferometric observations at 17 GHz of several small X-ray flares are presented along with soft X-ray observations of preflare active regions to show that a large mass increase accompanies the formation of an X-ray hot region in the corona. The total amount of energy contained in a hot coronal region is estimated, and a model is proposed in which a significant amount of the hot matter is supplied to the corona from the chromosphere during each flare. According to this model, energy produced by some coronal instability is transported by thermal conduction to the chromosphere, where dense gas is heated and subsequently expands into the corona. It is shown that impulsive heating of the chromosphere by nonthermal electrons cannot be the energy source of this model because the total energy supplied to the hot region during the gradual phase must be much greater than that supplied during the impulsive phase.

Ohki, K.↗
NASA NTRS1992

Status of solar wind modeling from the transition region outwards

In recent years, solar wind modeling has to some extent undergone a shift of emphasis, from attempts to produce high-speed streams far from the sun, to investigations of what conditions must exist in the solar corona and transition region in order to produce the observed conditions in both the solar wind and low corona. Thus there has been an increased awareness that the solar wind should really be treated as part of the solar atmosphere, and that the problems associated with heating and accelerating the solar wind should be treated in concert with the coronal (and perhaps chromospheric) heating problems. We will discuss several models which take this point of view but we will place particular emphasis on some outstanding problems, viz. the mass flux problem, some puzzling recent IPS data, our persistent difficulties with electron heat conduction, and observational uncertainties about the coronal and transition region boundary conditions which should be put into the models. We will conclude by suggesting some possible alternatives for future models.

Hollweg, J. V.↗
NASA NTRS1977

Flare heating by energetic nonthermal electrons

A model has been developed which includes, in a self-consistent fashion, the energy degradation and attendant bremsstrahlung emission of a high-energy electron beam, heating of a flare plasma by the beam, and subsequent cooling by thermal conduction and radiation. To assist in the interpretation of experimental diagnostics, the model also characterizes the radiative behavior of iron ions present in the flare plasma. Electron deposition is described by the Fokker-Planck equation for an initial power-law particle distribution. Results are presented for the bremsstrahlung radiation emitted by the incident beam as it impinges on the disturbed atmosphere. A comparison is made between direct beam heating and thermal-conduction heating of the flare plasma. The radiation emitted by several selected spectral lines of Fe ions is determined as a function of time during and after deposition

Davis, J.↗
NASA NTRS1975

Heat transfer tests on a 0.01-scale Rockwell configuration 3 space shuttle orbiter and tank (37-OT) in the Calspan 48-inch hypersonic shock tunnel (OH12/IH21), volume 1

Model information and data are presented from wind tunnel tests conducted on 0.01-scale models of the space shuttle orbiter and external tank. These tests were conducted in a hypersonic shock tunnel to determine heating rates on ascent and reentry configurations at various Reynolds numbers, Mach numbers, and angles of attack.

Kotch, M.↗
DOE OSTI2021

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2012

Prevention of Over-Pressurization During Combustion in a Sealed Chamber

The combustion of flammable material in a sealed chamber invariably leads to an initial pressure rise in the volume. The pressure rise is due to the increase in the total number of gaseous moles (condensed fuel plus chamber oxygen combining to form gaseous carbon dioxide and water vapor) and, most importantly, the temperature rise of the gas in the chamber. Though the rise in temperature and pressure would reduce with time after flame extinguishment due to the absorption of heat by the walls and contents of the sealed spacecraft, the initial pressure rise from a fire, if large enough, could lead to a vehicle over-pressure and the release of gas through the pressure relief valve. This paper presents a simple lumped-parameter model of the pressure rise in a sealed chamber resulting from the heat release during combustion. The transient model considers the increase in gaseous moles due to combustion, and heat transfer to the chamber walls by convection and radiation and to the fuel-sample holder by conduction, as a function of the burning rate of the material. The results of the model are compared to the pressure rise in an experimental chamber during flame spread tests as well as to the pressure falloff after flame extinguishment. The experiments involve flame spread over thin solid fuel samples. Estimates of the heat release rate profiles for input to the model come from the assumed stoichiometric burning of the fuel along with the observed flame spread behavior. The sensitivity of the model to predict maximum chamber pressure is determined with respect to the uncertainties in input parameters. Model predictions are also presented for the pressure profile anticipated in the Fire Safety-1 experiment, a material flammability and fire safety experiment proposed for the European Space Agency (ESA) Automated Transfer Vehicle (ATV). Computations are done for a range of scenarios including various initial pressures and sample sizes. Based on these results, various mitigation approaches are suggested to prevent vehicle over-pressurization and help guide the definition of the space experiment.

Gokoglu, Suleyman A.↗
DOE OSTI2023

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose–Einstein statistics into an “effective temperature” in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

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

Investigation of Effectiveness of Air-Heating a Hollow Steel Propeller for Protection Against Icing. 1: Unpartitioned Blades

An investigation to determine the effectiveness of icing protection afforded by air-heating hollow steel unpartitioned propeller blades has been conducted In the NACA Cleveland icing research tunnel. The propeller used was a production model modified with blade shank and tip openings to permit internal passage of heated air. Blade-surface and heated-air temperatures were obtained and photographic observations of Ice formations were made with variations In icing intensity and heating rate to the blades. For the conditions of Icing to which the propeller was subjected, it was found that adequate ice protection was afforded with a heating rate of 40 1 000 Btu per hour per blade. With less than 40,000 Btu per hour per blade, ice protection failed because of significant ice accretions on the leading edge. The chordwise distribution of heat was unsatisfactory with most of the available heat dissipated well back of the leading edge on both the thrust and camber face's instead of at the leading edge where it was most needed. A low utilization of available heat for icing protection is indicated by a beat-exchanger effectiveness of approximately 47 percent.

Mulholland, Donald R.↗
NASA NTRS1975

Induced anisotropy of thermal conductivity of polymer solids under large strains

A simple analytical form of induced anisotropy of heat conductivity of initially isotropic polymer solids results from employing the simplified theory of the three-chain model of the non-Gaussian network. The analytical form appears to be valid up to a stretch ratio of 2.65, which is the limit of existing experimental data. The effect of induced anisotropy on the temperature distribution, due to the large deformations, is illustrated for a highly expanded spherical shell and a cylindrical tube under a steady-state heat flow using the derived analytical form of the strain-dependent heat conductivity.

Peng, S. T. J.↗
DOE OSTI2021

Predicting Thermal Excursions During in situ Oxidative Regeneration of Packed Bed Catalytic Fast Pyrolysis Catalyst

Ex situ catalytic fast pyrolysis (CFP) uses a secondary reactor to upgrade biomass pyrolysis vapors to stabilized CFP oils with reduced oxygen content. In one configuration, the secondary reactor is operated as a packed-bed swing reactor system which allows coke-deactivated beds to be decarbonized in situ while other beds remain online for vapor upgrading. In situ decarbonization must be done carefully to avoid irreversible deactivation and/or physical degradation of catalyst pellets. Given that packed bed reactors are well known to have poor heat transfer characteristics, this is a critical issue impacting scaleability and commercial viability of the technology. To predict thermal excursions during regeneration, finite element computational models have been built to assist in scaling up oxidative decarbonization of a Pt/TiO2 CFP catalyst (0.5 mm spheres) from a bench scale packed bed with 100 g of catalyst to a pilot scale packed bed with 2 kg of catalyst and internal cooling tubes. Based on transient measurements of outlet temperature and effluent CO2 concentration, and using an assumed coke profile and activation energy, this paper demonstrates that specific combinations of effective thermal conductivity and wall heat transfer coefficient can fit bench scale oxidative regeneration data equally well. For the upscaled 2 kg bed, four bench-scale “best fit” parameter pairs give different predictions for location and magnitude of thermal excursions, with the maximum computed bed temperature gradients ranging from 30 °C cm-1 to as high as 3000 °C cm-1. The larger the fraction of heat removal by conduction through the cooling tubes, the greater the differences between the parameter pairs. The modelling results presented in this paper cast doubt on the industrial viability of the proposed combination of catalyst, bed and regeneration process, and point to the need for alternate reactor designs. However, there is considerable uncertainly in some of the key model parameters. The reliability of model predictions can be increased by adding more temperature measurements at key bed locations, testing additional variations in process conditions, performing careful bed dissections to determine the true coke profile, and perhaps most importantly, directly measuring the effective thermal conductivity of the catalyst pellets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the Volume of Fluid (VOF) and Sharp Interface multiphase 2D axisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. Turbulence modeling; turbulence damping at the interface; and using constant vs. temperature dependent fluid properties were shown to have the most profound influence on predicted tank pressures and temperatures. Current study indicates that including phase change at the interface into the computational model is crucial for accurate prediction of the tank self-pressurization process. The effect of accommodation coefficient was also studied. Tank pressure values predicted by the VOF, and Sharp Interface models are within 4% of the experimental ones.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-Phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the VOF and Sharp Interface multiphase 2Daxisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. These parameters include: the effect of turbulence modeling and turbulence damping at the interface; the effect of constant vs. temperature dependent fluid properties; the effect of accommodation coefficient; the effect of modeling phase change at the interface and the effect of the fill level.

Computational Fluid Dynamics↗
DOE OSTI2024

Thermally anisotropic building envelope for thermal management: finite element model calibration using field evaluation data

The thermally anisotropic building envelope (TABE) is an active building envelope that redistributes thermal loads in response to weather conditions and building energy demand. Conductive layers throughout the TABE distribute low-grade heat among hydronic loops, altering heat flow direction and intensity. Finite element models of TABE roof and wall panels were developed and calibrated using field evaluation data. The calibration results showed that heat flux differences between the experimental data and finite element models averaged –0.42% and 3.57%, with a maximum mean square error of 1.78 and 3.96 for roof and wall panels, respectively. A reduction in heat flux from the environment to the building living space over the entire testing period (weeks in July/August) was found to be 85% for roof panels and 335% (load reversed) for wall panels. Finally, these results indicate TABE can effectively harness low-grade thermal energy sources to achieve high energy efficiency and promote demand-side management.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1989

Comparison of the measured and predicted response of the Earth Radiation Budget Experiment active cavity radiometer during solar observations

The Earth Radiation Budget Experiment consists of an array of radiometric instruments placed in earth orbit by the National Aeronautics and Space Administration to monitor the longwave and visible components of the earth's radiation budget. Presented is a dynamic electrothermal model of the active cavity radiometer used to measure the earth's total radiative exitance. Radiative exchange is modeled using the Monte Carlo method and transient conduction is treated using the finite element method. Also included is the feedback circuit which controls electrical substitution heating of the cavity. The model is shown to accurately predict the dynamic response of the instrument during solar calibration.

Mahan, J. R.↗
NASA NTRS2001

Afterbody Heating Characteristics of a Proposed Mars Sample Return Orbiter

Aeroheating wind-tunnel tests were conducted on a 0.028 scale model of an orbiter concept considered for a possible Mars sample return mission. The primary experimental objectives were to characterize hypersonic near wake closure and determine if shear layer impingement would occur on the proposed orbiter afterbody at incidence angles necessary for a Martian aerocapture maneuver. Global heat transfer mappings, surface streamline patterns, and shock shapes were obtained in the NASA Langley 20-inch Mach 6 Air and CF4 Tunnels for post-normal shock Reynolds numbers (based on forebody diameter) ranging from 1,400 to 415,000, angles of attack ranging from -5 to 10 degrees at 0, 3, and 6 deg sideslip, and normal-shock density, ratios of 5 and 12. Laminar, transitional, and turbulent shear layer impingement on the cylindrical afterbody was inferred from the measurements and resulted in a localized heating maximum that ranged from 40 to 75% of the reference forebody stagnation point heating. Comparison of laminar heating prediction to experimental measurement along the orbiter afterbody highlight grid alignment challenges associated with numerical simulation of three-dimensional separated wake flows.

Horvath, Thomas J.↗