Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Conduction”

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 235 records · Page 13

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗

Development of a laminar boundary layer behind a suction point

A theoretical investigation is made of the development of a laminar boundary layer behind a suction slot that is assumed to cut off part of the boundary layer without exerting any sink effect. The development, which is approximate, is based on the heat conduction equation. The heat conduction equation enters the analysis through a linearization of the Prandtl-Mises form of the boundary-layer equation.

EQUATIONS - BOUNDARY LAYER↗

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗

Thermal Vacuum Testing of a Novel Loop Heat Pipe Design for the Swift BAT Instrument

An advanced thermal control system for the Burst Alert Telescope on the Swift satellite has been designed and an engineering test unit (ETU) has been built and tested in a thermal vacuum chamber. The ETU assembly consists of a propylene loop heat pipe, two constant conductance heat pipes, a variable conductance heat pipe (VCHP), which is used for rough temperature control of the system, and a radiator. The entire assembly was tested in a thermal vacuum chamber at NASA/GSFC in early 2002. Tests were performed with thermal mass to represent the instrument and with electrical resistance heaters providing the heat to be transferred. Start-up and heat transfer of over 300 W was demonstrated with both steady and variable condenser sink temperatures. Radiator sink temperatures ranged from a high of approximately 273 K, to a low of approximately 83 K, and the system was held at a constant operating temperature of 278 K throughout most of the testing. A novel LHP temperature control methodology using both temperature-controlled electrical resistance heaters and a small VCHP was demonstrated. This paper describes the system and the tests performed and includes a discussion of the test results.

Ottenstein, Laura↗

Cooling device featuring thermoelectric and diamond materials for temperature control of heat-dissipating devices

A cooling device for lowering the temperature of a heat-dissipating device. The cooling device includes a heat-conducting substrate (composed, e.g., of diamond or another high thermal conductivity material) disposed in thermal contact with the heat-dissipating device. During operation, heat flows from the heat-dissipating device into the heat-conducting substrate, where it is spread out over a relatively large area. A thermoelectric cooling material (e.g., a Bi.sub.2 Te.sub.3 -based film or other thermoelectric material) is placed in thermal contact with the heat-conducting substrate. Application of electrical power to the thermoelectric material drives the thermoelectric material to pump heat into a second heat-conducting substrate which, in turn, is attached to a heat sink.

Vandersande, Ian W.↗

Investigation of Freeze and Thaw Cycles of a Gas-Charged Heat Pipe

The traditional constant conductance heat pipes (CCHPs) currently used on most spacecraft run the risk of bursting the pipe when the working fluid is frozen and later thawed. One method to avoid pipe bursting is to use a gas-charged heat pipe (GCHP) that can sustain repeated freeze/thaw cycles. The construction of the GCHP is similar to that of the traditional CCHP except that a small amount of non-condensable gas (NCG) is introduced and a small length is added to the CCHP condenser to serve as the NCG reservoir. During the normal operation, the NCG is mostly confined to the reservoir, and the GCHP functions as a passive variable conductance heat pipe (VCHP). When the liquid begins to freeze in the condenser section, the NCG will expand to fill the central core of the heat pipe, and ice will be formed only in the grooves located on the inner surface of the heat pipe in a controlled fashion. The ice will not bridge the diameter of the heat pipe, thus avoiding the risk of pipe bursting during freeze/thaw cycles. A GCHP using ammonia as the working fluid was fabricated and then tested inside a thermal vacuum chamber. The GCHP demonstrated a heat transport capability of more than 200W at 298K as designed. Twenty-seven freeze/thaw cycles were conducted under various conditions where the evaporator temperature ranged from 163K to 253K and the condenser/reservoir temperatures ranged from 123K to 173K. In all tests, the GCHP restarted without any problem with heat loads between 10W and 100W. No performance degradation was noticed after 27 freeze/thaw cycles. The ability of the GCHP to sustain repeated freeze/thaw cycles was thus successfully demonstrated.

Ku, Jentung↗

Field-dependent nonelectronic contributions to thermal conductivity in a metallic ferromagnet with low Gilbert damping

Heat conduction in metals is typically dominated by electron transport since electrons carry both charge and heat. In magnetic metals magnons, or spin waves, excitations of the magnetic order can be used to transport information. Heat conduction via magnons has been previously shown mostly for insulating magnets with low Gilbert damping and resulting long spin-wave lifetimes where conduction electrons cannot contribute. Here we show that thin films of properly optimized metallic ferromagnetic (FM) alloys show significant nonelectronic contributions to heat conduction, which furthermore depend on the direction of an applied magnetic field. These measurements are enabled by micromachined thermal isolation platforms optimized for thermal conductivity measurements of thin-film systems. Electrical conductivity measurements on exactly the same samples allow application of the Wiedemann-Franz relation, which shows large nonelectronic contributions to thermal conductivity for the cobalt-iron alloy with 25% Co. This composition has been shown to have exceptionally low damping for a metallic FM. The thermal conductivity of a 75-nm-thick film of the 25% Co alloy changes by more than 20% at some temperatures, while a reference sample with 50% cobalt that has much higher damping shows no field-direction dependence. Our measurements indicate that applied magnetic fields alter the magnon lifetimes in these films and that these magnons contribute to thermal conductivity in this metallic magnetic alloy with low Gilbert damping.

