Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective temperature model”

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

Operational Parameter Database for Molten Salt Thermal Energy Storage Tank Modeling

The second generation of concentrated solar power (CSP) plants is characterized by the use of a central receiver (either cavity or external), two molten nitrate salt tanks (60 wt.% NaNO 3 and 40 wt.% KNO 3 ), and a steam Rankine power-generation cycle connected through a primary heat exchanger. Molten salt thermal energy storage (TES) tanks have been widely deployed in commercial CSP plants worldwide and have been essential for increasing plant dispatchability and capacity factor, while also reducing the levelized cost of electricity (LCOE). These systems enable energy storage at the gigawatt-hour scale, typically providing 6 to 17 hours of storage duration. Despite being a commercial technology, the multiple failures observed after only a few months or years of operation in plants around the world demonstrate the technology's relative infancy and highlight the need for further research to improve its reliability. The National Laboratory of the Rockies (NLR), in collaboration with industry partners and academic and research institutions, has been leading multiple projects funded by the U.S. Department of Energy (DOE). These projects focus on addressing molten salt tank failures by improving tank design and welding fabrication practices, evaluating new alloys and weld fillers, and providing guidelines for tank commissioning and safe operation. In particular, this report presents modeling results on the effect of key tank operation parameters during 60 minutes of operation, including the mass flow rate and temperature of the salt inflow, tank salt inventory temperature, and inventory level for a representative molten salt tank design. These results form a database of tank operation behaviors that captures the effects of each specific parameter during charging, charging/discharging, and discharging processes.

14 SOLAR ENERGY↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Real-fluid behavior in rapid compression machines: Does it matter?

Rapid compression machines (RCMs) have been extensively used to quantify fuel autoignition chemistry and validate chemical kinetic models at high-pressure conditions. Historically, the analyses of experimental and modeling RCM autoignition data have been conducted based on the adiabatic core hypothesis with ideal gas assumption, where real-fluid behavior has been completely overlooked, though this might be significant at common RCM test conditions. Here, this work presents a first-of-its-kind study that addresses two significant but overlooked questions for autoignition studies within RCMs in the fundamental combustion community: (i) experiment-wise, can unaccounted-for real-fluid behavior in RCMs affect the interpretation and analysis of RCM experimental data? and (ii) simulation-wise, can unaccounted-for real-fluid behavior in RCMs affect RCM autoignition modeling and the validation of chemical kinetic models? To this end, theories for real-fluid isentropic change are newly proposed and derived based on high-order Virial EoS, and are further incorporated into an effective-volume real-fluid autoignition modeling framework newly developed for RCMs. With detailed analyses, the strong real-fluid behavior in representative RCM tests is confirmed, which can greatly influence the interpretation of RCM autoignition experiments, particularly the determination of end-of-compression temperature and evolution of the adiabatic core in the reaction chamber. Furthermore, real-fluid RCM modeling results reveal that considerable error can be introduced into simulating RCM autoignition experiments when following the community-wide accepted effective-volume approach by assuming ideal-gas behavior, which can be as high as 64% in the simulated ignition delay time at compressed pressure of 125 bar and lead to contradictory validation results of chemical kinetic models. Therefore, we recommend the community to adopt frameworks with real-fluid behavior fully accounted for (e.g., the one developed in this study) to analyze and simulate past and future RCM experiments, so as to avoid misinterpretation of RCM autoignition experiments and eliminate the potential errors that can be introduced into the simulation results with the existing RCM modeling frameworks.

High-order Virial equation of state↗
DOE OSTI2025

E3SM‐GCAM: A Synchronously Coupled Human Component in the E3SM Earth System Model Enables Novel Human‐Earth Feedback Research

Abstract Modeling human‐environment feedbacks is critical for assessing the effectiveness of climate change mitigation and adaptation strategies under a changing climate. The Energy Exascale Earth System Model (E3SM) now includes a human component, with the Global Change Analysis Model (GCAM) at its core, that is synchronously coupled with the land and atmosphere components through the E3SM coupling software. Terrestrial productivity is passed from E3SM to GCAM to make climate‐responsive land use and CO 2 emission projections for the next 5‐year period, which are interpolated and passed to E3SM annually. Key variables affected by the incorporation of these feedbacks include land use/cover change, crop prices, terrestrial carbon, local surface temperature, and climate extremes. Regional differences are more pronounced than global differences because the effects are driven primarily by differences in land use. This novel system enables a new type of scenario development and provides a powerful modeling framework that facilitates the addition of other feedbacks between these models. This system has the potential to explore how human responses to climate change impacts in a variety of sectors, including heating/cooling energy demand, water management, and energy production, may alter emissions trajectories and Earth system changes. Plain Language Summary Earth system models help us understand how humans are changing the climate. Currently, these models do not include human systems, so predetermined greenhouse gas, aerosol, and land use change data are input to these models. These data do not reflect human responses to changes projected by Earth system models. We have added a human component to an Earth system model to represent human responses to environmental change and calculate corresponding greenhouse gas and land use change data instead of using predetermined data. Including this human component changes projections of land use, land carbon storage, and regional climate. Key Points We have incorporated a novel, advanceable human component in an Earth system model to simulate human‐Earth feedbacks Including terrestrial productivity feedbacks from the Earth to the human systems affects land change, crop prices, carbon, and climate Regional effects of including terrestrial productivity feedbacks are greater than global effects because land change is the main driver

