Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal simulation”

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 181 records · Page 10

Annual Report on NCSP Technical Support task in BNL during FY24

This work reports on the investigations done during FY 24 in collaboration with two 2024 Summer undergraduate interns through the DOE SULI and BNL SURP programs. These efforts were partially funded by the Nuclear Criticality Safety Program through the Technical Support Succession Plan task. The first project, developed by Ian Snider, focused on testing the impact of uncertainties in thermal cross sections for multiple materials in nuclear reactor simulations. The second project aimed implement concurrent into BNL’s machine-learning code to correct spin mis-assignments in neutron resonances, the Bayesian Resonance Reclassifier (BRR).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamic Pressure Drag on Rising Buoyant Thermals in a Neutrally Stable Environment

Abstract This study examines dynamic pressure drag on rising dry buoyant thermals. A theoretical expression for drag coefficient C d as a function of several other nondimensional parameters governing thermal dynamics is derived based on combining the thermal momentum budget with the similarity theory of Scorer. Using values for these nondimensional parameters from previous studies, the theory suggests drag on thermals is small relative to that on solid spheres in laminar or turbulent flow. Two sets of numerical simulations of thermals in an unstratified, neutrally stable environment using an LES configuration of the Cloud Model 1 (CM1) are analyzed. One set has a relatively low effective Reynolds number R e and the other has an order-of-magnitude-higher R e ; these produce laminar and turbulent resolved-scale flows, respectively. Consistent with the theoretical C d , the magnitude of drag is small in all simulations. However, whereas the laminar thermals have C d ≈ 0.01, the turbulent thermals have weakly negative drag ( C d ≈ −0.1). This difference is explained by the laminar thermals having near vertical symmetry but the turbulent thermals exhibiting considerable vertical asymmetry of their azimuthally averaged flows. In the laminar thermals, buoyancy rapidly becomes concentrated around the main centers of rotation located along the horizontal central axis, leading to expansion of thermals via baroclinic vorticity generation but doing little to break vertical symmetry of the flow. Vertical asymmetry of the azimuthally averaged flow of turbulent thermals is attributed mainly to small-scale resolved eddies that are concentrated in the upper part of the thermals.

Meteorology & Atmospheric Sciences↗

GOTHIC Aerosol Source Depletion Studies

Pacific Northwest National Laboratory has continued work to develop an aerosol-laden flow modeling capability with the Generation of Thermal Hydraulic Information in Containment (GOTHIC TM ) computer code to perform simulations for thermal hydraulic conditions and aerosol transport and deposition in spent fuel casks. This report describes our recent work to expand our model that was originally developed in 2019 to allow for thermal characterization, carrier gas flow determination, and tracking of particulate behavior throughout the entire canister volume. This was achieved through conversion of decay heat source in all fuel tubes within the model from heaters to thermal conductors in GOTHIC as well as remeshing the internal volume of the canister. These model improvements allowed for simulations of base case scenarios for comparison of similar efforts by our collaborators with alternative tools at Sandia and Oak Ridge. Several characteristics of the GOTHIC code were elucidated by this effort and were documented for consideration in how the code is used for this effort moving forward. For instance, treatment of particle size distributions, initial spatial distributions, mesh gradients, and minimum volume concentrations were identified as important issues for consideration. Ultimately, the current version of the code is capable of tracking temperatures, flow rates, and particle behavior throughout the canister internal volume. This report presents predicted depletion times and preferred deposition patterns within the canister. Additionally, this report documents code performance and insights for further model development and use for comparisons to other codes and for predicting system behavior in experiments being performed and planned by our collaborators.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Reactive Transport Simulations of High-Tempertature Geologic Thermal Energy Storage (GeoTES) in Deep Saline Formations - I/O Files

Simulation input and output files, post-processed figures and excel tables, and tecplot layout files for generating figures. These simulations were run with TOUGHREACT V4.12 by Lawrence Berkeley National Laboratory in 2021. This work was completed as part of the geologic thermal energy storage (GeoTES) research project reported in the final report for Phase I of this work, which is linked below.

