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 451 records · Page 25

Wind-US Code Physical Modeling Improvements to Complement Hypersonic Testing and Evaluation

This report gives an overview of physical modeling enhancements to the Wind-US flow solver which were made to improve the capabilities for simulation of hypersonic flows and the reliability of computations to complement hypersonic testing. The improvements include advanced turbulence models, a bypass transition model, a conjugate (or closely coupled to vehicle structure) conduction-convection heat transfer capability, and an upgraded high-speed combustion solver. A Mach 5 shock-wave boundary layer interaction problem is used to investigate the benefits of k- s and k-w based explicit algebraic stress turbulence models relative to linear two-equation models. The bypass transition model is validated using data from experiments for incompressible boundary layers and a Mach 7.9 cone flow. The conjugate heat transfer method is validated for a test case involving reacting H2-O2 rocket exhaust over cooled calorimeter panels. A dual-mode scramjet configuration is investigated using both a simplified 1-step kinetics mechanism and an 8-step mechanism. Additionally, variations in the turbulent Prandtl and Schmidt numbers are considered for this scramjet configuration.

Georgiadis, Nicholas J.↗

Optimized Heat Interception for Cryogen Tank Support Structure

We consider means for using the cooling available in boil-off gas to intercept heat conducted through the support structure of a cryogen tank. A one-dimensional model of the structure coupled to a gas stream gives an analytical expression for heat leak in terms of flow rate for temperature independent properties and laminar flow. A numerical model has been developed for heat transfer on a thin cylindrical tube with an attached vent line. The model is used to determine the vent path layout that will minimize heat flow into the cryogen tank. The results are useful for a number of applications, but the one of interest in this study is the minimization of the boil-off in large cryopropellant tanks in low Earth and low lunar orbit.

Canavan, E. R.↗

Optimized Heat Interception for Cryogen Tank Support

We consider means for using the cooling available in boil-off gas to intercept heat conducted through the support structure of a cryogen tank. A one-dimensional model of the structure coupled to a gas stream gives an analytical expression for heat leak in terms of flow rate for temperature independent-properties and laminar flow. A numerical model has been developed for heat transfer on a thin cylindrical tube with an attached vent line. The model is used to determine the vent path layout that will minimize heat flow into the cryogen tank. The results are useful for a number of applications, but the one of interest in this study is the minimization of the boil-off in large cryopropellant tanks in low Earth and low lunar orbit.

Canavan, Edgar R.↗

Analytic Theory of Titans Schumann Resonance: Constraints on Ionospheric Conductivity and Buried Water Ocean

This study presents an approximate model for the atypical Schumann resonance in Titan's atmosphere that accounts for the observations of electromagnetic waves and the measurements of atmospheric conductivity performed with the Huygens Atmospheric Structure and Permittivity, Wave and Altimetry (HASI-PWA) instrumentation during the descent of the Huygens Probe through Titan's atmosphere in January 2005. After many years of thorough analyses of the collected data, several arguments enable us to claim that the Extremely Low Frequency (ELF) wave observed at around 36 Hz displays all the characteristics of the second harmonic of a Schumann resonance. On Earth, this phenomenon is well known to be triggered by lightning activity. Given the lack of evidence of any thunderstorm activity on Titan, we proposed in early works a model based on an alternative powering mechanism involving the electric current sheets induced in Titan's ionosphere by the Saturn's magnetospheric plasma flow. The present study is a further step in improving the initial model and corroborating our preliminary assessments. We first develop an analytic theory of the guided modes that appear to be the most suitable for sustaining Schumann resonances in Titan's atmosphere. We then introduce the characteristics of the Huygens electric field measurements in the equations, in order to constrain the physical parameters of the resonating cavity. The latter is assumed to be made of different structures distributed between an upper boundary, presumably made of a succession of thin ionized layers of stratospheric aerosols spread up to 150 km and a lower quasi-perfect conductive surface hidden beneath the non-conductive ground. The inner reflecting boundary is proposed to be a buried water-ammonia ocean lying at a likely depth of 55-80 km below a dielectric icy crust. Such estimate is found to comply with models suggesting that the internal heat could be transferred upwards by thermal conduction of the crust, while convective processes cannot be ruled out.

TITAN↗

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters for an initial thermal non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two codes currently used at Glenn Research Center for Stirling modeling are Fluent and CFD-ACE. The codes porous-media models are equilibrium models, which assume solid matrix and fluid are in thermal equilibrium. This is believed to be a poor assumption for Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, use non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location and time during the cycle. Experimentally based information was used to define: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity was also estimated. Determination of model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Converter (TDC), which uses a random-fiber regenerator matrix. Emphasis is on use of available data to define empirical parameters needed in a thermal non-equilibrium porous media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates.

