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 325 records · Page 18

Electronic density response of warm dense matter

Matter at extreme temperatures and pressures—commonly known as warm dense matter (WDM)—is ubiquitous throughout our Universe and occurs in astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion and is realized in the laboratory using different techniques. A particularly important property for the understanding of WDM is given by its electronic density response to an external perturbation. Such response properties are probed in x-ray Thomson scattering (XRTS) experiments and are central for the theoretical description of WDM. In this work, we give an overview of a number of recent developments in this field. To this end, we summarize the relevant theoretical background, covering the regime of linear response theory and nonlinear effects, the fully dynamic response and its static, time-independent limit, and the connection between density response properties and imaginary-time correlation functions (ITCF). In addition, we introduce the most important numerical simulation techniques, including path-integral Monte Carlo simulations and different thermal density functional theory (DFT) approaches. From a practical perspective, we present a variety of simulation results for different density response properties, covering the archetypal model of the uniform electron gas and realistic WDM systems such as hydrogen. Moreover, we show how the concept of ITCFs can be used to infer the temperature from XRTS measurements of arbitrary complex systems without the need for any models or approximations. Finally, we outline a strategy for future developments based on the close interplay between simulations and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Role of dislocations on martensitic transformation temperatures and microstructure: A molecular dynamics study

Microstructure and defects strongly affect martensitic transformations in metallic alloys. Significant progress has been made in understanding the atomic-level processes that control the role of grain boundaries and precipitates in these solid-to-solid phase transformations. Yet, the role of dislocations and their structures on martensitic transformation temperature and the resulting microstructure remains unclear. Therefore, we used large-scale molecular dynamics simulations to study the forward and reverse transformation of a martensitic material modeled after Ni63Al37 under cyclic thermal loading. The simulations reveal that dislocations in the austenite phase act as one-dimensional seeds for the martensite phase, which is present at temperatures significantly above the martensite start value. We find a reduction in the dislocation density during cyclic thermal loading, which results in the increase in martensite and austenite transition temperatures, in agreement with experiments. Importantly, we extracted a critical martensitic nuclei size for developing stable domains and found that relatively low dislocation densities are needed to grow independent martensitic variants resulting in a multi-domain structure.

Physics↗

Benchmark of Neutron Thermalization in Graphite Using a Pulsed Slowing-Down-Time Experiment

A benchmark has been developed using a pulsed slowing-down-time experiment to isolate the thermalization process in graphite. The experiment was conducted at the Oak Ridge Electron Linear Accelerator facility at Oak Ridge National Laboratory, and it measured the time spectrum of neutrons leaking from a graphite pile during slowing down and thermalization within graphite. Simulations of the benchmark experiment were performed using the MCNP6.1 Monte Carlo code and the ENDF/B-VII.1 and ENDF/B-VIII.0 cross-section databases. The benchmark provides a time spectrum (i.e., time-dependent counts in a detector) that allows for validation of the graphite thermal scattering libraries (TSLs). The impact on the simulations using a suite of graphite TSLs was compared with the experimental results. Given the density of nuclear graphite, the TSL corresponding to graphite with 30% porosity, as implemented in ENDF/B-VIII.0, was found to most accurately represent the measured time spectrum corresponding to the thermal energy range with an average deviation of ±1.7%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ocean cavity regime shift reversed West Antarctic grounding line retreat in the late Holocene

Recent geologic and modeled evidence suggests that the grounding line of the Siple Coast of the West Antarctic Ice Sheet (WAIS) retreated hundreds of kilometers beyond its present position in the middle to late Holocene and readvanced within the past 1.7 ka. This grounding line reversal has been attributed to both changing rates of isostatic rebound and regional climate change. Here, we test these two hypotheses using a proxy-informed ensemble of ice sheet model simulations with varying ocean thermal forcing, global glacioisostatic adjustment (GIA) model simulations, and coupled ice sheet-GIA simulations that consider the interactions between these processes. Our results indicate that a warm to cold ocean cavity regime shift is the most likely cause of this grounding line reversal, but that GIA influences the rate of ice sheet response to oceanic changes. This implies that the grounding line here is sensitive to future changes in sub-ice shelf ocean circulation.

