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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Cryosphere and climate

This chapter will discuss two main issues related to the cryosphere and climate. One is the effect of sea ice and salinity gradients on ocean circulation, and in particular the possible role of sea ice transport on the ocean conveyer belt. The other is the effect of the cryosphere on climate, and in particular in high-latitude warming under increased CO2. In understanding the role of the cryosphere in both cases, it is useful to elucidate two types of toy sea ice models. Neither of these represents reality, but both are useful for illustrating the archetypal features of sea ice that control much of its large-scale behavior. The first model is a simple slab thermodynamic sea ice model as presented by Thorndike. In this model there are no dynamical effects and the thickness of ice is determined by surface heat budget and oceanic heat flux considerations, with the thickness of the ice critically affecting the effective conductivity whereby heat is transferred from the bottom ice boundary to the upper ice boundary. In this model all of the sea ice characteristics are controlled by the vertical heat fluxes from the atmosphere and ocean into the ice. The thickness is controlled by the ice's becoming an effective insulator as it thickens, thus reducing conductive heat loss to the atmosphere. A second model emphasizes the effects of dynamics. It considers the ice pack to be a collection of floes moving in response to synoptic wind fields and ocean currents. These motions create semipermanent leads (open areas) over which ice can grow rapidly.

Hibler, William D., III↗
NASA NTRS1984

Heat transfer and oil flow studies on a single-stage-to-orbit control-configured winged entry vehicle

Results are presented for oil flow and phase change paint heat transfer tests conducted on a 0.006 scale model of a proposed single stage to orbit control configured vehicle. The data were taken at angles of attack up to 40 deg at a free stream Mach number of 10 for Reynolds numbers based on model length of 0.5 x 10 to the 6th power, 1.0 x 10 to the 6th power and 2.0 x 10 to the 6th power. The magnitude and distribution of heating are characterized in terms of angle of attack and Reynolds number aided by an analysis of the flow data which are used to suggest the presence of various three dimensional flow structures that produce the observed heating patterns. Of particular interest are streak heating patterns that result in high localized heat transfer rates on the wing windward surface at low to moderate angles of attack. These streaks are caused by the bow-shock/wing-shock interaction and formation of the wing-shock. Embedded vorticity was found to be associated with these interactions.

Helms, V. T., III↗
NASA NTRS1992

An experimental study of heat transfer in a large-scale turbine rotor passage

An experimental study of the heat transfer distribution in a turbine rotor passage was conducted in a large-scale, ambient temperature, rotating turbine model. Heat transfer was measured for both the full-span suction and pressure surfaces of the airfoil as well as for the hub endwall surface. The objective of this program was to document the effects of flow three-dimensionality on the heat transfer in a rotating blade row (vs a stationary cascade). Of particular interest were the effects of the hub and tip secondary flows, tip leakage and the leading-edge horseshoe vortex system. The effect of surface roughness on the passage heat transfer was also investigated. Midspan results are compared with both smooth-wall and rough-wall finite-difference two-dimensional heat transfer predictions. Contour maps of Stanton number for both the rotor airfoil and endwall surfaces revealed numerous regions of high heat transfer produced by the three-dimensional flows within the rotor passage. Of particular importance are regions of local enhancement (as much as 100 percent over midspan values) produced on the airfoil suction surface by the secondary flows and tip-leakage vortices and on the hub endwall by the leading-edge horseshoe vortex system.

Blair, Michael F.↗
DOE OSTI2025

COMSOL Results for the Nominal Steady-State Operation of the Proposed 95-MW LEU Silicide Core for HFIR Conversion

