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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8

Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental characterization and potential energy savings of insulated cladding for U.S. residential buildings

This study evaluates a novel insulated cladding composed of foamed cement and vinyl siding that is designed as a retrofit product to enhance thermal performance in U.S. residential buildings. Experimental measurements showed that foamed cement exhibits a thermal conductivity of 0.0334–0.0365 W/m-K (R-value of 4.0–4.3 °F⋅ft 2 ⋅h/BTU), which is significantly higher than that of conventional fiber cement boards. Although mechanical testing confirmed the material’s brittleness and low strength, it also indicated its suitability for nonstructural insulation applications. Water vapor permeability testing demonstrated effective moisture resistance with the polymer coating. To evaluate the insulated cladding’s energy performance, we used the U.S. Department of Energy’s prototype single-family building model and simulated across 50 locations spanning 16 International Energy Conservation Code climate zones. The R-7 (R-value of 7 °F⋅ft 2 ⋅h/BTU per inch) retrofit cladding resulted in heating energy savings up to 275 therms (Fairbanks, AK) and cooling energy savings up to 1,547 kWh (Phoenix, AZ). The energy cost savings varied by region, but the highest annual cost savings ($\$$434/year) were observed in Santa Maria, CA, and the highest annual cost savings percentage (37 %) were noted in Monterey, CA. On a national scale, the insulated retrofit provided an average annual energy cost reduction of 16 %, highlighting the cladding system’s broad applicability and financial viability.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Thermal properties of U-Mo alloys irradiated under high fission power density

A variety of physical and thermal property measurements made on irradiated U-Mo alloy monolithic fuel samples with a Zr diffusion barrier and clad in Al alloy 6061 are reported. These properties form the basis for additional investigations in thermal behavior of U-Mo fuel under high performance research reactor irradiation conditions. Measurements were performed on samples harvested from fuel plates irradiated under high fission power density and surface heat flux and complement previous measurements on samples subjected to moderate fission power density and surface heat flux. Specific heat capacity of irradiated monolithic U-Mo increases with increasing temperature and fission density, but carries large uncertainties associated with thickness measurements used to determine the mass fraction of the layer constituents for extraction of the U-Mo specific heat capacity, specifically the very thin and nonuniform Zr diffusion barrier. Density of irradiated monolithic U-Mo decreases with increasing fission density resulting from fission gas swelling and porosity formation. Thermal diffusivity of irradiated monolithic U-Mo increased with increasing temperature and decreased with increasing fission density. Thermal conductivity of irradiated monolithic U-Mo increased with increasing temperature, and in general, appeared to decrease and become less sensitive to temperature with increasing fission density. Thermal conductivity measurements compared to a semi-empirical model agreed within ±30%. It is suspected that the main driver in deviation between model calculations and measurements is the fission gas swelling determination, which does not account for fission power density or irradiation temperature.

uranium (U), molybdenum, nuclear fuel, thermal con↗
DOE OSTI2021

Deep-learning-enabled Bayesian inference of fuel magnetization in magnetized liner inertial fusion

Fuel magnetization in magneto-inertial fusion (MIF) experiments improves charged burn product confinement, reducing requirements on fuel areal density and pressure to achieve self-heating. By elongating the path length of 1.01 MeV tritons produced in a pure deuterium fusion plasma, magnetization enhances the probability for deuterium–tritium reactions producing 11.8-17.1 MeV neutrons. Nuclear diagnostics thus enable a sensitive probe of magnetization. Characterization of magnetization, including uncertainty quantification, is crucial for understanding the physics governing target performance in MIF platforms, such as magnetized liner inertial fusion (MagLIF) experiments conducted at Sandia National Laboratories, Z-facility. We demonstrate a deep-learned surrogate of a physics-based model of nuclear measurements. A single model evaluation is reduced from O ( 10 – 100 ) (10–100) CPU hours on a high-performance computing cluster down to O ( 10 – 100 ) (10) ms on a laptop. This enables a Bayesian inference of magnetization, rigorously accounting for uncertainties from surrogate modeling and noisy nuclear measurements. The approach is validated by testing on synthetic data and comparing with a previous study. We analyze a series of MagLIF experiments systematically varying preheat, resulting in the first ever systematic experimental study of magnetic confinement properties of the fuel plasma as a function of fundamental inputs on any neutron-producing MIF platform. We demonstrate that magnetization decreases from B R ∼ 0.2 ~0.5 to B R ∼ 0.2 ~0.2 MG cm as laser preheat energy deposited increases from E preheat ∼ 460 ~460 J to E preheat ∼ 460 ~1.4 kJ. This trend is consistent with 2D LASNEX simulations showing Nernst advection of the magnetic field out of the hot fuel and diffusion into the target liner.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fiberboard Properties and Degradation Rates for Storage of 9975 Shipping Packages in KAC (2021 Status Report)