58 GEOSCIENCES↗

Engineering Thermal Transport across Layered Graphene–MoS 2 Superlattices

Layering two-dimensional van der Waals materials provides a high degree of control over atomic placement, which could enable tailoring of vibrational spectra and heat flow at the sub-nanometer scale. In this work, using spatially resolved ultrafast thermoreflectance and spectroscopy, we uncover the design rules governing cross-plane heat transport in superlattices assembled from monolayers of graphene (G) and MoS 2 (M). Using a combinatorial experimental approach, we probe nine different stacking sequences, G, GG, MG, GGG, GMG, GGMG, GMGG, GMMG, and GMGMG, and identify the effects of vibrational mismatch, interlayer adhesion, and junction asymmetry on thermal transport. Pure G sequences display evidence of quasi-ballistic transport, whereas adding even a single M layer strongly disrupts heat conduction. The experimental data are described well by molecular dynamics simulations, which include thermal expansion, accounting for the effect of finite temperature on the interlayer spacing. The simulations show that an increase of ~2.4% in the layer separation of GMGMG, relative to its value at 300 K, can lead to a doubling of the thermal resistance. Using these design rules, we experimentally demonstrate a five-layer GMGMG superlattice “thermal metamaterial” with an ultralow effective cross-plane thermal conductivity comparable to that of air.

2D materials↗

Thermal Modelling of Advanced Test Reactor Fuel in a Generalized Dry Storage System

Star-CCM+, a computational fluid dynamics (CFD) software was used to conduct modeling and simulation of the thermal performance of a dry storage configuration consisting of Department of Energy Standardized Canisters (DOESCs) loaded with aluminum-clad spent nuclear fuel (ASNF.) The configuration includes nine DOESCs loaded with Advanced Test Reactor ASNF contained within a stainless steel overcanister centered in a ventilated, concrete overpack. The simulations were used to estimate the maximum temperatures reached by backfill gases inside the overcanister and DOESCs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Building envelope anomaly characterization and simulation using drone time-lapse thermography

Defects in building envelopes deteriorate over time without being visible to the human eye, while significantly impacting energy performance due to unaccounted heat transfer. Defects can be characterized in the infrared (IR) spectrum. However, IR readings are typically recorded at singular points in time, when in several cases anomalies can only be revealed at specific times of the day, possibly in different seasons of the year. This paper presents a novel workflow for 3D envelope defect characterization and modeling using aerial time-lapse IR data collection using drones. A comprehensive envelope thermal profile is developed for a case study building employing the photogrammetry software Agisoft Photoscan, which generates temporal IR inspections of building skins using multiple thermography orthomosaics. Point-cloud data is then translated into a CAD model and thermal zones for whole Building Energy Modeling (BEM) using Honeybee as a frontend to EnergyPlus to showcase the potential of inclusion of detailed 4D data. Envelope contributions in this case study’s anomalies showed heat losses of 6447.6 kWh, and Energy Use Intensity (EUI) differences of ~2 kWh/m 2 /year from the baseline. Finally, why there is currently little translation of this work in BEM software is discussed, while identifying limitations and future research in the employment of time-lapse thermography using drones for more accurate building envelope inspection and modeling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Computational Fluid Dynamics Modeling of Solar Thermal Dry Reforming of Methane in a Parabolic Trough

Computational fluid dynamics simulations of solar-thermal dry reforming of methane using a parabolic trough configuration were performed. Parametric simulations of different combinations of gas flow rate, receiver tube emissivity, and geometric concentration ratio were conducted to determine configurations that could achieve the required catalyst temperatures of at least 700 °C to achieve high conversion of CH4 and CO2 to H2 and CO. Results showed that the concentration ratio of the parabolic trough collector had to be increased from ~70 to ~120 and the receiver-tube emissivity had to be reduced to ~0.2 to achieve bulk average catalyst temperatures of greater than 700 °C. Lower gas flow rates also reduced enthalpic heat losses and increased catalyst temperatures.

