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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗

DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗

DOE OSTI2025

A novel ignition model for low velocity impact of heterogeneous explosives based on interacting hot spots

While numerous studies have focused on the ignition of explosives occurring in high velocity impact and the associated shock-to-detonation transition, there has been growing interest in developing computational models focused on low-velocity impact regimes. A predictive low-velocity impact ignition model will be important for analyzing high explosive safety and potential accident scenarios. This work introduces a novel ignition model based on the concept of thermally interacting hot spots to simulate low velocity impacted heterogeneous explosives where observed ignition times are on the order of milliseconds. The model asserts that relevant hot spots are micron-sized, the typical separation between neighboring hot spots is on the order of a hundred microns, and that neighbors interact thermally through heat conduction across the interstitial region between them. To achieve tractable numerical solutions, hot spots are assumed to form a periodic array as opposed to the highly irregular positioning in an actual explosive. This idealization allows a single two hotspot system to characterize the ignition process. Consequently, the model is referred to as the two hot spot Frank-Kamenetskii ignition model. In the present study, hot spots are modeled as constant heat sources terms, but this can be extended to include grain-scale phenomena like frictional heating of micron-sized growing cracks that are confined under high pressure. Because the micron-sized features are below the scale that can be efficiently resolved at a systems level, an efficient subscale scheme based on the Method of Weighted Residuals (MWR) is used to efficiently solve the equations. In conclusion, we carry out numerical examples and analytic predictions illustrating the accuracy and the functioning of the model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

COMSOL Models and Simulation Results for Closed-Loop Geothermal Systems with Thermal Enhancements

The Closed-Loop Geothermal Working Group and the National Renewable Energy Laboratory (NREL) have studied the performance of disc-shaped and linear thermal enhancements (conductive material introduced around a wellbore) in closed-loop geothermal systems. COMSOL models were developed to investigate the increase in heat production in systems with thermal enhancements compared to those without, across a range of scenarios. This dataset includes two example COMSOL models, a summary of the simulation results, and a link to an accompanying paper. Studied are idealized versions of two thermal enhancement designs: (1) existing or newly created fractures surrounding the wellbore, filled with highly conductive material, and (2) a 'fishbone structure,' where short wellbores are drilled radially outward from the main wellbore and filled with highly conductive material.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Self-consistent integrated modeling of combined hybrid discharge-laser produced plasma devices for extreme ultraviolet metrology

Discharge- and laser-produced plasma (DLPP) devices are being used as light sources for extreme ultraviolet (EUV) generation. A key challenge for both, DPP and LPP, is achieving sufficient brightness to support the throughput requirements of nanometrology tools. To simulate the environment of a hybrid DLPP device and optimize the EUV output, we have developed an integrated HEIGHTS-DLPP computer simulation package. The package integrates simulation of two evolving plasmas (DPP and LPP) and includes modeling of a set of integrated self-consistent processes: external power source and plasma energy balance, plasma resistive magnetohydrodynamics, plasma heat conduction, detailed radiation transport (RT), and laser absorption and refraction. We simulated and optimized DLPP devices using Xe gas as a target material. We synchronized the external circuit parameters, chamber gas parameters, and laser beam temporal and spatial profiles to achieve maximum EUV output. The full 3D Monte Carlo scheme was integrated for detailed RT and EUV output calculations in Xe using more than 3600 spectral groups. The modeling results are in good agreement with Julich Forschungszentrum experimental data. Theoretical models, developed and integrated into the HEIGHTS package, showed wide capabilities and flexibility. In conclusion, the models and package can be used for optimization of the experimental parameters and settings, investigation of DLPP devices with complex design, analyzing the impact of integrated spatial effects and working timeline arrangement on the final EUV output, and EUV source size, shape, and angular distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

High-Temperature Latent Heat Storage System using Transportable Heat Pipes for Versatile Integration with Emerging Microreactors

