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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Influence of Cooling Channel Geometry on the Thermal Response in Silicon Nitride Plates Studied

Engine manufacturers are continually attempting to improve the performance and efficiency of internal combustion engines. Usually they raise the operating temperature or reduce the cooling air requirement for the hot section turbine components. However, the success of these attempts depends on finding materials that are lightweight, are strong, and can withstand high temperatures. Ceramics are among the top candidate materials considered for such harsh applications. They hold low-density, high-temperature strength, and thermal conductivity, and they are undergoing investigation as potential materials for replacing nickel-base alloys and superalloys that are currently used for engine hot-section components. Ceramic structures can withstand higher operating temperatures and a harsh combustion environment. In addition, their low densities relative to metals help reduce component mass. The long-term objectives of the High Temperature Propulsion Components (HOTPC) Project are to develop manufacturing technology, thermal and environmental barrier coatings (TBC/EBC), and the analytical modeling capability to predict thermomechanical stresses in minimally cooled silicon nitride turbine nozzle vanes under simulated engine conditions. Two- and three-dimensional finite element analyses with TBC were conducted at the NASA Glenn Research Center. Nondestructive evaluation was used to determine processing defects. The study included conducting preliminary parametric analytical runs of heat transfer and stress analyses under steady-state conditions to demonstrate the feasibility of using cooled Si3N4 parts for turbine applications. The influence of cooling-channel shapes (such as circular, square, and ascending-order cooling channels) on cooling efficiency and thermal stresses was investigated. Temperature distributions were generated for all cases considered under both cooling and no-cooling conditions, with air being the cooling medium. The table shows the magnitude of the maximum and minimum temperature obtained for the plates under air cooling. Each channel's cross-sectional shape delivered a different temperature; however, the two-dimensional analyses for circular and square or equal-side rectangular holes produced close results. Moreover, the model of the panel with ascending order cooling channels experienced the lowest temperature. A difference of near 260 C was found among the three cooling-hole configurations investigated. The ascending-order cooling channels arrangement showed superior performance by attaining the lowest temperature (1077 C) in comparison to the circular (1379 C) and square (1343 C) channels for the same cooling-hole size. This indicates that the panel with ascending-order cooling channels is the most suitable configuration regardless of the complexity involved in its manufacture. More details pertaining to this study are reported.

Abdul-Aziz, Ali↗
NASA NTRS2021

Goddard Latent Heating Retrieval Algorithm for TRMM and GPM

The Goddard Convective-Stratiform Heating (CSH) algorithm has been used to retrieve latent heating (LH) associated with clouds and cloud systems in support of the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) mission. The CSH algorithm requires the use of a cloud-resolving model (CRM) to simulate LH profiles to build look-up tables (LUTs). In this paper we describe the current V6 CSH and its differences/similarities versus the previous V5 CSH. Long-term CRM simulations are conducted to identify the impact of CRM resolution and the convective-stratiform separation method on LH structure/profiles. The TRMM and GPM Combined algorithm derived surface rainfall rates and their associated precipitation properties are the input to the CSH algorithm. Since there is no direct measurement of LH structure, the performance of the CSH V6 algorithm is examined by comparing its vertically-integrated heating (or equivalent surface rainfall rate) against the surface rainfall rate derived from the TRMM/GPM Combined algorithm. The results suggest the CSH three-month and zonal mean equivalent surface rain rates are in good agreement with the Combined rates over the ITCZ region; the agreement is best over ocean. CSH three-month and zonal mean equivalent surface rain rates are larger than the Combined rates over land in both the tropics and sub-tropics. CSH three-month mean equivalent surface rainfall rates also have local differences with the Combined rates that can be smoothed by area averaging to larger horizontal resolutions (from the CSH standard grid of 0.25° × 0.25° to 0.5° × 0.5° or 1.0° × 1.0°). CSH equivalent surface rainfall rates have more light rain rates but less larger rain rates compared to the GPM Combined surface rain rates.