Thermal, mechanical, and physical properties have been measured on fiberboard samples following accelerated aging for up to 16 years. The aging environments include elevated temperature up to 250 ºF (the maximum allowed service temperature for fiberboard in 9975 packages) and elevated humidity. Accelerated aging results have been analyzed and used to build aging models. Correlations relating several properties (thermal conductivity, energy absorption, weight, dimensions, and density) to their rate of change in potential storage environments have been developed. Combined with an estimate of the actual conditions the fiberboard experiences in K-Area Complex (KAC), these models allow development of service life predictions. KAC completed calculations that supported extending the service life of 9975 packages with 3013 and non3013 containers in storage from 15 years to 20 years. These calculations addressed the potential for degraded fiberboard properties following 20 years in storage using degradation models developed previously. The updated degradation models developed in this report have been compared to the assumptions in these calculations, and the current results for fiberboard dimensions, density, radial thermal conductivity, and energy absorption remain consistent with those calculations. The axial thermal conductivity degradation rate slightly exceeds calculations estimate of 0.5% per year between 130 and 160 °F. Additional data will continue to be collected to permit future refinements to the models and assumptions. Some of the predicted degradation rates presented in this report are extreme. However, these relate to environments that do not exist within KAC or would be postulated only as upset conditions that would not likely persist for an extended period. For a typical package stored in KAC with approximately 10-12 watts internal heat load or less, and ambient temperatures below 90 ºF, the fiberboard assemblies experience storage conditions less severe than any of the aging environments. Fiberboard in conforming packages with lower internal heat loads should experience little or no degradation and is expected to provide a service life beyond the currently approved 20 year storage period. Packages with higher internal heat loads may not continue to perform their required safety functions beyond 20 years. Ultimately, service life will be determined by the cumulative effect of degradation from all the conditions these packages might encounter. The results and model predictions presented in this report are applicable to 9975 packages with cane or softwood fiberboard overpack assemblies. These degradation models do not address the effects of nonconforming conditions such as the presence of excess moisture and mold or beetle infestations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Fiberboard Properties and Degradation Rates for Storage of 9975 Shipping Package in KAC (2022 Status Report)

Thermal, mechanical, and physical properties have been measured on fiberboard samples following accelerated aging for up to 17 years. The aging environments include elevated temperature up to 250 °F (the maximum allowed service temperature for fiberboard in 9975 packages) and elevated humidity. Accelerated aging results have been analyzed and used to build aging models. Correlations relating several properties (thermal conductivity, energy absorption, weight, dimensions, and density) to their rate of change in potential storage environments have been developed. Combined with an estimate of the actual conditions the fiberboard experiences in K-Area Complex (KAC), these models allow development of service life predictions. KAC completed calculations that supported extending the service life of 9975 packages with 3013 containers in storage from 20 years to 40 years. These calculations addressed the potential for degraded fiberboard properties following 40 years in storage using degradation models developed previously. Due to differences in containers, 9975 packages with non-3013 containers are only approved for 20 years of storage. The updated degradation models developed in this report have been compared to the assumptions in those calculations, and the current results for fiberboard dimensions, density, radial thermal conductivity, and energy absorption remain consistent with those calculations. The axial thermal conductivity degradation rate slightly exceeds calculations estimate of 0.5% per year between 130 and 160 °F. Additional data will continue to be collected to permit future refinements to the models and assumptions. Some of the predicted degradation rates presented in this report are extreme. However, these relate to environments that do not exist within KAC or would be postulated only as upset conditions that would not likely persist for an extended period. For a typical package stored in KAC with approximately 10-12 watts internal heat load or less, and ambient temperatures below 90 °F, the fiberboard assemblies experience storage conditions less severe than any of the aging environments. Fiberboard in conforming packages with lower internal heat loads should experience little or no degradation and is expected to provide a service life beyond the currently approved 40 year storage period. Packages with higher internal heat loads may not continue to perform their required safety functions beyond 40 years. Ultimately, service life will be determined by the cumulative effect of degradation from all the conditions these packages might encounter. The results and model predictions presented in this report are applicable to 9975 packages with cane or softwood fiberboard overpack assemblies. These degradation models do not address the effects of non-conforming conditions such as the presence of excess moisture and mold or beetle infestations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigation of Point-Contact Strategies for CFD Simulations of Pebble-Bed Reactor Cores

