Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermal modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 649 records · Page 36

Release on the Virtual Test Bed of an MSRE thermal hydraulics model

The nuclear industry is taking leaps in innovations with companies seeking a sustainable energy future through advanced nuclear reactors. The \gls{DOE}’s \gls{neams} program seeks to substantiate and bolster the deployment of advanced reactors through flexible multifidelity, multiphysics simulations of advanced nuclear reactors. Applications like SAM for one-dimensional systems thermal-hydraulics, and Pronghorn for multidimensional coarse mesh thermal-hydraulics, are geared to support innovations in industry by facilitating design, optimization, and licensing of advanced nuclear reactors. Coupling systems thermal-hydraulics and computational fluid dynamics codes can be difficult as the pressure coupling converges slowly; however, it is important to obtain the desired accuracy in each part of the primary loop. The authors of this model created an \gls{odc}~\cite{Mau23} approach to coupling SAM and Pronghorn. Leveraging this coupling technique, a \gls{msre} model was developed and released to the \gls{neams}/\gls{nric} \gls{vtb}. The \gls{msre} was chosen to be modeled because of the wealth of experimental data available and because of the strong physics coupling between the core and primary circuit \cite{doi:10.13182/NT8-2-118}. This document contextualizes the history of the \gls{msre}, describes the thermal hydraulics models used, and detail the implementation of multidimensional thermal-hydraulics and system codes based on the \gls{odc} method~\cite{Penn} for the \gls{msre} model. Finally, this document presents how other modelers could apply the SAM-Pronghorn \gls{odc} for other advanced reactor models. Current item is the set of slides for ANS Winter 23. The release of the model and the ANS summary have already been approved

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Ablation Modeling of Ares-I Upper State Thermal Protection System Using Thermal Desktop

The thermal protection system (TPS) for the Ares-I Upper Stage will be based on Space Transportation System External Tank (ET) and Solid Rocket Booster (SRB) heritage materials. These TPS materials were qualified via hot gas testing that simulated ascent and re-entry aerothermodynamic convective heating environments. From this data, the recession rates due to ablation were characterized and used in thermal modeling for sizing the thickness required to maintain structural substrate temperatures. At Marshall Space Flight Center (MSFC), the in-house code ABL is currently used to predict TPS ablation and substrate temperatures as a FORTRAN application integrated within SINDA/G. This paper describes a comparison of the new ablation utility in Thermal Desktop and SINDA/FLUINT with the heritage ABL code and empirical test data which serves as the validation of the Thermal Desktop software for use on the design of the Ares-I Upper Stage project.

Sharp, John R.↗

Enabling Long-Duration Lunar Equatorial Operations With Thermal Wadi Infrastructure

Long duration missions on the Moon s equator must survive lunar nights. With 350 hr of cryogenic temperatures, lunar nights present a challenge to robotic survival. Insulation is imperfect, so it is not possible to passively contain enough heat to stay warm through the night. Components that enable mobility, environmental sensing and solar power generation must be exposed, and they leak heat. Small, lightweight rovers cannot store enough energy to warm components throughout the night without some external source of heat or power. Thermal wadis, however, can act as external heat sources to keep robots warm through the lunar night. Electrical power can also be provided to rovers during the night from batteries stored in the ground beside wadis. Buried batteries can be warmed by the wadi s heat. Results from analysis of the interaction between a rover and a wadi are presented. A detailed three-dimensional (3D) thermal model and an easily configurable two-dimensional (2D) thermal model are used for analysis.

Jones, Heather L.↗

Development of an Emulation-simulation Thermal Control Model for Space Station Application

Analysis techniques to evaluate the effects of changing thermal loads and the methods utilized to control temperature distributions in the orbital space station are essential. Analysis techniques including a user-friendly computer program, were developed which should prove useful in thermal design and system analysis of the the space station. The program uses a data base and user input to compute costs, sizes, and power requirements for individual components and complete systems. User input consists of selecting mission parameters, selecting thermal acquisition configurations, transport systems and distances, and thermal rejection configurations.

Colwell, G. T.↗

Development of an emulation-simulation thermal control model for space station application

An improved capability for comparing various techniques for thermal management in the Space Station was developed. Current planning for the orbiting space station calls for a dual keel configuration. The thermal control system (TCS) for the space is composed of a central TCS and internal thermal control systems for the modules, as well as service facilities and attached payloads. The internal TCS may be attacted to the central TCS through a thermal bus. The central TCS is composed of a main transport system which collects waste thermal energy from each of the modules and transports it through coolant lines to the main rejection system. The waste heat loads in the modules arise from electrical and electronic equipment as well as metabolic loads in the manned modules. Several candidate technologies are being considered for acquiring the waste heat loads, for transporting the thermal energy between the acquisition and rejection systems, and for rejecting the waste heat to space. The analysis techniques described were developed for use in evaluating reliability, weights, costs, volumes, and power requirements for configurations using different candidates and different mission parameters.