Ho, Clifford↗

Simulations of the fault current limiting operation of a long-length REBCO CORC ® superconducting cable cooled by helium gas

Conductor-on-round-core (CORC®) cables composed of rare-earth-barium-copper-oxide high-temperature superconducting (HTS) tapes are of great interest for power transmission applications due to their many advantages such as high power density, light weight, and low loss. Closed circulation loops of cryogenic helium gas can be used to cool HTS cables down to low temperatures to significantly improve their current-carrying capacity. Coupled circuit-electromagnetic-thermal finite element simulations implemented in the COMSOL Multiphysics package were developed, validated, and then used for simulating the fault current limiting (FCL) performance and the cooling processes of an 8-layer CORC® cable cooled with a flow of cryogenic helium gas. In the simulations, the temperature dependence of the electrical and thermal properties of all component materials is implemented for improved accuracy. To overcome computational challenges caused by the considerable difference in geometrical scales (i.e. few-µm-thick HTS layers versus 10 m-long HTS cable), the model is divided into two separate simulations. The first simulation is performed on the transverse cross-section of the cable to calculate the electric field, heating power and temperature rise in each component of a CORC® cable during FCL operation. The heating power calculated in the first simulation is transferred to the second model to simulate the cooling of a 10 m-long cable after the fault is cleared. The effect of the helium gas flow rate on the cooling process is also investigated to develop strategic approaches for optimizing cooling systems for HTS cables with FCL capability. The simulations indicated that a 40 ms fault with a voltage drop of 20 V m −1 along the cable can result in a temperature increase from 60 K to about 165 K inside the cable, and it takes about 500 s to cool the cable back to nearly 60 K with a flow of cold helium gas at a rate of 5 g s −1 .

24 POWER TRANSMISSION AND DISTRIBUTION↗

Design, fabrication, and testing of a water-cooled tungsten-foam high-heat flux plasma-facing module

Here, we present a comprehensive study detailing the design, fabrication and testing a new modular “finger” component, featuring a hexagonal tungsten foam armor tile bonded to the top surface of the water-cooled finger unit. The mechanical, thermal, and particle transport properties of the tungsten foam armor were determined for a wide range of porosity and cell size, and then combined in an equivalent solid material with “effective” properties suitable for large-scale simulations. A series of 3-D multiphysics simulations for coupled thermal and mechanical behavior of the “finger” module are used to optimize the design before fabrication. The simulation results show that, for the expected experimental conditions, a safe range of heat flux and foam porosity is possible for both water-cooling and helium-cooling. Based on these multiphysics simulations, an optimum finger cooling module has been fabricated at Ultramet, Inc. A series of experiments were designed to study the thermomechanical performance of the “finger” unit under steady-state and transient thermal loading at surface heat flux up to 16 MW/m 2 in the High Energy Flux Test facility (HEFTY) at UCLA. Testing conditions included low-power, steady-state plasma exposure at 6 MW/m 2 , followed by severe high heat flux cycles at 16 MW/m 2 . Results of these tests are explained through comparisons with multiphysics simulations. To eliminate the effects of thermomechanical damage carried out on a module on its thermo-mechanical response on a subsequent test, three different yet identically designed and fabricated finger modules were tested. The results show that the current design is capable of cooling the walls of fusion devices at surface heat flux up to 6 MW/m 2 without failure, and that the failure heat flux limit is below 16 MW/m 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Simulation Framework for Analyzing the Impact of Stochastic Occupant Behaviors on Demand Flexibility in Typical Commercial Buildings