Latent heating retrieval↗
NASA NTRS1995

Flow in serpentine coolant passages with trip strips

Under the subject contract, an effort is being conducted at Scientific Research Associates, Inc. (SRA) to obtain flow field measurements in the coolant passage of a rotating turbine blade with ribbed walls, both in the stationary and rotating frames. The data obtained will be used for validation of computational tools and assessment of turbine blade cooling strategies. The configuration of the turbine blade passage model is given, and the measuring plane locations are given. The model has a four-pass passage with three 180 turns. This geometry was chosen to allow analyses of the velocity measurements corresponding to the heat transfer results obtained by Wagner. Two passes of the passage have a rectangular cross-section of 1.0 in x 0.5 in. Another two passes have a square cross-section of 0.5 in x 0.5 in. Trips with a streamwise pitch to trip height (P/e) = 5 and trip height to coolant passage width (e/Z) = 0.1, were machined along the leading and trailing walls. These dimensions are typical of those used in turbine blade coolant passages. The trips on these walls are staggered by the half-pitch. The trips are skewed at +/- 45 deg, and this allows the effect of trip orientation to be examined. Experiments will be conducted with flow entering the model through the 1.0 in x 0.5 in rectangular passage (Configuration C) and the 0.5 in x 0. 5 in square passage (Configuration D) to examine the effect of passage aspect ratio. Velocity measurements were obtained with a Reynolds number (Re) of 25,000, based on the hydraulic diameter of and bulk mean velocity in the half inch square passage. The coordinate system used in presenting the results for configurations C and D, respectively, is shown. The first, second and third passes of the passage will be referred to as the first, second and third passages, respectively, in later discussion. Streamwise distance (x) from the entrance is normalized by the hydraulic diameter (D). Vertical (y) and tangential (z) distances are normalized by the half passage height (H) and width (Z), respectively. The x coordinate and U component are positive in the streamwise direction. The y coordinate and V component are positive against gravity. The z coordinate and W component are positive in the direction of rotation. The velocities are normalized by the bulk mean velocity (Ub) of 3.44 m/s based on the half-inch square passage. The contours of the 1.0 in x 0.5 in and 0.5 in x 0.5 in passages were evaluated from 11 x 30 and 9 x 30 measurement grids, respectively.

Tse, D. G.-N.↗
NASA NTRS1994

A realistic model of evaporation for a liquid droplet

An intuitive delineation along with dimensional considerations and experimental evidences are presented to show that in a general case, the evaporation of a liquid droplet undergoes three regimes through the process. Initially, the heat transfer inside the evaporating droplets is conduction controlled; then, in the second stage, convective heat transfer may take over; finally, the convections subside, and the process returns to conduction controlled mode.

Chai, An-Ti↗
NASA NTRS2001

Augmentation of Stagnation Region Heat Transfer Due to Turbulence from a DLN Can Combustor

Heat transfer measurements have been made in the stagnation region of a flat plate with a circular leading edge. Electrically heated aluminum strips placed symmetrically about the leading edge stagnation region were used to measure spanwise-averaged heat transfer coefficients. The maximum Reynolds number obtained, based on leading edge diameter, was about 100,000. The model was immersed in the flow field downstream of an approximately half-scale model of a can-type combustor from a low NO(x), ground-based power-generating turbine. The tests were conducted with room temperature air; no fuel was added. Room air flowed into the combustor through six vane-type fuel/air swirlers. The combustor can contained no dilution holes. The fuel/air swirlers all swirled the incoming airflow in a counterclockwise direction (facing downstream). A five-hole probe flow field survey in the plane of the model stagnation point showed the flow was one big vortex with flow angles up to 36 deg at the outer edges of the rectangular test section. Hot-wire measurements showed test section flow had very high levels of turbulence, around 28.5%, and had a relatively large axial-length scale-to-leading edge diameter ratio of 0.5. X-wire measurements showed the turbulence to be nearly isotropic. Stagnation heat transfer augmentation over laminar levels was around 77% and was about 14% higher than predicted by a previously developed correlation for isotropic grid-generated turbulence.