Hartley, James G.↗

Residual Stresses Modeled in Thermal Barrier Coatings

Thermal barrier coating (TBC) applications continue to increase as the need for greater engine efficiency in aircraft and land-based gas turbines increases. However, durability and reliability issues limit the benefits that can be derived from TBC's. A thorough understanding of the mechanisms that cause TBC failure is a key to increasing, as well as predicting, TBC durability. Oxidation of the bond coat has been repeatedly identified as one of the major factors affecting the durability of the ceramic top coat during service. However, the mechanisms by which oxidation facilitates TBC failure are poorly understood and require further characterization. In addition, researchers have suspected that other bond coat and top coat factors might influence TBC thermal fatigue life, both separately and through interactions with the mechanism of oxidation. These other factors include the bond coat coefficient of thermal expansion, the bond coat roughness, and the creep behavior of both the ceramic and bond coat layers. Although it is difficult to design an experiment to examine these factors unambiguously, it is possible to design a computer modeling "experiment" to examine the action and interaction of these factors, as well as to determine failure drivers for TBC's. Previous computer models have examined some of these factors separately to determine their effect on coating residual stresses, but none have examined all the factors concurrently. The purpose of this research, which was performed at DCT, Inc., in contract with the NASA Lewis Research Center, was to develop an inclusive finite element model to characterize the effects of oxidation on the residual stresses within the TBC system during thermal cycling as well as to examine the interaction of oxidation with the other factors affecting TBC life. The plasma sprayed, two-layer thermal barrier coating that was modeled incorporated a superalloy substrate, a NiCrAlY bond coat, and a ZrO2-8 wt % Y2O3 ceramic top coat. We examined the effect on stress during burner rig thermal cycling of the following independent variables: creep in the bond coat and top coat, oxidation, bond coat coefficient of thermal expansion, number of thermal cycles, and interfacial roughness. All these factors were suspected of influencing TBC failure. The model showed that all the material properties studied had a significant effect on the coating's residual stresses if the interface of the bond coat was rough. Bond coat expansion, bond coat oxidation, and bond coat creep had the highest effect on coating stresses, and these were highly interactive. The model also showed that the mechanism of stress generation during thermal cycling changed with the number of thermal cycles. Bond coat and top coat creep dominated stress generation during early thermal cycles, greatly increasing delamination stresses at the peaks of the bond coat. Therefore, creep is the prime driver for delamination cracking early in life, but cracking is limited to the bond coat peak region. Oxidation of the bond coat, on the other hand, tended to dominate stress generation during later cycles by greatly increasing delamination stresses over bond coat valleys. These results indicate that oxidation is the driver for the continued cracking necessary to cause ceramic layer spallation.

Freborg, A. M.↗

Thermal Analysis of Potted Litz Wire for High-Power-Density Aerospace Electric Machines

Increasing the power density and efficiency of electric machines (motors and generators) is integral to bringing Electrified Aircraft (EA) to commercial realization. To that end an effort to create a High Efficiency Megawatt Motor (HEMM) with a goal of exceeding 98% efficiency and 1.46 MW of power has been undertaken at the NASA Glenn Research Center. Of the motor components the resistive losses in the stator windings are by far the largest contributor (34%) to total motor loss. The challenge is the linear relationship between resistivity and temperature, making machine operation sensitive to temperature increases. In order to accurately predict the thermal behavior of the stator the thermal conductivity of the Litz wire-potting-electrical insulation system must be known. Unfortunately, this multi material system has a wide range of thermal conductivities (0.1 W/m-K – 400 W/m-K) and a high anisotropy (axial vs transverse) making the prediction of the transverse thermal conductivity an in turn the hot spot temperatures in the windings is difficult. In order to do this a device that simulates the thermal environment found in the HEMM stator was designed. This device is not unlike the motorettes (little motors) that are described in IEEE standards for testing electrical insulation lifetimes or other electric motor testing. However, because the HEMM motor design includes significant rotor electrical and thermal considerations the term motorette was not deemed appropriate. Instead statorette (or little stator) was adopted as the term for this test device. This paper discussed the design, thermal heat conjugate analysis (thermal model), manufacturing and testing of HEMM's statorette. Analysis of the results is done by thermal resistance network model and micro thermal model and is compared to analytical predictions of thermal conductivity of the insulated and potted Litz wire system.