E3SM↗
DOE OSTI2021

Global gyrokinetic simulation of neoclassical ambipolar electric field and its effects on microturbulence in W7-X stellarator

Global neoclassical simulations of a model equilibrium of the W7-X stellarator find an ambipolar electric field with either an ion root or an electron root solution depending on the temperature ratio between electrons and ions. The ambipolar electric field is then used as an equilibrium field in the turbulence simulations of ion temperature gradient (ITG) instability. The shear of the ambipolar electric field has modest effects on the ITG linear instability, nonlinear saturation amplitude, and turbulent transport in the ion root case. However, in the electron root case, the ambipolar electric field significantly reduces the linear ITG growth rate, turbulence intensity, and radial correlation length, resulting in an ion heat conductivity comparable to the neoclassical transport level in the strong shear region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Fuel Flow Modeling under Thermal Fluidics Considerations within Operational Domains of Molten Salt Reactors

Molten Salt Reactors (MSRs) offer significant versatility supported by their design flexibility and safety advantages due to high temperatures, design customization capabilities, strong negative reactivity feedback effects, and salt chemistry characteristics. Different control schedules and different mass flow rates in the MSR primary systems impact neutronics and thermal hydraulics coupling. The result is two-fold. First, it leads to adaptability in various applications when the customization of a MSR unit is desired. Second, it calls for the quantification of dynamics and safety characteristics to assess the operational domain and make sure reactor stability and inherent safety characteristics are maintained. This paper provides an in-depth transient analysis of the operational implications of MSRs having heat-generating fuel and heat-transporting salt mixed forming liquid fuel salts contained within primary systems. A coupling model was developed to simulate fuel flow in molten salt reactors. A zero-dimensional reactor kinetics model was used with a one-dimensional heat transfer model in the reactor core. The temperature reactivity feedbacks of the fuel salt resulting from the thermal-hydraulics model were used in the reactor kinetics model to complete the coupled code. The model was applied to simulate the Molten Salt Breeder Reactor (MSBR) response at steady state and transients. The developed model was shown to be a suitable tool for the dynamic analysis of molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

An international study on THM modelling of the full-scale heater experiment at Mont Terri laboratory

We present results from an international model comparison study of the Full-Scale Emplacement (FE) experiment in Opalinus Clay at the Mont Terri Laboratory, Switzerland. Based on a provided parameter set the teams decided which parameters they adopted for their models, whether they considered the excavation and the ventilation phase in addition to the heating phase and if they included technical features like the shotcrete or the EDZ. The teams were able to reproduce the measured parameters temperature, relative humidity and pore pressure. The modelled results for temperature agree very closely between the teams especially in the sensors in Opalinus Clay. All teams were able to reproduce the redistribution of water in the bentonite backfill due to heating. The evolution of the relative humidity showed similar trends with differences in the intensity of the dry out effect. To model the pore pressure evolution is more complex because it comprises the full interaction of the coupled THM processes. The spread between the pore pressure modelled by the teams was larger, with some teams overestimating the pressure increase due to heating and some teams overestimating the extent of drainage. The agreement of modelled results with measurements improves with larger distance to the heater. We conclude that the EDZ and the shotcrete potentially influence the behaviour of the rock causing higher differences closer to the heater. Further research is needed to better implement those influences into the models. Based on the calibrated models, the future evolution of temperature, relative humidity and pore pressure was predicted over the next 10 years following a change of the heat power applied in 2023 and 2024. Again, the predicted temperatures agree very closely between the teams. Most teams do not expect an increase in relative humidity during the next 10 years after the initial dry-out.

58 GEOSCIENCES↗
DOE OSTI2022

Unraveling the role of EGR olefins at advanced combustion conditions in the presence of nitric oxide: Ethylene, propene and isobutene