36 MATERIALS SCIENCE↗

Cool Flame Quenching

Cool flame quenching distances are generally presumed to be larger than those associated with hot flames, because the quenching distance scales with the inverse of the flame propagation speed, and cool flame propagation speeds are often times slower than those associated with hot flames. To date, this presumption has never been put to a rigorous test, because unstirred, non-isothermal cool flame studies on Earth are complicated by natural convection. Moreover, the critical Peclet number (Pe) for quenching of cool flames has never been established and may not be the same as that associated with wall quenching due to conduction heat loss in hot flames, Pe approx. = 40-60. The objectives of this ground-based study are to: (1) better understand the role of conduction heat loss and species diffusion on cool flame quenching (i.e., Lewis number effects), (2) determine cool flame quenching distances (i.e, critical Peclet number, Pe) for different experimental parameters and vessel surface pretreatments, and (3) understand the mechanisms that govern the quenching distances in premixtures that support cool flames as well as hot flames induced by spark-ignition. Objective (3) poses a unique fire safety hazard if conditions exist where cool flame quenching distances are smaller than those associated with hot flames. For example, a significant, yet unexplored risk, can occur if a multi-stage ignition (a cool flame that transitions to a hot flame) occurs in a vessel size that is smaller than that associated with the hot quenching distance. To accomplish the above objectives, a variety of hydrocarbon-air mixtures will be tested in a static reactor at elevated temperature in the laboratory (1g). In addition, reactions with chemical induction times that are sufficiently short will be tested aboard NASA's KC-135 microgravity (mu-g) aircraft. The mu-g results will be compared to a numerical model that includes species diffusion, heat conduction, and a skeletal kinetic mechanism, following the work on diffusion-controlled cool flames by Fairlie et,al., 2000.

Pearlman, Howard↗

Deployable Heat Pipe Radiator

A 1.2- by 1.8-m variable conductance heat pipe radiator was designed, built, and tested. The radiator has deployment capability and can passively control Freon-21 fluid loop temperatures under varying loads and environments. It consists of six grooved variable conductance heat pipes attached to a 0.032-in. aluminum panel. Heat is supplied to the radiator via a fluid header or a single-fluid flexible heat pipe header. The heat pipe header is an artery design that has a flexible section capable of bending up to 90 degrees. Radiator loads as high as 850 watts were successfully tested. Over a load variation of 200 watts, the outlet temperature of the Freon-21 fluid varied by 7 F. An alternate control system was also investigated which used a variable conductance heat pipe header attached to the heat pipe radiator panel.

Edelstein, F.↗

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Description and orbit data of variable-conductance heat-pipe system for the communications technology satellite

A variable-conductance heat-pipe system (VCHPS) with methanol as the working fluid and a nitrogen and helium mixture as the control gas was used for the thermal control of a 200 W RF traveling wave tube of the Communication Technology Satellite. Three stainless steel heat pipes (one redundant) and an aluminum radiator were designed to transfer 196 watts for an evaporator temperature of 50 C. The system has operated for three years with no noticeable change in performance. On four occasions the heat pipes apparently deprimed. A short time after reducing the tube power, the heat pipes reprimed and the system continued to operate normally. The description, qualification testing, and orbit data of the VCHPS are presented.

Gedeon, L.↗

A dynamic thermal algal growth model for pilot-scale open-channel raceways

Here this study describes the makeup of and results from a simulation consisting of a lumped thermal model integrated with a dynamic algae growth model to simulate the microalgae productivity of an open-channel raceway cultivation system. The thermal model considers the dynamic effects of weather, light absorption, convective heat transfer, radiation heat transfer, conductive heat transfer, thermal capacitance, and water control strategies. The dynamic algae growth model solves a set of ordinary differential equations consisting of growth functions dependent on incident radiation, temperature, nutrient availability, basal metabolism, and losses due to dark- and photo-respiration. Relative errors in the predicted ash-free dry weight of Nannochloropsis oceanica are 12.5%, -6.1%, and 4.4% for three separate replicates of cultivation cycles (~4.2 m 2 ) performed at AzCATI during the ATP 3 Unified Field Studies in Fall, Spring, and Summer. This research demonstrates that thermal modeling is an essential contributor to the validation of microalgae growth models.

09 BIOMASS FUELS↗

Insights into Mechanisms of Screw Feeder Plugging by Heated Pine Forestry Residues Using in-situ and Correlative Microscopy