Motor↗

Mars 2020 Mobility Actuator Thermal Testing and Model Correlation

This paper describes the thermal testing and model correlation of a mobility actuator planned for use on the Mars 2020 Rover. The mobility actuator is identical to those which have been successfully flown and operated on the Mars Science Laboratory (MSL) Curiosity rover since its successful landing on Mars in August of 2012. The actuator consists of a motor, brake, and encoder paired with a four stage planetary gear box. In this thermal test, the actuator was instrumented with a number of thermocouples on both the interior and exterior of the gearbox. Heaters and a cold plate were used to generate thermal gradients across the actuator in vacuum, low pressure GN2, and low pressure CO2 environments in an effort to correlate a thermal model and develop a better understanding of how heat flows through the mechanism. This testing resulted in the successful model correlation of a simplified thermal model, and yielded important insights regarding the conductance of ball bearings and gear-to-gear contact. Ball bearing thermal conductance in a low pressure environment can be estimated by using correlations for vacuum thermal conductance along with a multiplier to account for increased gas conduction, and gear-to-gear conductance can be estimated by accounting for gas and grease conduction between two gears.

Novak, Keith↗

Modeling of Thermal Barrier Coatings

The project examined the effectiveness of studying the creep behavior of thermal barrier coating system through the use of a general purpose, large strain finite element program, NIKE2D. Constitutive models implemented in this code were applied to simulate thermal-elastic and creep behavior. Four separate ceramic-bond coat interface geometries were examined in combination with a variety of constitutive models and material properties. The reason for focusing attention on the ceramic-bond coat interface is that prior studies have shown that cracking occurs in the ceramic near interface features which act as stress concentration points. The model conditions examined include: (1) two bond coat coefficient of thermal expansion curves; (2) the creep coefficient and creep exponent of the bond coat for steady state creep; (3) the interface geometry; and (4) the material model employed to represent the bond coat, ceramic, and superalloy base.

Ferguson, B. L.↗

Lunar Surface Crater Thermal Effects on Lander Radiator Performance

Lunar surface craters smaller than the spatial resolution of surface meshes used in typical Lunar surface thermal models (10 to 60 meters per pixel) may impact the accuracy of thermal model extrema predictions. The goal of this study is to investigate the thermal sensitivity of representative lander systems with realistic thermal surface orientations in bare and cratered terrain environments at relevant Artemis mission locations. This thermal analysis task investigates the impact of lunar surface craters on lander radiator performance by comparing heat rejection capability results between bare and cratered terrain environments. This study examines external body-mounted lander radiator thermal performance across varying lander heights (5.5m, 20m, 50m) and radiator orientations (horizontal, 45º tilted, and vertical) at two representative Artemis mission latitudes (-89.5ºS and -82.5ºS), spanning from the Shackleton Connecting Ridge to Mons Mouton Plateau.

Lunar Surface↗

Global isostatic geoid anomalies for plate and boundary layer models of the lithosphere

Commonly used one dimensional geoid models predict that the isostatic geoid anomaly over old ocean basins for the boundary layer thermal model of the lithosphere is a factor of two greater than that for the plate model. Calculations presented, using the spherical analogues of the plate and boundary layer thermal models, show that for the actual global distribution of plate ages, one dimensional models are not accurate and a spherical, fully three dimensional treatment is necessary. The maximum difference in geoid heights predicted for the two models is only about two meters. The thermal structure of old lithosphere is unlikely to be resolvable using global geoid anomalies. Stripping the effects of plate aging and a hypothetical uniform, 35 km, isostatically-compensated continental crust from the observed geoid emphasizes that the largest-amplitude geoid anomaly is the geoid low of almost 120 m over West Antarctica, a factor of two greater than the low of 60 m over Ceylon.

Hager, B. H.↗

Mineralogy and Thermal Properties of V-Type Asteroid 956 Elisa: Evidence for Diogenitic Material from the Spitzer IRS (5-35 Micrometers) Spectrum