Tew, Roy C.↗

Numerical Evaluation of Effective Thermal Conductivity of PCM with Metal Foam Incorporating Buoyancy Effects for Thermal Energy Storage

The thermal energy storage (TES) system has the capability to efficiently preserve thermal energy directly derived from the energy source, minimizing any conversion losses. Especially latent heat storage offers distinct advantages, including a substantial increase in energy storage density and minimization of temperature fluctuations within the plants. However, the phase change material (PCM) employed in latent heat storage has low thermal conductivity. Consequently, various studies are being conducted to enhance heat transfer. One approach to enhance heat transfer involves utilizing metal foam to maximize the heat transfer area. However, modeling metal foam with its intricate structure is a challenging task in numerical analysis. For this reason, ongoing research focuses on simplifying the modeling of metal foam. Nevertheless, fully encompassing all the characteristics of actual metal foam proves to be a challenging task for the simplified analytical model. The objective of this paper is to interpret the simple lattice metal foam analysis model from the perspective of behavior induced by buoyancy. When comparing the analysis results of solid PCM and liquid PCM with the same thermal conductivity under changes in porosity and gravity direction, we conducted an analysis to discern the trends in effective thermal conductivity that are overestimated due to convection. In the analysis, a constant heat flux of 10 kW and a constant surface boundary condition of 350 K were applied, and a sensitivity study regarding the mesh was conducted. The results indicate that, from the perspective of gravity in the simple lattice model, the solid analysis yields an effective thermal conductivity 29-47% higher compared to the liquid analysis. Additionally, as porosity increases, there is an observed increase of 24-33% in effective thermal conductivity.

25 ENERGY STORAGE↗

Multiphysics simulation of TRISO fuel compacts and the effects of homogenization on silver release predictions

This work studies the impact of explicit and homogenized modeling approaches on the multiphysics simulation of TRistructural ISOtropic (TRISO) fuel compacts in prismatic High Temperature Gas Reactors (HTGRs) and the silver release predictions. TRISO fuel compacts exhibit complex double heterogeneity that significantly affects heat conduction, neutron transport, and fission products release. In this work, we use Cardinal, a multiphysics tool based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to couple neutron transport and heat conduction. OpenMC is used for neutron transport and the MOOSE heat transfer module is used for the heat conduction. Then, we use BISON for the silver release predictions. Two different modeling approaches—explicit modeling of individual TRISO particles and homogenized representation using effective thermal properties—are compared at high TRISO packing fractions (20% and 40%) across varying power densities. Results demonstrate that homogenization significantly under-predicts peak temperatures, for the case of high power/TRISO, there is a difference of 83.73 K in the maximum temperature. Additionally, homogenization underestimates the silver release predictions compared to explicit modeling. Specifically, the temperature and power density differences lead to significant differences in silver release predictions. This work demonstrates the importance of accurately modeling heterogeneity of the TRISO particles to reliably predict fission product release and assess reactor safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Coal-derived conductive pavement for winter de-icing: prototype, modeling, and simulation

Existing pavement de-icing methods result in high installation and maintenance costs, traffic delays, excessive weariness and corrosion, and negative environmental and safety impacts. To overcome these challenges, this paper presents an innovative pathway to designing and constructing smart self-heating pavements using a low-cost coal-char bearing asphalt material. The conductive asphalt incorporates coal char, i.e., a key byproduct of the coal pyrolysis process, into the Stone Mastic Asphalt (SMA) mixture. This asphalt containing coal-derived solid carbon exhibits highly tailorable electrical conductivity, satisfactory mechanical and thermophysical properties, and superior cost-efficiency as compared to other conductive pavement materials. The de-icing performance was also demonstrated by laboratory experiments on a bench-scale prototype. Furthermore, an efficient thermal network model was developed and validated by experiments to investigate the transient thermal behavior and energy performance of the Ohmic heating pavement system. Furthermore, the whole-year energy simulation case studies were conducted on a bridge pavement with an annual energy use of 24.6–1444.8 kWh/m 2 , showcasing its potential in field applications across cool humid, cold humid, and subarctic/arctic climate zones.

42 ENGINEERING↗

Multiphase Modeling in a Parallel Plate Fluidized Bed Receiver for Concentrating Solar Power

