Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heating”

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,153 records · Page 64

Sensitivity and vulnerability to summer heat extremes in major cities of the United States

Many cities are experiencing increases in extreme heat because of global temperature rise combined with the urban heat island effect. The heterogeneity of urban morphology also leads to fine-scale variability in potential for heat exposure. Yet, how this rise in temperature and local variability together impacts urban residents differently at exposure-relevant scales is still not clear. Here we map the Universal Thermal Climate Index, a more complete indicator of human heat stress at an unprecedentedly fine spatial resolution (1 m), for 14 major cities in the United States using urban microclimate modeling. We examined the different heat exposure levels across different socioeconomic and racial/ethnic groups in these cities, finding that income level is most consistently associated with heat stress. We further conducted scenario simulations for a hypothetical 1 °C increase of air temperature in all cities. Results show that a 1 °C increase would have a substantial impact on human heat stress, with impacts that differ across cities. The results of this study can help us better evaluate the impact of extreme heat on urban residents at decision-relevant scales.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Heat Exchanger Technology for Integrated Energy System: Methodology and Preliminary Application

This report proposes a novel heat exchanger evaluation method, based on integrated QFD-AHP technique, to support strategic optimal selection of heat exchanger (HX) technology for integrated energy systems. The term “optimal selection” is defined as selecting a product (i.e., HX) that best meets various aspects of end user demand, including both technical and economic requirements. Compared to previous HX evaluation methods, the proposed approach allows decision makers to explicitly take into account the practical demand of customers (i.e., HX end users) which may vary depending on the HX application target. In addition, the proposed evaluation method can be used not only for high-level decisions to select HX type, but also for low-level decisions such as the selection of the key components including HX material or heat transfer fluid. The focus of this study is to describe the new HX evaluation method and to illustrate the overall HX evaluation and selection processes using the three nuclear heat application cases that represent the thermal integration between advanced nuclear reactors and industrial processes. As a preliminary study, three different types of HX, i.e., shell-and-tube heat exchanger, printed circuit heat exchanger, and heat pipe heat exchanger, were compared.

42 ENGINEERING↗

Robust Heat-Flux Sensors for Coal-Fired Boiler Extreme Environments

In this project, robust heat-flux measurement systems were developed. The heat-flux sensors utilize thermoelectric effects to directly transduce the heat-flux inputs to analog electrical voltage signals. They were constructed from dedicated materials that can withstand temperatures of at least 1000°C and maintain adequate performance at these conditions for prolonged periods of time. The proposed approaches took into account numerous considerations, including system cost, sensor head resilience, sensor footprint, data accuracy, response time, and maintenance requirements. Through modern thermoelectric materials design, methodical materials selection and rigorous testing in materials characterization labs and medium-scale fire research facilities, we have demonstrated functioning laboratory prototypes, upon which one could base industrial heat-flux sensing platforms capable of operating in the challenging high-temperature, corrosive environments of the boilers of coal-fired power plants. A distributed sensor array for heat-flux measurements throughout the furnace water-wall, the superheater area and the economizer coils can provide critical data for the power plant control systems to increase efficiency, improve safety and reduce down times. For example, the combined heat-flux sensor/control systems can contribute to the optimization of burner and boiler operations under flexible loads, the optimization of heat-exchange conditions and overall reduction of heat rate and emissions, the prediction of imminent overheating conditions, and the optimization of the soot-blowing protocols.

20 FOSSIL-FUELED POWER PLANTS↗

Comparing four heat-inducible promoters in stably transformed sugarcane regarding spatial and temporal control of transgene expression reveals candidates to drive stem-preferred transgene expression