Engineering design studies are being performed to determine the feasibility of converting the High Flux Isotope Reactor (HFIR) from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel at Oak Ridge National Laboratory. This activity is sponsored by the Office of Reactor Conversion and Uranium Supply (ORCUS) under the auspices of the US Department of Energy National Nuclear Security Administration’s Office of Material Management and Minimization. HFIR is a very high flux, pressurized, light water–cooled and moderated, flux trap–type research reactor with a core made of involute shaped U 3 O 8 /Al cermet fuel plates and coolant channels. HFIR currently operates at a thermal power of 85 MW and supports key national and international missions in neutron scattering, isotope production, materials/fuels irradiation, neutron activation analysis, gamma irradiation, and neutrino research. Advanced multiphysics computational fluid dynamics models have been developed in the COMSOL Multiphysics software to simulate the steady-state operating conditions for the proposed low-and high-density LEU U 3 Si 2 -Al (uranium silicide dispersion) fuel designs. The COMSOL models for HFIR inner and outer fuel element models incorporate various essential inputs and physics such as spatially dependent nuclear heat deposition, multilayer heat conduction, conjugate heat transfer, turbulent flows (using Reynolds-averaged Navier Stokes turbulence models), structural mechanics (thermal–structural interactions and fuel swelling), and oxide layer build-up. This report presents the best-estimate thermal hydraulics results for the low- and high-density optimized silicide LEU core designs at 95 MW steady-state nominal operation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1983

A static model of chromospheric heating in solar flares

The response of the solar chromosphere to flare processes, namely nonthermal electrons, thermal conduction, and coronal pressure, is modeled. Finite difference methods employing linearization and iteration are used in obtaining simultaneous solutions to the equations of steady-state energy balance, hydrostatic equilibrium, radiative transfer, and atomic statistical equilibrium. The atmospheric response is assumed to be confined to one dimension by a strong vertical magnetic field. A solution is obtained to the radiative transfer equation for the most important optically thick transitions of hydrogen, magnesium, and calcium. The theoretical atmospheres discussed here are seen as elucidating the role of various physical processes in establishing the structure of flare chromospheres. At low coronal pressures, conduction is found to be more important than nonthermal electrons in establishing the position of the transition region. Only thermal conduction can adequately account for the chromospheric evaporation in compact flares. Of the mechanisms considered, only nonthermal electrons bring about significant heating below the flare transition region.

Ricchiazzi, P. J.↗
NASA NTRS1976

Results of phase change paint heat transfer tests utilizing 0.040 scale 50% forebody models (No. 82-0) of the Rockwell International space shuttle orbiter in AEDC VKF hypersonic tunnel B (test OH54A)

Results of aerodynamic heating tests conducted in October 1974 on a space shuttle orbiter model using the phase change paint technique are presented. The model was a 0.040 scale representation of the forward 50 percent of the orbiter. Surface roughness effects on boundary layer transition were investigated. Roughness was simulated by using steel balls varying in diameter from 0 (no balls) to 0.039 inch with 0.040 inch wide by 0.080 inch deep gaps. A nominal Mach number of 8 was tested with Reynolds number varying from 0.75 through 3.5 million per foot. Angle of attack was varied from 20 deg to 40 deg.

Dye, W. H.↗
DOE OSTI2021

Validation of theory-based models for the control of plasma currents in W7-X divertor plasmas

A theory-based model for the control of plasma currents for steady-state operation in W7-X is proposed and intended for model-based plasma control. The conceptual outline implies the strength of physics-based models: it offer approaches applicable to future conditions of fusion devices or next-step machines. The application at extrapolated settings is related to the validity range of the theory model. Therefore, the predictive power of theory-based control models could be larger than for data-driven approaches and limitations can be predicted from the validity range for the prediction of bootstrap currents in W7-X. The model predicts the L/R response when density or heating power is changed. The model is based on neoclassical bootstrap current calculations and validated for different discharge conditions. While the model was found to be broadly applicable for conducted electron-cyclotron-heated discharges in W7-X, limits were found for cases when the polarization of the electron cyclotron heating was changed from X2 to O2-heating. The validity assessment attempts to quantify the potential of the derived model for model-based control in the operational space (density, heating power) of W7-X.

neoclassical modelling↗
NASA NTRS1981

Finite element modeling of nonisothermal polymer flows

A finite element formulation designed to simulate polymer melt flows in which both conductive and convective heat transfer are important is described, and the numerical model is illustrated by means of computer experiments using extruder drag flow and entry flow as trial problems. Fluid incompressibility is enforced by a penalty treatment of the element pressures, and the thermal convective transport is modeled by conventional Galerkin and optimal upwind treatments.

