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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Modeling of Melt Growth During Carbothermal Processing of Lunar Regolith

The carbothermal processing of lunar regolith has been proposed as a means to produce carbon monoxide and ultimately oxygen to support human exploration of the moon. In this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. Carbon monoxide is further processed in other reactors downstream to ultimately produce oxygen. The amount of oxygen produced crucially depends on the amount of regolith that is molten. In this paper we develop a model of the heat transfer in carbothermal processing. Regolith in a suitable container is heated by a heat flux at its surface such as by continuously shining a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and its temperature rises until it attains the melting point. The energy from the heat flux is then used for the latent heat necessary to change phase from solid to liquid, after which the temperature continues to rise. Thus a small melt pool appears under the heated zone shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat and supplies the energy required to melt more of the regolith, and the size of the molten zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the melt is balanced by radiative losses from the surface. In this paper, we model the melting and the growth of the melt zone with time in a bed of regolith when a portion of its surface is subjected to a constant heat flux. The heat flux is assumed to impinge on a circular area. Our model is based on an axisymmetric three-dimensional variation of the temperature field in the domain. Heat transfer occurs only by conduction, and effects of convective heat transport are assumed negligible. Radiative heat loss from the surface of the melt and the regolith to the surroundings is permitted. We perform numerical computations to determine the shape and the mass of the melt at steady state and its time evolution. We first neglect the volume change upon melting, and subsequently perform calculations including it. Predictions from our model are compared to test data to determine the effective thermal conductivities of the regolith and the melt that are compatible with the data

Balasubramaniam, R.↗
NASA NTRS1994

On the vertical thermal structure of Io's atmosphere

A radiative-thermal conduction model for the vertical thermal structure of Io's atmosphere is developed with solar heating by SO2 absorption in UV and near-IR bands and non-local thermodynamic equilibrium (LTE) cooling by SO2 nu(sub 1), nu(sub 2), nu(sub 3), vibrational bands and rotational lines. The model predicts the existence of a mesopause in Io's atmosphere when the surface pressure exceeds approximately 10 nbar. The radiative time constant for establishing a mesosphere/mesopause on Io is only approximately 20 min, whereas the thermospheric radiative time constant is about 1 hr. These time constants are significantly shorter than the diurnal time scale and competitive with dynamic time scales. In the thermosphere when solar UV heating dominates, the asymptotic thermospheric temperature is approximately 270 K, only 140 K greater than the surface temperature because at high altiudes non-LTE cooling by SO2 rotation lines exceeds cooling in the nu(sub 2) virbrational band. Solar-heating-only modles are incapable of generating warm enough atmospheres to satisfy the observational inferences from UV and especially millimeter-wave meausrements. Joule heating driven by the penetration of Jupiter's corotational electric field into Io's conducting ionosphere is demonstrated to be the dominant heating mechanism in the subnanobar regions of Io's atmosphere with tempertures ranging from 150 to 1000 K as a function of decreasing pressure from 1 to 0.1 nbar, The asymoptotic thermospheric temperature can attain a value as high as 1800 K.

Strobel, Darrell F.↗
NASA NTRS1993

Observational and modeling studies of heat, moisture, precipitation, and global-scale circulation patterns

The research sponsored by this grant is a continuation and an extension of the work conducted under a previous contract, 'South Pacific Convergence Zone and Global-Scale Circulations'. In the prior work, we conducted a detailed investigation of the South Pacific convergence zone (SPCZ), and documented many of its significant features and characteristics. We also conducted studies of its interaction with global-scale circulation features through the use of both observational and modeling studies. The latter was accomplished toward the end of the contract when Dr. James Hurrell, then a Ph.D. candidate, successfully ported the NASA GLA general circulation model (GCM) to Purdue University. In our present grant, we have expanded our previous research to include studies of other convectively-driven circulation systems in the tropics besides the SPCZ. Furthermore, we have continued to examine the relationship between these convective systems and global-scale circulation patterns. Our recent research efforts have focused on three objectives: (1) determining the periodicity of large-scale bands of organized convection in the tropics, primarily synoptic to intraseasonal time scales in the Southern Hemisphere; (2) examining the relative importance of tropical versus mid-latitude forcing for Southern Hemisphere summertime subtropical jets, particularly over the Pacific Ocean; and (3) estimating tropical precipitation, especially over oceans, using observational and budget methods. A summary list of our most significant accomplishments in the past year is given.