Small heat shock protein (sHSP) promoters contain cis-regulatory elements that facilitate transcription in response to heat stress, making them valuable tools for functional studies through controlled gene expression and the precise regulation of gene-editing tools or morphogenic regulators. To evaluate their utility, GUS reporter gene expression driven by four plant-sourced HSP promoters (p GmHSP17.5 , p HvHSP17 , p ZmHSP17.7 , and p ZmHSP26 ) was compared across various tissues of stably transformed sugarcane before and after heat treatment. At 22°C, all promoters showed minimal activity in leaves and roots, although p ZmHSP17.7 and p HvHSP17 displayed moderate expression in stems. Following heat treatment, all promoters exhibited their highest activity in stems, followed by leaves and roots. In stem tissues, p GmHSP17.5 displayed heat-induced uidA expression comparable to the constitutive p ZmUbi promoter. Notably, heat-induced reporter gene activity in stem middle sections of single-copy transgenic lines containing p ZmHSP17.7 , p HvHSP17 , or p ZmHSP26 exceeded p ZmUbi -derived uidA activity by 9.7-fold, 3.8-fold, and 3.0-fold, respectively, with 346- to 3,672-fold induction compared to control conditions. Most promoters showed peak expression in the middle sections of the stem, while p HvHSP17 was the most active in the stem apices. Histochemical analysis revealed that p ZmHSP17.7 and p HvHSP17 were active in both parenchyma cells and vascular bundles within sugarcane stems. Among leaf tissues, mature leaves exhibited greater expression than senescing or immature leaves, while root activity remained consistently minimal across all promoters. Temperature-course experiments identified distinct activation thresholds: 34°C–36°C for p ZmHSP17.7 , 36°C for p ZmHSP26 , 36°C–38°C for p HvHSP17 , and 40°C–42°C for p GmHSP17.5 . Drought stress also induced reporter gene transcription in stems under HSP promoters, although with lower fold induction than heat treatment. These findings provide valuable tools for gene function studies and biotechnology applications, including heat stress tolerance research, controlled transgene expression in metabolic engineering, precision gene editing, and developmental biology studies.

60 APPLIED LIFE SCIENCES↗

Design of Phase-Change Thermal Storage Device in a Heat Pump for Building Electric Peak Load Shaving: Preprint

Replacing carbon-intensive fossil fuel heating systems with electric heat pumps powered by renewables is a promising approach to decarbonize the building sector. However, one of the technical barriers of this approach is the large scale of heat demand, which will put excessive stress on the electricity grid. Integrating thermal energy storage (TES) into the heating systems can help alleviate this problem, by shifting thermal load and thus shaving peaks in the building electric load. Therefore, it is critical to understand how to design a thermal storage device in a heat pump for peak load shaving. In this study, we developed a numerical model for a cascaded vapor compression heat pump system integrating a phase change thermal storage device. This novel system can control the net thermal charging and discharging rate of the TES independently from the building's thermal load, which allows for precise control of electric power use. In the current study, we controlled the system to shave building electric peak load during a cold winter morning, and to charge TES during the relatively warm afternoon while still providing space heating. We used the model to evaluate the system performance with and without peak shaving. We also investigated the effect of PCM transition Tt on the peak reduction and electric energy saving potentials. The results show that peak shaving scheme effectively reduces the peak electric power consumption during a predefined discharge time window. When comparing to the no shaving case, for PCM Tt = 10 degrees C, the peak electric load reduction is 23.5%. When comparing to an air-source heat pump with back up electric heater, as Tt increases from 0 to 20 degrees C, the peak reduction increases from 46.1% to 50.9%. Integrating PCM Tt = 10 degrees C with peak shaving leads to the 45.5% of electric energy saving, which is the highest among the three transition temperatures.

electric energy saving↗

Numerical Investigation of High Delta T Sensible Storage Integrated CO2 Heat Pump: Preprint

To assist building heating electrification, this paper numerically investigates a load flexible heat pump system for commercial buildings. The system consists of a CO2 vapor compression cycle, a sensible thermal storage tank, and an air handling unit. The thermal storage medium is inexpensive, non-toxic and stable anti-freeze solution (30% potassium acetate). The air handing unit has an indoor coil and a ventilation coil. The system can be used to manage building electric load. During peak hours, the heat pump is off and the hot solution water is discharged from the tank to heat up the indoor air and ventilation air. During the hour of charge, the heat pump delivers hot solution water to the tank and to the air. The tank can also stand by while the heat pump provides space heating directly. We selected a medium sized office building located in Minnesota as the representative building and used EnergyPlus to obtain its 24 hour load data. We designed three storage tank volumes assuming 50 degrees C, 65 degrees C and 80 degrees C tank temperatures to independently provide the building load for 4 hours in the morning. The higher the tank temperature, the smaller the required volume, and thus higher energy density. The effective energy density is 78 with an 80 degrees C tank, and 40 kWhth/m3 with 50 degrees C. We simulated the tank integrated heat pump performance subjected to the 24-hour building load profile and ambient data. The baseline is the same system without storage tank. There was a trade-off between the storage energy density and the charging COP. The charge hour COP was 2.77 to charge the tank to 80 degrees C, and 3.01 to 50 degrees C. The proposed system could shift building load from the peak hours (8:00 - 12:00) to off-business hour (23:00 - 7:00+1). It eliminated 100% compressor electricity use during the peak hours, and avoided a peak electric power of 34 kW. The 65 degrees C tank saved 9.5 kWhe (4%) considering all day operation, which was the best balance between energy density and the system operation efficiency among the three options.

