Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Conduction”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

From Serendipity to Exploration: Synthesis and Properties of PtSiSb

PtSiSb is a compound with covalent Si–Sb bonds, which are relatively rare instances in solid-state materials. The original report on PtSiSb indicates synthetic challenges for the isolation of single-phase samples, as well as potential discrepancies in the predicted electronic structure and experimental transport properties. In this work, we considered PtSiSb as a case study with the objective of demonstrating the importance of synthesis as a bridge between theoretical predictions and experimental properties. Using guidance from in situ powder X-ray diffraction, the synthesis of single-phase polycrystalline samples was developed, allowing further experimental characterization of transport and catalytic properties. Thermal conductivity and heat capacity demonstrated anomalous behavior, indicating an unconventional phonon scattering mechanism in PtSiSb. Charge transport measurements indicated semimetallic behavior. Hydrogen evolution reaction (HER) tests in acidic electrolyte revealed that the PtSiSb cathode needs an overpotential η 10 = 120 mV to achieve a current density j = −10 mA/cm 2 . This is lower activity than elemental Pt but comparable to known and active transition-metal phosphide HER catalysts, such as CoP, with η 10 = 93 mV when measured at identical conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Stability Criteria for Self-Propagating Reaction Waves in Co/Al Multilayers

The propagation of self-sustained formation reactions in sputter-deposited Co/Al multilayers is known to exhibit a design-dependent instability. Multilayers having thin bilayers (<55 nm period) exhibit stable propagating waves, whereas those with a larger period react unstably. The specific two-dimensional (2D) instability observed involves the transverse propagation of a band in front of a stalled front commonly referred to as a “spin band.” Previous finite-element studies have shown that these instabilities are thermodynamically driven by the forward conduction of heat away from the flame front. However, the magnitude of that loss is inherently tied to the bilayer design in traditional bimetallic multilayers, which couples any proposed stability criteria to a varying critical diffusion distance. Here, this work utilizes a recently developed class of materials known as “inert-mediated reactive multilayers” to decouple the thermodynamic and kinetic contributions to propagating wave stability by reducing the stored chemical energy density in normally stable bilayer designs. By depositing an inert product phase (B2-CoAl) within the mid-plane of Co and Al reactant layers, spin instabilities arise as a function of both diluted volume and critical diffusion distance. From there, a stability criterion is determined for Co/Al multilayers based on enthalpy loss from the reaction zone, and its physical significance is explored.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Agreement of Analytical and Simulation-Based Estimates of the Required Land Depth in Climate Models

Previous analytical and simulation-based analyses suggest that deeper land surface models are needed to realistically simulate the terrestrial thermal state in climate models, with implications for land-atmosphere interactions. Analytical approaches mainly focused on the subsurface propagation of harmonics such as the annual temperature signal, and a direct comparison with climate-change model output has been elusive. Here, this study addresses the propagation of a harmonic pulse fitted to represent the timescale and amplitude of anthropogenic warming. Its comparison to land model simulations with stepwise increased bottom boundary depth leads to an agreement between the simulation-based and analytical frameworks for long-term climate trends. Any depth increase gradually decreases the relative error in the subsurface thermodynamics, and a minimum depth of 170 m is recommended to simulate the ground climate adequately. The approach provides an accurate estimate of the required land-model depth for climate-change simulations and assesses the relative bias in insufficiently deep land models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Implementation and Evaluation of a Unified Turbulence Parameterization Throughout the Canopy and Roughness Sublayer in Noah‐MP Snow Simulations

Abstract The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

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

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Boundary layer description of directional polymer crystallisation