This study numerically investigated the effects of various contact strategies on the thermal hydraulic behavior within a structured bed of 100 explicitly modeled pebbles. Four contact strategies and two thermal hydraulic conditions were considered. The strategies to avoid contact singularities include decreasing the pebble diameter, increasing the pebble diameter, bridging the pebble surfaces near the contact region, and capping the pebble surfaces near the contact region. One strategy, Strategy 3a, which involves bridging with a cylinder equal to 10% of the pebble diameter, was selected as the baseline strategy because it addressed the contact singularity while minimizing the geometric changes that affect the bed porosity. The two thermal hydraulic conditions were full-power operation (Case 1) and pressurized loss of forced cooling or PLOFC (Case 2). Simulations of the conjugate heat transfer within the structured bed were performed using the Reynolds-averaged Navier–Stokes approach with the realizable k-ϵ turbulence model and two-layer all y + wall treatment. The thermal-fluid quantities of interest were compared between the contact strategies for each case. In Case 1, the hydraulic behavior was sensitive to the contact strategy, with large differences in the pressure drop (30%) and volume-average velocity (4%). The thermal behavior was not sensitive, with less than a 0.5% difference across the strategies. To better understand the separate effects of each heat transfer mode, Case 2 was divided into the following subcases: conduction (Case 2a); conduction/radiation (Case 2b); and conduction/radiation/convection (Case 2c). Case 2a represents an early phase of the PLOFC transient. Case 2b represents an intermediate phase of the PLOFC transient, with the pebble temperatures sufficiently high for the radiative heat transfer to be non-negligible. Case 2c represents a late phase of the PLOFC transient after the establishment of the natural circulation of the heat transfer fluid. For Case 2, large differences in the contact strategy were observed only in Case 2a with only conduction. The difference in the maximum pebble temperature was 23% in Case 2a, 2% in Case 2b, and 0.3% in Case 2c.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗
DOE OSTI2021

Energy-saving potential evaluation for primary schools with occupant-centric controls

Recent studies demonstrated that there is significant energy-saving potential for primary schools, which heating, ventilation, and air-conditioning (HVAC) systems with occupant-centric control (OCC) is an excellent candidate to save energy. However, such an energy impact has yet to be evaluated for different climate zones as well as different energy code versions, but is critical for technology transfer and deployment. Therefore, this paper conducts comprehensive evaluation on the energy-saving potentials for the primary schools with two advanced OCC strategies: presence-based and counting-based. In this work, ninety-six building energy models with stochastic behavior of occupants are developed and simulated, which consist of two building energy code versions, 16 climate zones, and baseline case (without OCC) and two advanced cases (with OCC). To evaluate the energy saving potentials for OCC, primary schools in the U.S. are used as an example. The results show that there is significant energy-saving potential for primary schools by considering OCC strategies, especially the counting-based case. The energy-saving potential is up to 10.2% for presence-based OCC and 12.41% for counting-based OCC. Furthermore, both climate and code version have a significant impact on energy savings from OCC strategies. The energy-saving potentials vary from 1.79% to 12.41% for different climates and code versions. This evaluation can also contribute to quantify the nationwide energy saving potential for other countries.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Evaluation of an Improved Convective Triggering Function: Observational Evidence and SCM Tests

