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 145 records · Page 8

Estimation of soil moisture from diurnal surface temperature observations

A coupled heat and moisture balance model was used to determine the thermal inertia of a grass covered top soil under different meteorological conditions. Relations between thermal inertia and soil moisture were established using the De Vries models for thermal conductivity and heat capacity to relate soil moisture and thermal inertia as a function of soil type. A sensitivity study of the surface roughness length and thermal inertia on diurnal surface temperature shows the necessity of focusing on the night time surface temperature rather than on the day time surface temperature, in order to estimate the soil moisture content of the top soil.

Vandegriend, A. A.↗

Snow thermal conductivity and conductive flux in the Central Arctic: Estimates from observations and implications for models

During the Arctic winter, the conductive heat flux through the sea ice and snow balances the radiative and turbulent heat fluxes at the surface. Snow on sea ice is a thermal insulator that reduces the magnitude of the conductive flux. The thermal conductivity of snow, that is, how readily energy is conducted, is known to vary significantly in time and space from observations, but most forecast and climate models use a constant value. This work begins with a demonstration of the importance of snow thermal conductivity in a regional coupled forecast model. Varying snow thermal conductivity impacts the magnitudes of all surface fluxes, not just conduction, and their responses to atmospheric forcing. Given the importance of snow thermal conductivity in models, we use observations from sea ice mass balance buoys installed during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition to derive the profiles of thermal conductivity, density, and conductive flux. From 13 sites, median snow thermal conductivity ranges from 0.33 W m −1 K −1 to 0.47 W m −1 K −1 with a median from all data of 0.39 W m −1 K −1 from October to February. In terms of surface energy budget closure, estimated conductive fluxes are generally smaller than the net atmospheric flux by as much as 20 W m −2 , but the average residual during winter is −6 W m −2 , which is within the uncertainties. The spatial variability of conductive heat flux is highest during clear and cold time periods. Higher surface temperature, which often occurs during cloudy conditions, and thicker snowpacks reduce temporal and spatial variability. These relationships are compared between observations and the coupled forecast model, emphasizing both the importance and challenge of describing thermodynamic parameters of snow cover for modeling the Arctic as a coupled system.

54 ENVIRONMENTAL SCIENCES↗

The 0.040-scale space shuttle orbiter base heating model tests in the Lewis Research Center space power facility

Space shuttle base heating tests were conducted using a 0.040-scale model in the Plum Brook Space Power Facility of The NASA Lewis Research Center. The tests measured heat transfer rates, pressure distributions, and gas recovery temperatures on the orbiter vehicle 2A base configuration resulting from engine plume impingement. One hundred and sixty-eight hydrogen-oxygen engine firings were made at simulated flight altitudes ranging from 120,000 to 360,000 feet.

Dezelick, R. A.↗

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The effects of inlet turbulence and rotor/stator interactions on the aerodynamics and heat transfer of a large-scale rotating turbine model. Part 4: Aerodynamic data tabulation

A combined experimental and analytical program was conducted to examine the effects of inlet turbulence and airfoil heat transfer. The experimental portion of the study was conducted in a large-scale (approx. 5X engine), ambient temperature, rotating turbine model configured in both single-stage and stage-and-a-half arrangements. Heat transfer measurements were obtained using low-conductivity airfoils with miniature thermocouples welded to a thin, electrically heated surface skin. Heat transfer data were acquired for various combinations of low or high inlet turbulence intensity, flow coefficient, first stator-rotor axial spacing, Reynolds number and relative circumferential position of the first and second stators. Aerodynamic measurements obtained include distributions of the mean and fluctuating velocities at the turbine inlet and, for each airfoil row, midspan airfoil surface pressures and circumferential distributions of the downstream steady state pressures and fluctuating velocities. Results include airfoil heat transfer predictions produced using existing 2-D boundary layer computation schemes and an examination of solutions of the unsteady boundary layer equations.

Dring, R. P.↗

Experiments on thermoacoustic convection heat transfer in gravity and zero-gravity environments

The results of an experimental study of thermoacoustic convection (TAC) heat transfer in gravity and zero-gravity environments are presented. The experimental apparatus consisted of a cylinder containing air as the compressible fluid. The enclosed air was heated electrically at the top surface which consisted of a thin high-resistance steel foil connected to a power source. Thermocouples were used to measure the transient temperature of the air on the axis of the cylinder and the heated surface in the both zero-gravity and gravity environments. The zero-gravity tests were performed in the Zero-Gravity Drop Tower Facility of NASA-Lewis Research Center. The experimental results were corrected for the error due to radiation absorption by the thermocouples. A conduction-only numerical heat transfer model was developed to compute the transient air temperature in the cylindrical geometry. The results were compared to the experimental data to determine the significance of the thermoacoustic convection heat transfer mechanism. It is observed that the rate of heat transfer to the air measured during the experiments is consistently higher than that obtained by the conduction-only solution indicating a significant presence of the TAC heat transfer. Further experiments are planned to measure directly (1) the radiative heat transfer contribution to the rise in the air temperature, and (2) the air pressure oscillations within the cylinder that are responsible for the convective heat transfer mode.