With the growing interest in the development of micro nuclear reactors and their application to microgrid, high-temperature thermal energy storage (TES) has also been receiving great attention as a technology to reinforce the energy utilization and profitability of this class of reactors. Given the importance and research needs of the TES technology, a research project is underway at Idaho National Laboratory (INL) to develop high-temperature TES that can be easily interfaced with emerging micro nuclear reactors. Currently under development is the Heat pipe-Integrated Thermal Battery (HITB), a novel latent heat storage design that enables flexible and robust integration with any design of currently considered micro nuclear reactors. The development of detailed design based on the proposed TES concept, HITB, requires not only development and validation of modeling and simulation capabilities, but also research that spans thermal performance analysis, component- and system-scale design optimization, and experimental demonstrations including high-temperature feasibility testing. The current research efforts (and plans) are divided into the following four areas: (i) Thermal and structural modeling and analysis at both system- and component-scales for evaluating performance and conducting design optimization, (ii) Small-scale experimental demonstrations for the key components of the proposed HITB design, (iii) Modelica-based dynamic model development for developing a control strategy under dynamic operation within integrated energy systems configuration, (iv) Integrated experimental testing and demonstration. This paper introduces the overall research efforts that have been and are being conducted to develop and demonstrate the novel design of the high-temperature latent heat storage device, HITB. The key design features, ongoing research efforts including the development of modeling and simulation capabilities, small-scale demonstrations of the HITB’s key functionalities prior to scaled-up experimental demonstration, are discussed. Then, this paper concludes with a remark on the future research plans for this project.

25 ENERGY STORAGE↗

DOE OSTI2022

Analytical benchmark solution for 1-D neutron transport coupled with thermal conduction and material expansion

In this work we present fully analytical solutions for a class of finite, homogeneous, 1-D slab benchmark problems with nonlinear temperature feedback effects. The proposed class of benchmarks include multiplicative 1-D neutron transport (limited to quasistatic S{sub 2} with μ = ±1) coupled with thermal conduction, convection, Doppler broadening, and expansion effects along the length of the slab. This class of benchmark models, along with the corresponding analytical solutions, are valuable for validating multiphysics analysis frameworks that support coupled neutronics/thermal/structural calculations. Analytical solutions for the benchmarks are obtained by introducing an ansatz that the equilibrium flux and temperature distributions in the slab have the same shape. Specific values for the total microscopic cross section, σ{sub t,0}, and conductive heat transfer coefficient, h, that satisfy the assumed ansatz are then determined. A discussion of the procedure for generating benchmark models and analytical solutions is provided, along with numerical results for an example set of model parameters and thoughts on practical applications of the benchmarks for multiphysics code validation. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Non-Contact Mass Density and Thermal Conductivity Measurements of Organic Thin Films Using Frequency–Domain Thermoreflectance

Thin-film organic materials are broadly used to study amorphous stabilization of active pharmaceuticals, control explosive detonation phenomena, and introduce insulation in novel thermal barriers. Their synthesis, however, introduces defects and thickness variations that warrant careful characterization of local thermophysical properties such as thermal conductivity and mass density. In this study, wide bandwidth (200 Hz to 20 MHz) frequency–domain thermoreflectance (FDTR) is demonstrated to simultaneously extract the thermal conductivity and mass density of 1 μm physical vapor-deposited indomethacin films on Si and SiO 2 substrates, as well as 10 and 100 μm films on Si. By assuming a bulk specific heat capacity, mass densities are determined with FDTR measurements of volumetric heat capacity and are in good agreement with the literature, as well as models based upon a dependence on porosity and the kinetic theory for phonons. Lastly, it is found that for broad-band FDTR measurements, insulating substrates provide improved fidelity for the extraction of thermal conductivity and volumetric heat capacity in organic thin films. Overall, this work demonstrates the potential for FDTR as a non-contact method to determine microscale mass density variations across the surface and thickness of organic thin films.

36 MATERIALS SCIENCE↗

DOE OSTI2021