15 GEOTHERMAL ENERGY↗

Concept Descriptions for the VTR Rabbit System and Driver Fuel Test Assemblies

Two of the experiment vehicles being developed for the Versatile Test Reactor (VTR) are presented here. The first is a rabbit system that will enable rapid insertion of small test capsules into the high fast flux of the VTR core for relatively short durations. The rabbit concept development includes the construction/demonstration of a near-full-scale system in a deep-water pool to demonstrate functionality, development of a concept of operations and initial procedures, and validation of thermal-hydraulic modeling. In addition, modeling efforts are underway to simulate the thermal and neutronic environment of a rabbit capsule. The second type of experiment vehicle presented here is a driver fuel test assembly for inserting fuel and materials tests into the core by replacing a driver fuel assembly. Here, a novel design for dismountable test assemblies is proposed for the VTR.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Review—Thermal Safety Management in Li-Ion Batteries: Current Issues and Perspectives

Approaches for thermal management of lithium-ion (Li-ion) batteries do not always keep pace with advances in energy storage and power delivering capabilities. Root-cause analysis and empirical evidence indicate that thermal runaway (TR) in cells and cell-to-cell thermal propagation are due to adverse changes in physical and chemical characteristics internal to the cell. However, industry widely uses battery management systems (BMS) originally designed for aqueous-based batteries to manage Li-ion batteries. Even the "best" BMS that monitor both voltage and outside-surface temperature of each cell are not capable of preventing TR or TR propagation, because voltage and surface-mounted temperature sensors do not track fast-emerging adverse events inside a cell. Most BMS typically include a few thermistors mounted on select cells to monitor their surface temperature. Technology to track intra-cell changes that are TR precursors is becoming available. Simultaneously, the complex pathways resulting in cell-to-cell TR propagation are being successfully modelled and mapped. Innovative solutions to prevent TR and thermal propagation are being advanced. These include modern BMS for rapid monitoring the internal health of each individual cell and physical as well as chemical methods to reduce the deleterious effects of rapid cell-to-cell heat and material transport in case of TR.

25 ENERGY STORAGE↗

Phase Field Modeling of Chemical Reaction Related Damage Evolution in Environmental Barrier Coatings

The advent of next-generation engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the thermally grown oxides (TGO) and phase transformation in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and porosity evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

fast-diffusion path↗

Project Report: High Frequency Axion LDRD - A Multiport Approach to Thermal Noise and Scattering Parameter Simulation of Cryogenic Experiments

In this report, a simple algorithm for detailed system-level thermal noise analysis is developed, demonstrated, and verified. This method uses noise-wave theory and noise covariance matrices to cascade noise and scattering parameters of multiport devices at different temperatures. This method addresses the effects of return loss, multiport isolation/coupling, and static temperature differentials between components, and will work in cases where the noise temperature is at or near the quantum noise limit. An ideal multidevice network will first be demonstrated to show that this method’s results are consistent with the Friis cascade when component parameters such as return loss and isolation are ideal. Following the ideal multidevice example, a cryogenic experiment is conducted to demonstrate that the proposed simulation method is successful when real data are used.

42 ENGINEERING↗

Challenges and approaches to interpretive modeling of boundary plasma and neutral transport in a closed, pumped divertor