We report the role of EGR (exhaust gas recirculation) olefinic constituents at advanced combustion conditions in the presence of nitric oxide is unraveled in this study through experimental and modeling efforts using a twin-piston rapid compression machine operating at a stochiometric fuel loading with 20% EGR by mass, pressures of 20 and 40 bar, and temperatures from 680 to 950 K. Five different levels of olefin addition, focusing on ethylene, propene and isobutene, with a fixed amount of NO at 70 ppm are doped into test mixtures of PACE-20, a multi-component gasoline surrogate, where olefin addition effects are characterized through changes in ignition times and heat release rates. Experiments indicate that all three EGR olefins inhibit autoignition reactivity and low-temperature heat release at Tc < 850 K, with isobutene exhibiting the greatest impact, while at Tc > 850 K, low ethylene and propene additions promote reactivity. A recently updated chemical kinetic model, with detailed gasoline/NO x interacting and olefin/NO x interacting chemistry incorporated, is adopted to simulate the experiments. Simulation results are somewhat inconsistent with the experiments, where the model captures the inhibiting effects of all olefins on first-stage ignition reactivity, while consistently predicting a promoting effect on main ignition reactivity. Sensitivity and rate of production analyses reveal that adding olefins greatly alters the role of the consuming pathways for the olefins and their primary derivatives at the initial stage of the oxidation process, particularly with the presence of NO, where the olefins and their derivatives interact with both NO x species such as NO 2 and other species such as OH and HO 2 . The olefin/NO x interactions are particularly pronounced with propene and isobutene addition, where these lead to increased ignition reactivity by facilitating NO production and are mostly responsible for the disagreement between the model and experiments. Further investigations of olefin interacting chemistry, particularly those with NOx species, are needed for chemistry models to accurately predict the complicated effects of EGR.

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

On the nonlinear temperature dependence of residual stresses in solid oxide fuel cells

The residual stresses in a dense yttria-stabilized zirconia (YSZ) layer bonded to a porous NiO-YSZ substrate were measured by X-ray diffraction in air between 25°C and 900°C. In this work, it was found that the residual stresses on the outer surface of the YSZ layer were compressive, that their magnitude increased with decreasing temperature, and that the temperature dependence was nonlinear particularly at temperatures below 600°C. An elastic model that incorporates micromechanics and laminate theory was formulated to successfully predicts the nonlinearity of the temperature dependence of the residual stresses by accounting for the effect of the order-disorder transition experienced by YSZ around 600°C and the magnetic transition of NiO at the Néel temperature around 265°C, on the properties of these materials. The implications of the results for the reliability of solid oxide fuel cells are discussed.

25 ENERGY STORAGE↗
DOE OSTI2025

Interactions Between Climate Mean and Variability Drive Future Agroecosystem Vulnerability

ABSTRACT Agriculture is crucial for global food supply and dominates the Earth's land surface. It is unknown, however, how slow but relentless changes in climate mean state, versus random extreme conditions arising from changing variability , will affect agroecosystems' carbon fluxes, energy fluxes, and crop production. We used an advanced weather generator to partition changes in mean climate state versus variability for both temperature and precipitation, producing forcing data to drive factorial‐design simulations of US Midwest agricultural regions in the Energy Exascale Earth System Model. We found that an increase in temperature mean lowers stored carbon, plant productivity, and crop yield, and tends to convert agroecosystems from a carbon sink to a source, as expected; it also can cause local to regional cooling in the earth system model through its effects on the Bowen Ratio. The combined effect of mean and variability changes on carbon fluxes and pools was nonlinear, that is, greater than each individual case. For instance, gross primary production reduces by 9%, 1%, and 13% due to change in mean temperature, change in temperature variability, and change in both temperature mean and variability, respectively. Overall, the scenario with change in both temperature and precipitation means leads to the largest reduction in carbon fluxes (−16% gross primary production), carbon pools (−35% vegetation carbon), and crop yields (−33% and −22% median reduction in yield for corn and soybean, respectively). By unambiguously parsing the effects of changing climate mean versus variability and quantifying their nonadditive impacts, this study lays a foundation for more robust understanding and prediction of agroecosystems' vulnerability to 21st‐century climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Future precipitation increase constrained by climatological pattern of cloud effect