As one of the primary users of the electric grid, buildings and building equipment, including heating, ventilation, and air conditioning (HVAC) systems, can be leveraged to provide the flexible demand needed to balance the grid. Typical strategies to achieve demand flexibility are to reduce electricity use during peak or critical periods by shutting down equipment or relaxing system setpoints, which will inevitably impact the occupants’ comfort. When occupants feel uncomfortable, they may take actions to regain their comfort, and some of those actions (such as turning on a personal fan) may have a negative impact on meeting the demand response goal. Therefore, it is important to incorporate occupant behaviors into the assessment ofthe building demand flexibility potential. In this study, a simulation framework that includes simulation of zone thermal loads, an HVAC system, and occupant behaviors, was developed to investigate the impact of occupant behaviors on demand flexibility. A case study was conducted using a small office model from the U.S. Department of Energy (DOE) Commercial Prototype Building Models to simulate the building envelope and zone loads. An agent-based occupant thermal behavior model was adapted to forecast occupants’ thermal comfort and their resulting thermal behaviors. An artificial neural network (ANN) based airflow model trained from a computational fluid dynamics (CFD) model of the zone was adopted to better predict the ambient environment of each occupant. An air-source heat pump simulation model that was calibrated from a real two-stage air-source heat pump system was used as the HVAC system. A typical load shedding event during peak hours was studied. Repeated simulations were conducted to capture the stochastic effects of occupant behaviors. The interplay between the demand flexibility, occupant comfort and behavior were analyzed by evaluating key performance indicators, including the energy use, occupant discomfort duration, and occupant behavior duration during the peak period. The results suggest that this framework can be used to analyze typical commercial buildings and their HVAC systems in terms ofdemand flexibility potential under the impact of occupant behaviors.

Chen, Zhelun↗

Simulations of Nb$_3$Sn Layer RF Field Limits Due to Thermal Impedance

Nb3Sn performance in RF fields is limited to fields far below its superheating critical field. Poor thermal conductivity of Nb3Sn has been speculated to be the reason behind this limit. In order to better understand the contribution of Nb3Sn thermal conductivity to its RF performance, we simulated numerically with Matlab program*(based on SRIMP and HEAT codes) thelimiting fields under different realistic conditions. Our simulations indicate that limiting fields observed presently in the experiments with RF fields cannot be explained by the thermal impedance of Nb3Sn alone. The results change significantly in the presence of higher losses due to extrinsic mechanisms.

43 PARTICLE ACCELERATORS↗

Defects vibrations engineering for enhancing interfacial thermal transport in polymer composites

To push upper boundaries of thermal conductivity in polymer composites, understanding of thermal transport mechanisms is crucial. Despite extensive simulations, systematic experimental investigation on thermal transport in polymer composites is limited. To better understand thermal transport processes, we design polymer composites with perfect fillers (graphite) and defective fillers (graphite oxide), using polyvinyl alcohol (PVA) as a matrix model. Measured thermal conductivities of ~1.38 ± 0.22 W m -1 K -1 in PVA/defective filler composites is higher than those of ~0.86 ± 0.21 W m -1 K -1 in PVA/perfect filler composites, while measured thermal conductivities in defective fillers are lower than those of perfect fillers. We identify how thermal transport occurs across heterogeneous interfaces. Thermal transport measurements, neutron scattering, quantum mechanical modeling, and molecular dynamics simulations reveal that vibrational coupling between PVA and defective fillers at PVA/filler interfaces enhances thermal conductivity, suggesting that defects in polymer composites improve thermal transport by promoting this vibrational coupling.

42 ENGINEERING↗

Simulations of SITI Cookoff Experiments Carried Out with Different Lots of PBX 9502

A pressure dependent cookoff model for PBX 9502 was developed by Hobbs’ et. al. PBX 9502 is composed of 95% by mass triaminotrinitrobenzene (TATB) and a 5% by mass chlorotrifluoroethylene/vinylidine fluoride binder. The objective in this study is to implement this cookoff model in Aria to simulate Sandia Instrumented Thermal Ignition (SITI) experiments that were carried out with different manufacturing lots of PBX 9502. The SITI design consists of solid cylinders (1" diameter × 1" height) of insensitive high explosive (IHE) confined by a cylindrical aluminum case. An electric heater is wrapped around the lateral surface of the case. This heater produces a temperature heating ramp on the outer surface of the case. Internal thermocouples measure the IHE temperature rise from the center to locations close to the IHE-aluminum interface. The energetic material is heated until thermal ignition occurs. Pressure is measured with a static pressure transducer installed on top of the confinement case. Two–dimensional axisymmetric heat conduction finite element models were implemented to simulate these experiments using four options of the PBX 9502 cookoff model. In addition, the predictive ability of this thermal decomposition model is evaluated using Latin Hypercube Sampling (LHS) techniques.