This study provides a strong observational support for a recently developed convective triggering function that uses the large-scale dynamic convective available potential energy (dCAPE) as a constraint combined with an unrestricted air parcel launch level (ULL) to relax the unrealistic strong coupling of convection to surface heating and capture nocturnal elevated convection. Both case study and statistical analysis are conducted using the observations collected from the Department of Energy (DOE)'s Atmospheric Radiation Measurement (ARM) program at its Southern Great Plains (SGP) and Manaus (MAO) sites. They show that dCAPE has a much stronger correlation with precipitation than CAPE and ULL is essential to detect elevated convection above boundary layer under both midlatitude and tropical conditions and for both afternoon and nighttime deep convection regimes. Sensitivity tests with the single-column model (SCM) of DOE's Energy Exascale and Earth System Model (E3SM) indicate that the role of dCAPE in suppressing daytime convection is more effective for tropical convection than midlatitude convection. Even though the dCAPE can suppress the overestimated convection, ULL plays a much bigger role in improving the diurnal cycle of precipitation than dCAPE. It not only helps capture nocturnal elevated convection, but also significantly removes the spurious morning precipitation seen in the default model, due to the release of unstable energy at night. However, the use of ULL has led to an overestimation of light-to-moderate precipitation (1-10 mm/day) due to more convection being triggered above boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Enhancing Turbulent Mixing and Microphysical Uniformity in a Tall Convection‐Cloud Chamber Through Idealized Heterogeneity of Boundaries

A large convection cloud chamber has been proposed for exploring aerosol–cloud–drizzle interactions under well‐controlled turbulent conditions. Recent theoretical and numerical studies suggest that a convection cloud chamber with two heated and two cooled sidewalls can significantly enhance the liquid water content and thus benefit drizzle initiation. However, a chamber with such a sidewall configuration develops stable stratification and extremely weak turbulence therein. In this study, we conduct large‐eddy simulations of a tall convection chamber with five different sidewall configurations consisting of alternating warm and cold patches. For each configuration, the total surface area of warm patches equals that of cold patches, resulting in the same expected cloud‐free supersaturation based on a flux budget model. Results show that changing the sidewall configuration, while keeping all other factors constant, can substantially enhance turbulent mixing and improve the uniformity of thermodynamic and cloud microphysical properties in the bulk region of the chamber. In addition, turbulence strength is positively correlated with liquid water content and negatively correlated with cloud droplet number concentration, consistent with theoretical predictions. Our results highlight the advantage of building a large cloud chamber using modular patches with individually controllable temperature and humidity to achieve well‐mixed conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Divertor Plasma Detachment Control Neural Network

DivControlNN is a state-of-the-art software tool that leverages advanced machine learning techniques to predict and control divertor plasma behavior in fusion reactors. Plasma, a highly energetic and electrically charged gas, requires meticulous management to protect reactor components and maintain optimal energy production. Conventional simulation methods, although extremely detailed, typically demand extensive computational time-making them unsuitable for real-time control scenarios. DivControlNN addresses this challenge by learning from tens of thousands of high-fidelity simulations, thereby creating a rapid surrogate model that can deliver near-instantaneous predictions. At the core of its functionality is a sophisticated technique known as latent space mapping, which condenses complex, high-dimensional plasma data into a compact, lower-dimensional representation. This streamlined representation enables the system to quickly forecast essential plasma properties and determine the precise conditions required for effective detachment. Detachment is a crucial process in which the plasma is cooled before reaching the divertor plates, thereby reducing heat loads and mitigating material erosion. In recent experiments conducted on the KSTAR tokamak in South Korea, DivControlNN successfully guided the detachment process without any fine-tuning-even when applied to a new tungsten divertor configuration. By achieving a computational speed-up of over one hundred million times compared to traditional simulation methods while maintaining low prediction errors, DivControlNN stands to significantly enhance real-time control and diagnostic capabilities in future fusion reactors. This breakthrough paves the way for safer, more reliable reactor operation and represents a major advancement toward realizing fusion energy as a practical, sustainable, and clean power source.