Parang, Masood↗

Development of Solar Wind Model Driven by Empirical Heat Flux and Pressure Terms

We are developing a time stationary self-consistent 2D MHD model of the solar corona and solar wind as suggested by Sittler et al. (2003). Sittler & Guhathakurta (1999) developed a semiempirical steady state model (SG model) of the solar wind in a multipole 3-streamer structure, with the model constrained by Skylab observations. Guhathakurta et al. (2006) presented a more recent version of their initial work. Sittler et al. (2003) modified the SG model by investigating time dependent MHD, ad hoc heating term with heat conduction and empirical heating solutions. Next step of development of 2D MHD models was performed by Sittler & Ofman (2006). They derived effective temperature and effective heat flux from the data-driven SG model and fit smooth analytical functions to be used in MHD calculations. Improvements of the Sittler & Ofman (2006) results now show a convergence of the 3-streamer topology into a single equatorial streamer at altitudes > 2 R(sub S). This is a new result and shows we are now able to reproduce observations of an equatorially confined streamer belt. In order to allow our solutions to be applied to more general applications, we extend that model by using magnetogram data and PFSS model as a boundary condition. Initial results were presented by Selwa et al. (2008). We choose solar minimum magnetogram data since during solar maximum the boundary conditions are more complex and the coronal magnetic field may not be described correctly by PFSS model. As the first step we studied the simplest 2D MHD case with variable heat conduction, and with empirical heat input combined with empirical momentum addition for the fast solar wind. We use realistic magnetic field data based on NSO/GONG data, and plan to extend the study to 3D. This study represents the first attempt of fully self-consistent realistic model based on real data and including semi-empirical heat flux and semi-empirical effective pressure terms.

Sittler, Edward C., Jr.↗

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗

The role of thermal contact resistance in pyrotechnic ignition

This paper describes a pyrotechnic ignition model based on transient heat conduction from a heated bridgewire to a pyrotechnic that is placed in contact with it. The boundary condition used at the interface was a thermal contact conductance estimated at 31,200 W per sq m-K between the wire and the pyrotechnic. Ignition was assumed to occur when a 2.5-micron layer of pyrotechnic next to the bridgewire reached a critical ignition temperature. The times to ignition predicted by this model for constant current firings were in good agreement with experimentally observed times to fire at 3.5- and 5-ampere current levels and ambient temperatures from 144 K to 366 K.

Sernas, V.↗

Uncertainty Analysis of Coaxial Thermocouple Calorimeters Used in Arc Jets

Recent introduction of Coaxial Thermocouple type calorimeters into the NASA Ames arc jet facilities has inspired an analysis of 2D conduction effects internal to this type of calorimeter. Lateral conduction effects violate the 1D finite slab inverse analysis which is typically used to deduce the heat transfer to such calorimeters. The spherical shaped nose associated with most calorimeters (rather than flat) leads to a bias error that over-estimates the stagnation heating. Non-uniform heating on the face of spherically shaped calorimeters leads to conduction losses to the colder rim of the calorimeter which causes an underestimate of the stagnation heating. These two effects come into play at different times of the calorimeter's exposure to the arc jet, so they do not cancel. The spherical body effects come into play in the early stages of exposure, while the non-uniform heating effect becomes most severe at the later stages of exposure. The bias associated with spherical effects can be avoided by rewriting the 1D finite slab inverse analysis code to solve for 1D conduction in spherical coordinates. However, reducing the bias error associated non-uniform heating requires a somewhat ad hoc modification to the 1D finite element inverse analysis.

Isotherms↗

Oceanic lithosphere and asthenosphere: The thermal and mechanical structure

A coupled thermal and mechanical solid state model of the oceanic lithosphere and asthenosphere is presented. The model includes vertical conduction of heat with a temperature dependent thermal conductivity, horizontal and vertical advection of heat, viscous dissipation or shear heating, and linear or nonlinear deformation mechanisms with temperature and pressure dependent constitutive relations between shear stress and strain rate. A constant horizontal velocity u sub 0 and temperature t sub 0 at the surface and zero horizontal velocity and constant temperature t sub infinity at great depth are required. In addition to numerical values of the thermal and mechanical properties of the medium, only the values of u sub 0, t sub 0 and t sub infinity are specified. The model determines the depth and age dependent temperature horizontal and vertical velocity, and viscosity structures of the lithosphere and asthenosphere. In particular, ocean floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of the age of the ocean floor.