Vincent, Dayton G.↗
NASA NTRS1980

Acoustic behavior of fibrous bulk materials

A semiempirical model of the acoustic behavior of fibrously constructed bulk materials of Hersh and Walker (1979) is generalized to take into account the filtration or removal of particles by fibrous mats and heat conduction between the material fibers and the surrounding fluid. Equations governing the propagation and attenuation of sound waves in a fibrous material are derived on the basis of a solution of the Navier Stokes equation for momentum conservation and a one-dimensional model of heat transfer between the sound field and the fibers. A comparison of the two propagation constants and material impedance specified by the model and experimental measurements for Kevlar 29 indicates the accuracy of the model over a wide range of sound frequencies, material porosities and specimen thickness. It is also found that heat transfer effects are relatively unimportant, while the attenuation constants and material characteristic impedance are a function of fiber orientation relative to sound field propagation direction.

Hersh, A. S.↗
DOE OSTI2024

Evaluating Heat Extraction Performance of Closed-Loop Geothermal Systems with Thermally Conductive Enhancements in Conduction-Only Reservoirs: Preprint

We investigated the impact of disc-shaped and linear thermal enhancements on performance of closed-loop geothermal systems using steady-state and transient COMSOL numerical simulations. Thermal enhancements refer to thermally conductive material introduced in the rocks surrounding a wellbore to compensate for the relatively low rock thermal conductivity and increase heat extraction. Materials proposed include composite cements with thermally conductive additives such as metals or graphite to obtain overall cement thermal conductivities of 10 to 100x the rock thermal conductivity. Our simulation approach was modeling a subsegment (e.g., 100 m long) of a wellbore with idealized (e.g., perfect disc-shaped) thermal enhancements, and we report results as relative increase in heat extraction with respect to systems without thermal enhancements. We did not evaluate technical feasibility or cost of installing thermal enhancements. Simulations indicate that for 5-mm thick disc-shaped thermal enhancements with 5-m radius, repeated every 1 m along the wellbore and with thermal conductivity of 100x the rock thermal conductivity, the thermal output increases roughly 20% with respect to a closed-loop system without thermal enhancements. For fishbone structure thermal enhancements, we estimate a roughly 7% increase in thermal output for 5-m long radially outward pipes with half the radius of the main wellbore, repeated every 2.5 m, with thermal conductivity of 100x the rock thermal conductivity.

advanced geothermal systems↗
DOE OSTI2025

FY-25 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the system level modeling using RELAP5-3D of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY25. This year’s work focuses on a new tank configuration where the inlet of the tank was lowered in elevation by 45”. The stability boundaries of the NSTF are thoroughly studied and stability maps are constructed based on the stability and the oscillation patterns of the system. Five distinct operational modes are identified, namely single-phase liquid, uniform double peak oscillations, uniform sinusoidal oscillations, stable two-phase flow, and non-uniform oscillations. Next, the riser inlet throttling case of experimental test Run-104 is simulated with the RELAP5 model where good agreement is obtained between the model and the experimental data. The simulation also highlights the effects of backflow of water from the tank to the upper region of the chimney. Additionally, the decay heat removal test of Run-99 is simulated with the RELAP5 model. Comparison is carried out between this run and a similar run with the mid-tank inlet of Run-74 performed in FY22. With the lower tank inlet, the RELAP5 model is able to predict the experimental data more accurately than the previous mid tank inlet configuration. The discrepancy in model prediction accuracy highlights the non-symmetrical spatial effects in the tank that would otherwise be more easily captured with higher fidelity models. Lastly, two exploratory studies are conducted to investigate the behaviors of the NSTF when 1) heating is provided to the downcomer and 2) a bypass channel is added between the horizontal chimney section to the downcomer.

42 ENGINEERING↗
DOE OSTI2026

Thermal-Hydraulic-Mechanical Modeling and Simulation of Sodium-Potassium–Cooled MARVEL Microreactor Core