An efficient numerical model for predicting residual stress and strain in parts manufactured by laser powder bed fusion

Abstract Computational modeling of additively manufactured structures plays an increasingly important role in product design and optimization. For laser powder bed fusion processes, the accurate modeling of stress and distortion requires large amount of computational cost due to very localized heat input and evolving complex geometries. The current study takes advantage of a graphics processing unit accelerated explicit finite element analysis code and approximated heat conduction analysis to predict the macroscopic thermo-mechanical behavior in laser selective melting. Adjacent layers and tracks were lumped to reduce the number of time steps and elements in the finite element model. The effects of track and layer grouping on prediction accuracy and solution efficiency are investigated to provide a guidance for a cost-effective simulation. Thin-wall builds from Inconel alloy 625 (IN625) powders were simulated by applying the developed modeling approach to get the detailed residual stress and distortion at a computational speed 50 times higher than conventional approach. Under repeated heating and cooling cycles, a high tensile stress was produced near surfaces of a build due to a larger shrinkage on surface than that in central area. It is also shown that horizontal stresses concentrate near the root and top layers of the IN625 build. The predicted residual elastic strain distribution was validated by the experimental measurement using x-ray synchrotron diffraction.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Low-temperature nanoscale heat transport in a gadolinium iron garnet heterostructure probed by ultrafast x-ray diffraction

Time-resolved x-ray diffraction has been used to measure the low-temperature thermal transport properties of a Pt/Gd3Fe5O12//Gd3Ga5O12 metal/oxide heterostructure relevant to applications in spin caloritronics. A pulsed femtosecond optical signal produces a rapid temperature rise in the Pt layer, followed by heat transport into the Gd3Fe5O12 (GdIG) thin film and the Gd3Ga5O12 (GGG) substrate. The time dependence of x-ray diffraction from the GdIG layer was tracked using an accelerator-based femtosecond x-ray source. The ultrafast diffraction measurements probed the intensity of the GdIG (1 −1 2) x-ray reflection in a grazing-incidence x-ray diffraction geometry. The comparison of the variation of the diffracted x-ray intensity with a model including heat transport and the temperature dependence of the GdIG lattice parameter allows the thermal conductance of the Pt/GdIG and GdIG//GGG interfaces to be determined. Complementary synchrotron x-ray diffraction studies of the low-temperature thermal expansion properties of the GdIG layer provide a precise calibration of the temperature dependence of the GdIG lattice parameter. The interfacial thermal conductance of the Pt/GdIG and GdIG//GGG interfaces determined from the time-resolved diffraction study is of the same order of magnitude as previous reports for metal/oxide and epitaxial dielectric interfaces. The thermal parameters of the Pt/GdIG//GGG heterostructure will aid in the design and implementation of thermal transport devices and nanostructures.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Comparative analysis of thermal management systems in electric vehicles at extreme weather conditions: Case study on Nissan Leaf 2019 Plus, Chevrolet Bolt 2020 and Tesla Model 3 2020