An experimental discharge from the DIII-D tokamak is modeled using the SOLPS-ITER code suite and compared against measurements in the pumped and relatively closed upper divertor. Uncertainties of boundary plasma simulations are identified by attempting to match code inputs to experimental conditions, including iteratively solving transport coefficients to match upstream experimental profiles using varying quantities of core particle flux, different pumping models, and various assumptions of ion thermal transport. Simulated boundary conditions for particle injection at the core interface are shown to be relevant to the plasma solution at the divertor targets, even if upstream transport is modified so that plasma profiles are comparatively similar, although seperatrix density is not held constant. When upstream plasma profiles are matched to experimental measurements by varying diffusive transport coefficients, using either poloidally symmetric or ballooning structure, the model finds a majority of injected energy being transported radially off the computational domain, in conflict with experimental radiated power measurements and heat flux measurements at the divertor target. Imposing a maximum thermal diffusivity or radially shifting the experimental separatrix location of the fitted profiles to increase power conducted to the targets by increasing the upstream electron temperature does not significantly modify this result. Including a thermalizing plenum volume in the simulation domain is shown to maintain the experimental volumetric pumping rate without knowing the neutral energy distribution incident on the pump duct a priori. By modifying transport parameters to match different assumptions for ion temperature, downstream neutral pressure changes by more than a factor of two, suggesting that attention to ion thermal transport may be a critical parameter for simulations to accurately resolve recycling and neutral transport, particularly in a closed divertor geometry. In addition to quantifying various modeling uncertainties, this work motivates both further experimental study and modeling improvements to improve predictive capabilities.

divertor↗

Stabilizing NMC 811 Li-Ion Battery Cathode through a Rapid Coprecipitation Process

Improving the capacity of LiNi x Mn y Co 1-x-y O 2 (NMC) by increasing the Ni content results in materials with poor cycling performance and reduced thermal stability. These problems can be alleviated by creating concentration gradients in the materials, but this requires more expensive batch processing. Herein, a simple process is described that can create a coating that penetrates into the precursor particles, creating unique gradients in the precursor material. The coated material demonstrates improved cycling performance and thermal stability as compared to uncoated materials. Additionally, rapid heating to simulate a thermal runaway event while observing these materials with X-rays gives a deeper understanding of how the coating affects the decomposition process.

25 ENERGY STORAGE↗

Enabling Load Following Capability in the Transatomic Power MSR

This project is dedicated towards designing a fuel processing system that enables liquid-fueled molten salt reactors (MSR) to load follow by removing the dissolved xenon in the fuel salt. As one of the Gen-IV nuclear reactor concepts, the molten salt reactor receives increasing development interests in the recent years. One distinguishing feature of the liquid-fueled molten salt reactor is its improved ability to operate in a load-following mode by including the unique online fission product removal system. Load-following means that the reactor changes its power output based on the demand on the grid. Most of the current operating nuclear reactors have limited load-following ability and operate as the base load on the grid. Due to the rapid increase of solar energy, the requirement on load-following capacity is significantly increased because of the varying power output of the solar panels, yet the traditional load-following capacity is expected to decrease as the decarbonization of the grid continues. Therefore, the ability to perform load-following operation for the nuclear reactors will greatly enhance the resilience of the grid and make nuclear energy more economically competitive. This load-following feature is included in many commercial molten salt reactor designs, such as the designs by Transatomic Power, Terrestrial Energy, and Flibe Energy. Unfortunately, detailed analysis of the fuel processing system for commercial scale MSRs is still lacking, as well as how the fuel processing quantitively impacts the load-following operation. Moreover, experimental data for many of the underlying physics of fuel processing is limited. This project aims to pave the way for the fuel processing technology to advance to the commercial stage by performing combined experimental and simulation research. During the project period, four interconnected aspects of the development of the fuel processing system in liquid-fueled molten salt reactors are investigated. These aspects are the simulation and analysis of the fission product removal system, the fuel cycle simulation, the coupled reactor neutronics and thermal hydraulics transient simulation, and the gaseous fission product removal experiment. Multiphase CFD simulations are performed for components of the processing systems, and simplified air-water experiments are carried out to provide validation data. It is concluded that the CFD simulation can satisfactorily predict the system level performance of the components, and engineering models are constructed based on this success. Fuel cycle analysis is performed for two representative MSR design, the MSBR and the Transatomic Power MSR. Open-source code SaltProc is developed to incorporate the unique fuel processing system of the MSRs. It is concluded that the removal of xenon is essential for load-following operation in thermal spectrum MSR and Molten Salt Breeder Reactor. For the Transatomic Power MSR, the xenon poisoning effect is negligible due to its relatively fast neutron spectrum, though the overall fuel cycle economics still benefits from the removal of xenon. Coupled reactor neutronics and thermal hydraulics transient simulation is performed specifically for the Transatomic Power MSR. It is concluded that the reactor core design could perform power ramping fast enough to satisfy load-following operation. Combining the findings from each aspect, it is concluded that the load-following operation of a thermal neutron MSR is dependent upon the removal of xenon, which could be achieved for a commercial sized reactor using continuous inert gas sparging in a separate system with reasonable dimensions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analysis of Conduction Cooling Strategies for Wire Arc Additive Manufacturing