42 ENGINEERING↗

Effect of Internal Structure on Warpage in a Large-Scale Additive Manufacturing Process With Bio-Derived Composites

Extrusion-based large-format additive manufacturing (LFAM) often results in unintended deformation and failures due to thermal residual stress between layers. This LFAM involves a hot molten polymer deposited on a previously deposited and cooled layer, generating a temperature mismatch between the layers. The temperature discrepancy causes a thermal contraction mismatch, generating thermal residual stress. Due to the residual stress, printed structures experience warpage and delamination Therefore, understanding the heat distribution behavior is essential to preventing undesired deformation and failures. The goal of this study is to investigate the effect of infill patterns on heat distribution and deformation in an LFAM system. Wood fiber-reinforced polylactic acid was used as a sustainable and low-cost reinforcement material. A numerical model was developed to predict the temperature and deformation field with thermo-mechanical properties, which were measured by a dynamic mechanical analysis and a thermal expansion test. Simulations were performed on a box geometry with three different infill patterns. The simulation model was verified with temperature data gathered from an infrared camera. The results show good agreement between predicted and measured temperature profiles. The developed simulation model was applied to a roof tray with different infill patterns for a case study. We found that the different thermal mass distributions resulting from the various infill patterns affected the heat distribution and deformation. These findings contribute to a better understanding of the thermal and mechanical behavior of LFAM.

Jo, Eonyeon↗

Effect of resistive wall on thermal quench in JET disruptions

In this paper, experimental data, simulations and theory are presented of a JET tokamak thermal quench. The emphasis is on the timescale of the bulk plasma thermal energy loss. The simulations suggest that the thermal energy loss is caused by a resistive wall tearing mode, and experimental data is consistent with this conclusion. The time scale of the thermal quench is the inverse of the mode growth rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A mesoscopic digital twin that bridges length and time scales for control of additively manufactured metal microstructures

We present our recent development of an integrated mesoscale digital twin (DT) framework for relating processing conditions, microstructures, and mechanical responses of additively manufactured (AM) metals. In particular, focusing on the laser powder bed fusion technique, we describe how individual modeling and simulation capabilities are coupled to investigate and control AM microstructural features at multiple length and time scales. We review our prior case studies that demonstrate the integrated modeling schemes, in which high-fidelity melt pool dynamics simulations provide accurate local thermal profiles and histories to subsequent AM microstructure simulations. We also report our new mechanical response modeling results for predicted AM microstructures. In addition, we illustrate how our DT framework has been validated through modeling–experiment integration, as well as how it has been practically utilized to guide and analyze AM experiments. Finally, we share our perspectives on future directions of further development of the DT framework for more efficient, accurate predictions and wider ranges of applications.

36 MATERIALS SCIENCE↗

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage (CRADA Report)

NREL and Element 16 collaborated on sulfur thermal energy storage modeling using NREL’s high performance computing (HPC) resources to assist its application in industrial processes. Industrial process heat (IPH) accounts for ~70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately, 1500 TWht (~60% Terawatt hour thermal) of IPH demand is in the temperature range of 100-300°C. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. This HPC4EI project facilitated Element 16’s development of low-cost and high-impact molten sulfur TES for dispatchable IPH. The development of a high-fidelity model validated by experimental data and HPC simulations enabled the successful resolution of the complex interplay between fluid dynamics and heat transfer processes during transient operation of sulfur TES, overcoming the numerical challenges posed by the non-linear temperature-dependent physical properties of sulfur. The project helped accelerate Element 16’s molten sulfur TES product design and support its broad applications.

25 ENERGY STORAGE↗