CO2 heat pump↗

Heating effects on jack pine pyrogenic organic matter properties from a pyrocosm study in 2022

This dataset contains data associated with the preprint “Fire removes preexisting pyrogenic organic matter from the ecosystem through the mechanisms of both direct combustion and increasing mineralizability” (Luo et al., 2025b), which is the complementary study to the published paper “Reburning pyrogenic organic matter: a laboratory method for dosing dynamic heat fluxes from above” (Luo et al., 2025a). We designed a full-factorial experiment with different burial depths of jack pine (Pinus banksiana Lamb) pyrogenic organic matter (PyOM) (Surface, 1 cm, and 5 cm) and different heat-flux profiles (High, Low, and Control) to examine how subsequent fires affect the properties of preexisting PyOM. We measured total carbon (C), pH, dissolved organic carbon (DOC), dissolved inorganic carbon (DIC), and mineralized C (as CO₂-C, from a 12-week incubation).We found that high heat flux and/or surface placement resulted in substantial direct C losses through combustion. Intermediate heat exposure produced both combustion losses and increases in DOC and mineralizability, which may have complex long-term implications: an increased dissolved fraction of PyOM may promote downward transport into mineral soils and potentially contribute to deeper, longer-term C storage, but it may also make PyOM more susceptible to microbial decomposition. Under the lowest heat flux and deepest burial, most PyOM was retained, and changes in DOC and C mineralization were minimal. Finally, PyOM pH, an important chemical property, decreased under low-temperature heating but increased under higher temperatures.We uploaded pH data for all samples (“pH_of_all_samples.csv”); pH and temperature-related data (peak temperature and degree hours) for samples in High and Low heat-flux treatments (“pH_vs_peakT_and_degree_hours_only_for_heated_samples.csv”); total C data (“CN_pct_C_stock_C_loss_in_samples.csv”); DOC and DIC data (“doc_dic.csv”); and mineralized C (CO₂-C) data (“CO2-C_all_original.csv”). Additional details can be found in the Methods & Sampling section.All datasets uploaded to ESS-DIVE are clearly labeled, cleaned, and include both raw and derived data, ready for reuse in other analyses. All analysis code and raw datasets are also available on GitHub: https://github.com/MengmengLuo/Fire-removes-preexisting-pyrogenic-organic-matter-from-the-ecosystem.

54 ENVIRONMENTAL SCIENCES↗

The Effect of Slag Heat Treatment on Ingot Quality during Cold Start Electroslag Remelting

Heat treatments of a commercial 40CaF 2 –30CaO–30Al 2 O 3 slag are explored to study their effect on quality of ingots produced using cold start electroslag remelting (ESR). Reducing heat treatment (or roasting) time or temperature for slag degassing can lead to energy savings as well as longer lifetimes for various components. Heat treatments in vacuum with pressure monitoring and air are investigated along with several temperatures, heating rates, and/or holding times. A research scale ESR furnace is used with steel electrodes. Significant differences in ingot quality, such as the occurrence of pores and a rough sidewall near the bottom of the ingot produced using slag heat treated in air at 580 °C, are observed. Further, adjustments in the heat treatment temperature for heat treatment in air eventually lead to an ingot quality comparable to that obtained using slag degassing in a controlled atmosphere at higher temperatures. Using differential thermal analysis, it is found that moisture is primarily removed from the slag at ≈460 °C and 700 °C. A safe slag roasting temperature is concluded to be 750 °C. Improper slag heat treatment leads to hydrogen concentrations from 8 to 15 ppm in about one‐quarter of the ingot volume.

36 MATERIALS SCIENCE↗

Active Divertor Heat Flux Control using Impurity Powder Dropper

Divertor plasma-facing components (PFCs) in a tokamak are typically designed to withstand average steady-state heat loads of about 5–10 MW/m 2 , a limit that applies to both solid and liquid lithium (LL) PFCs. Exceeding these design values can result in surface damage to tungsten PFCs or excessive lithium (Li) evaporation in liquid lithium divertor (LLD) PFCs. Since exceeding the divertor heat load limits has serious consequences, it is therefore prudent to develop a tool to reduce the divertor heat load and bring the heat load to within the design limit without affecting the plasma performance. Active low Z impurity injection such as Li has been suggested as a potential solution to mitigate excess heat flux as suggested previously, given that non-coronal radiation can be quite large ~ 20–30 MJ per mole of injected Li. Li is considered desirable for reducing the edge neutral recycling helping to improve plasma energy confinement. In this paper, we model the Impurity Power Dropper (IPD) to investigate its potential of divertor heat flux control. The IPD is typically located at the top of the tokamak device and uses a vertical drift tube of a few meters. In the 2 m drift tube case, the IPD powder is accelerated to ~ 6 m/sec before reaching the plasma with the upper divertor configuration, matching the condition for the in-board side pellet injection case. By modeling the IPD geometry we determined the IPD powder deposition profile, and thus the non-coronal radiation and ionization profiles in time as well. From the enhanced radiation power loss, it is therefore possible to reduce the divertor heat load using the divertor simulation code. In conclusion, the IPD divertor heat flux control can be tested in the facilities with IPD including ST-40, DIII-D, EAST, WEST and NSTX-U.