The U.S. Department of Energy's Microreactor Program, with Idaho National Laboratory's development of a nuclear microreactor applications test bed named MARVEL, aims to support R&D for the deployment of small, transportable reactors across civilian, industrial, and defense sectors. The MARVEL microreactor, an 85-kWth thermal fission reactor, incorporates TRIGA nuclear fuel and a sodium-potassium eutectic as its primary coolant, designed for safety and efficiency, with natural circulation eliminating the risk of critical heat flux conditions. The reliance on natural circulation for primary cooling means the reactor avoids using fuel spacers to minimize core pressure drop, which could disrupt the primary coolant's natural flow. However, the reactor core’s tight P/D ratio of 1.05, in the absence of fuel spacers, could pose a risk of fuel rod contact and increased peak cladding temperatures. To ensure the reactor safety, this study conducted computational modeling and simulations to investigate the reactor's thermal-hydraulic-mechanical characteristics, including the reactor core heat transfer coefficients, the potential for rod-to-rod contact, and assessed its impact on peak cladding temperature and overall reactor safety. The computational analyses of the MARVEL microreactor core revealed that the thermal deformation of fuel rods under worst-case scenario may lead to the fuel rod contact, but the peak cladding temperatures will remain significantly lower than the safety criteria, ensuring the safety operation reactor without fuel spacers under normal operating conditions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
NASA NTRS1994

Multispectral observations of chromospheric evaporation in the 1991 November 15 X-class solar flare

We analyze simultaneous H(alpha) images and spectra (from Mees Solar Observatory), and soft and hard X-ray images and spectra (from YOHKOH) during the early phase of an X1.5/3B flare. We investigate the morphological relationship between chromospheric downflows, coronal upflows, and particle precipitation sites, and the energetic relationship between conductive heating, nonthermal particle heating, and the chromospheric response. We find that the observations consistently fit the chromospheric evaporation model. In particular, we demonstrate that the observed upflowing coronal and downflowing chromospheric plasma components originate in the same locations, and we show that our unique set of optical and X-ray observations can clearly distinguish between conductively driven and electron beam driven evaporation.

Wulser, Jean-Pierre↗
DOE OSTI2024

Modeling the Impacts of Material Properties on Oscillatory Neuron Behavior

In this study, neuromorphic computing, which mimics the functions of biological brains, offers improvements in both latency and energy efficiency over typical von Neumann computing architectures. Spiking neural networks can be especially power-efficient because they encode information temporally and can use more sparse electrical inputs. Here, we study the design of volatile memristors (variable resistors with memory) for neuronal devices, with particular consideration toward the feasibility of all-on-chip oscillation using built-in capacitance. We use circuit simulations to model the behavior of oscillator neurons with a range of realistic material properties. We find that energy inputs increase with insulating-phase resistivity, thermal conductivity, and device aspect ratio. However, we also find that the minimum capacitance needed for oscillation decreases with increasing insulating-phase resistivity, which opposes the constraints for power efficiency. Based on published data on NbO 2 , VO 2 , and EuNiO 3 , we find that existing materials can be engineered for all-on-chip spiking using their parasitic capacitance.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hydro-mechanical behavior of heated bentonite buffer for geologic disposal of high-level radioactive waste: A bench-scale X-ray computed tomography investigation

The behavior of heated bentonite buffer is critical for the security and long-term performance of a geological repository for high-level radioactive waste (HLW). While laboratory column experiments have been conducted to investigate compacted bentonite and coupled THMC (thermal-hydro-mechanical and chemical) processes for a moderate temperature range of up to 100 °C, data for a higher temperature range are limited. Understanding bentonite behavior and coupled THMC processes under higher temperatures (e.g., up to 200 °C) could allow for a more economic repository design and would expand the data and knowledge base for more reliable modeling. In this study, a bench-scale experiment was conducted in a compacted bentonite column experiencing both heating up to 200 °C in the center and hydration from a sand-clay boundary surrounding the column. During the experiment run for 1.5 years, frequent X-ray computed tomography (CT) scanning of bentonite provided insights into the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. After the experiment, a comprehensive post-dismantling characterization of bentonite samples was conducted. Results showed that the bentonite hydration was axi-symmetrical despite the initial heterogeneity due to packing, confirming the ability of bentonite to seal fast flow/transport paths. Further, compared to a non-heated control experiment, the heated column showed greater CT density variations along the radial distance, indicating that homogenization of bentonite might be more difficult if a temperature gradient is maintained in the repository. Precipitation of an anhydrite layer occurred in the inner hot zone, pointing to potential concerns about salt precipitation causing canister corrosion. Ultimately, the experiments provided a high-resolution window into the strongly dynamic and coupled behavior of bentonite exposed to heating, hydration and swelling, which will be valuable for improving modeling of coupled processes, especially for the early state of a HLW repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2016

Impacts of Snow Darkening by Absorbing Aerosols on Eurasian Climate

