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 739 records · Page 41

Thermal coupling-decoupling mechanism of heat transfer across van der Waals interfaces in n-eicosane

Progress toward improving thermal transport properties of n-eicosane, a promising phase change material for thermal energy storage near human body temperature, is hindered by lack of a general theoretical framework to model thermal interfacial conductance (TIC) of the van der Waals (vdW) force bonded molecular interfaces, which have little compositional/structural mismatch. Combining molecular dynamics simulations at temperatures up to the melting point and a time-domain analysis of interfacial heat currents in individual molecules, we unveil that (1) the heat flux across a molecular vdW interface has an alternating coupling-decoupling pattern at the single molecule level due to thermal fluctuations, and (2) the net heat transfer during the thermally coupled periods is orders of magnitude larger than that during the thermally decoupled periods. Furthermore, we demonstrate that a vdW interface's TIC correlates almost linearly with two atomistic-scale coupling-decoupling parameters: αc and σc, that represent the duration and strength of the thermal coupling, respectively. This thermal coupling-decoupling mechanism proposes a new paradigm that explicitly accounts for the effects of structural flexibility and dynamical fluctuation on heat transfer across molecular interfaces, and will likely provide new insights into modeling thermal transport properties in a wide-ranging soft materials, including self-assembled monolayers, paraffins or lipids.

42 ENGINEERING↗

Heat transfer model for moving packed-bed particle-to-sCO 2 heat exchangers integrated with metal foams

Particle-to-supercritical carbon dioxide (sCO 2 ) heat exchangers (HXs) play a vital role in coupling heat transfer fluid (HTF) from high-temperature thermal receivers to power cycle working fluids (WF). Heat transfer enhancement is essential for adopting particle-based moving packed-bed heat exchangers (MPBHXs) in next-generation thermal energy storage (TES) systems, as MPBHXs usually exhibit low particle bed-to-wall heat transfer coefficients. High-porosity metal foams have shown their effectiveness in heat transfer enhancement. Here, this work presents a continuum heat transfer model to demonstrate the heat transfer enhancement in MPBHXs when the particle bed channel is filled with high porosity metal foams compared to open channel MPBHXs. The presence of metal foams increases the effective thermal conductivity in the particle channel and enhances the interstitial heat transfer coefficient between the moving particle bed and the stationary metal foams. The present model considers coupled two-dimensional (2D) heat transfer in the particle channel with metal foams and 1D heat transfer in the sCO 2 channel and the dividing wall. The temperature profiles for the particle bed, metal foam, dividing wall, and sCO 2 stream, as well as the local heat flux profiles between those, are studied in detail for various foam porosities, from which the particle bed-to-wall heat transfer coefficient, overall heat transfer coefficient, and total heat transfer rate for the MPBHX are determined. The effects of major MPBHX design and operating parameters on the particle bed-to-wall heat transfer coefficient and overall heat exchange capacity are thoroughly examined; thus, the MPBHX performance improvement with metal foams is quantified for each case. The present heat transfer model can provide valuable insights into the metal-foam MPBHX design and optimization, scale-up, and operating parameters selection.

42 ENGINEERING↗

An analytical approach to thermal modeling of Bridgman-type crystal growth. I - One-dimensional analysis

A relatively simple one-dimensional thermal model of the Bridgman growth process has been developed which is applicable to the growth of small diameter samples with conductivities similar to those of metallic alloys. The heat flow in a translating rod is analyzed in a way that is applicable to Biot numbers less than unity. The model accommodates an adiabatic zone, different heat transfer coefficients in the hot and cold zones, and changes in sample material properties associated with phase change. The analysis is applied to several simplified cases. The effect of the rod's motion is studied in a three-zone furnace for a rod sufficiently long that end effects can be neglected; end effects are then investigated for a motionless rod. Finally, the addition of a fourth zone, an independently controlled booster heater between the main heater and the adiabatic zone, is evaluated for its ability to increase the gradient in the sample at the melt interface and to control the position of the interface.

Naumann, R. J.↗

Liquid metal micro heat pipes for space radiator applications