We present the thermal infrared (5-35 micrometer) spectrum of 956 Elisa as measured by the Spitzer Infrared Spectrograph ("IRS"; Houck,1.R. et .11. [20041. Astrophys, 1. SuppL 154, 18-24) together with new ground-based lightcurve data and near-IR spectra. From the visible lightcurve photometry, we determine a rotation period of 16.494 +/- 0.001 h, identify the rotational phase of the Spitzer observations, and estimate the visible absolute magnitude (Hv) at that rotational phase to be 12.58 +/- 0.04. From radiometric analysis of the thermal flux spectrum, we find that at the time of observation 956 Elisa had a projected radius of 5.3 +/- 0.4 km with a visible albedo pv = 0.142+/- 0.022, significantly lower than that of the prototype V-type asteroid, 4 Vesta. (This corresponds to a radius of 5.2 +/- 0.4 km at lightcurve mean.) Analysis with the standard thermal model (STM) results in a sub-solar temperature of 292.3 +/- 2.8 K and beaming parameter eta = 1.16 +/- 0.05. Thermophysical modeling places a lower limit of 20 J m(exp -2)K(exp -1)s(exp -1/2) on the thermal inertia of the asteroid's surface layer (if the surface is very smooth) but more likely values fall between 30 and 150 J m(exp -2)K(exp -1)s(exp -1/2) depending on the sense of rotation. The emissivity spectrum, calculated by dividing the measured thermal flux spectrum by the modeled thermal continuum, exhibits mineralogically interpretable spectral features within the 9-12 micrometer reststrahlen band, the 15-16.5 micrometer Si-O-Si stretching region, and the 16-25 micrometer reststrahlen region that are consistent with pyroxene of diogenitic composition: extant diogenitic pyroxenes fall within the narrow compositional range W0(sub 2+/-1)En(sub 74+/-2)Fs(sub 24+/-1). Spectral deconvolution of the 9-12 micrometer reststrahlen features indicates that up to approximately 20% olivine may also be present, suggesting an olivine-diogenite-like mineralogy. The mid-IR spectrum is inconsistent with non-cumulate eucrite as the major component on the surface of 956 Elisa, although cumulate eucrite material may be present at abundances lower than that of the diogenite component. Analysis of new near-IR spectra of 956 Elisa with the Modified Gaussian Model (MGM; Sunshine, J,M., Pieters, C.M., Pratt, S.F.[1990]. J. Geophys. Res, 95 (May), 6955-6966) results in two pyroxene compositions: 75% magnesian low-Ca pyroxene and 25% high-Ca pyroxene. High-Ca pyroxene is not evident in the mid-IR data, but may belong to a component that is underrepresented in the mid-IR spectrum either because of its spatial distribution on the asteroid or because of its particle size. High-Ca pyroxenes that occur as exsolution lamellae may also be more evident spectrally in the NIR than in the mid-IR. In any case, we find that the mid-IR spectrum of 956 Elisa is dominated by emission from material of diogenite- like composition, which has very rarely been observed among asteroids.

Lim, Lucy F.↗

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗

Multiplicity dependence of f 0 (980) production in pp collisions at $\mathbf{\sqrt{s}}$ = 13 TeV

The dependence of f 0 (980) production on the final-state charged-particle multiplicity is reported for proton–proton (pp) collisions at the centre-of-mass energy, $\mathbf{\sqrt{s}}$ = 13 TeV. The production of (980) is measured with the ALICE detector via the decay channel in a midrapidity region of |y| < 0.5. The evolution of the integrated yields and mean transverse momentum of f 0 (980) as a function of charged-particle multiplicity measured in pp at $\mathbf{\sqrt{s}}$ = 13 TeV follows the trends observed in pp at $\mathbf{\sqrt{s}}$ = 5.02 TeV and in proton–lead (p–Pb) collisions at$\mathbf{\sqrt{s_{NN}}}$ = 5.02 TeV. Particle yield ratios of f 0 (980) to π ± and K* (892) are found to decrease with increasing charged-particle multiplicity. These particle ratios are compared with calculations from the canonical statistical thermal model as a function of charged-particle multiplicity. The thermal model calculations provide a better description of the decreasing trend of particle ratios when no strange or antistrange quark composition for f 0 (980) is assumed, which suggests that the data do not support significant hidden strangeness in the f 0 (980).

Abualrob, I. J. [University of Houston] (ORCID:000↗

Accelerated Fuel Qualification of Fast Modular Reactor Fuel in a Thermal Reactor: Modeling and Simulation Paired with Irradiation Testing

The accelerated fuel qualification (AFQ) methodology is applied by simulating accelerated fuel tests of the General Atomics Electromagnetic Systems’ fuel system for its 44-MW(electric) gas-cooled, fast-spectrum fast modular reactor (FMR). This fuel is comprised of UO 2 pellets in SiGA® cladding, a silicon carbide ceramic matrix composite. Fast reactors, like the FMR, offer many benefits, including high fuel utilization and flexibility, but may require a lengthy material design process if tests are performed using fast neutron irradiation alone. A thermal neutron irradiation can instead be used to rapidly test how well key components of the current material models extend to high burnup. Thermal neutrons produce a different radial power distribution within the pin than fast neutrons. However, the temperature and burnup values for the two neutron types are comparable, and the differences between the simulated fuel responses are relatively small, demonstrating the weak sensitivity of the physics-based fuel model calculations on the neutron type and the irradiation rate. Furthermore, the deformation of the SiGA cladding saturates after about 1 displacement per atom for both neutron spectra. In an accelerated fuel test, the irradiation time required to reach the target fuel burnup can be reduced by a factor of 3 by using a small rodlet with a 45% smaller pellet diameter while maintaining the same linear power. Therefore, the time for data collection up to high burnup can be significantly reduced while maintaining the same temperature profile, which largely determines the material response. Tests of fuel rodlets of standard and compact size will be carried out in the Idaho National Laboratory’s Advanced Test Reactor (ATR), including full size and compact rodlets with varying gap sizes. By applying physics-based mechanistic modeling and simulation in accordance with the AFQ methodology, this type of compact rodlet testing in a thermal test reactor captures the necessary phenomena to test fuel material models up to high burnup and to simulate the expected impact of fast neutron radiation on the fuel in FMR operations. Finally, this approach to testing fast reactor fuels in existing thermal test reactors, paired with advanced physics-based mechanistic modeling and simulation, is expected to be applicable to a range of advanced fuels and will decrease the overall fuel qualification timeframe from decades to years.

Advanced test reactor (ATR)↗

GPHS Module Pre-Impact test Temperature Prediction with ANSYS [PowerPoint]

The purpose of the project was to: Develop a finite element thermal model of the General Purpose Heat Source (GPHS) Module loaded with four fueled clads; Perform steady state analysis to validate the thermal model via energy balance; Perform transient thermal analysis to predict clad temperatures at various times; and, Transient model results will be used to inform module soak time needed to reach a desired clad temperature prior to impact.

42 ENGINEERING↗

Predicting Maximum Temperatures of a Li-ion Battery on a Simulated Flight Profile using a Model-based Prognostics

One of the challenges in using Li-ion packs in aeronautics is their safety, and thermal runaway (TR) is a major concern. The current engineering solutions to prevent a Li-ion pack from a catastrophic TR require additional mass and volume to isolate cells. The excess mass could be reduced by improving detection and, thus, preventing a TR event. One of the possible early warning indicators of a TR is crossing a threshold temperature. We have developed an approach, based on the Unscented Kalman Filter (UKF), to predict the likelihood of reaching the threshold temperature for simulated flight profiles. The current battery prognostics algorithms for aerospace predict state-of-charge (SOC) and end-of-life (EOL) [1]. We extended this two-level algorithm to predict the maximum temperature during discharge. The amount of heat generated in a cell depends on factors such as cell chemistry, cell packaging, total cycles, operating temperature, and abuse history [2]. Our semi-empirical thermal model depends on three phenomenological parameters which account for those factors. In addition, a two-parameter reduced-order model is developed to predict the temperature rise for short bursts of “random-walk” (RW) discharge current sequence, which simulates a flight's current-loading profile. The performance of these models on different datasets and types of current loading will be presented. To predict the maximal temperatures for future cycles we must estimate the evolution of thermal parameters as the batteries age. It is found that the parameters of the 3-parametric thermal model cannot be estimated only from the RW data. To address the issue, we will present two alternative approaches: i) expanding the datasets to include discharge profiles beyond RWs; ii) model reduction to a two-parametric model. The two approaches will be illustrated by an application to the cycling data from a commercial LG 18650 cell. References: 1. M. Daigle, C.S. Kulkarni, End-of-discharge and End-of-life Prediction in Lithium-ion Batteries with Electrochemistry-based Aging Models, in: AIAA Infotech @ Aerospace, American Institute of Aeronautics and Astronautics, San Diego, California, USA, 2016. 2. M. Börner, et. al, Correlation of aging and thermal stability of commercial 18650-type lithium ion batteries, Journal of Power Sources. 342 (2017) 382–392.

Thermal runaway↗

Models of thermal relaxation mechanics for particle simulation methods

A comparative study is conducted of the simulation of molecular vibrational relaxation in rarefied flows by means of existing models and a novel phenomenological microscopic model. This improved technique iterates between translation-rotation and rotation-vibration exchanges; equilibrium is in this way promoted if the model for the latter process is compatible with quantized oscillators. The iteration-equipartition model retains computational simplicity and promotes thermal equilibrium, even in the case of nondegenerate anharmonic quantized oscillators and multispecies gas mixtures.

Haas, Brian L.↗