With the surge in electric vehicle (EV) adoption and the need for extended driving ranges, optimizing energy efficiency, particularly through thermal management, is critical, especially in extreme weather. Managing the substantial energy needed for cabin climate control and battery temperature regulation can increase energy demands by over 50 %, severely limiting range. This study conducts a comparative analysis of thermal management systems (TMS) in three popular EV vehicles, 2020 Chevrolet Bolt, 2019 Nissan Leaf Plus, and 2020 Tesla Model 3, evaluating their distinct TMS configurations and performance under varied weather conditions. Using both numerical simulations and experimental data collected on a controlled test bench at Argonne National Laboratory, we assess how TMS architecture and operational modes influence energy consumption and range. A comprehensive TMS model was developed, integrating cabin and battery thermal sub-models in the Autonomie software platform, to simulate temperature fluctuations and range impacts. Cabin climate was modeled using a mono-zonal approach, while battery cell temperature distribution was estimated through a 2D nodal structure. Each vehicle's distinct TMS setup was evaluated: the Chevrolet Bolt and Tesla Model 3 use a dual evaporator vapor compression cycle with a PTC heater for the cabin and a coolant loop for battery thermal management; the Nissan Leaf Plus employs a heat pump with a PTC heater for the cabin and air-cooling for the battery. Tests conducted at ambient temperatures of 35°C, 22°C, -7°C, and -18°C reveal significant differences in energy use and range reduction across both configurations and conditions. At 35°C, the Tesla Model 3, Chevrolet Bolt, and Nissan Leaf Plus have a range reduction of 8%, 9%, and 13%, respectively, due to air conditioning. In winter, heating technology is paramount; at -7°C, the Nissan Leaf's heat pump configuration achieves a lower range reduction (19.3%) compared to the Tesla and Chevrolet Bolt PTC heaters, which reduce range by 28.3% and 31%, respectively. Further, this study provides valuable insights for automotive engineers, EV technology researchers, and thermal management system designers aiming to enhance electric vehicle performance by understanding how different weather conditions and TMS architectures impact energy consumption and driving range.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Validation and Verification for INL Modelica-based TEDS models Via Experimental Results

This report provides an overview on the verification and validation (V&V) of the Thermal Energy Distribution System (TEDS) model developed in the Modelica process modeling ecosystem using experimental data. Model development has led to the creation of a dynamic process model of the experimental TEDS facility housed within the Energy Systems Laboratory (ESL) at Idaho National Laboratory (INL). The model was then used during the preconstruction phase of the experimental effort to inform experimental design (e.g., insulation requirements, bypass line placement, expected performance of components) and to test innovative control schemes prior to the initial operation. The TEDS model developed in Modelica includes the primary components of the TEDS experimental unit: a 200kW Chromalox heater; a single-tank packed-bed thermal energy storage system filled with 0.125-inch alumina (Al2O3) beads; an ethylene-glycol-to-Therminol-66 heat exchanger; system piping; five control valves; and all associated temperature, pressure, and volumetric flow sensors. Using the Institute of Electrical and Electronics Engineers (IEEE) V&V methodologies, considered the gold standard in the engineering field, the model was verified using a combination of static analysis, spatial convergence, and regression tests. Then using dynamic time warping (DTW) initial runs to validate and tune the TEDS model versus the experiment were conducted. This tuning method was accomplished using the INL Risk Analysis Virtual ENvironment (RAVEN) software package. Tuning is required to account for physical phenomena that are less understood within the empirical heat transfer correlations. Through the commencement of this work, a systems-level model of TEDS with associated control systems, sensors, piping diameters, and component capabilities has been created. This model was utilized in the pre-experimental phase to inform system design, insulation thicknesses, and potential control schemes to operate the system effectively and safely. Then, initial experimental startup and operational data were used to demonstrate the validation and tuning methodology. This process demonstrates the classical two-step approach of a model informing experimental design followed by the experiment validation and tuning the model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗

DOE OSTI2020

Construction of a Code Verification Matrix for Heat Conduction With Finite Element Code Applications

When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.

42 ENGINEERING↗

DOE OSTI2024

High-Pressure Melting Experiments of Fe 3 C and a Thermodynamic Model of Fe-C Liquids for the Earth's Core

Melting experiments of Fe 3 C were conducted to 85 GPa in laser-heated diamond anvil cells with in situ X-ray diffraction and post-experiment textural observation. From the determined pressure-temperature conditions of the melting curve for Fe 3 C, together with literature data on the melting point of diamond and eutectic point of the system Fe-Fe 3 C/Fe 7 C 3 under high pressures, we established a self-consistent thermodynamic model for high-pressure melting of the system Fe-C including the mixing parameters for liquids. The results show that mixing of Fe and C liquids is negatively nonideal from 1 bar to the pressure at the center of the Earth. The departure from ideal mixing becomes progressively larger with increasing pressure, which leads to greatly stabilized liquids under core pressures. The modeled carbon content in eutectic melts under core pressures is 3.3–4.4 wt%. From the Gibbs free energy, we derived an internally consistent parameters for Fe-C outer cores which included the crystallizing points at their bottoms, isentropic thermal profiles, and densities and longitudinal seismic wave speeds (Vp). While the addition of carbon in excess of the eutectic melt composition effectively reduces the density of iron liquid, the Vp of iron liquid is not greatly changed. Therefore, the low density and high Vp of PREM relative to pure iron cannot be reconciled by an Fe-C liquid. Therefore, the Earth's core cannot be approximated by the system Fe-C and should include another light element.