Micromachining is a chemical means of etching three-dimensional structures, typically in single-crystalline silicon. These techniques are leading toward what is coming to be referred to as MEMS (micro electro mechanical systems), where in addition to the ordinary two dimensional (planar) microelectronics, it is possible to build three-dimensional micromotors, electrically-actuated microvalves, hydraulic systems, and much more on the same microchip. These techniques become possible because of differential etching rates of various crystallographic planes and materials used for semiconductor microfabrication. The University of Cincinnati group in collaboration with NASA Lewis formed micro heat pipes in silicon by the above techniques. Work is ongoing at a modest level, but several essential bonding and packaging techniques have been recently developed. Currently, we have constructed and filled water/silicon micro heat pipes. Preliminary thermal tests of arrays of 125 micro heat pipes etched in a 1 inch x 1 inch x 250 micron silicon wafer have been completed. These pipes are instrumented with extremely small P-N junctions to measure their effective conductivity and their maximum operating power. A relatively simple one-dimensional model has been developed in order to predict micro heat pipes' operating characteristics. This information can be used to optimize micro heat pipe design with respect to length, hydraulic diameter, and number of pipes. Work is progressing on the fabrication of liquid-metal micro heat pipes. In order to be compatible with liquid metal (sodium or potassium), the inside of the micro heat pipes will be coated with a refractory metal (such as tungsten, molybdenum, or titanium).

Gerner, F. M.↗

Development of a Computational Chemical Vapor Deposition Model: Applications to Indium Nitride and Dicyanovinylaniline

A computational chemical vapor deposition (CVD) model is presented, that couples chemical reaction mechanisms with fluid dynamic simulations for vapor deposition experiments. The chemical properties of the systems under investigation are evaluated using quantum, molecular and statistical mechanics models. The fluid dynamic computations are performed using the CFD-ACE program, which can simulate multispecies transport, heat and mass transfer, gas phase chemistry, chemistry of adsorbed species, pulsed reactant flow and variable gravity conditions. Two experimental setups are being studied, in order to fabricate films of: (a) indium nitride (InN) from the gas or surface phase reaction of trimethylindium and ammonia; and (b) 4-(1,1)dicyanovinyl-dimethylaminoaniline (DCVA) by vapor deposition. Modeling of these setups requires knowledge of three groups of properties: thermodynamic properties (heat capacity), transport properties (diffusion, viscosity, and thermal conductivity), and kinetic properties (rate constants for all possible elementary chemical reactions). These properties are evaluated using computational methods whenever experimental data is not available for the species or for the elementary reactions. The chemical vapor deposition model is applied to InN and DCVA. Several possible InN mechanisms are proposed and analyzed. The CVD model simulations of InN show that the deposition rate of InN is more efficient when pulsing chemistry is used under conditions of high pressure and microgravity. An analysis of the chemical properties of DCVA show that DCVA dimers may form under certain conditions of physical vapor transport. CVD simulations of the DCVA system suggest that deposition of the DCVA dimer may play a small role in the film and crystal growth processes.

Cardelino, Carlos↗

Modeling the seasonal cycle of CO2 on Mars: A fit to the Viking lander pressure curves

We have constructed a more accurate Mars thermal model, similar to the one used by Leighton and Murray in 1966, which solves radiative, conductive, and latent heat balance at the surface as well as the one-dimensional heat conduction equation for 40 layers to a depth of 15 meters every 1/36 of a Martian day. The planet is divided into 42 latitude bands with a resolution of two degrees near the poles and five degrees at lower latitudes, with elevations relative to the 6.1 mbar reference areoid. This estimate of the Martian zonally averaged topography was derived primarily from radio occultations. We show that a realistic one-dimensional thermal model is able to reproduce the VL1 pressure curve reasonably well without having to invoke complicated atmospheric effects such as dust storms and polar hoods. Although these factors may cause our deduced values for each model parameter to differ from its true value, we believe that this simple model can be used as a platform to study many aspects of the Martian CO2 cycle over seasonal, interannual, and long-term climate timescales.

Wood, S. E.↗

Thermal math model analysis of FRSI test article subjected to cold soak and entry environments

A multi-objective test program was conducted at the NASA/JSC Radiant Heat Test Facility in which an aluminum skin/stringer test panel insulated with FRSI (Flexible Reusable Surface Insulation) was subjected to 24 simulated Space Shuttle Orbiter ascent/entry heating cycles with a cold soak in between in the 10th and 20th cycles. A two-dimensional thermal math model was developed and utilized to predict the thermal performance of the FRSI. Results are presented which indicate that the modeling techniques and property values have been proven adequate in predicting peak structure temperatures and entry thermal responses from both an ambient and cold soak condition of an FRSI covered aluminum structure.

Gallegos, J. J.↗

Effective thermal conductivity of composites with fibre-matrix debonding