Roylance, D.↗
NASA NTRS1996

Analytical Investigation of a Reflux Boiler

A thermal model of a single Ultralight Fabric Reflux Tube (UFRT) was constructed and tested against data for an array of such tubes tested in the NASA-JSC facility. Modifications to the single fin model were necessary to accommodate the change in radiation shape factors due to adjacent tubes. There was good agreement between the test data and data generated for the same cases by the thermal model. The thermal model was also used to generate single and linear array data for the lunar environment (the primary difference between the test and lunar data was due to lunar gravity). The model was also used to optimize the linear spacing of the reflux tubes in an array. The optimal spacing of the tubes was recommended to be about 5 tube diameters based on maximizing the heat transfer per unit mass. The model also showed that the thermal conductivity of the Nextel fabric was the major limitation to the heat transfer. This led to a suggestion that the feasibility of jacketing the Nextel fiber bundles with copper strands be investigated. This jacketing arrangement was estimated to be able to double the thermal conductivity of the fabric at a volume concentration of about 12-14%. Doubling the thermal conductivity of the fabric would double the amount of heat transferred at the same steam saturation temperature.

Simon, William E.↗
NASA NTRS1990

Lower solar chromosphere-corona transition region. I - Theoretical models with small temperature gradients

A study of transition region models including the effects of classical thermal conduction, heating, and radiative cooling is carried out with attention directed toward the problem of understanding the observed emission in the lower transition region. It is found that the observationally inferred emission measure curve implies a near-balance between heating and radiative cooling in the lower transition region, and that the presence of strong hydrogen Ly-alpha cooling leads to the existence of singularities in the solutions of the force balance and energy balance equations when such a near-balance between heating and cooling is assumed. These singularities place strong constraints on the nature of viable models of the lower transition region and must be considered when Ly-alpha cooling is important. Previously suggested explanations of the observed emission from the lower transition region are considered in the context of the results of the present study, and conditions for the applicability of these suggested explanations are discussed.

Woods, D. Tod↗
NASA NTRS1984

Atmospheric heating of meteorites: Results from nuclear track studies

A quantitative model to estimate the degree of annealing of nuclear tracks in mineral grains subjected to a variable temperature history was proposed. This model is applied to study the track annealing records in different meteorites resulting from their atmospheric heating. Scale lengths were measured of complete and partial track annealing, delta X sub 1 and delta X sub 2, respectively. In mineral grain close to fusion crust in about a dozen meteorites. Values of delta X sub 1 and delta X sub 2 depend on extent and duration of heating during atmospheric transit and hence on meteorite entry parameters. To estimate track annealing, the temperature history during atmospheric heating at different distances from the crusted surface of the meteorite is obtained by solving heat conduction equation in conjunction with meteorite entry model, and use of the annealing model to evaluate the degree of annealing of tracks. It is shown that the measured values of delta X sub 1 and delta X sub 2 in three of the meteorites studied are consistent with values using preatmospheric mass, entry velocity and entry angle of these meteorites.

Jha, R.↗
NASA NTRS2004

Boiling Heat Transfer Mechanisms in Earth and Low Gravity: Boundary Condition and Heater Aspect Ratio Effects

Boiling is a complex phenomenon where hydrodynamics, heat transfer, mass transfer, and interfacial phenomena are tightly interwoven. An understanding of boiling and critical heat flux in microgravity environments is of importance to space based hardware and processes such as heat exchange, cryogenic fuel storage and transportation, electronic cooling, and material processing due to the large amounts of heat that can be removed with relatively little increase in temperature. Although research in this area has been performed in the past four decades, the mechanisms by which heat is removed from surfaces in microgravity are still unclear. Recently, time and space resolved heat transfer data were obtained in both earth and low gravity environments using an array of microheaters varying in size between 100 microns to 700 microns. These heaters were operated in both constant temperature as well as constant heat flux mode. Heat transfer under nucleating bubbles in earth gravity were directly measured using a microheater array with 100 m resolution operated in constant temperature mode with low and high subcooled bulk liquid along with images from below and from the side. The individual bubble departure diameter and energy transfer were larger with low subcooling but the departure frequency increased at high subcooling, resulting in higher overall heat transfer. The bubble growth for both subcoolings was primarily due to energy transfer from the superheated liquid layer relatively little was due to wall heat transfer during the bubble growth process. Oscillating bubbles and sliding bubbles were also observed in highly subcooled boiling. Transient conduction and/or microconvection was the dominant heat transfer mechanism in the above cases. A transient conduction model was developed and compared with the experimental data with good agreement. Data was also obtained with the heater array operated in a constant heat flux mode and measuring the temperature distribution across the array during boiling. The instantaneous heat transfer into the substrate was numerically determined and subtracted from the supplied heat to obtain the wall to liquid heat flux.

