SEARCH · Search NASA
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.
❑ Introduction ❑COMSOL thermal simulations for present conditions - Results - Validation ❑COMSOL thermal simulations for PIP-II conditions ❑Proposed changes - New busbars section - Results ❑Peak temperature tables ❑Conclusions
This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.
Increase of primary beam power for neutrino beam-lines leads to a reduced lifespan for production targets. New concepts for robust targets are emerging from the field of High Power Targetry (HPT); one idea being investigated by the HPT R&D Group at Fermilab is an electrospun nanofiber target. As part of their evaluation, samples with different densities were sent to the HiRadMat facility at CERN for thermal shock tests. The samples with the higher density, irradiated under a high intensity beam pulse, exhibit major damage at the impact site whereas those with the lower density show no apparent damage. The exact cause of this failure was unclear at the time. In this paper, we present the results of multiphysics simulations of the thermal shock experienced by the nanofiber targets that suggest the failure originates from the reduced permeability of the high density sample to air flow. The air present in the porous target expands due to heating from the beam, but is unable to flow freely in the high density sample, resulting in a larger back pressure that blows apart the nanofiber mat. We close with a discussion on how to further validate this hypothesis.
Increase of primary beam power for neutrino beam-lines leads to a reduced lifespan for production targets. The field of High Power Targetry (HPT) is generating new concepts to meet the need for robust targets. One idea being investigated by the HPT Research and Development Group at Fermilab is an electrospun nanofiber target. As part of their evaluation, samples with different densities were sent to the HiRadMat facility at CERN for thermal shock tests. The samples with the higher density, irradiated under a high intensity beam pulse, exhibit major damages at the impact site whereas those with the lower density show no apparent damages. The exact cause of this failure was unclear at the time. In this paper, we present the results of multiphysics simulations of the thermal shock experienced by the nanofiber targets that suggest the failure originates from the reduced permeability of the high density sample to airflow. The air present in the porous target expands due to heating from the beam, but is unable to flow freely in the high density sample, resulting in a larger back pressure that blows apart the mat. We close with a discussion on how to further validate this hypothesis.
Three-dimensional (3D) characterization provides opportunities for understanding processing-structure relationships in additively manufactured (AM) materials. Bulk samples of Inconel 718 were fabricated via electron beam melting (EBM) in order to study microstructural development as a function of energy input and beam scan strategy. TriBeam tomography of bulk Inconel 718 microstructures built under steady-state growth conditions reveals the sensitivity of microstructure formation and evolution to machine process parameters. Here, samples manufactured using a narrow range of energy input per unit build area result in varied grain morphologies and crystallographic textures. Using TRUCHAS, a thermal simulation software, the thermal history of bulk scan strategies was predicted, and combined with a calibrated microstructure-processing map to accurately predict bulk grain morphologies. The solidification parameters and the 3D measured nucleation density are used to predict the transition between columnar and equiaxed grain morphologies, providing a process map to guide AM parameter choices to locally control as-printed microstructure. A two-dimensional metric for characterizing bulk grain morphology was also found to agree well with predictions from the process map calibrated by 3D data. Combined with 3D tomography and thermal modelling, the physics of structure development were understood at a new level of detail with respect to the competing processes of grain nucleation and epitaxial growth.
The purpose of this paper is to accurately evaluate the performance and overall efficiency of an electrical machine driving a propeller (as a part of a hybrid or electric propulsion system) during the flight mission profile using electromagnetic-thermal simulations. The study focuses on a 250 kW, 5000 rpm, 1.25 kHz, 36/30 slot/pole surface permanent magnet (SPM) machine with additively manufactured (AM) hollow conductors and integrated heat pipes (HP). The machine drive parameters (frequency, and current) are tuned to maximize the efficiency for a given flight cycle, while providing the highest specific power (SP) during takeoff. Afterward using Multiphysics simulations, the transient thermal behavior of the machine is studied and temperature profile for different parts of the machine is extracted. Finally, steady-state temperatures during cruising are estimated and overall efficiency is re-evaluated.
Nano-size confinement induces many intriguing non-Fourier heat conduction phenomena, such as nonlinear temperature gradients, temperature jumps near the contacts, and size-dependent thermal conductivity. Over the past decades, these effects have been studied and interpreted by nonequilibrium molecular dynamics (NEMD) and phonon Boltzmann transport equation (BTE) simulations separately, but no theory that unifies these two methods has ever been established. In this work, we unify these methods using a quantitative mode-level comparison and demonstrate that they are equivalent for various thermostats. We show that different thermostats result in different non-Fourier thermal transport characteristics due to the different mode-level phonon excitations inside the thermostats, which explains the different size-dependent thermal conductivities calculated using different reservoirs, even though they give the same bulk thermal conductivity. Specifically, the Langevin thermostat behaves like a thermalizing boundary in phonon BTE and provides mode-level thermal-equilibrium phonon outlets, while the Nose-Hoover chain thermostat and velocity rescaling method behave like biased reservoirs, which provide a spatially uniform heat generation and mode-level nonequilibrium phonon outlets. These findings explain why different experimental measurement methods can yield different size-dependent thermal conductivity. They also indicate that the thermal conductivity of materials can be tuned for various applications by specifically designing thermostats. Finally, the unification of NEMD and phonon BTE will largely facilitate the study of thermal transport in complex systems in the future by, e.g., replacing computationally unaffordable first-principles NEMD simulations with computationally less expensive spectral BTE simulations.
A round bobbin test has been widely used to measure critical current Ic of practical superconductors including Nb–Ti, Nb3Sn, and Bi-2212 round strands at 4.2 K and in high magnetic fields but rarely used for rare earth-barium-copper-oxide (REBCO) coated conductor tapes. Here, we applies this method to REBCO tapes and test their Ic at 77 K, self-field, and 4.2 K with a magnetic field of 14 T applied parallel to the tape. The sample carries very high current densities; it is 1.3 m long and has a dense array of voltage taps that allows probing localized thermal runaways and correlation with Ic variations along the length. At 77 K, the Ic values of all sections were measured and found to be rather uniform. However, at 4.2 K, we found that a 3 cm section has 90% of the overall voltage drop and it went into thermal runaways upon further increasing current. Thus, at 4.2 K, we were prevented by localized thermal runaways and unable to determine Ic of most of the REBCO conductor sections; this behavior is different from Bi-2212 and Bi-2223 for which Ic of all turns of a round barrel sample can be determined. We simulate REBCO tapes during current ramping to a thermal runaway and analyze impacts of various factors including cooling, conductor n-values, and thickness of stabilizer on the thermal runaway current Iq. The combined experiments and simulations identify n-value as a crucial factor that leads to the observed differences between Bi-2212 and REBCO. Together with cooling and the standard deviation of Ic variations along the length, the n-value determines whether Iq is above or below the average Ic of a long piece of conductor and their ratio Iq/Ic,ave, greatly affecting quench detection and thermal runaway prevention strategies of Bi-2212 and REBCO magnets.
Developing new materials to be applied in extreme environments is an opportunity and a challenge for the future. High entropy alloys are new materials that seem promising approaches to work in nuclear fusion reactors. In this work, FeTaTiVW high entropy alloys were developed and characterized with Molecular Dynamic and Hybrid Molecular Dynamic Monte Carlo simulations. The simulation results show that phase separation originates a lower potential energy per atom and a high level of segregation compared to those of a uniform solid solution. Moreover, the experimental diffractogram of the milled powder shows the formation of a body-centred cubic-type structure and the presence of TiO2. In addition, the microstructure of the consolidated material evidenced three phases: W-rich, Ti-rich, and a phase with all the elements. This phase separation observed in the microstructure agrees with the Hybrid Molecular Dynamic Monte Carlo simulation. Moreover, the consolidated material’s thermal conductivity and specific heat are almost constant from 25 °C to 1000 °C, and linear expansion increases with increasing temperature. On the other hand, specific heat and thermal expansion values are in between CuCrZr and W values (materials chosen for the reactor walls). The FeTaTiVW high entropy alloy evidences a ductile behaviour at 1000 °C. Therefore, the promising thermal properties of this system can be attributed to the multiple phases and systems with different compositions of the same elements, which is exciting for future developments.
The ability to do thermal simulations for entire additive manufacturing builds is a key computational problem facing the additive manufacturing community; however, complex numerical models considering multiple physical phenomena currently do not have the capacity for simulations at this scale. To this end, conduction only analytic models offer a viable approach due to the massive drop in computational expense. In this work, we extend an existing implementation which uses a governing equation which can be evaluated at any point in space and time. This implementation already utilizes OpenMP with a spatial decompositions scheme stemming from a melt pool tracking algorithm. Furthermore, we then combine this with a parallel in time (PinT) approach to make the problem highly parallelizable. The new scheme, which uses MPI for internode communication and OpenMP for intranode communication, is shown to scale very well across multiple computational nodes. This approach results in the ability to simulate the 3D solidification conditions for entire layers of additively manufactured parts in minutes making part scale thermal simulations more practical.
Dynamic Energy Transport and Integration Lab (DETAIL) at Idaho National Laboratory is to support experimental demonstration and validation research on Nuclear-Renewable Hybrid Energy System [1]. The Thermal Energy Distribution System (TEDS) is a thermal-hydraulic flow loop in DETAIL with its own dedicated control system to support the integration of co-located multiple experimental systems, where a packed-bed thermal energy storage (TES) is installed as a thermal buffer and storage unit. Among various TES options, the packed-bed TES is adopted in TEDS because it offers a low-cost single-tank thermal storage option compared to the traditional two-tank TES. However, since the TES tank is filled with granular fillers having different thermophysical properties from those of TES tank wall, there is a potential thermo-mechanical issue to be carefully addressed like thermal ratcheting which may pose a significant design concern for the packed-bed TES tank. Thermal ratcheting is a thermomechanical process caused by the repeated rearrangement of granular filler inside a TES tank during continuous thermal cycling operation of the packed-bed TES system. If the thermally induced stress exceeds the yield strength of a TES tank wall, catastrophic consequences may happen like rupture of the TES tank. Thus, it is important to understand the phenomenon to ensure the robust operation of the packed-bed TES tank. Given that thermal ratcheting is caused by complex interaction of thermal transport in the porous bed and solid mechanics, the accurate prediction of transient thermal behavior of the packed-bed TES tank, which is the focus of this paper, is critical to the reliable thermal ratcheting analysis. This paper discusses the numerical modeling, simulation, and validation studies that are ongoing at INL to investigate the transient thermal behavior of the packed-bed TES. Of particular concern is the transient thermal process occurring in the packed-bed TES unit that is operated in conjunction with the INL TEDS. The main goal of this research is three-fold: (i) provide preliminary insights into the transient thermal behavior of the TEDS TES tank, (ii) support the thermal measurement and validation plan for TEDS experiment, and (iii) provide transient thermal boundary conditions to support the reliable thermal ratcheting analysis of the TEDS TES tank. For the transient thermal modeling and analysis, a CFD model was developed, and the validity of the modeling approach was examined via comparing the numerical simulation results with the experimental data obtained from various design characteristics of packed-bed TES tanks. Then, the present modeling method was applied for the transient thermal analysis of the TEDS TES tank, and the results are discussed along with the potential improvement of data acquisition strategy for the future TEDS experiments for more precise validation study.
The fluoride-salt-cooled high-temperature reactor (FHR) is one type of the advanced reactors and has been attracting great interest from the research institutes and commercial companies in the recent years. However, currently the technology is relatively immature. To facilitate the design and safety analysis of FHRs, Kairos Power published the generic FHR (gFHR) benchmark. In this paper, a hybrid method combining the stochastic code and deterministic code is developed to simulate the gFHR benchmark. Serpent 2 is employed as a few-group cross section generator and the cross sections are applied to the finite difference neutron diffusion code AGREE. A good consistency with the gFHR benchmark is achieved. The agreement between the Serpent 2 results and the AGREE results shows that the hybrid method is applicable to FHRs. The thermal-hydraulics are coupled to neutronics for the steady-state calculations and good agreement between AGREE and SAM is achieved. (authors)
Reactor core materials exposed to rapid thermal transients during accidents can experience significant changes in mechanical properties as irradiation hardness is recovered. Understanding material performance during transients using computer modeling requires accurate descriptions of mechanical property evolution with time and temperature. A simple experimental capability was developed to expose irradiated Zircaloy to a simulated thermal transient using immersion in molten tin (at 675°C) to rapidly heat and hold before rapidly cooling with a water quench. Multiple heating-cooling cycles using a nonirradiated tensile specimen showed simplicity and viability. Based on preliminary test results, a robust system was established for routine evaluation of simulated thermal accident scenarios in the Low Activation Materials Design and Analysis (LAMDA) laboratory at Oak Ridge National Laboratory (ORNL). This paper summarizes system development and nonirradiated Zircaloy commissioning test results. Future experiments will conduct tensile tests on irradiated Zircaloy after various rapid temperature profiles.
Within the European SFR (Sodium-cooled Fast Reactor) - Safety Measures Assessment and Research Tools (ESFR-SMART) project, steady-state coupled simulation of the ESFR core has been performed using several core analysis packages, with the objective of quantifying the coupling effect. Focus is on the fuel Doppler effect and coolant expansion effect. Standalone neutronics calculations in TRACE/PARCS (PSI), DYN3D (HZDR) and WIMS (Jacobs) showed superb agreement with the reference Serpent power distribution (root mean square - rms discrepancies of 1.3%, 1.5% and 0.7% respectively). Results for COUNTHER (CIEMAT) were also in reasonable agreement, but with a somewhat higher discrepancy of 3.7%. Temperature distributions from thermal-hydraulic calculations were also compared and are found to be in good agreement. The effect of Doppler and coolant density feedback on core power distribution was predicted by TRACE/PARCS, DYN3D and WIMS to be between 0.4% and 0.8% rms difference in assembly powers. Reactivity coefficients for perturbations in the inlet temperature, flow rate and core power were shown to be negative for these three codes, with values of roughly -0.5 pcm/C. degree, -0.3 pcm/C. degree and -3.5 pcm/C. degree respectively. A preliminary investigation of differential thermal expansion effects indicates that this may have a significant effect on core power distribution of a few %, greater than anticipated a priori and may warrant inclusion in coupled core analysis to ensure the accurate calculation of power distributions. (authors)
The electrodes and solid-state electrolytes in protonic ceramic electrochemical cells (PCECs) experience significant lattice expansions when exposed to high steam concentrations at elevated temperatures. In this paper, phonon calculations based on a new machine learning potential (MLP) are employed to elucidate the volume expansions of the proton-conducting PrNi x Co 1-x O 3-δ (PNC) lattices, manifested under a combined influence of oxygen vacancies (V$^{\cdot\cdot}_O$ ) and proton uptake (OH$^{\cdot}_O$ ) in the bulk at varying Ni/Co occupancies. It is revealed that the Ni/Co occupancy contributes to thermal and chemical expansions differently, where thermal expansions are related to Co occupancy. In contrast, chemical expansions are more closely associated with the Ni occupancy. Both V$^{\cdot\cdot}_O$ and OH$^{\cdot}_O$ lead to higher thermal expansions when compared to the pristine PNC. The temperature increase will negatively impact the hydration-induced chemical expansions. For combined thermal and chemical expansions, it is predicted that the strategies that boost the PCEC's electrochemical performance may harm the electrode–electrolyte interfacial stability, when the Ni occupancy is high, due to severe chemical expansions. Mitigating chemical expansions of the Ni-abundant PNC will benefit the interfacial stability. Finally, the presented computational methods for phonon calculations, based on emerging machine learning interatomic potential techniques are anticipated to have a lasting impact on future PCEC development.
Not Available
The prototypical NAI-DMAC TADF emitter is investigated in the gas phase- and high-packing fraction limits at finite temperature, by combining first principles molecular dynamics with a quantum thermostat to account for nuclear quantum effects (NQE).
Abstract not provided.