VanFossen, G. James↗
NASA NTRS1997

Experimental and Computational Studies of the Flow Over a Sting Mounted Planetary Probe Configuration

This paper summarizes the results of a series of experimental studies in the LENS shock tunnel and computations with DSMC and Navier Stokes codes which have been made to examine the aerothermal and flowfield characteristics of the flow over a sting-supported planetary probe configuration in hypervelocity air and nitrogen flows. The experimental program was conducted in the LENS hypervelocity shock tunnel at total enthalpies of 5and 10 MJkg for a range of reservoir pressure conditions from 70 to 500 bars. Heat transfer and pressure measurements were made on the front and rear face of the probe and along the supporting sting. High-speed and single shot schlieren photography were also employed to examine the flow over the model and the time to establish the flow in the base recirculation region. Predictions of the flowfield characteristics and the distributions of heat transfer and pressure were made with DSMC codes for rarefied flow conditions and with the Navier-Stokes solvers for the higher pressure conditions where the flows were assumed to be laminar. Analysis of the time history records from the heat transfer and pressure instrumentation on the face of the probe and in the base region indicated that the base flow was fully established in under 4 milliseconds from flow initiation or between 35 and 50 flow lengths based on base height. The measurements made in three different tunnel entries with two models of identical geometries but with different instrumentation packages, one prepared by NASA Langley and the second prepared by CUBRC, demonstrated good agreement between heat transfer measurements made with two different types of thin film and coaxial gage instrumentation. The measurements of heat transfer and pressure to the front face of the probe were in good agreement with theoretical predictions from both the DSMC and Navier Stokes codes. For the measurements made in low density flows, computations with the DSMC code were found to compare well with the pressure and heat transfer measurements on the sting, although the computed heat transfer rates in the recirculation region did not exhibit the same characteristics as the measurements. For the 10MJkg and 500 bar reservoir match point condition, the measurements and heat transfer along the sting from the first group of studies were in agreement with the Navier Stokes solutions for laminar conditions. A similar set of measurements made in later tests where the model was moved to a slightly different position in the test section indicated that the boundary layer in the reattachment compression region was close to transition or transitional where small changes in the test environment can result in larger than laminar heating rates. The maximum heating coefficients on the sting observed in the present studies was a small fraction of similar measurements obtained at nominally the same conditions in the HEG shock tunnel, where it is possible for transition to occur in the base flow, and in the low enthalpy studies conducted in the NASA Langley high Reynolds number Mach 10 tunnel where the base flow was shown to be turbulent. While the hybrid Navier- StokedDMSC calculations by Gochberg et al. (Reference 1) suggested that employing the Navier- Stokes calculations for the entire flowfield could be seriously in error in the base region for the 10 MJkg, 500 bar test case, similar calculations performed by Cornell, presented here, do not.

Holden, Michael S.↗
NASA NTRS1984

Megaregolith thickness, heat flow, and the bulk composition of the Moon

Models developed to assess the effects of megaregolith on lunar thermal evolution are discussed. It is confirmed that the two sites where lunar heat flow was measured are probably unrepresentative, with heat flows about 25% higher than regional averages, due to focussing of heat flow towards regions with thin megaregolith. Numerous lines of evidence indicate that the megaregolith is generally 2 to 3 km thick under highlands (which cover about 83% of the total lunar surface), and 1 km thick under maria. In most models, megaregolith thickness is assumed to be roughly 6x greater over highlands than over maria. Based on sparse data for porosity among lunar rock types, and the correlation between thermal conductivity and porosity, it is assumed that megaregolith conductivity is roughly 20 kiloerg s(-1)cm(-1)K(-1), and bedrock conductivity is roughly 7x greater. It is also found that insulation by megaregolith exacerbates the problem of reconciling modest temperatures inferred for the (present) matle with a high rate of heat production; an upper limit of 30 ng/g for the bulk-Moon U content can be derived from this constraint alone.