The fractional increase in global mean precipitation (?P ¯/P ¯) is a first-order measure of the hydrological cycle intensification under anthropogenic warming. However, ?P ¯/P ¯ varies by a factor of more than three among model projections, hindering credible assessments of the associated climate impacts. The uncertainty in ?P ¯/P ¯ stems from uncertainty in both hydrological sensitivity (global mean precipitation increase per unit warming) and climate sensitivity (global mean temperature increase per forcing). Here, by investigating hydrological and climate sensitivities in a unified surface-energy-balance perspective, we find that both sensitivities are significantly correlated with surface shortwave cloud feedback, which is further linked to the climatological pattern of cloud shortwave effect. The observed pattern of cloud effect thus constrains both sensitivities and consequently constrains ?P ¯/P ¯. The 5%-95% uncertainty range of ?P ¯/P ¯ from 1979-2005 to 2080-2100 under the high-emission (moderate-emission) scenario is constrained from 6.34±3.53% (4.19±2.28%) in the raw ensemble-model projection to 7.03±2.59% (4.63±1.71%). The constraint thus suggests a higher most-likely ?P ¯/P ¯ and reduces the uncertainty by ~25%, providing valuable information for impact assessments.

Zhou, Wenyu↗
DOE OSTI2021

Oxidation of an iso-paraffinic alcohol-to-jet fuel and n-heptane mixture: an experimental and modelling study

Oxidation of ATJ/n-heptane blends was studied over a wide range of test conditions, using single pulse shock tubes. Test conditions were designed to study the effect of pressure (4 and 50 bar), and fuel loading (~100-1400 ppm) on the oxidation of the blends across a wide range of temperatures (800-1300 K). These effects were observed by measuring concentrations of intermediate species formed using in-line Gas Chromatography (GC) and GCxGC TOF-MS. Results showed that increasing the initial fuel load does not have a significant impact on the results. However, increasing the pressure shifts the mixture reactivity to a lower temperature by about 150 K and causes the fuel to be oxidized instead of decomposing pyrolytically as at lower pressures. Additional experiments for pure ATJ and pure n-heptane were performed at conditions matching the ATJ/n-heptane 50 bar experiments to analyze the differences between the pure and mixed fuels. Speciation measurements were compared against predictions from a detailed kinetic model. The ability of the kinetic model to capture the effects of varying experimental test conditions on the evolution of intermediate species is discussed and kinetic analyses have been conducted to identify the important reaction pathways.

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

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Constraints on Λ CDM extensions from the SPT-3G 2018 $EE$ and $TE$ power spectra

Here, we present constraints on extensions to the Λ CDM cosmological model from measurements of the E-mode polarization autopower spectrum and the temperature-E-mode cross-power spectrum of the cosmic microwave background (CMB) made using 2018 SPT-3G data. The extensions considered vary the primordial helium abundance, the effective number of relativistic degrees of freedom, the sum of neutrino masses, the relativistic energy density and mass of a sterile neutrino, and the mean spatial curvature. We do not find clear evidence for any of these extensions, from either the SPT-3G 2018 dataset alone or in combination with baryon acoustic oscillation and Planck data. None of these model extensions significantly relax the tension between Hubble-constant, H 0 , constraints from the CMB and from distance-ladder measurements using Cepheids and supernovae. The addition of the SPT-3G 2018 data to Planck reduces the square-root of the determinants of the parameter covariance matrices by factors of 1.3–2.0 across these models, signaling a substantial reduction in the allowed parameter volume. We also explore CMB-based constraints on H 0 from combined SPT, Planck, and ACT DR4 datasets. While individual experiments see some indications of different H 0 values between the TT, TE, and EE spectra, the combined H 0 constraints are consistent between the three spectra. For the full combined datasets, we report H 0 = 67.49 ± 0.53 km s -1 Mpc -1 , which is the tightest constraint on H0 from CMB power spectra to date and in 4.1σ tension with the most precise distance-ladder-based measurement of H 0 . The SPT-3G survey is planned to continue through at least 2023, with existing maps of combined 2019 and 2020 data already having ~ 3.5 x lower noise than the maps used in this analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A new dynamic zOnal model with air-diffuser (DOMA) - Application to thermal comfort prediction

A new Dynamic zOnal Model with Air-diffuser (DOMA) was developed. Several case studies were investigated and tested to evaluate and validate this program using measurement data. This new model was integrated into a TRaNsient SYstems Simulation program library and coupled with the multi-zone thermal model. The DOMA/TRNSYS coupled model was then used to predict room temperature distribution over an entire day of a single-zone building. The results show that increasing the heating outputs of the electric floor system, for example, from 75 to 200 W/m 2 , would not effectively improve the indoor thermal comfort, since the thermostat will reach the set point first and then turn off the system before the room gets enough heat and reach a comfortable level. This indicates the importance of selecting an appropriate location and set point for the thermostat when using a floor heating system. This potential thermal comfort issue can only be identified through the two-node model with a dynamic zonal model rather than the conventional PMV model, which thus suggests that for optimizing indoor thermal comfort of a building equipped with a time-sensitive control strategy and/or HVAC system, the TSENS results obtained from the two-node model integrated with DOMA are more appropriate than PMVs.