Increasing concerns about global carbon emissions due to fossil fuel use has created interest in renewable fuel alternatives, such as lignocellulosic biomass. With an estimated global biomass production probability ten-fold higher than projected need, functionalizing and scaling biomass conversion processes, like pyrolysis, is a promising component of a diverse sustainable energy program (Kan, Strezov, & Evans, 2016). Pyrolysis reactors employed in biomass conversion often encounter the known problem of screw feeder plugging, halting work and reducing efficiency of biorefineries. Although the screw feeder is not heated directly, conductive heat transfer from the fluidized bed pyrolysis reactor creates a temperature gradient, where particles experience temperatures between 22°C and 500°C from the time they enter the feeder to being emptied into the reactor (Digman, Joo, & Kim, 2009; Nieminen & Folke, 1982). The impact of unintentional heating on this highly variable feedstock may result in particle morphology modification, surface texture transformation, and the evolution of viscous compounds and may contribute to screw feeder plugging (Haas, Nimlos, & Donohoe, 2009; Ingram et al., 2008). We hypothesized that particles of various anatomical fractions (including needles, branches, bark, cambium, and whitewood) will experience dramatic changes during this heating regime, which will be characterized by decreased particle size, increased surface roughness, degradation of key structural macromolecules, and generation of viscous compounds that would result in increased particle friction and cohesion, contributing to screw feeder plugging.

biofuels↗

Uranium Density, Thermal Conductivity, Specific Heat, and Thermal Diffusivity

Thermophysical properties for pure uranium are collected here as a common and maintained resource with a focus on casting. Source material details will be provided for decisions on fitting emphasis. The four material properties represented in this document were chosen based upon the definition of thermal diffusivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

COMSOL Results for the Nominal Steady-State Operation of the Proposed 95-MW LEU Silicide Core for HFIR Conversion

Engineering design studies are being performed to determine the feasibility of converting the High Flux Isotope Reactor (HFIR) from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel at Oak Ridge National Laboratory. This activity is sponsored by the Office of Reactor Conversion and Uranium Supply (ORCUS) under the auspices of the US Department of Energy National Nuclear Security Administration’s Office of Material Management and Minimization. HFIR is a very high flux, pressurized, light water–cooled and moderated, flux trap–type research reactor with a core made of involute shaped U 3 O 8 /Al cermet fuel plates and coolant channels. HFIR currently operates at a thermal power of 85 MW and supports key national and international missions in neutron scattering, isotope production, materials/fuels irradiation, neutron activation analysis, gamma irradiation, and neutrino research. Advanced multiphysics computational fluid dynamics models have been developed in the COMSOL Multiphysics software to simulate the steady-state operating conditions for the proposed low-and high-density LEU U 3 Si 2 -Al (uranium silicide dispersion) fuel designs. The COMSOL models for HFIR inner and outer fuel element models incorporate various essential inputs and physics such as spatially dependent nuclear heat deposition, multilayer heat conduction, conjugate heat transfer, turbulent flows (using Reynolds-averaged Navier Stokes turbulence models), structural mechanics (thermal–structural interactions and fuel swelling), and oxide layer build-up. This report presents the best-estimate thermal hydraulics results for the low- and high-density optimized silicide LEU core designs at 95 MW steady-state nominal operation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Time-Dependent Structure of the Upper Atmosphere

The physical properties of the upper atmosphere are determined mainly by heat conduction, the heat sources and the barometric law. An analysis of the integro-differential equation which describes these physical processes has been carried out. It is found that heating of the thermosphere due to absorption of the solar extreme ultraviolet (EUV) radiation alone cannot explain the observed diurnal variation of density and temperature, since it would yield a maximum of these properties at about 17h local time, instead of 14h where it is observed. Secondly, if the EUV flux is adjusted to give the observed average temperature, then the diurnal variation in density would be much too large compared with the observed amplitude. Thirdly, it would require an extremely high efficiency for the conversion of EUV radiation into heat, if we compare the required flux with Hinteregger's measurements of the EUV flux. Thus, it is necessary to have another heat source in addition to the heating due to absorption of EUV radiation. If an additional beat source is used, which has a maximum at about 9h local time and a flux of 1 erg cm−1 sec−1, a time-dependent model of the upper atmosphere is obtained that is in good agreement with the observed densities. There is evidence that this additional heat source derives its energy ultimately from the solar corpuscular radiation. In this paper we present the results of calculations for a model in the equatorial and temperature zones of the earth, for those times when the average solar activity corresponds to a solar radiation flux of 200×110−22 Wm−2 (cps)−1 at 10.7-cm wavelength. The physical properties (temperature, density, pressure, scale height, mean molecular weight and the number densities of N2, O2, O, He and H) are given as a function of local time and for the altitudes between 120 km and 2050 km.

ATMOSPHERE↗

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.↗