Nearly fifty years ago Lovinger and Gryte suggested that the directional crystallization of a polymer was analogous to the quiescent isothermal crystallization experiment but at a supercooling where the crystal growth velocity was equal to the velocity of the moving front. Our experiments showed that this equivalence holds in a detailed manner at low directional velocities. To understand the underlying physics of these situations, we modeled the motion of a crystallization front in a liquid where the left side boundary is suddenly lowered below the melting point (Stefan's problem) but with the modification that the crystallization kinetics follow a version of the Avrami model. Our numerical results surprisingly showed that the results of the polymer analog track with the Stefan results which were derived for a simple liquid that crystallizes completely at its melting point; in particular, the position of the crystal growth-front evolved with time exactly as in the Stefan problem. The numerical solution also showed that the temperature in the immediate vicinity of the growth-front decreased with increasing front velocity, which is in line with Lovinger and Gryte's ansatz. To provide a clear theoretical understanding of these numerical results we derive a boundary layer solution to the governing coupled differential equations of the polymer problem. The analytical results are in agreement with our observations from experiments and numerical computations but show that this equivalence between the small molecule and polymer analog only holds in the limit where the crystallization enthalpy is much larger than the rate at which heat is conducted away in the polymer. In particular, in the context of the temperature profile, the enthalpy generated by the crystallisation process which is spread out over a narrow spatial region can be approximated as a point source whose location and temperature correspond to the Lovinger–Gryte ansatz.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Thermal Design Consideration of Medium Voltage High Frequency Transformers

This paper describes a thermal design using epoxy potting for high frequency transformer (HFT) for medium voltage high power electronic applications, namely emerging solid-state transformers. To address HFT’s thermal stress problem, thermal performance of the HFT with air and epoxy potting cooling method are both simulated by magnetothermal-fluid coupling method (MTFCM). The effectiveness of the simulation analysis is verified in a 100 kW DC-DC resonant converter. Influencing factors including distance between the high frequency transformer and the metal enclosure, epoxy thermal conductivity and heat sink placements are simulated and analyzed. Finally, an optimal cooling structure is proposed for future high power HFTs.

Wang, Haoming↗
DOE OSTI2023

Dating submarine landslides using the transient response of gas hydrate stability

Abstract Submarine landslides are prevalent on the modern-day seafloor, yet an elusive problem is constraining the timing of past slope failure. We present a novel age-dating technique based on perturbations to underlying gas hydrate stability caused by slide-impacted seafloor changes. Using three-dimensional (3-D) seismic data, we mapped an irregular bottom simulating reflection (BSR) underneath a submarine landslide in the Orca Basin, Gulf of Mexico. The irregular BSR mimics the pre-slide seafloor geometry rather than the modern bathymetry. Therefore, we suggest that the gas hydrate stability zone (GHSZ) is still adjusting to the post-slide sediment temperature. We applied transient conductive heat-flow modeling to constrain the response of the GHSZ to the slope failure, which yielded a most likely age of ca. 8 ka, demonstrating that gas hydrate can respond to landslides even on multimillennial time scales. We further provide a generalized analytical solution that can be used to remotely date submarine slides in the absence of traditional dating techniques.

Geology↗
DOE OSTI2025

Moltensaltpropnet

MoltenSaltPropnet is a physics-informed machine learning framework that aims to predict the thermophysical properties of molten fluoride and chloride salt mixtures, which are crucial for the design and safety of Generation IV molten salt reactors. The code processes data from the Molten-Salt Thermal Properties Database (MSTDB-TP) and the Janz compendium, converting critically evaluated correlations into fast, differentiable surrogate models for density, viscosity, thermal conductivity, and heat capacity across 448 distinct salt systems. The implementation consists of several key components: 1. Data Curation: The code parses and cleans the raw data, normalizing elemental mole fractions and extracting relevant regression coefficients for various thermophysical properties. 2. Feature Engineering: It generates fixed-length numerical descriptors that encapsulate the composition and temperature, incorporating polynomial interaction terms and dimensionality-reduction techniques to optimize model performance. 3. Coefficient Learning: Four different machine learning architectures are employed: a deep residual network (ResNet), a Kolmogorov–Arnold network (KAN), a sparsity-inducing neural network (SNN), and classical regression models. Each model learns to predict coefficients that define the temperature-dependent correlations for the thermophysical properties. 4. Property Reconstruction: The predicted coefficients are used to compute temperature-dependent property values, ensuring positivity and monotonic trends through a composite loss function that enforces physical constraints. 5. User Interface: An open-source web application enables users to filter the database, train task-specific models, and visualize the results, allowing for rapid exploration of candidate salt mixtures. MoltenSaltPropnet bridges the gap between limited experimental data and high-fidelity reactor simulations, providing a powerful tool for researchers in the field of molten salt reactors and advanced nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2023

Low-order aquifer thermal energy storage model for geothermal system simulation

This paper presents a low-order aquifer thermal energy storage (ATES) model for simulation of combined subsurface and above-surface energy systems. The model is included in the Modelica IBPSA Library, which is a free open-source library with basic models for building and district energy and control systems. The model uses a lumped-component method, in which the transient conductive-convective heat and mass transfer equation is radially discretized. To verify the accuracy of the model, we present an inter-model comparison from a simulation test suite. Results show that the Modelica ATES model is in good agreement, with a normalized mean bias error for yearly variation of aquifer temperatures of 1.6×10−2 and 9×10−5 at 1 m and 10 m distance from the well.

Maccarini, Alessandro↗
DOE OSTI2020

Thermophysical Property Measurements on Salt Mixture

Thermal conductivity, specific heat, viscosity, and density are critical properties for understanding thermal hydraulics of heat transfer fluids (HTFs) like molten salts. These properties are important for design and safe operation of any HTF loop. These characterizations are required to benchmark computational chemical models that are used to predict properties, but existent data are limited. These measurements are extremely difficult to make due to the operating temperature range and corrosive properties of the salt.Density measurements are needed to predict volumetric behavior. Thermal expansion mismatches can cause mechanical failure of structural materials and important components that are in contact with the molten salt. All of these properties are related to heat transfer in HTFs. The heat generated in a reactor must be high enough to derive power; however it should not be so hot that it cannot be safely removed via the HTF at an appropriate rate. It is important to understand the thermal hydraulics required to operate at the desired power generation level and to avoid situations in which thermal shifts can damage components,cause failure of structural materials, and/or lead to accidents. Understanding these thermal properties can aid in component and process design and can contribute to optimized operational and safety parameters.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Benefit FOA FY2015 - Solid State Magnetocaloric Air Conditioner (Final Report)

The active magnetic refrigeration (AMR) cycle holds great promise to reduce energy consumption of space cooling and heating. It is estimated that it could achieve up to 25% improvement in energy efficiency over conventional vapor compression systems. Prototypes have been built and demonstrated the working principles and energy benefits of the technology. These prototypes however were of low specific cooling power, defined as cooling power obtainable from unit mass of magnetocaloric materials (MCM). Specific cooling power is inversely proportional to the cost of the system. Therefore, for the magnetocaloric technology to be commercially viable, the specific cooling power must be increased. To increase the specific cooling power, for any MCM, the frequency of operation must increase. However, the maximum frequency that an AMR can operate at is limited by the residence time that is required to transfer heat in and out of the regenerator during the cold and got blows. This time length is dependent on the thermal conductivity of the heat transfer medium. The higher the thermal conductivity of the heat transfer medium, the faster heat can be transferred between the medium and the MCM.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Multi-fidelity thermal modeling of laser powder bed additive manufacturing

Laser powder bed fusion (LPBF) Additive manufacturing (AM) has attracted interest as an agile method of building production metal parts to reduce design-build-test cycle times for systems. However, predicting part performance is difficult due to inherent process variabilities. This makes qualification challenging. Computational process models have attempted to address some of these challenges, including mesoscale, full physics models and reduced fidelity conduction models. The goal of this work is credible multi-fidelity modeling of the LPBF process by investigating methods for estimating the error between models of two different fidelities. Two methods of error estimation are investigated, adjoint-based error estimation and Bayesian calibration. Adjoint-based error estimation is found to effectively bounding the error between the two models, but with very conservative bounds, making predictions highly uncertain. Bayesian parameter calibration applied to conduction model heat source parameters is found to effectively bound the observed error between the models for melt pool morphology quantities of interest. However, the calibrations do not effectively bound the error in heat distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Boson Fermion Nucleus Polynuclear Structure: Monograph #11

The boson fermion nucleus (BFN) structure indicates the possibility of stable, nonradioactive polynuclear structures bound by the strong nuclear force. In the process of making polynuclear structures, BFN self-ordering on the nuclear level as a stable structure resembles the self-ordering of carbon atoms at the atomic level into buckminsterfullerene, graphene sheet, or graphene nanotubes. This natural assembly derives from the modulation of electron orbital collapse, which decreases coulomb repulsion between atoms and allows a natural reordering among nuclear structures by modulating the permeability of free space. Polynuclear material would be unearthly in tensile strength, thermal and electrical conductivity, and heat of disassociation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

An Overview of the Molten Salt Thermal Properties Database–Thermophysical, Version 2.1.1 (MSTDB-TP v.2.1.1)

The current status of the Molten Salt Thermodynamic Database- Thermal Physical (MSTDB-TP) is reported. Building off of MSTDB-TP and through the direction of the Roadmap for thermal property measurement of Molten Salt Reactor systems, MSTDB-TP 2.1.1 has now been released containing a total 448 salt entries (data from 140+independent studies) compared to the original commit of 62, containing thermophysical properties including density, viscosity, thermal conductivity, and heat capacity. Along with the new database release, further advancement of Saline, a C++application programming interface (API), and a graphical user interface (GUI) has facilitated increased user/developer interaction with the database. Furthermore, estimation techniques, first principles calculations (Ab-Initio) and interpolation/extrapolation methods (Redlich-Kister/Muggianu), have shown great promise in the future of thermophysical property determination for filling out compositional spaces and investigating experimentally difficult salts (hazardous/expensive). This report describes the database composition, development and advancement of database tools, and the strategy of advancing and implementing estimation data into future iterations of the database for MSTDB-TP.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Thermogravimetric Analysis (TGA) for Carbon Fiber and Glass Fiber Epoxy Composites and their Constituents

In this work, thermogravimetric analysis (TGA) was performed on samples of a carbon fiber epoxy composite, a glass fiber epoxy composite, and a mixed carbon fiber/glass fiber epoxy composite, as well on each constituent material (polymer epoxy, carbon fibers and glass fibers). TGA was conducted for heating rates from 1-20 C/min with purified purge gases of nitrogen and dry air. For the fiberglass composite, we find that ~70% of the material remains after heating in air to 1200 C. For the carbon fiber epoxy composite, we observe greater mass loss as the carbon fibers can oxidize, leaving little material by the end of the test. The mixed composite, which has a 2:1 ratio of glass fibers to carbon fibers, experienced a total mass loss between the two other composites. By determining the relationship between the thermal decomposition of a composite material and its constituent materials, we can predict the fire behavior of novel composites during the material design phase.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Refining NEAMS MARVEL Reactor Model Accuracy

This report outlines the progress of Idaho National Laboratory (INL) in improving the fidelity of the multiphysics model of the Microreactor Applications Research Validation and Evaluation (MARVEL) reactor developed in FY24. The multiphysics model firstly developed in FY24 leverages three single-physics models coupled via the MOOSE’s MultiApps and Transfer systems. The first single-physics model, that functions as the driver application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). The second single-physics model uses BISON to compute the solid temperature in the reactor. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. In this report we detail the following improvement and enhancements: • An updated mesh leveraging the latest feature in the reactor module (see Section 3.1.).

42 - ENGINEERING↗