Construction & Building Technology↗
DOE OSTI2021

TRANSIENT THERMAL MODELING OF THE HIGH BURNUP DEMONSTRATION RESEARCH PROJECT CASK USING STAR-CCM+ AND COBRA-SFS

The Department of Energy in collaboration with the Electric Power Research Institute is in the process of conducting the High Burnup Demonstration Project. Where the objective is to characterize the performance of high-burnup fuel in long term storage. As part of this demonstration, a TN-32B dry storage cask was instrumented and loaded with spent nuclear fuel at North Anna Nuclear Generating Station in November 2017. The project cask provides a unique opportunity to gain information on spent fuel and cask performance from an in-service operational system. The cask was instrumented with thermocouples inside of the fuel assemblies, then loaded and dried using normal procedures. After the drying process a thermal soak period where the cask was left indoors was used to obtain steady temperatures for model comparison along with surface temperatures. This paper details thermal modeling validation work that was done to model both the steady state and transient cases. Two modeling tools were used to predict temperatures in the cask. The general purpose CFD and heat transfer code STAR-CCM+ was used with both a detailed pin-by-pin model was used along with a more efficient k-effective simplification. The other code is COBRA-SFS, a purpose-built detailed thermal modeling tool developed and maintained at PNNL. Results from all the modeling tools were compared blind to the test data. Each model utilized general design information and compared reasonably well to the blind steady state data using convection and surface temperature boundary conditions. Post-test adjustments were later made to better reflect the “as built” conditions of the cask. Based on lessons learned from the steady state analysis the vacuum drying process was modeled in both codes. Utilizing the measured loading conditions the transient models were able to compare very well with the measured data. Overall the project showed an ability to model spent fuel storage conditions very well and future work is planned to generalize the methodology used for vacuum drying.

Thermal Analysis, Spent Fuel, Dry Storage↗
DOE OSTI2022

Development and Validation of Resistance-Capacitance Model for Phase Change Material Embedded in Porous Media

A simple and effective way to enhance PCM thermal conductivity is by embedding porous materials with high thermal conductivities. In this work, CFD was used to evaluate the melting of PCM embedded in metal foam supplied with constant heat flux. Good agreement between experimental and numerical results has been found with a mean PCM local temperature deviation of 1.62K and 2.97K in two locations. The CFD simulation time was 56 hours for 6.5 hours of flow time. A computationally efficient resistance capacitance-based model (RCM) was developed to improve simulation speed. The proposed RCM exhibited mean local temperature deviations of 1.86K and 2.96K when compared against the experimental data for two thermocouple positions, however, it was approximately 2.5x10 5 times faster than CFD. The charging time deviation was less than 1% compared to CFD, and no significant change in thermal energy stored was observed since the temperature differences were small enough to result in insignificant sensible heat deviation. A simulation cost-effectiveness index (CEI) considering accuracy and computational affordability was introduced; the CEI of the RCM solver is 10 3 times better compared to CFD for the prediction of local PCM temperature.

25 ENERGY STORAGE↗
DOE OSTI2020

Lower length scale informed improvements to Bison U-Pu-Zr fuel swelling model

Due to renewed interest in metallic fuels from the U-(Pu)-Zr material system, significant improvements have recently been made to BISON's capability to simulate metallic fuel, particularly in the area of swelling. In this report, lower length scale simulations that have been conducted to inform the engineering-scale models in BISON during FY20 are described. For the high-temperature regions of the fuel that are dominated by the $\gamma$ phase of U-(Pu)-Zr, two new models for swelling have implemented in BISON. Phase-field simulations were conducted to calculate parameters for these and other BISON models that control when gaseous swelling ceases and fission gas release begins. For the low-temperature regions of the fuel dominated by the phase of uranium, microstructurally resolved simulations of pore growth and crystal deformation were performed to understand the effect of irradiation-generated crystal shape changes. Conformally meshed microstructures with pores were simulated with a crystal plasticity model and with a burnup eigenstrain model. Plastic deformation will not cause pore size to increase; burnup eigenstrain (irradiation-induced crystal shape changes) will increase the volume of existing porosity and cause it to become anisotropic in shape. Additional model development work for early-stage fission gas cluster/bubble formation in the $\gamma$ phase were performed, by coupling with defect clustering within cluster dynamics framework. Sensitivity analysis were conducted to identify the critical thermo-kinetic parameters needed for rigorous assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