Debonding of the fiber-matrix interface is a major cause for the degradation of the mechanical properties and the loss of thermal conductivity of fiber-reinforced composites. This paper discusses two analytical approaches for modeling the thermal conduction problem of composites. One is based on the concept of modeling the thermal barrier by an equivalent heat transfer coefficient at the fiber-matrix interface, as described by Hasselman and Johnson (1987) and Benveniste and Miloh (1986). The other approach, suggested by Hatta and Taya (1986), is by treating a composite with debonded interface as a coated-fiber composite. The major advantage of the latter aproach is that the thickness of the fiber coating can be realistically modeled depending upon the extent of degradation of the composite with the thermal conductivity of the coating as that of air.

Fadale, T. D.↗

Temperature and density structure of the corona and inner solar wind

A combined theoretical and observational approach is used to construct models of the inner corona applicable to quiet open-field regions and coronal holes. These models, extending from the top of the transition region at about 1.003 solar radii out to 3 solar radii, are calculated by employing parameters that describe the coronal geometry and heating. Thermal conductivity, heating, radiative losses, and solar-wind convection are included in the full energy-balance equation, which is integrated to obtain density and temperature models; the unknown heating term is parameterized in terms of the mechanical flux incident at the base and the characteristic scale height for the dissipation of this flux; and boundary conditions at the base are fixed by EUV-derived models of the transition region. Two representative models, one for a coronal hole and one for a quiet region, are examined in detail as illustrative examples of the technique used. It is noted that the models obtained predict most of the pertinent observations reasonably well but significantly overestimate the radio emission at 80 and 160 MHz and significantly underestimate the intensity of the Fe XV line at 284 A, especially in coronal holes.

Kopp, R. A.↗

Building Model Calibration: Validation of GridLAB-D Thermal Dynamics Modeling

This report investigates how well GridLAB-D’s house model characterizes the thermal dynamics of buildings given the overpredicted diurnal electric load swing observed in the Distribution System Operation with Transactive (DSO+T) study. This study seeks to validate GridLAB-D’s house model by calibrating it to data from the well-instrumented Pacific Northwest National Laboratory Lab Homes in Richland, WA. The datasets chosen included multiple pre-cooling and pre-heating testing periods where the indoor air temperature was allowed to float over a multi-hour length of time to mimic diurnal behavior. The multi-season calibrations were evaluated by comparing the heating/cooling electric power, indoor air temperature, and the rise/decay time during temperature float periods with Lab Homes data. The default GridLAB-D assumptions for latent load fraction, air heat capacity, mass heat capacity, window-to-wall-ratio, overall envelope conductance (assumed as NORMAL thermal integrity level), and solar heat gain coefficient were compared with the calibrated model to confirm when the default assumptions were adequate and the impact of calibration on the accuracy of modeling the thermal dynamics of homes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Controls on the CO2 seasonal cycle

Surface pressure measurement performed by the Viking landers show substantial variations in pressure on seasonal timescales that are characterized by two local minima and two local maxima. These variations have widely been attributed to the seasonal condensation and sublimation of CO2 in the two polar regions. It has been somewhat of a surprise that the amplitude of the minimum and maximum that is dominated by the CO2 cycle in the north was much weaker than the corresponding amplitude of the south-dominated extrema. Another surprise was that the seasonal pressure cycle during years 2 and 3 of the Viking mission was so similar to that for year 1, despite the occurrence of two global dust storms during year 1 and none during years 2 and 3. An energy balance model that incorporates dynamical factors from general circulation model (GCM) runs in which the atmospheric dust opacity and seasonal date were systematically varied was used to model the observed seasonal pressure variations. The energy balance takes account of the following processes in determining the rates of CO2 condensation and sublimation at each longitudinal and latitudinal grid point: solar radiation, infrared radiation from the atmosphere and surface, subsurface heat conduction, and atmospheric heat advection. Condensation rates are calculated both at the surface and in the atmosphere. In addition, the energy balance model also incorporates information from the GCM runs on seasonal redistribution of surface pressure across the globe. Estimates of surface temperature of the seasonal CO2 caps were used to define the infrared radiative losses from the seasonal polar caps. The seasonal pressure variations measured at the Viking lander sites were closely reproduced.

Pollack, J. B.↗

Combustion hot section technology