Active lithium injection↗

Technoeconomic optimization of superalloy supercritical CO 2 microtube shell-and-tube-heat exchangers

Heat exchangers are critical components of aerospace systems and improve efficiency of operation by providing waste heat recovery. Supercritical CO has emerged as a promising candidate working fluid due to its potential in increasing the power density of heat exchangers. In this paper, a generalized costing model is developed to estimate the capital costs incurred to manufacture microtube shell-and-tube heat exchangers. This model is utilized in conjunction with an accurate and efficient 2D numerical shell-and-tube heat exchanger performance prediction model to conduct optimization studies with two key objectives – minimization of cost and maximization of heat exchanger power density – on supercritical CO microtube heat exchangers utilizing superalloy Haynes 282 as the solid material. A methodology is then demonstrated to optimize these heat exchangers using Particle Swarm Optimization for aerospace applications and highly compact and cost-effective optimal designs with power density around 20 and cost per conductance less than 5 $ ∙ K/W are obtained.

42 ENGINEERING↗

Thermal hydraulic characteristics of liquid sodium and supercritical CO 2 in a diffusion bonded heat exchanger

Compact heat exchangers are being examined for use in advanced, high-temperature nuclear or concentrated solar power plants to improve thermal efficiencies and reduce capital costs. At the time of this work, no experimental investigation of liquid sodium in the popular diffusion-bonded heat exchanger design had been performed. Herein this paper describes the thermal-hydraulic and overall operations performance observed during testing of a 316/L stainless steel “zig-zag” microchannel heat exchanger operating between a (0.5 kg/s) sodium loop and a (0.3 kg/s) recuperated supercritical CO 2 system with outlet conditions up to 550 °C at 16 MPa. Observed pressure drop is presented in the form of an effective friction factor correlation which aligns closely with recent experimental data. The thermal performance of this heat exchanger, measured by its overall conductance under various flow conditions, was decomposed to examine the independent convection behaviors of both the CO 2 and sodium streams. These results are presented in the form of Nusselt correlations fitted to the experimental data. The sodium heat transfer correlation aligns closely with the widely accepted Lubarsky-Kaufman correlation, while the CO 2 heat transfer correlation lies between previous correlations developed by computational and experimental methods. This was the first long-duration use of sodium in a microchannels (~1.8 mm hydraulic diameter) heat exchanger. The use of an oxide control system enabled trouble-free operation of the heat exchanger for over 150 h of high-temperature operation. This first-of-a-kind demonstration using liquid sodium and supercritical carbon dioxide shows great promise for further application of this technology in advanced nuclear or solar power plants.

42 ENGINEERING↗

Refractory-based thermal energy storage for industrial process heat: one-dimensional modeling, control, and optimization

The variable and weather-dependent output of wind and solar power plants present a substantial challenge for planning and operating electricity-systems, particularly in the absence of cost-effective and dispatchable energy storage technologies. This study investigates a high-temperature, electrically heated, refractory-based thermal energy storage (RTES) system that stores electrical energy as sensible heat in dense ceramic bricks over the 950–1800 °C range. The stored heat can be discharged as a controlled hot-gas stream for industrial heating, fuel substitution in high-temperature processes, or electricity generation. The main novelty is a comprehensive modelling, control, mapping, and optimization framework that integrates one-dimensional transient gas–solid heat transfer, fan-assisted discharge, bypass-flow regulation, reheating logic, fan-power evaluation, insulation-loss assessment, and genetic-algorithm-based design optimization. The model uses feedback from outlet temperature and delivered power to regulate discharge, while a two-stage genetic algorithm optimizes brick-channel geometry, gas-flow operation, and multilayer insulation thicknesses. Storage capacities below 50 MWh and discharge powers of 5–30 MW are analyzed to evaluate hold time, thermal delivery, fan-power penalty, heat loss, state-of-charge evolution, and indicative capital cost. Results demonstrate that optimized and well-insulated refractory-based thermal energy storage units can provide stable, efficient, and repeatable heat delivery over multiple discharge cycles. The generated performance and cost maps support modular refractory thermal energy storage as a practical option for large-scale integration of wind and solar generation and for high-temperature industrial process heat.