Kim, Jungho↗
NASA NTRS1983

Constraints on rift thermal processes from heat flow and uplift

The implications of heat flow data available from five major Cenozoic continental rift systems for the processes of continental rifting are discussed, and simple thermal models of lithospheric thinning which predict uplift are used to further constrain the thermal processes in the lithosphere during rifting. Compilations of the heat flow data are summarized and the salient results of these compilations are briefly discussed. The uplift predictions of the slow and rapid thinning models, in which thinning is assumed to occur at a respectively slower and faster rate than heat can be conducted into the lithosphere, are presented. Comparison of uplift rates with model results indicates that the lithosphere is in a state between the two models. While uplift is predicted to continue after thinning has ceased due to thermal relaxation of the lithosphere, the rapid thinning model is always predicted to apply to surface heat flow, and an anomaly in this flow is not predicted to develop until after thinning has stopped.

Morgan, P.↗
DOE OSTI2024

Thermal Conductivity Degradation in High Burnup U-Pu-Zr Fuel

Recent advancements in the characterization of irradiated U-Pu-Zr fuels have revealed complexities that challenge existing understanding of constituent redistribution. Traditionally, models have proposed three concentric regions within the fuel, each characterized by distinctive phases and porosity. However, through detailed analysis of high burnup U-Pu-Zr, we discovered the presence of four distinct constituent redistribution regions. Particularly novel is the observation of significant Pu redistribution, a previously unreported phenomenon that necessitates a reevaluation of current models. This work aims to delve deeper into these findings, seeking to correlate mesoscale measurements of thermal diffusivity and respective thermal conductivity with the phases present in each redistribution region. To achieve this objective, we employed mesoscale thermoreflectance methods using the unique, Idaho National Laboratory (INL) developed, Thermal Conductivity Microscope (TCM) at INL’s Irradiated Materials Characterization Laboratory. The TCM employs two tightly focused lasers: one for heating to generate periodic thermal waves in the substrate, and another spatially separated probe laser to detect changes in the optical reflectivity of the gold-coated substrate resulting from thermal wave diffusion. We conducted several thermal diffusivity measurements within each region of constituent redistribution of a U-19Pu-10Zr fuel pin cross section irradiated to 11 at. % burnup. The TCM measurement positions strategically aligned with transmission electron microscopy (TEM) lift-out locations previously collected from the fuel sample. Complementary microstructural analysis techniques such as optical and scanning electron microscopy (OM/SEM), electron probe microanalysis for chemical compositions, and TEM-based selective area electron diffraction (SAED) analysis for crystallographic insights into each phase were also utilized. This comprehensive approach allowed us to correlate local thermal diffusivity data with microstructural characteristics, enabling the computation of local thermal conductivity at each position. The significance of this contribution lies in its pioneering use of the TCM for ternary fuel mesoscale examination, shedding light on the previously overlooked effects of Pu redistribution on local thermal conductivity. By informing current models capturing constituent redistribution and heat transfer, our findings pave the way for more accurate predictions of metallic fuel performance. Moreover, this work sets the stage for future comparisons with similar TCM examinations on U-19Pu-10Zr fuels at ultra-low burnup, facilitating a comprehensive understanding of thermal property changes across different burnup levels. Ultimately, our study not only enriches our understanding of the thermophysical properties of individual redistribution regions within U-Pu-Zr fuel but also offers valuable insights for the design and operational parameters of proposed next-generation fast reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