Warren, P. H.↗
NASA NTRS1977

An implicit-iterative solution of the heat conduction equation with a radiation boundary condition

For the problem of predicting one-dimensional heat transfer between conducting and radiating mediums by an implicit finite difference method, four different formulations were used to approximate the surface radiation boundary condition while retaining an implicit formulation for the interior temperature nodes. These formulations are an explicit boundary condition, a linearized boundary condition, an iterative boundary condition, and a semi-iterative boundary method. The results of these methods in predicting surface temperature on the space shuttle orbiter thermal protection system model under a variety of heating rates were compared. The iterative technique caused the surface temperature to be bounded at each step. While the linearized and explicit methods were generally more efficient, the iterative and semi-iterative techniques provided a realistic surface temperature response without requiring step size control techniques.

Williams, S. D.↗
NASA NTRS1977

A design handbook for phase change thermal control and energy storage devices

Comprehensive survey is given of the thermal aspects of phase change material devices. Fundamental mechanisms of heat transfer within the phase change device are discussed. Performance in zero-g and one-g fields are examined as it relates to such a device. Computer models for phase change materials, with metal fillers, undergoing conductive and convective processes are detailed. Using these models, extensive parametric data are presented for a hypothetical configuration with a rectangular phase change housing, using straight fins as the filler, and paraffin as the phase change material. These data are generated over a range of realistic sizes, material properties, and thermal boundary conditions. A number of illustrative examples are given to demonstrate use of the parametric data. Also, a complete listing of phase change material property data are reproduced herein as an aid to the reader.

Humphries, W. R.↗
NASA NTRS1984

Solar Pond Fluid Dynamics and Heat Transfer

The primary objective of the solar pond research was to obtain an indepth understanding of solar pond fluid dynamics and heat transfer. The key product was the development of a validated one-dimensional computer model with the capability to accurately predict time-dependent solar pond temperature, salinities, and interface motions. Laboratory scale flow visualization experiments were conducted to better understand layer motion. Two laboratory small-scale ponds and a large-scale outdoor solar pond were designed and built to provide quantitative data. This data provided a basis for validating the model and enhancing the understanding of pond dynamic behavior.

Jones, G. F.↗
NASA NTRS1989

Coolant side heat transfer with rotation: User manual for 3D-TEACH with rotation

This program solves the governing transport equations in Reynolds average form for the flow of a 3-D, steady state, viscous, heat conducting, multiple species, single phase, Newtonian fluid with combustion. The governing partial differential equations are solved in physical variables in either a Cartesian or cylindrical coordinate system. The effects of rotation on the momentum and enthalpy calculations modeled in Cartesian coordinates are examined. The flow of the fluid should be confined and subsonic with a maximum Mach number no larger than 0.5. This manual describes the operating procedures and input details for executing a 3D-TEACH computation.

Syed, S. A.↗
NASA NTRS1991

Extinction of a Free Methanol Droplet in Microgravity

The combustion of isolated, unsupported methanol droplets under microgravity conditions has been studied in the NASA-Lewis 2.2-sec drop tower. The substitution of O2/He mixtures for O2/N2 mixtures results in a substantial increase in the oxygen index required for ignition, the rate of gasification, and the diameter at extinction. As a result, droplets with initial size larger than 1 mm can be produced, deployed, ignited, and burned to extinction, all under microgravity conditions in this facility. A comprehensive time dependent numerical model incorporating detailed multicomponent molecular transport mechanisms and complex elementary chemical kinetics was developed and applied to simulate the entire transient droplet combustion processes (ignition/quasi-steady burning/extinction). Model outputs show good agreement with the experimental measurements, without parameter adjustment. Analysis of modeling results shows that the thermal conductivity of the gas phase under combustion environments increased by a factor of two when the 50-percent O2/50 percent He mixture is used instead of air. The increased thermal conductivity accelerates the heat transfer from the flame front to the droplet surface, resulting in a factor of two increase in the burning rate.