Xu, Xueqiao [Lawrence Livermore National Laborator↗
DOE OSTI2024

Design, Fabrication and Testing of Surveillance Test Articles for MSR Materials Degradation Management

This report details the design, fabrication, and testing of surveillance test articles aimed at assessing material damage in reactor-relevant environments for effective degradation management. Two types of surveillance test articles with reduced sizes were developed based on design algorithms and finite element modeling: welded design and interlocking design. A furnace heating setup was adopted to apply multiple thermal cyclic loading profiles on the test articles with a temperature range of 500°C - 700°C, while the strain response was monitored using a digital image correlation technique. The testing results demonstrated the successful capturing of expected strain range for welded design while machining tolerance should be improved to engage strain coupling in the interlocking design. Mid-term (500 hours) and long-term (1500 hours) cyclic tests were conducted on welded test articles. A constant strain range of ~0.6% was observed at the specimen with testing under 500 hours, while a gradual decrease of strain at specimen was observed after 500 hours. Non-destructive evaluation through X-ray computed tomography confirmed the microcracks in the welds at specimen-driver joints after cyclic test that caused the strain change. Creep testing of the specimen after long-term cyclic test revealed a short creep life than expected. A multi-profile cyclic test was also conducted on a test article and demonstrated consistent strain response under different temperature ramp rates. The report also briefly discussed the challenges and future research efforts to advance test article development for material surveillance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Estimating soil thermal inertia profiles from the passive equilibration of a temperature probe: Supporting data and code

This dataset contains measurements of soil thermal properties and organic layer thickness at three field sites on the Seward Peninsula near Nome, Alaska, USA. Code is also provided for estimation of soil thermal properties from temperature data, including a forward heat transfer mode, a library of simulated temperature data, and an inversion scheme to retrieve thermal properties from temperature measurements (Lamb et al., 2025). This dataset and code is published in support of the manuscript “Estimating Soil Thermal Inertia Profiles from the Passive Equilibration of a Temperature Probe” by Lamb et al., 2025, in which an experimental thermal measurements technique is developed. The dataset contains measurements of soil thermal properties and organic layer thickness obtained with an industry standard thermal properties analyzer and direct observations, as well as data generated using the experimental technique to estimate thermal properties and soil structure from temperature observations. All data is reported in .csv format, and all code is written in Python 3.7 using standard packages. This research was performed as a part of the NGEE Arctic project, which aims to advance model predictions of arctic carbon cycle responses to a changing climate over the 21st Century.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

RELAP5-3D validation studies based on the High Temperature Test facility

In the spring and summer of 2019, experiments were conducted at the High Temperature Test Facility (HTTF) that form the basis of an upcoming high-temperature gas-cooled reactor (HTGR) thermal hydraulics (T/H) benchmark. HTTF is an integral effects test facility for HTGR T/H modeling validation. This paper presents RELAP5-3D models of two of those experiments: PG-27, a pressurized conduction cooldown (PCC); and PG-29, a depressurized conduction cooldown (DCC). These models used the RELAP5-3D model of HTTF originally developed by Paul Bayless as a starting point. The sensitivity analysis and uncertainty quantification code, RAVEN was used to perform calibration studies for the steady-state portion of PG-27. Here we developed four PG-27 calibrations based on steady-state conditions. These calibrations all used an effective thermal conductivity equal to 36 % of the measured thermal conductivity, but they differed with respect to the frictional pressure drops and radial conduction models. These models all captured the trends in steady-state temperature distributions and transient temperature behavior well. All four calibrations show room for improvement in predicting the transient temperature rise. The smallest error in temperature rise during the transient was a 21 % underprediction, and the largest was a 48 % underprediction. The errors in transient temperature rise are largely a result of a mismatch in power density between the RELAP5-3D model and the experiment due to the location of active heater rods along the boundary between heat structures in the model. The best of these calibrations was applied to PG-29 to model the DCC. Once again, temperatures during the transient were underpredicted but trends in temperature were captured. The RELAP5-3D model captured trends in the data but could not reproduce measured temperatures exactly. This result is not attributed to deficiencies in the experimental data or to RELAP5–3D itself. Rather, this result likely arises due to the some of the assumptions and decisions made when the RELAP5-3D model was first developed, prior to the execution of HTTF experiments. An agreement in prediction of temperature trends but challenges reproducing HTTF temperatures within measurement uncertainty is consistent with previous analyses of HTTF in the literature. Future RELAP5-3D validation activities centered around HTTF may be able to provide greater insight into the code’s capabilities for HTGR modeling with a more finely nodalized model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Designing ITER motional Stark effect line shift (MSE-LS) spectrometers

As a part of ITER beam aided diagnostics, the design of Motional Stark Effect (MSE) diagnostic observing the emission from the Balmer-α line is underway. The physics of Stark splitting shows that the Stark manifold is polarization dependent, and the energy splitting results in a line shift proportional to the electric field. Due to the challenges of maintaining the calibration of the plasma facing mirrors in ITER, the conventional MSE polarimetry measurement technique is replaced with a spectral approach that is deemed more favorable in the ITER environment. The MSE line shift (LS) diagnostic is designed to quantify the Lorentz electric field magnitude by measuring the Stark manifold using visible spectroscopy. In the presence of large magnetic fields and high energy heating beams of 1 MeV, the expected Stark splitting is much larger than in typical devices. The MSE-LS design has unique challenges requiring careful consideration and modeling of its viewing geometry and photon budget. The MSE-LS approach on ITER is promising but has stringent demands on the allowable errors for the statistical and systematic fitting uncertainties. In this study, a full system model and numerical simulations of data for each sightline are completed. For a range of optical transmission fractions, photon noise analysis is conducted to determine the statistical uncertainties. This provides guidance on the spectrometer throughput, dispersion at the detector, optics, and other design choices. A conceptual design of a high throughput spectrometer with a volume phase transmission grating is presented.

Instruments & Instrumentation↗
DOE OSTI2025

Measuring the thermal conductivity of hydrogels with a bidirectional 3w method

Hydrogels are soft, water-absorbing polymer materials with diverse applications in biomedicine and agriculture. Recently, hydrogels have been proposed to encapsulate water-soluble phase change materials which store energy in their latent heat of solidification. In these applications, the thermal conductivity of these materials affects their performance. Few methods exist for measuring the thermal conductivity of small quantities of hydrogels. Here, we describe an implementation of the bidirectional 3w technique to measure the thermal conductivity of hydrogels with particular attention to their moisture content. Our implementation of the technique can probe sample volumes as little as ~20 mL and yields the thermal conductivity without requiring fitting of additional thermal parameters. We numerically simulate 3w sensor designs with frequency-domain 3-D models to quantify and reduce errors introduced by the choice of substrate and insulation layer thickness. Frequencies in the ~1−20 Hz range yield less error for the materials considered here. We verify our setup with measurements on water and report values for polyacrylamide and poly(2-acrylamido-2-methylpropane sulfonic acid) (PAMPS) hydrogels. Our swollen hydrogels exhibited thermal conductivities nearly equivalent to water, 0.6 W m-1 K-1, and we estimate thermal conductivities of 0.43 and 0.42 W m-1 K-1 for neat polyacrylamide and PAMPS, respectively. Finally, we estimate an error of ±7%, consistent with other 3ω methods, with the largest error coming from the sensor calibration. We find our implementation of the bidirectional 3w method gives reasonable results and can be employed for prototyping soft materials relevant for thermal storage.

3-omega, thermal conductivity, hydrogel, moisture ↗
DOE OSTI2023

Spark Channel Dynamics of Electrostatic Discharges

When two differently-charged objects are brought in close proximity to each other, the resulting high electric fields can cause electron avalanche breakdown of the air gap separating the objects, a process known as electrostatic discharge (ESD). If enough initial charge is stored on the objects, the electrical breakdown can proceed to ionize the air to such a degree that a highly conductive filament of plasma forms in the gap, known as a spark channel. The spark electrically bridges the air gap, resulting in a rapid pulse of current that neutralizes the charge difference. The current pulse produces significant heating of the gas in the spark, resulting in dissociation, ionization, thermal radiation, and hydrodynamic expansion. ESD presents a hazard to electrically-sensitive devices, with consequences such as economic losses (e.g. damaged electronics) or unsafe response (e.g. unintended ignition of flammable gas mixtures, initiation of detonators, etc.). For this thesis, the ESD spark is taken to occur between two conducting electrodes, with the spark channel being axisymmetric in a cylindrical coordinate system centered on the channel. An RLC-type circuit is used for the discharge model of the ESD event. The spark is treated as a time-dependent resistance that is in series with a capacitance, an inductance, and (optionally) a load resistance representing a “victim” component under threat from the ESD event. The primary motivation of this work is to use a numerical hydrodynamic model to understand the energy dissipation and transport processes in the spark. The model consists of the compressible Euler equations of mass, momentum, and energy conservation together with an Eddington/P1 approximation for thermal radiation transport. To close the hydrodynamic system, an equation of state (EOS) was fitted from tabular data for air that accounts for the dissociation and ionization of air species. The hydrodynamic equations are solved using a conservative Lagrangian finite volume method. These partial differential equations are coupled to the circuit equations by calculation of the spark resistance via numerical integration of the electrical conductivity of the channel. Computational results are compared against experimental measurements of discharge current and radial density of the spark channel.

42 ENGINEERING↗