The deposition of absorbing aerosols on snow surfaces reduces snow-albedo and allows snowpack to absorb more sunlight. This so-called snow darkening effect (SDE) accelerates snow melting and leads to surface warming in spring. To examine the impact of SDE on weather and climate during late spring and early summer, two sets of NASA GEOS-5 model simulations with and without SDE are conducted. Results show that SDE-induced surface heating is particularly pronounced in Eurasian regions where significant depositions of dust transported from the North African deserts, and black carbon from biomass burning from Asia and Europe occur. In these regions, the surface heating due to SDE increases surface skin temperature by 3-6 degrees Kelvin near the snowline in spring. Surface energy budget analysis indicates that SDE-induced excess heating is associated with a large increase in surface evaporation, subsequently leading to a significant reduction in soil moisture, and increased risks of drought and heat waves in late spring to early summer. Overall, we find that rainfall deficit combined with SDE-induced dry soil in spring provide favorable condition for summertime heat waves over large regions of Eurasia. Increased frequency of summer heat waves with SDE and the region of maximum increase in heat-wave frequency are found along the snow line, providing evidence that early snowmelt by SDE may increase the risks of extreme summer heat wave. Our results suggest that climate models that do not include SDE may significantly underestimate the effect of global warming over extra-tropical continental regions.

GEOS-5↗
DOE OSTI2022

Computational Methodology to Simulate Pyrometallurgical Processes in a Secondary Lead Furnace

Pyrometallurgical recovery of nonferrous metals involves a combination of thermally intensive transformations during exothermic gas-phase reactions, endothermic decomposition of solid charge, and melting of simpler solids. In the recovery of secondary lead, simultaneous thermal effects in a reverberatory-style furnace cause a melt pool to accumulate at the bottom, with lighter solids (slag) floating above and gaseous products from decomposition of the charge diffusing through the gas–slag interface. Species from oxy-fuel combustion of natural gas, species profiles from smelting reactions, and the formation of a melt pool consisting primarily of lead are simulated via a time-averaged formulation. Predictions of outflow are compared with preset inflow profiles to ensure conservation of mass. Thermal profiles for solid, liquid, and gas phases are presented by species. A novel method is implemented to model the latent heat of fusion using a heterogeneous chemical reaction. The simulation is conducted in Simcenter STAR-CCM+ v. 16.02.009-R8.

Rao, Vivek↗
NASA NTRS1987

The loss and depth of CO2 ice in comet nuclei

The chemical differentiation of a homogeneous, initially unmantled cometary nucleus whose composition encompasses water ice, unclathrated CO2 ice, and silicate dust, is presently simulated by an analytical model to test whether the selective sublimation of any free CO2 ice in the new comet will produce a surface layer of water ice and dust covering the undifferentiated core. The model encompasses the effects of nucleus rotation, arbitrary rotatiuon axis orientation, latitude, nucleus heat conduction, and CO2 gas outflow restriction by the water ice and dust layer.

Fanale, Fraser P.↗
NASA NTRS1997

Hierarchical Parallelism in Finite Difference Analysis of Heat Conduction

Based on the concept of hierarchical parallelism, this research effort resulted in highly efficient parallel solution strategies for very large scale heat conduction problems. Overall, the method of hierarchical parallelism involves the partitioning of thermal models into several substructured levels wherein an optimal balance into various associated bandwidths is achieved. The details are described in this report. Overall, the report is organized into two parts. Part 1 describes the parallel modelling methodology and associated multilevel direct, iterative and mixed solution schemes. Part 2 establishes both the formal and computational properties of the scheme.

Padovan, Joseph↗
NASA NTRS2002

Graphite Ablation and Thermal Response Simulation Under Arc-Jet Flow Conditions

The Two-dimensional Implicit Thermal Response and Ablation program, TITAN, was developed and integrated with a Navier-Stokes solver, GIANTS, for multidimensional ablation and shape change simulation of thermal protection systems in hypersonic flow environments. The governing equations in both codes are demoralized using the same finite-volume approximation with a general body-fitted coordinate system. Time-dependent solutions are achieved by an implicit time marching technique using Gauess-Siedel line relaxation with alternating sweeps. As the first part of a code validation study, this paper compares TITAN-GIANTS predictions with thermal response and recession data obtained from arc-jet tests recently conducted in the Interaction Heating Facility (IHF) at NASA Ames Research Center. The test models are graphite sphere-cones. Graphite was selected as a test material to minimize the uncertainties from material properties. Recession and thermal response data were obtained from two separate arc-jet test series. The first series was at a heat flux where graphite ablation is mainly due to sublimation, and the second series was at a relatively low heat flux where recession is the result of diffusion-controlled oxidation. Ablation and thermal response solutions for both sets of conditions, as calculated by TITAN-GIANTS, are presented and discussed in detail. Predicted shape change and temperature histories generally agree well with the data obtained from the arc-jet tests.

