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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin↗
DOE OSTI2023

Testing of an Optical Fiber--Based Gamma Thermometer in the High Flux Isotope Reactor Gamma Irradiation Facility

This report describes the design, thermal modeling, and gamma irradiation testing of an optical fiber–based gamma thermometer (OFBGT), which was irradiated in the High Flux Isotope Reactor (HFIR) Gamma Irradiation Facility (GIF). OFBGTs are a promising technology for application in nuclear reactors because they can provide a distributed measurement of gamma ray heating rate, unlike thermocouple-based gamma thermometers, which are fixed in-core sensors that can be used in boiling water reactors to calibrate local power range monitors. OFBGTs measure gamma ray heating rate by measuring the temperature difference between a pair of optical fibers; one fiber is in thermal contact with a heat sink (usually the reactor coolant), and the other is in thermal contact with a thermally isolated mass. The device can be calibrated with a heating wire within the thermal mass. The OFBGT that was designed and fabricated at Oak Ridge National Laboratory can measure distributed gamma ray heating rate over an effective measurement length of 61 cm, and the outer diameter of the sensor is 12.7 mm, giving the prototypical sensor design a relatively small footprint. The sensor housing is backfilled with Ar to ensure a well-predicted thermal response that is not affected by humidity or chemical interactions during operation. For calibration, the sensor design uses a Ni–Cr wire, which can be supplied with currents from 0 to 1 A to capture the wide range of potential gamma ray heating rates expected in HFIR’s spent fuel elements. The OFBGT was thermally modeled analytically and numerically; both models account for temperature-dependent thermal conductivities of the materials and show good agreement. The thermal response of the sensor inside spent HFIR fuel elements was simulated for times up to 1 year after discharge of the fuel element. Out-of-pile open-air tests indicated that the steady-state response of the sensor matches modeled results within experimental uncertainty. Calibration tests were performed using electrical heating in the HFIR spent fuel pool, above the fuel elements, to establish a relationship between the difference in spectral shift measured by optical fibers located inside and outside the OFBGT and the applied electrical heating. Subsequently, the OFBGT was placed within the the spent fuel element from HFIR cycle 501 to measure the spatial profile of the gamma heating rates. Results showed good agreement between the theoretical and measured gamma dose rate profiles, with maximum deviations of ~10% or less.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2010

The Prediction of Broadband Shock-Associated Noise from Dualstream and Rectangular Jets Using RANS CFD

Supersonic jets operating off-design produce broadband shock-associated noise. Broadband shock-associated noise is characterized by multiple broadband peaks in the far-field and is often the dominant source of noise towards the sideline and upstream direction relative to the jet axis. It is due to large scale coherent turbulence structures in the jet shear layers interacting with the shock cell structure. A broadband shock-associated noise model recently developed by the authors predicts this noise component from solutions to the Reynolds averaged Navier-Stokes equations using a two-equation turbulence model. The broadband shock-associated noise model is applied to dualstream and rectangular nozzles operating supersonically, heated, and off-design. The dualstream jet broadband shock-associated noise predictions are conducted for cases when the core jet is supersonic and the fan jet is subsonic, the core jet is subsonic and the fan jet is supersonic, and when both jet streams operate supersonically. Rectangular jet predictions are shown for a convergent-divergent nozzle operating both over- and under-expanded for cold and heated conditions. The original model implementation has been heavily modified to make accurate predictions for the dualstream jets. It is also argued that for over-expanded jets the oblique shock wave attached to the nozzle lip contributes little to broadband shock-associated noise. All predictions are compared with experiments.

Miller, Steven A.↗
NASA NTRS1972

The stable auroral red arc of October 31 - November 1, 1968 and its interaction with the neutral atmosphere

Satellite observations of electron temperature and the topside electron density structure are used to calculate the red arc according to the thermal conduction model. In this model, energy from the magnetosphere flows in the electron gas along geomagnetic field lines into the ionosphere. This energy heats the ambient F-region electrons sufficiently to excite the oxygen atoms to the O I(1D) level by collisional impact giving rise to the 6300 A emission characteristic of the arc. For the 31 October/1 November 1968 red arc, the calculated emission rate, geographical position, and horizontal extent of the red arc are in agreement with the photometric data obtained by airglow observatories. Almost all of the energy conducted into the red arc is ultimately transferred to the neutral gas through elastic and inelastic collisions. This energy drives a large thermal cell and the circulation extends the influence of the arc thousands of kilometers beyond the region of direct heat input. The calculated neutral gas temperature response to electron heating within the arc is small.

