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 55 records · Page 3

Modeling a Sodium Heat Pipe Experiment at SPHERE Using Sockeye

The Single Primary Heat Extraction and Rejection Emulator (SPHERE) facility at Idaho Na- tional Laboratory was recently utilized to generate data for the startup and steady operation of a high-performance, sodium heat pipe over the course of 1000 hours, as a test of detrimental, long-term effects of heat pipe operation. The setup consists of a single, sodium heat pipe enclosed in a stainless-steel vacuum chamber, heated radiatively via a cylindrical ceramic fiber heater configuration and cooled via a water-cooled calorimeter. Measurements include temperatures at several axial locations along the outer surface of the heat pipe, the power provided to the heaters, and the heat removal rate of the calorimeter. In this work, this data is utilized to validate heat pipe models in the heat pipe application Sockeye, which is based upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Sockeye provides various heat pipe models at an engineering scale appropriate for the multiphysics simulation of microreactors, which may feature several hundred heat pipes. This work details models of this experiment in SPHERE using various heat pipe models with Sockeye, including heat conduction-based models and compressible flow models of the heat pipe interior. These models are compared to the experimental data to assess the accuracy of several aspects of heat pipe modeling, including frozen startup, the effect of non-condensable gases, and the coupling of the heat pipe to its environment.

97 - MATHEMATICS AND COMPUTING↗

Thermal modeling of phase change solidification in thermal control devices including natural convection effects

Natural convection effects in phase change thermal control devices were studied. A mathematical model was developed to evaluate natural convection effects in a phase change test cell undergoing solidification. Although natural convection effects are minimized in flight spacecraft, all phase change devices are ground tested. The mathematical approach to the problem was to first develop a transient two-dimensional conduction heat transfer model for the solidification of a normal paraffin of finite geometry. Next, a transient two-dimensional model was developed for the solidification of the same paraffin by a combined conduction-natural-convection heat transfer model. Throughout the study, n-hexadecane (n-C16H34) was used as the phase-change material in both the theoretical and the experimental work. The models were based on the transient two-dimensional finite difference solutions of the energy, continuity, and momentum equations.

Ukanwa, A. O.↗

Thermal Hydraulics Analysis of a Divertor Monoblock Using SALAMANDER

Divertors are critical components in magnetic confinement fusion devices. One of the Divertor's crictial roles is absorbing the highest heat flux from the plasma. The Divertor accomplishes this through the use of divertor monoblocks and cooling channels. This work presents a thermal hydraulic analysis of a divertor monoblock using the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering and Research (SALAMANDER). SALAMANDER is an open source tool developed at Idaho National Laboratory for conducting multiphysics and multscale analysis of Plasma Facing Components (PFCs). For the monoblock, solid heat conduction is modeled in the block with a constant heat flux of 1E7 W/m^2 applied to the top of the block to represent peak heat flux from the plasma. The cooling channel (which is pressurized water at a Reynolds number of ~1,000,000) is modeled using a k-epsilon turbulence model with standard wall functions. The heat transfer between the solid and fluid domain is modeled using the Dittus-Boelter correlation for the convective heat transfer coefficient. This work aims to show the importance of a multiphysics modeling approach when performing simulations on PFCs.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling Heat Pipes with Non-Condensable Gases

The report showcases advancements in modeling heat pipes, considering both scenarios with and without non-condensable gases. To achieve this, two distinct modeling approaches are juxtaposed. Firstly, the effective conduction model, implemented in the MOOSE-based code Sockeye, is employed. Secondly, a first-of-a-kind two-phase Euler-Euler Computational Fluid Dynamics (CFD) model is developed using the STAR-CCM+ code. Both models undergo validation against experimental data, acknowledging the inherent uncertainties associated with each modeling assumption. Interestingly, the non-tuned CFD model surpasses the performance of the calibrated conduction model for heat pipes operating under both conditions: with and without non-condensable gases. It's worth noting, however, that the CFD models entail significantly longer runtimes compared to the conduction models. Nevertheless, the insights garnered from the CFD model shed invaluable light on the intricate operational dynamics of heat pipes. Future work involves broadening the validation scope of these models and continuing their development to enhance their utility as robust tools for heat pipe design and operational support.

42 ENGINEERING↗

Response Of Ceramic Insulation To Aerothermodynamic Heating

Thermal conductivities determined from engineering model used to predict accurately thermal responses of multicomponent insulation on spacecraft. Types of insulation studied include such composite parts as alumina-enhanced thermal barrier backed by low-density fibrous refractory material. Internal responses of various types of multicomponent insulation to aerothermodynamic heating calculated with help of computer program that implements engineering model for thermal conductivity. Transfer of heat by internal radiation plays major role in responses of fibrous materials at high temperatures and incorporated into model and program. Used in future spacecraft and exposed during reentry and aerobraking maneuvers in outer atmosphere to environments even harsher than that of Space Shuttle.