Metal additive manufacturing (AM) processing consists of numerous parameters which take time to optimize for various geometries. One aspect of the metal AM process that continues to be explored is the control of thermal energy accumulation during component manufacturing due to the melting and solidification of the feedstock. Excessive energy accumulation causes thermal failure of the component while minimal energy accumulation causes lack of fusion with the build plate or previous layer. The ability to simulate the thermal response of an AM component can increase research efficiency by reducing the time to optimize thermal energy accumulation. This paper presents an effective implementation of finite element analysis to determine the thermal response of a wire arc additive manufactured component with various build plate sizes and cooling methods including, integral build plate cooling, oversized build plates with passive cooling, and non-integral build plate cooling. The use of integral build plate cooling channels was shown to decrease the interpass temperature at the conclusion of the build process by 55% and build plate temperature by 96% compared to the conventionally deposited sample with 20 second dwell time. The use of a tall build plate with passive cooling was shown to reduce the interpass temperature by 32% as compared to the conventionally deposited sample with 20 second dwell time. Each cooling strategy evaluated decreased the interpass temperature within a range of 20–55% which enables higher deposition rates and decreased dwell times during depositions. The cooling strategies are designed to be implemented in a hybrid or retrofit AM platform to mitigate concerns of the thermal input from the additive process having detrimental effects on the precision of the machining process. This paper shows that accurate simulations of all strategies can be used to accurately predict the thermal response of the various strategies discussed. These cooling strategies will allow for increased deposition rates with comparable interpass temperature and decreased dwell time, increasing deposition efficiency. This model and these simulations are verified by experimental results. It is concluded that passive strategies, such as the over-sized tall build plate, can be used when liquid coolant in the AM environment could negatively affect the deposition process. Active cooling strategies, such as the integral build plate cooling could be used if low thermal conductivity materials are deposited or higher material deposition rates are desired. This paper discusses the use of active and passive cooling used during AM and shows how a simulation model can be used to make design choices for cooling strategies. The model also enables verification of select critical process parameters such as dwell times for a desired interpass temperature.

Heinrich, Lauren↗

Thermal conductivity of 3D-printed block-copolymer-inspired structures

This study primarily focuses on examining the impact that geometric structure has on thermal conductivity of multi-phase constructs in different 3D-printed poly(lactic acid), PLA, samples. The investigated structures are inspired by morphologies formed by diblock copolymers: lamellae, hexagonally packed cylinders, and gyroid. This research also investigates how volume percentage and material combination influence the thermal conductivity of these structures. Further, the samples can be tailored to simulate various thermal management structures observed in practical applications, such as thermal interface materials in electronic devices. Thermal conductivity ratio is controlled using air, the least conductive material at 0.026 W/(m K), PLA at 0.136 W/(m K), and thermal paste at 5.11 W/(m K). Different models were tested against thermal conductivity measurements in order to capture the effect of material type (PLA-Air versus PLA-Thermal Paste), volume percentage, structure, and orientation. Simple, effective medium models were good predictions of thermal conductivity in lamellar structures, but it was necessary to develop models for conduction through cylindrical and gyroid structures. Finally, all results were normalized to find a universal model that is independent of structure and material. This approach provides a simple method to predict how to reduce or enhance transport properties and heat management capabilities of 3D printed objects.