58 GEOSCIENCES↗

DOE OSTI2023

Numerical investigations to identify environmental factors for field-scale reactive transport of pathogens at riverbank filtration sites

While induced bank filtration is a proven method for facilitating sustainable drinking water production, it is at risk from surface water contaminations (e.g., pathogens). Induced bank filtration and pathogen transport in groundwater have been studied extensively. However, long-term studies that consider real-world conditions are missing. These conditions include seasonal changes to environmental conditions and waterworks operations. Therefore, to analyze the effect of seasonal changes on the transport of human pathogenic viruses and their indicators in induced bank filtration, concentrations of adenoviruses and pathogen indicators were monitored over 16 months at an active bank filtration plant at the Rhine River, in Düsseldorf (Germany). Based on this data, a 2D groundwater model was created in PFLOTRAN that simulated flow, heat transport, conservative transport of chloride and the resulting electrical conductivity, reactive transport of oxygen and nitrate, and colloid-based transport of coliforms, somatic coliphages, and adenoviruses. The results show that reduced travel time was the key factor determining periods with a low removal of coliforms and somatic coliphages in the aquifer. Travel time was controlled by river level variations during rainy seasons, and the waterworks extraction rates during dry seasons. Further, for adenovirus transport, travel distance in the subsurface appeared to be the key factor, while travel time had no significant impact. Coliform removal increased when the colmation layer permeability decreased, while coliphage and adenovirus removal was unaffected by the colmation layer permeability. Seasonal changes in temperature and oxygen content did not significantly impact the removal of coliphages and adenoviruses in groundwater. Denitrifying conditions correlated with a lowered coliform removal, but the modelling could not establish a connection between denitrifying conditions and coliform removal. Our study showed that removal of pathogens and pathogen indicators at induced bank filtration plants varies greatly in time and space (e.g., for coliforms from 1 to 4 log-levels at 20m travel distance), and that adenovirus transport differs considerably from transport of coliforms and somatic coliphages.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Assessing thermal resilience of an assisted living facility during heat waves and cold snaps with power outages

Extreme hot and cold weather events are becoming more frequent, intense, and longer due to climate change. When these events occur coincidentally with power outages, the resulting extreme indoor temperatures pose a severe health hazard for occupants. This study conducted a holistic modeling and analysis of an assisted living facility, where senior residents live, to assess its thermal resilience performance under a six-day heat wave in 2015 and a three-day cold snap in 2021 with power outages. Impacts of 13 energy efficiency measures on thermal resilience and backup power capacity of the facility were evaluated. Three thermal resilience metrics: the SET (standard effective temperature) degree-hours, the Heat Index, and the Hours of Safety, were used and calculated from the EnergyPlus simulation models. Furthermore, major findings are: (1) the facility would suffer from extreme temperatures during the cold and hot events without a power supply, not meeting the passive survivability requirements; (2) most passive envelope measures improve thermal resilience for both hot and cold events, but making the building envelope airtight results in conflicting performance between the hot and cold events; (3) natural ventilation is an effective measure to mitigate summer indoor overheating; and (4) the energy efficiency package can reduce backup power capacity by 19% for the three-day cold snap. It is recommended that building technologies and design strategies be evaluated to consider co-benefits of energy use, thermal resilience, and backup power needs through building energy codes or policies for existing and new buildings, which are transitioning for decarbonization and climate resilience.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