gLaSDI: Parametric physics-informed greedy latent space dynamics identification

A parametric adaptive physics-informed greedy Latent Space Dynamics Identification (gLaSDI) method is proposed for accurate, efficient, and robust data-driven reduced-order modeling of high-dimensional nonlinear dynamical systems. In the proposed gLaSDI framework, an autoencoder discovers intrinsic nonlinear latent representations of high-dimensional data, while dynamics identification (DI) models capture local latent-space dynamics. Here, an interactive training algorithm is adopted for the autoencoder and local DI models, which enables identification of simple latent-space dynamics and enhances accuracy and efficiency of data-driven reduced-order modeling. To maximize and accelerate the exploration of the parameter space for the optimal model performance, an adaptive greedy sampling algorithm integrated with a physics-informed residual-based error indicator and random-subset evaluation is introduced to search for the optimal training samples on the fly. Further, to exploit local latent-space dynamics captured by the local DI models for an improved modeling accuracy with a minimum number of local DI models in the parameter space, a -nearest neighbor convex interpolation scheme is employed. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including Burgers equations, nonlinear heat conduction, and radial advection. The proposed adaptive greedy sampling outperforms the conventional predefined uniform sampling in terms of accuracy. Compared with the high-fidelity models, gLaSDI achieves 17 to 2,658× speed-up with 1 to 5% relative errors.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Transient thermal response of ablating bodies

A numerical study of transient thermal response of a blunt-nosed axisymmetric body made of Teflon is presented using a two-layer thermal model. It is shown that phase change and transverse heat conduction have a considerable effect on the internal temperature field. Comparison of the numerical results with experimental data shows that the single-layer thermal model does not predict the real feature of the thermal field, whereas the results of the two-layer thermal model agree reasonably well with the experiment.

Arai, N.↗
DOE OSTI2022

Quantifying Thermal Output of Energetic Materials (LDRD Final Report)

We present the results of an LDRD project, funded by the Nuclear Deterrence IA, to develop capabilities for quantitative assessment of pyrotechnic thermal output. The thermal battery igniter is used as an exemplar system. Experimental methodologies for thermal output evaluation are demonstrated here, which can help designers and engineers better specify pyrotechnic components , provide thermal output guidelines for new formulations, and generate new metrics for assessing component performance and margin given a known failure condition. A heat-transfer analysis confirms that the dominant mode of energy transfer from the pyrotechnic output plume to the heat pellet is conduction via deposition of hot titanium particles. A simple lumped-parameter model of titanium particle heat transfer and a detailed multi-phase model of deposition heat transfer are discussed. Pyrotechnic function, as defined by "go/no-go" standoff testing of a heat pellet, is correlated with experimentally measured igniter plume temperature, titanium metal particle temperature, and energy deposition. Three high-speed thermal diagnostics were developed for this task. A three-color imaging pyrometer, acquiring 100k images per second on three color channels, is deployed for measurement of titanium particle temperatures. Complimentary measurements of the overall igniter plume emission ("color") temperature were conducted using a transmission-grating spectrograph in line-imaging mode. Heat flux and energy deposition to a cold wall at the heat-pellet location were estimated using an eroding thermocouple probe, with a frequency response of ~5 kHz. Ultimate "go/no-go" function in the igniter/heat-pellet system was correlated with quantitative thermal metrics, in particular surface energy deposition and plume color temperature. Titanium metal-particle and plume color temperatures both experience an upper bound approximated by the 3245-K boiling point of TiO 2 . Average metal-particle temperatures remained nearly constant for all standoff distances at T = 2850 K, ± 300 K, while plume color temperature and heat flux decay with standoff—suggesting that heat-pellet failure results from a drop in metal-particle flux and not particle temperature. At 50% likelihood of heat-pellet failure, peak time-resolved plume color temperatures drop well below TiO 2 boiling to ~2000 - 2200 K, near the TiO 2 melting point. Estimates of peak heat flux decline from up to 1 GW/m 2 for near-field standoffs to below 320 MW/m 2 at 50% failure likelihood.

36 MATERIALS SCIENCE↗