Roble, R. G.↗
NASA NTRS1973

Role of convection in the moon.

Conduction solutions of the problem of the moon's thermal history can be unstable to solid state convection. To examine the role of solid convection, the stability of models in which the initial distribution of radioactive heat sources is homogeneous and the initial temperature profile is due to accretional heating is analyzed in detail. Growth rates of instabilities are compared with the appropriate conduction times in order to determine the effective viscosities for which solid convection is a dominant process. It is found that instability growth may not have been rapid enough to prevent the extensive partial melting predicted by conduction theory, but convection after resolidification is likely to have occurred, and thus the moon would have cooled faster than it would if it were purely conducting.

Cassen, P.↗
NASA NTRS2008

Cold Start of a Radiator Equipped with Titanium-Water Heat Pipes

Radiator panels utilizing titanium-water heat pipes are being considered for lunar applications. A traditional sandwich structure is envisioned where heat pipes are embedded between two high thermal conductivity face sheets. The heat pipe evaporators are to be thermally connected to the heat source through one or more manifolds containing coolant. Initial radiator operation on the lunar surface would likely follow a cold soak where the water in the heat pipes is purposely frozen. To achieve heat pipe operation, it will be necessary to thaw the heat pipes. One option is to allow the sunlight impinging on the surface at sunrise to achieve this goal. Testing was conducted in a thermal vacuum chamber to simulate the lunar sunrise and additional modeling was conducted to identify steady-state and transient response. It was found that sunlight impinging on the radiator surface at sunrise was insufficient to solely achieve the goal of thawing the water in the heat pipes. However, starting from a frozen condition was accomplished successfully by applying power to the evaporators. Start up in this fashion was demonstrated without evaporator dryout. Concern is raised over thawing thermosyphons, vertical heat pipes operating in a gravity field, with no wick in the condenser section. This paper presents the results of the simulated cold start study and identifies future work to support radiator panels equipped with titanium-water heat pipes.

Jaworske, Donald A.↗
DOE OSTI2026

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗
NASA NTRS2025

Metallurgical and Thermal Properties of 57Ni-40Ti-3Hf (wt%)

Selected properties of 57Ni-40Ti-3Hf (wt%) were measured to further develop this material for use in aerospace rolling-element bearings. The average grain size and microindentation hardness of the material at room temperature were 90 μm and 549 HV, respectively. From room temperature to 500 °C, electrical resistivity increased from 8.9×10–7 to 11.3×10–7 Ω⋅m, thermal expansion coefficients varied from 1.1×10–5 to 1.2×10–5 mm/mm⋅°C, specific heat increased from 0.44 to 0.53 W⋅sec/g⋅K, thermal diffusivity rose from approximately 0.003 to 0.005 cm2/sec and thermal conductivity increased from approximately 0.09 to 0.19 W/cm⋅K. The behavior of the material can be explained by well-known metallurgical theory. These results are expected to provide a basis for design and modeling of future aerospace components.

rule of mixtures↗
DOE OSTI2024

Heat transfer optimization of uo 2 -mo fuel using genetic algorithms

Two genetic algorithm (GA) methods were applied to thermal finite element models to optimize the heat transfer efficacy of a UO 2 -Mo composite fuel pellet with typical pressurized water reactor fuel geometry. Mo additions to UO 2 have been shown to increase the thermal conductivity, thus reducing centerline temperatures and temperature gradients. Previous studies evaluated uniformly dispersed Mo or continuous Mo internal geometries (e.g., fins, plates, discs) that were selected using engineering intuition. The current study uses two different implementations of the same GA to optimize Mo placement and minimize the fuel temperature with the only constraint being a maximum 10% Mo volume fraction. One approach superimposed Mo line elements onto the monolithic UO 2 pellet model, and the other converted entire UO 2 volume elements to Mo. The former method generated 1D heat transfer connections between nodes, whereas the latter method allowed for the formation of 3D structures. Features of the optimal fuel design produced by the GAs included dispersed Mo near the centerline that shifted the peak fuel temperature outward by 0.6 mm, Mo chains in the high-heat-flux region in the mid-to-outer radial zone, and a large continuous structure that spanned the full radius and height of the pellet and accounted for 87.7 % of the total Mo in the pellet. Analysis of this design indicates that the optimal Mo configuration is a balance between creating continuous heat transfer pathways and optimally dispersing Mo to minimize the heat transfer distance through UO 2 . This architecture ultimately produced an effective thermal conductivity of 11.3 W/m·K under the assumed boundary conditions. This result is higher than any previous values from the literature. In conclusion, potential fabrication methods and challenges are discussed in addition to the implications on fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Impact of heating and cooling loads on battery energy storage system sizing in extreme cold climates