A novel high temperature particle solar receiver is developed by using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a packed bed at the lower end, leaving the system from the bottom at a constant flow rate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built and a near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2K.

CSP↗

Preliminary Thermal-hydraulics and Selected Safety Analysis for Holos-Quad Reactor Design

This report presents thermal hydraulics and safety analyses of the 10 MWe Holos-Quad micro-reactor design developed by HolosGen LLC. These analyses were executed under steady-state and station black-out (SBO) transient conditions using the System Analysis Module (SAM), with a focus on ensuring peak fuel temperatures do not exceed the safety thresholds of 1,250°C for steady state normal operating conditions and 1,600°C for transient conditions. For steady-state temperature predictions, a coupled 1D fluid-to-3D solid heat conduction methodology was used. The steady-state simulations included a model of the central core region where power peaking occurs to predict the maximum steady-state fuel temperature, and a coarse-mesh full-core model to provide the initial condition for the consequent transient SBO decay heat removal simulation. The maximum fuel temperature predicted from the central core region simulations under the normal operating condition was 1,222°C, which is below the design limit 1,250°C. To simulate a SBO transient, a full core model was necessary because of the non-axial symmetric core design and heat loss from the outer core structures to the environment. To reduce the required computational costs, a novel scaled subassembly approach was used to model the full core domain. With the full core steady-state temperature distribution provided as the initial condition, the passive decay heat removal of the core was simulated with the scaled subassembly full core model under SBO transient conditions. This analysis showed that the peak fuel temperatures of the Holos-Quad core remained below their steady-state values during the transient, therefore much lower than the safety limit 1,600°C. This is mainly due to the large thermal inertia of the graphite matrix, relatively low power density, and the large surface-to-volume-ratio of the core. The preliminary analyses presented in this report confirm the inherent safety of the Holos-Quad micro-reactor thermal-hydraulics design, as the peak fuel temperatures under both normal and off-normal operational conditions remain within design and safety limits.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Cryogenic Selective Surfaces

There are many challenges involved in deep-space exploration, but several of these can be mitigated, or even solved, by the development of a coating that reflects most of the Sun’s energy, yet still provides far-infrared heat emission. Such a coating would allow non-heat-generating objects in space to reach cryogenic temperatures without using an active cooling system. This would benefit deep-space sensors that require low temperatures, such as the James Webb Telescope focal plane array. It would also allow the use of superconductors in deep space, which could lead to magnetic energy storage rings, lossless power delivery, or perhaps a large-volume magnetic shield against galactic cosmic radiation. However, perhaps the most significant enablement achieved from such a coating would be the long-term, deep space storage of cryogenic liquids, such as liquid oxygen (LOX). In our Phase I NIAC study, we realized that a combination of scattering particles and a silver backing could yield a highly effective, very broadband, reflector that could potentially reflect more than 99.9% of the Sun’s irradiant power. We developed a sophisticated model of this reflector and theoretically showed that cryogenic temperatures could be achieved in deep space at one astronomical unit (1 AU) from the Sun. We showed how this new reflector could minimize heat conduction into the cryogenic tanks by coating the tank support struts. We then modelled a strawman architecture for a mission to Mars, using a coated LOX tank, coated struts, and infrared shields, to show that with our new coating it would be possible to maintain liquid oxygen passively. As a result of this work a patent application was generated and a paper published in Optics Letters. Our Phase II NIAC study had two primary goals, to develop a rigid version of the cryogenic selective surface proposed in Phase I and to test its performance in a simulated deep space environment. During the first year of the project the work concentrated on developing rigid tiles of BaF2, leading to tiles as large as 4 inches in diameter that transmitted very little visible light. In addition, during the first year a simulated deep space environment was created using a vacuum chamber and cryocooler. Using this facility, we showed that our BaF2 tiles absorbed less than ¼% of 375 nm radiation, a significant milestone for the work. During the second year of the project, we continued to develop the BaF2 tiles and we put significant effort into the construction of a deep space environment where we could project simulated solar radiation onto a sample. In the spring of 2018, we conducted our first solar simulator test with BaF2 and saw about 3.6% absorption. This is better than the state-of-the-art, but disappointing since predictions were for much lower absorption. We, erroneously, attributed this absorption to water retention by the BaF2, and decided to change materials. We considered several oxides and settled on yttrium oxide (Y2O3) for further development, because it is broadband, lightweight, has high index, and is hydrophobic. In July 2018 we conducted our first test of a rigid tile of Y2O3 in the simulated deep space environment and saw significant absorption again. We then realized that the issue was not water, but mid-wave radiation passing through the tile and being absorbed by the temperature sensor and the varnish used to hold it in place. We wrapped the sensor in silver foil, re-ran the test, and saw much lower absorption; only 1.1%. We then re-ran the BaF2 tile and saw 1.4% absorption. These values are almost adequate to maintain LOX in deep space, but we suspect that there are still issues in our test apparatus; we suspect thermocouple wires may be absorbing radiation. Further, post-NIAC, testing will better determine the performance of our new solar reflector. In order to restrict the size of this report, we will only briefly describe topics that we have previously published, allowing us to devote more time to new material. So minimal material will be devoted to modeling the material and deep space cryogenic storage, while longer sections will cover our material development, simulated deep space testing, and new applications. The Launch Service Program (LSP) requested that we explore ways to use this new coating to maintain LOX in low Earth Orbit and that work is described. In addition, the Nuclear Thermal Propulsion (NTP) Program asked us to explore ways to reduce the heat load for liquid hydrogen, resulting in the development of a spray-on version of the coating that should significantly improve in-space multi-layer insulation performance.