Stewart, David A.↗

A Model of Thermal Conductivity for Planetary Soils: Theory for Cemented Soils - 2

A numerical model of heat conduction through particulate media made of spherical grains cemented by various bonding agents is presented. The pore-filling gas conductivity, volume fraction, and thermal conductivity of the cementing phase are tunable parameters. Cement fractions <0.001-0.01% in volume have small effects on the soil bulk thermal conductivity. A significant conductivity increase (factor 3-8) is observed for bond fractions of 0.01 to 1% in volume. In the 1 to 15% bond fraction domain, the conductivity increases continuously but less intensely (25-100% conductivity increase compared to a 1% bond system). Beyond 15% of cements, the conductivity increases vigorously and the bulk conductivity rapidly approaches that of bedrock. The composition of the cements (i.e. conductivity) has little influence on the bulk thermal inertia of the soil, especially if the volume of bond <10%. These results indicate that temperature measurements are sufficient to detect cemented soils and quantify the amount of cementing phase, but the mineralogical nature of the bonds and the typical grain size are unlikely to be determined from orbit. On Mars, a widespread surface unit characterized by a medium albedo (0.19-0.26) and medium/high thermal inertia (200-600 J s(0.5)/sq m/K) has long been hypothesized to be associated with a duricrust. The fraction of cement required to fit the thermal data is less than approx.1-5% by volume. This small amount of material is consistent with orbital observations, confirming that soil cementation is an important factor controlling the thermal inertia of the Martian surface

Piqueux, S.↗

Stanford Thermal Earth Model for the Conterminous United States

Provided here are various forms of the Stanford Thermal Earth Model, as well as the data and methods used for its creation. The predictions produced by this model were visualized in two-dimensional spatial maps across the modeled depths (0-7 km) for the conterminous United States. The thermal earth model is made available as an application programming interface (API) and as feature layers on ArcGIS, which are both provided via links below. A data-driven spatial interpolation algorithm based on physics-informed graph neural networks was used to develop these national temperature-at-depth maps. The model satisfied the three-dimensional heat conduction law by predicting subsurface temperature, surface heat flow, and rock thermal conductivity. Many physical quantities, including bottomhole temperature, depth, geographic coordinates, elevation, sediment thickness, magnetic anomaly, gravity anomaly, gamma-ray flux of radioactive elements, seismicity, and electric conductivity were used as model inputs. Surface heat flow, temperature, and thermal conductivity predictions were constructed for depths of 0-7 km at an interval of 1 km with spatial resolution of 18 km2 per grid cell. The model showed superior temperature, surface heat flow and thermal conductivity mean absolute errors of 4.8C, 8.1 mW/m2 and 0.07 W/(C-m), respectively..

15 GEOTHERMAL ENERGY↗

Extinguishment of a Diffusion Flame Over a PMMA Cylinder by Depressurization in Reduced-Gravity

Extinction of a diffusion flame burning over horizontal PMMA (Polymethyl methacrylate) cylinders in low-gravity was examined experimentally and via numerical simulations. Low-gravity conditions were obtained using the NASA Lewis Research Center's reduced-gravity aircraft. The effects of velocity and pressure on the visible flame were examined. The flammability of the burning solid was examined as a function of pressure and the solid-phase centerline temperature. As the solid temperature increased, the extinction pressure decreased, and with a centerline temperature of 525 K, the flame was sustained to 0.1 atmospheres before extinguishing. The numerical simulation iteratively coupled a two-dimensional quasi-steady, gas-phase model with a transient solid-phase model which included conductive heat transfer and surface regression. This model employed an energy balance at the gas/solid interface that included the energy conducted by the gas-phase to the gas/solid interface, Arrhenius pyrolysis kinetics, surface radiation, and the energy conducted into the solid. The ratio of the solid and gas-phase conductive fluxes Phi was a boundary condition for the gas-phase model at the solid-surface. Initial simulations modeled conditions similar to the low-gravity experiments and predicted low-pressure extinction limits consistent with the experimental limits. Other simulations examined the effects of velocity, depressurization rate and Phi on extinction.

Goldmeer, Jeffrey Scott↗

Heat Transfer in a Superelliptic Transition Duct