Efficient operation of battery energy storage systems requires that battery temperature remains within a specific range. Current techno-economic models neglect the parasitic loads heating and cooling operations have on these devices, assuming they operate at constant temperature. In this work, these effects are investigated considering the optimal sizing of battery energy storage systems when deployed in cold environments. Here, a peak shaving application is presented as a linear programming problem which is then formulated in the PYOMO optimization programming language. The building energy simulation software EnergyPlus is used to model the heating, ventilation, and air conditioning load of the battery energy storage system enclosure. Case studies are conducted for eight locations in the United States considering a nickel manganese cobalt oxide lithium ion battery type and whether the power conversion system is inside or outside the enclosure. The results show an increase of 42% to 300% in energy capacity size, 43% to 217% in power rating, and 43% to 296% increase in capital cost dependent on location. This analysis shows that the heating, ventilation, and air conditioning load can have a large impact on the optimal sizes and cost of a battery energy storage system and merit consideration in techno-economic studies.

25 ENERGY STORAGE↗
DOE OSTI2023

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

FY20 Verification of BISON Using Analytic and Manufactured Solutions

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution Modeling and Simulation (M&S) tools for modeling of Light Water Reactors (LWRs). This hub, Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including BISON, CTF, MPACT, and MAMBA. BISON is a fuel performance code which models the thermo-mechanical behavior of nuclear fuel using high performance M&S. It is capable of modeling traditional LWR fuel rods, fuel plates, and TRi-structural ISOtropic (TRISO) fuel particles. It can employ three-dimensional Cartesian, two-dimensional axisymmetric cylindrical, or one-dimensional radial spherical geometry. It includes empirical models for a large variety of fuel physics: temperature- and burnup-dependent thermal properties, fuel swelling and densification, fission gas production, cladding creep, fracture, cladding plasticity, and gap/plenum models. This document details a series of code verification test problems that are used to test BISON. These problems add confidence that the BISON code is a faithful representation of its underlying mathematical model. The suite of verification tests are mapped to the underlying conservation equations solved by the code: heat conduction, mechanics, and species conservation. Twenty-two problems are added for the heat conduction solution, two for the mechanics solution, and none for species conservation. Method of Manufactured Solutions (MMS) capability is demonstrated with three problems, and temperature drops across the fuel gap are tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1979

The effect of energetically produced O2/+/ on the ion temperatures of the Martian thermosphere

The heating rates and ion temperatures of the Martian thermosphere resulting from the thermalization of the energetic O2(+) ions produced by the reactions of the solar ionization products CO2(+) and O(+) with neutral particles are calculated and the effects of small magnetic fields on the ion thermal balance are investigated. The energy transfer and transport of energetic ions is modelled by solving the continuity equation for the ions over a series of finite energy cells yielding the equilibrium densities as a function of energy and altitude. Particles which are able to proceed upwards without collisions are dealt with separately from the continuity equation. It is shown that the thermalization of the energetic O2(+) ions can greatly increase ion temperatures above 200 km compared to those calculated for only ambient electron heating. Current solar wind interaction models predict that small horizontal magnetic fields act to restrict the ion thermal conductivity and to increase upper altitude ion temperatures. The combined effects of these processes provide a partial agreement with measurements made by Viking 1.

Rohrbaugh, R. P.↗
NASA NTRS1980

Semiempirical models of chromospheric flare regions

Homogeneous plane-parallel semiempirical flare model atmospheres which reproduce observations in lines and continua of H I, Si I, C I, Ca II, and Mg II have a thin transition zone at the top of the enhanced chromosphere, indicating a significant amount of heating from the zone to the temperature minimum level. The minimum temperature is located deeper and is higher than in the quiet-sun and active-region models. The results do not agree with the particle-heated theoretical models, and it is suggested that the models of Brown (1973) and Henoux and Nakagawa (1977, 1978) do not include an essential term for heat conduction in their energy balance equations. It is concluded that substantial Ly-alpha radiative heating occurs in the upper chromosphere resulting from the conductive energy flux in the transition zone where the Ly-alpha line cools the gas.

Machado, M. E.↗