Cho, S. Y.↗
NASA NTRS1993

The differential emission measure of nested hot and cool magnetic loops

The detailed thermal structure of the magnetized solar transition region, as measured by its differential emission measure (DEM(T)), is poorly known. Building on the fact that the solar surface is strongly magnetized and thereby structured, proposals have been made that envision a significant lower-temperature contribution to the energy balance from (ion) heat flux across an arcade of different temperature loops. In this paper, we describe a self-consistent 2D MHD simulation, which includes the full thermal effects of parallel stability and anisotropic conduction, of a nested-loop model of the thermal and magnetic structure of the transition region. We then demonstrate that the predicted DEM agrees with observations in the conceptually elusive T less than 10 exp 5 K regime.

Van Hoven, G.↗
NASA NTRS2000

Analysis of Aluminum-Nitride SOI for High-Temperature Electronics

We use numerical simulation to investigate the high-temperature (up to 500K) operation of SOI MOSFETs with Aluminum-Nitride (AIN) buried insulators, rather than the conventional silicon-dioxide (SiO2). Because the thermal conductivity of AIN is about 100 times that of SiO2, AIN SOI should greatly reduce the often severe self-heating problem of conventional SOI, making SOI potentially suitable for high-temperature applications. A detailed electrothermal transport model is used in the simulations, and solved with a PDE solver called PROPHET In this work, we compare the performance of AIN-based SOI with that of SiO2-based SOI and conventional MOSFETs. We find that AIN SOI does indeed remove the self-heating penalty of SOL However, several device design trade-offs remain, which our simulations highlight.

Biegel, Bryan A.↗
NASA NTRSDate not supplied

Design and Comparison of Lunar Sample Return Options

Returning lunar samples to earth is a high priority for NASA in order to better understand Lunar Geology, Lunar formation, and potential resource utilization. Core samples of frozen volatiles from permanently shadowed regions are sought after, but require intensive cryogenic cooling to preserve. Returning samples to Earth using the Orion capsule is restrictive for moderate amounts of mass/volume and thus alternate means of transport are investigated in this paper. The trade-space of potential options includes varying: sample mass amount, sample containment conditions, transfer path from Lunar surface to Earth surface, and various vehicle system trades. An integrated solver is applied with parametric estimators and known system heuristics to effectively estimate size, mass, and power of various subsystems. Modeled subsystems include sample containment, EDL technologies(heat shield, parachute, recovery), power, thermal, avionics, structures, and propulsive systems. Transfer path options include various trajectory options, transfer methods, staging points, and vehicle staging options. A Design of Experiments was conducted to explore the tradespace, to identify sensitivities, and to identify promising designs. Results are detailed for various stakeholder objectives including total mass, sample mass, science return, architecture impacts, and operational flexibility.

Lunar outpost↗
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions to the finite element problem, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
NASA NTRS1985

Megaregolith thickness, heat flow, and the bulk composition of the moon

Models developed to assess the effects of megaregolith on lunar thermal evolution are discussed. It is confirmed that the two sites where lunar heat flow was measured are probably unrepresentative, with heat flows about 25 percent higher than regional averages, due to focussing of heat flow towards regions with thin megaregolith. Numerous lines of evidence indicate that the megaregolith is generally 2 to 3 km thick under highlands (which cover about 83 percent of the total lunar surface), and 1 km thick under maria. In most models, megaregolith thickness is assumed to be roughly 6x greater over highlands than over maria. Based on sparse data for porosity among lunar rock types, and the correlation between thermal conductivity and porosity, it is assumed that megaregolith conductivity is roughly 20 kiloerg s(-1)cm(-1)K(-1), and bedrock conductivity is roughly 7x greater. It is also found that insulation by megaregolith exacerbates the problem of reconciling modest temperatures inferred for the (present) mantle with a high rate of heat production; an upper limit of 30 ng/g for the bulk-Moon U content can be derived from this constraint alone.

Rasmussen, K. L.↗