Schubert, G.↗

Improved Ionospheric Electrodynamic Models and Application to Calculating Joule Heating Rates

Improved techniques have been developed for empirical modeling of the high-latitude electric potentials and magnetic field aligned currents (FAC) as a function of the solar wind parameters. The FAC model is constructed using scalar magnetic Euler potentials, and functions as a twin to the electric potential model. The improved models have more accurate field values as well as more accurate boundary locations. Non-linear saturation effects in the solar wind-magnetosphere coupling are also better reproduced. The models are constructed using a hybrid technique, which has spherical harmonic functions only within a small area at the pole. At lower latitudes the potentials are constructed from multiple Fourier series functions of longitude, at discrete latitudinal steps. It is shown that the two models can be used together in order to calculate the total Poynting flux and Joule heating in the ionosphere. An additional model of the ionospheric conductivity is not required in order to obtain the ionospheric currents and Joule heating, as the conductivity variations as a function of the solar inclination are implicitly contained within the FAC model's data. The models outputs are shown for various input conditions, as well as compared with satellite measurements. The calculations of the total Joule heating are compared with results obtained by the inversion of ground-based magnetometer measurements. Like their predecessors, these empirical models should continue to be a useful research and forecast tools.

Weimer, D. R.↗

The lunar magma ocean - A transient lunar phenomenon

The time required for the solidification of a lunar magma ocean is considered. In the absence of a thermal boundary crust, heat loss is by radiation and the ocean will solidify in a few decades. However, only a few millimeters of crust would effectively limit radiation. The current investigation has the objective to present estimates of conduction boundary (crust) thickness based on several models of the mechanism of crustal formation. Three different models for estimating conduction boundary layer thickness are discussed. They are based on the formation of a chill margin, a plagioclase flotation layer, and a layer which is continually destroyed by meteorite impacts and regrows by conductive heat loss. Based on these models, the maximum solidification time for a 200 km thick ocean is about 100 million years. A more probable solidification time based on a plagioclase flotation crust model is about 60 million years.

Minear, J. W.↗

Measuring Titanium Heat Pipe Thermal Conductivity for CO2 Adsorption Systems

Current models of the Air-Cooled Temperature Swing Adsorption Compressor (AC-TSAC) being developed for long duration space exploration utilize cooling chambers and heat plates to adsorb and desorb contaminants -13 . To reduce resources required and lower the equivalent system mass (ESM), the study focuses on determining the applications of titanium heat pipe systems to operate both phases, allowing for future research into the development of a heat pipe chamber to hold the zeolite formulations.

multifunctional sorbent↗

The MOOSE Thermal Hydraulics Module

The Multiphysics Object-Oriented Simulation Environment (MOOSE) is an open-source object-oriented finite element framework written in C++ (Lindsay et al., 2022). The Thermal Hydraulics Module (THM) is an optional MOOSE physics module that provides capabilities for studying thermal hydraulic systems. Its core capability lies in assembling a network of coupled components, for instance, pipes, junctions, and valves. THM provides several new systems to MOOSE to enable and facilitate thermal hydraulic simulations, most notably the Components system, which provides a higher-level syntax to MOOSE’s lower-level objects. This system is extensible by the user, but the current library primarily includes components based on a one-dimensional, single-phase, variable-area, compressible flow model, as well as heat conduction.

42 ENGINEERING↗

Improved Particle Heat Transfer by way of Bimodal Particle Distributions for High Temperature Solar Thermal Energy

High temperature solar thermal facilities are looking to increase operating temperatures through novel heat transfer media, one such being solid particles. These particles operating at high temperatures will require transferring their thermal energy into another working fluid like supercritical carbon dioxide which can be used in advanced power cycles. Achieving high heat transfer between the particles and supercritical carbon dioxide is essential to high efficiency and low-cost operation. Therefore, optimizing the thermal conductivity of these particles is one potential way to ensure high performance. Traditionally, unimodal particle distributions have been employed in high temperature particle solar power plants. However, ambient temperature testing of bimodal particle distributions has revealed a superior thermal conductivity when compared to its unimodal counterpart at the same temperature. This data was obtained by certified, off-the-shelf instruments that can effectively simulate the conditions a particle would be exposed to in a high temperature solar thermal system. Data obtained in this way suggests that the increased thermal conductively imputed by a bimodal particle distribution is significant at working temperatures in solar facilities. Furthermore, the thermal conductivity of these bimodal particle distributions peaks when the best combination of large and small particles is applied. At high temperatures, binary particle distributions are compared to monodispersed distributions of larger particles where heat transfer is more prolific due to the increased surface radiation. Various thermal conductivity, porosity and heat exchanger models are explored in conjunction with data acquired up to 700 C.

Stout, Dallin (ORCID:0009000294586091)↗