Local heat transfer measurements were experimentally mapped using a transient liquid-crystal heat transfer technique on the surface of a circular-to-rectangular transition duct. The transition duct had a length-to-diameter ratio of 1.5 and an exit-plane aspect ratio of 3. The crosssectional geometry was defined by the equation of a superellipse. The cross-sectional area was the same at the inlet and exit but varied up to 15 percent higher through the transition. The duct was preheated to a uniform temperature (nominally 64 C) before allowing room temperature air to be suddenly drawn through it. As the surface cooled, the resulting isothermal contours on the duct surface were revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record was made of the surface temperature and time data for all points on the duct surfaces during each test. Using this surface temperature-time data together with the temperature of the air flowing through the model and the initial temperature of the model wall, the heat transfer coefficient was calculated by employing the classic one-dimensional, semi-infinite wall heat transfer conduction model. Test results are reported for inlet diameter-based Reynolds numbers ranging from 0.4x106 to 2.4x106 and two grid-generated freestream turbulence intensities of about 1 percent, which is typical of wind tunnels, and up to 16 percent, which may be more typical of real engine conditions.

Poinsatte, Philip↗

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↗

Heat pipe life and processing study

The merit of adding water to the reflux charge in chemically and solvent cleaned aluminum/slab wick/ammonia heat pipes was evaluated. The effect of gas in the performance of three heat pipe thermal control systems was found significant in simple heat pipes, less significant in a modified simple heat pipe model with a short wickless pipe section. Use of gas data for the worst and best heat pipes of the matrix in a variable conductance heat pipe model showed a 3 C increase in the source temperature at full on condition after 20 and 246 years, respectively.

Antoniuk, D.↗

Development of a conduction-based model for analyzing frozen startup of alkali-metal heat pipes

One key area of interest in heat pipe modeling/simulation is to analyze the startup behavior of the liquid-metal heat pipes (LMHPs) from a frozen state. This so-called ‘frozen startup’ process involves a complex set of nonlinear mass and heat transport phenomena, including phase transitions from solid to liquid and vapor, multiphase interactions, microporous wick flow, and compressible vapor dynamics. The complexity of these processes makes it challenging to simulate LMHP’s frozen startup using conventional numerical methods or commercial computational fluid dynamics (CFD) software. This paper presents a simplified conduction-based modeling approach that can provide practical insights into the entire LMHP frozen startup process, while alleviating the challenges of modeling its complex physics. The theoretical foundation and physical assumptions of the proposed model are based solely on heat-conduction equation, allowing for a more tractable simulation without sacrificing essential physical accuracy. The proposed model was implemented in a commercial CFD software, and its prediction was compared with the experimental data obtained from sodium heat-pipe startup experiments. The comparison highlights the proposed model's ability to capture the transient thermal behavior of LMHP during frozen startup. This study not only validates the conduction-based frozen startup modeling method but also shows its potential as a practical and efficient tool for understanding the startup performance of the LMHP systems.

Microreactor↗

Modeling a Sodium Heat Pipe Experiment at SPHERE Using Sockeye

The Single Primary Heat Extraction and Rejection Emulator (SPHERE) facility at Idaho National Laboratory was recently utilized to generate data for the startup and steady operation of a high-performance, sodium heat pipe over the course of 1,000 hours to test the detrimental, long-term effects of heat pipe operation. The setup consisted of a single, sodium heat pipe enclosed in a stainless-steel vacuum chamber, heated radiatively via a cylindrical ceramic-fiber heater configuration and cooled via a water-cooled calorimeter. Measurements included temperatures at several axial locations along the outer surface of the heat pipe, the power provided to the heaters, and the heat removal rate of the calorimeter. In this work, we use this data to validate heat pipe models in Sockeye, a heat pipe application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Sockeye provides various heat pipe models at an engineering scale appropriate for the multiphysics simulation of microreactors, which may feature several hundred heat pipes. This work details models of this experiment at SPHERE using various heat pipe models with Sockeye, including heat-conduction-based and compressible flow models of the heat pipe interior.

97 - MATHEMATICS AND COMPUTING↗

Reinterpretation of Mariner 9 IRIS data on the basis of a simulation of radiative-conductive convective transfer in the dust laden Martian atmosphere

Time dependent vertical distributions of atmospheric temperature and static stability were determined by a radiative-convective-conductive heat transfer model attuned to Mariner 9 IRIS radiance data. Of particular interest were conditions of both the dust-laden and dust-free atmosphere in the middle latitudes on Mars during the late S.H. summer season. The numerical model simulates at high spatial and temporal resolution (52 atmospheric and 30 subsurface levels; with a time-step of 7.5 min.) the heat transports in the ground-atmosphere system. The algorithm is based on the solution of the appropriate heating rate equation which includes radiative, molecular-conductive and convective heat transfer terms. Ground and atmosphere are coupled by an internal thermal boundary condition.

Pallman, A. J.↗