Chen, Y.-K.↗
NASA NTRS2012

Plasma Actuators for Turbomachinery Flow Control

This report is Part I of the final report of NASA Cooperative Agreement contract no. NNX07AC02A. The period of performance was January 1, 2007 to December 31, 2010. This report includes the project summary, a list of publications and reprints of the publications that appeared in archival journals. Part II of the final report includes a Ph.D. dissertation and is published separately as NASA/CR-2012-2172655. The research performed under this project was focused on the operation of surface dielectric barrier discharge (DBD) devices driven by high voltage, nanosecond scale pulses plus constant or time varying bias voltages. The main interest was in momentum production and the range of voltages applied eliminated significant heating effects. The approach was experimental supplemented by computational modeling. All the experiments were conducted at Princeton University. The project provided comprehensive understanding of the associated physical phenomena. Limitations on the performance of the devices for the generation of high velocity surface jets were established and various means for overcoming those limitations were proposed and tested. The major limitations included the maximum velocity limit of the jet due to electrical breakdown in air and across the dielectric, the occurrence of backward breakdown during the short pulse causing reverse thrust, the buildup of surface charge in the dielectric offsetting the forward driving potential of the bias voltage, and the interaction of the surface jet with the surface through viscous losses. It was also noted that the best performance occurred when the nanosecond pulse and the bias voltage were of opposite sign. Solutions include the development of partially conducting surface coatings, the development of a semiconductor diode inlaid surface material to suppress the backward breakdown. Extension to long discharge channels was studied and a new ozone imaging method developed for more quantitative determination of surface jet properties.

Miles, Richard, B↗
NASA NTRS2019

Heat Analysis Manager (HAM), a Thermal Desktop API Based Heat Map Generation Software

Thermal engineers often create custom heat maps to analyze their thermal model. However, generating a heat map is difficult because thermal simulation only readily provide attributes of simulated nodes such as temperature, capacitance, heat generation, and a network of conductances. Heat flow values are a derived quantity from the nodal attributes, and the data processing and management of heat flow between nodes quickly become difficult for large models. Deriving a network of heat flow values requires vast amount of calculations and data handling, heat map generation process generally suffers from processing speed, loss of accuracy, and/or presentation of data in a useful format. Heat Analysis Manager (HAM) is a Thermal Desktop (TD) based free multi-purpose tool developed to aid thermal engineers in analyzing their thermal model, including a heat map generation functionality. HAM’s heat map generator retains accuracy and fast processing speed by utilizing TD’s application programming interface (API) and built-in TD’s “Qflow from Results.” Furthermore, HAM’s heat map output is presented in an easily customizable format in Excel, allowing users to create various custom visual heat maps. A full description of how HAM utilizes TD’s API to create a customizable heat map is provided. A simple model demonstration is included along with step-by-step procedures on creating custom heat maps. HAM’s heat map result has been verified against TD’s and other heat map generation software, and verification methods are also included.

Thermal desktop↗
NASA NTRSDate not supplied

Harness Thermal Heat Loss Measurement for a Lunar Surface-Deployed LEMS Artemis III Payload

The Lunar Environment Monitoring Station (LEMS) is an autonomous, survive-the-lunar-night seismic suite to be deployed on the Lunar surface by the Artemis III crew and designed to operate continuously for two years. It will see the extremes of the Lunar south pole thermal environment where surface temperatures range from -200°C to +20°C and where nighttime duration is at least 354 hours. Given power and mass constraints, the thermal system is limited to 2.5 Watts of heat during the lunar night. The electrical harnessing named the Signal and Power Passthrough (SAPP) was designed to minimize heat loss while meeting power and signal integrity requirements. To mitigate risk due to uncertainty associated with the harnessing materials, routing, and tie-downs, a flight-like thermal conductance test was performed to measure the heat loss. The test methodology, results, and model correlation are presented.

TVAC↗