36 MATERIALS SCIENCE↗

Cabin Thermal Management Analysis for SuperTruck II Next-Generation Hybrid Electric Truck Design

In this article, we present a multistage, coupled thermal management simulation approach, informed by physical testing where available, to aid design decisions for PACCAR's SuperTruck II hybrid truck cabin concept. Focus areas include cabin insulation, battery sizing, and sleeper curtain position, as well as heating, ventilating, and air-conditioning (HVAC) component and accessory configurations, to maintain or improve thermal comfort while saving energy. The authors analyzed weather data and determined the national vehicle miles traveled weighted temperature and solar conditions for long-haul trucks. Example weather day profiles were selected to approximate the 5th and 95th percentile weighted conditions. A daylong drive cycle was developed to impose appropriate external wind conditions during rest and driving periods. Using the National Renewable Energy Laboratory's vehicle HVAC modeling and simulation tool VTCab, HVAC load design trade-off studies for the new truck geometry concept were completed. Parameters analyzed included effects of paint color, insulation, glass transmissivity, and curtain location. Simulation results helped with early design material selections for efficient cabin climate control. A detailed three-dimensional computer-aided engineering (CAE), computational fluid dynamics (CFD), radiation, and human physiology co-simulation, referred to in this article as CAE Thermal-CFD, was used to evaluate thermal comfort and energy impacts of diffuser configurations and air supply settings in driving and hoteling modes. Analysis revealed that it is more difficult to heat the cabin in hoteling mode during the winter than to cool the space in the summer. This seasonal load profile drives the requirement of additional energy storage for heating comfort. To determine the battery capacity requirement, multiday HVAC operation drive cycle simulations were then completed, showing that a 15-kWh battery would be required for HVAC operation during hoteling. Results helped reduce cabin thermal loads, determine component sizing requirements, and improve occupant comfort to save fuel and contribute to the economic viability of the hybrid system.

33 ADVANCED PROPULSION SYSTEMS↗

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this study, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E x B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical ion thermal transport simulated by NEO, account for 30-60% of the total experimental heat flux. In addition, the neoclassical electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne /L Te . Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for neoclassical particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of neoclassical particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. In conclusion, the implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-stage preconditioners for thermal–compositional–reactive flow in porous media

We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe conver-gence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Péclet number) to diffusion dominated (low Péclet number). The number of linear iterations is reduced by 40–85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.

42 ENGINEERING↗

Thermal-Hydraulic-Mechanical Modeling and Simulation of Sodium-Potassium–Cooled MARVEL Microreactor Core

The U.S. Department of Energy's Microreactor Program, with Idaho National Laboratory's development of a nuclear microreactor applications test bed named MARVEL, aims to support R&D for the deployment of small, transportable reactors across civilian, industrial, and defense sectors. The MARVEL microreactor, an 85-kWth thermal fission reactor, incorporates TRIGA nuclear fuel and a sodium-potassium eutectic as its primary coolant, designed for safety and efficiency, with natural circulation eliminating the risk of critical heat flux conditions. The reliance on natural circulation for primary cooling means the reactor avoids using fuel spacers to minimize core pressure drop, which could disrupt the primary coolant's natural flow. However, the reactor core’s tight P/D ratio of 1.05, in the absence of fuel spacers, could pose a risk of fuel rod contact and increased peak cladding temperatures. To ensure the reactor safety, this study conducted computational modeling and simulations to investigate the reactor's thermal-hydraulic-mechanical characteristics, including the reactor core heat transfer coefficients, the potential for rod-to-rod contact, and assessed its impact on peak cladding temperature and overall reactor safety. The computational analyses of the MARVEL microreactor core revealed that the thermal deformation of fuel rods under worst-case scenario may lead to the fuel rod contact, but the peak cladding temperatures will remain significantly lower than the safety criteria, ensuring the safety operation reactor without fuel spacers under normal operating conditions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