Robert C. Youngquist,↗

A comparison of laboratory measured temperatures with predictions for a spar/skin type aircraft structure

A typical spar/skin aircraft structure was heated nonuniformly in a laboratory and the resulting temperatures were measured. The heat transfer NASTRAN computer program was used to provide predictions. Calculated temperatures based on a thermal model with conduction, radiation, and convection features compared closely to measured spar temperatures. Results were obtained without the thermal conductivity, specific heat, or emissivity with temperature. All modes of heat transfer (conduction, radiation, and convection) show to affect the magnitude and distribution of structural temperatures.

Jenkins, J. M.↗

Theoretical study of the seasonal behavior of the global ionosphere at solar maximum

The seasonal behavior of the global ionosphere was studied using a time-dependent three-dimensional physical model (developed by Shunk and his coworkers) of the ionosphere at altitudes between 120 and 800 km. This model accounts for field-aligned diffusion, cross-field electrodynamic drifts both the equatorial region and at high latitudes, interhemispheric flow, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow, and local heating and cooling processes. The model studies were carried out for both June and December solstice conditions at solar maximum and for low geomagnetic activity. The ionospheric features predicted by the model agreed qualitatively with the available measurements.

Sojka, J. J.↗

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area: Preprint

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs drilled in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated including a single, double, and multi-lateral design. The subsurface within area is characterized by deep, hot (> 200degreesC) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lies within a target depth of 6 km where temperatures are estimated to approach 300degreesC. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down technoeconomic analysis model. The results show that the performance and cost optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60degreesC/km. By injecting 20 degreesC-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth can be achieved, resulting in an LCOE and LCOH of $136/MWh and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗

In situ temperature profile measurements with high-energy X-rays as a probe of optical floating zone crystal growth environment

The ability of optical floating zone (OFZ) furnaces to rapidly produce large single crystals of complex emerging materials has had a transformative effect on many scientific fields that require samples of this type. However, the crystal growth process within the OFZ furnace is not well understood owing to the challenges involved in monitoring the high-temperature crystal growth process. Novel beamline-compatible optical furnaces that approximate the inhomogeneous growth environment within an OFZ furnace have been fabricated and tested in high-energy synchrotron beamlines. In this paper, it is demonstrated that temperature profiles can be effectively extracted from powder diffraction data collected on polycrystalline ceramic rods heated at their tip. Furthermore, these measured temperature profiles can be accurately reproduced using a heat-transfer model that accounts for solid-state thermal conduction, partial sample lamp power absorption, convective air cooling and radiative cooling, allowing key thermal parameters such as thermal conductivity to be extracted from experimental data.

36 MATERIALS SCIENCE↗

The structure and heating of the chromosphere-corona transition region

The structure and heating (or energy balance) of the transition region and the role of the transition region in the structure and heating of the solar atmosphere as a whole are investigated. The features of the structure of the atmosphere and radiative energy losses of the atmosphere are summarized. A static, planar model of the solar temperature which has a temperature profile representative of the actual solar atmosphere is considered. Then a static, planar model of the transition region which is heated by thermal conduction from the corona and cooled by radiative losses is developed. A general conclusion is that the temperature profile of the transition region and lower corona results primarily from the energy balance of the corona, while the number density is determined by the energy balance of the transition region.

Moore, R. L.↗