25 ENERGY STORAGE↗

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↗

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The reservoir network: A new network topology for district heating and cooling

Thermal district networks are effective solutions to substitute fossil fuels with renewable energy sources for heating and cooling. Moreover, thermal networking of buildings allows energy efficiency to be further increased. The waste heat from cooling can be reused for heating in thermal district systems. Because of bidirectional energy flows between prosumers, thermal networks require new hydraulic concepts. In this work, we present a novel network topology for simultaneous heating and cooling: the reservoir network. The reservoir network is robust in operation due to hydraulic decoupling of transfer stations, integrates heat sources and heat sinks at various temperature levels and is flexible in terms of network expansion. Furthermore, we used Modelica simulations to compare the new single-pipe reservoir network to a basecase double-pipe network, taking yearly demand profiles of different building types for heating and cooling. The electric energy consumed by the heat pumps and circulations pumps differs between the reservoir and base case networks by less than 1%. However, if the reservoir network is operated with constant instead of variable mass flow rate, the total electrical consumption can increase by 48% compared to the base case. As with any other network topology, the design and control of such networks is crucial to achieving energy efficient operation. Investment costs for piping and trenching depend on the district layout and dimensioning of the network. If a ring layout is applied in a district, the reservoir network with its single-pipe configuration is more economical than other topologies. For a linear layout, the piping costs are slightly higher for the reservoir network than for the base case because of larger pipe diameters.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A first principles framework to predict the transient performance of latent heat thermal energy storage

Thermal energy storage (TES) is increasingly recognized as an essential component of efficient Combined Heat and Power (CHP), Concentrated Solar Power (CSP), Heating Ventilation and Air Conditioning (HVAC), and refrigeration as it reduces peak demand while helping to manage intermittent availability of energy (e.g., from solar or wind). Latent Heat Thermal Energy Storage (LHTES) is a viable option because of its high energy storage density. Parametric analysis of LHTES in terms of dimensionless numbers is highly desired as a tool to model LHTES systems. One approach is to develop a model equation so as to minimize the error between the model and data obtained from experiments or simulations. While this approach can produce an accurate correlation applicable within the range of data used for its creation, it does not provide physical understanding of the rate-limiting process controlling the transient behavior of the device. In this paper we present an alternative approach whereby the potential rate-limiting processes are identified from first principles and then the key process is determined as a function of time as a LHTES device is charged. For example, in a simple geometry, the melt-fraction can be expected to vary linearly in time if the heat transfer rate is limited by natural convection of the phase changing material and we show it scales with the PCM Grashof number as $Gr^1_p$ and PCM Prandtl number as $Pr_p^{(1/3)}$. On the other hand, if surface area of solid PCM limits the heat transfer rate, the melt fraction increases asymptotically to reach full melting. The existence of these linear and asymptotic regions and the $Gr^1_pP r^{1/3}_p$ shape of the melt fraction curve is verified using our database of 64 simulations. Of practical importance in designing LHTES devices is the melt fraction at which the heat transfer rate ceases to be limited by convection, after which the heat storage rate deteriorates. For our geometry, this is found to be about 90%. This test case of our methodology shows the value of our approach, that predicting heat storage rate based on the rate-limiting physical phenomenon as a function of time is an effective approach to modeling LHTES devices.

25 ENERGY STORAGE↗

Development of a high-temperature Inconel 625 heat exchanger by model design and binder jetting additive manufacturing

A nickel-based Inconel 625 superalloy heat exchanger for high-temperature applications was developed via binder jetting additive manufacturing. The material properties were characterized first on printed and sintered parts. Two sintering temperatures were used to investigate the effects of the temperature on densification and microstructure. For channel geometry design, the heat transfer capabilities of the heat exchanger were optimized for the cross-section geometry of fluid flow channels in a counterflow configuration, and a stress analysis was conducted to investigate the effects of channel geometry. After headers were incorporated into the heat exchanger as a one-piece component, a prototype was printed. A complete depowdering was achieved via compressed air blowing, and sintering was performed with the developed profile. Ultimately, a one-piece heat exchanger with a nearly full density was obtained. Heat transfer tests were performed on this unit, and the results were compared with those from the simulations. In this study, systematic processes were developed for an additively manufactured Inconel-based heat exchanger for high-pressure, high-temperature heat transfer applications.

36 MATERIALS SCIENCE↗

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↗