The overall objective of the Turbine Engine Hot Section Technology Combustion Project is to develop and verify improved and more accurate analysis methods for increasing the ability to design with confidence the combustion system for advanced aircraft turbine engines. The analysis methods developed will be generically applicable to combustion systems and not restricted to one specific engine or manufacturer. This project's approach was to first assess and evaluate existing combustor aerothermal analysis models by means of a contracted effort initiated during FY 1982. This evaluation effort has assessed and quantified known models' strengths and deficiencies. During FY 1984 the Aerothermal Modeling Program, Phase 2 will be initiated, which is expected to have contracted model development efforts in the areas of improved numerical methods for turbulent viscous flows, flow interactions, and fuel spray flow foekd interactions. A Phase 3 effort is planned to address remaining model deficiencies. The primary inhouse effort in this area will be the determination of high pressure flame radiation characteristics in a full annular combustor. This experiment will be conducted in the NASA LeRC High Pressure Facility with the results compiled into a comprehensive flame radiation and liner heat flux model.

Ercegovic, D. B.↗

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of non-thermal electrons and x-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick HI, CaII and MgII transitions using the non-LTE prescription in the program MULTI. Hydrogen ionization rates from x-ray photo-ionization and collisional ionization by non-thermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring.

Hawley, Suzanne L.↗

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane-parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of nonthermal electrons and X-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick H I, Ca II and Mg II transitions using the non-local thermodynamic equilibrium (non-LTE) prescription in the program MULTI. Hydrogen ionization rates from X-ray photoionization and collisional ionization by nonthermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in the pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring. We have computed the chromospheric radiation that results from a range of coronal heating rates, with particular emphasis on the widely observed diagnostic H(alpha). Our conclusion is that the H(alpha) fluxes and profiles actually observed in flares can only be produced under conditions of a low-pressure corona with strong beam heating. Therefore we suggest that H(alpha) in flares is produced primarily at the footprints of newly heated loops where significant evaporation has not yet occurred. As a single loop evolves in time, no matter how strong the heating rate may become, the H(alpha) flux will diminish as the corona becomes denser and hence more effective at stopping the beam. This prediction leads to several observable consequences regarding the spatial and temporal signatures of the X-ray and H(alpha) radiation during flares.

Hawley, Suzanne L.↗

Continuous stabilization of polyacrylonitrile (PAN) - carbon nanotube (CNT) fibers by Joule heating

We report stabilization is one of the important processes used to convert polyacrylonitrile (PAN) fibers into carbon fiber and is usually driven by energy-intensive and time-consuming convective heating. A continuous and energy-efficient alternative is desirable. Here, a continuous Joule heating process has been used for the first time to stabilize PAN fibers containing carbon nanotubes (CNT). Single-component PAN/CNT fibers and bi-component fibers with core-sheath geometry were fabricated through the gel spinning method. Joule heating was conducted continuously on a 200 filament bundle for single-component fibers and on bundles with up to 6000 filaments for the bi-component fibers. A model of heat transfer and thermodynamics occurring during Joule heating is developed and utilized to understand the reductions in stabilization time and energy consumed when compared to convective heating. The minimum energy needed for stabilization through Joule heating is estimated to be less than 1% of the energy used through convective heating.

42 ENGINEERING↗

Site specific porosity-thermal performance correlations in neutron irradiated U-10Zr fuel

In this work, we present a site-specific three-dimensional analysis of porosity evolution in neutron-irradiated U–10Zr metallic fuel using high-resolution synchrotron X-ray tomography. Focused ion beam cubic lift-outs from four radial positions—fuel center, middle, edge, and fuel-cladding interaction (FCCI) zones—were imaged, reconstructed, and segmented to quantify pore volume, density, morphology, and pore connectivity. Porosity increased modestly from 5.5% at the center to 10.5% at the edge, yet pore number density increased by over two orders of magnitude in the fuel portion near the FCCI interface (from 4.8 x 10 4 to 6.0 x 10 6 pores/mm 3 ). Morphological classification revealed a progression from small spherical pores at the center to equiaxed and tortuous networks at the periphery, with enhanced orientation along the radial direction. Connectivity and permeability analysis reveal that the FCCI region maintains dense, highly interconnected pores despite a moderate volume fraction, enabling rapid fission gas transport and lanthanide migration. Effective thermal conductivity models incorporating sodium logging confirm that these site-specific pore features critically influenced heat transport during reactor irradiation. These findings demonstrate that pore topology—not just porosity fraction—controls thermal performance and thermal pathways in U–10Zr fuels, impacting fuel thermal performance in a reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