Search NASASearch

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.

1,242 records · Page 3

Spacecraft thermal balance testing using infrared sources

A thermal balance test (controlled flux intensity) on a simple black dummy spacecraft using IR lamps was performed and evaluated, the latter being aimed specifically at thermal mathematical model (TMM) verification. For reference purposes the model was also subjected to a solar simulation test (SST). The results show that the temperature distributions measured during IR testing for two different model attitudes under steady state conditions are reproducible with a TMM. The TMM test data correlation is not as accurate for IRT as for SST. Using the standard deviation of the temperature difference distribution (analysis minus test) the SST data correlation is better by a factor of 1.8 to 2.5. The lower figure applies to the measured and the higher to the computer-generated IR flux intensity distribution. Techniques of lamp power control are presented. A continuing work program is described which is aimed at quantifying the differences between solar simulation and infrared techniques for a model representing the thermal radiating surfaces of a large communications spacecraft.

G. B. T. Tan

Thermomechanical Modeling of Woven Materials With Particle-Based, Explicit-Fiber Simulations

Fiber-based materials are extensively used to protect spacecraft during entry. Insulative fibers, often in a fiber network or woven, provide rigidity, strength, and control of material anisotropy and density. Woven thermal protection materials, such as ADEPT (Adaptable, Deployable Entry and Placement Technology), 3D-MAT (3-Dimensional Multifunctional Ablative Thermal Protection), and 3MDCP (3D Woven Mid-Density Carbon Phenolic), enable missions with stronger and denser materials for entry profiles with high shear and heat flux. Vulnerabilities to woven thermal protection materials include manufacturing-induced material property variation, and impact from micrometeoroids. Simulating woven materials under these conditions require models that can resolve hierarchal structures, thermomechanical behavior, and failure. To address this, we simulate weave thermal conduction and mechanical deformation. We simulate the full weave with a coarse-grained yarn model is presented. The model combines a validated, high-resolution single 3MDCP yarn model and phenolic resin model. Instead of modeling every fiber, each yarn ply with order 10, instead of order 1000, fibers. The discrete element bonded particle model (DEM-BPM) of fibers captures the thermal and mechanical behavior within and between fibers. We study the proportion of heat transfer and stress via the contact network, fiber bonds, and overall weave geometry.

bonded particle

Thermal Design and Thermal Vacuum Testing of the StarBurst Instrument

The StarBurst Multimessenger Pioneer is a small satellite mission serving as a wide-field gamma-ray observatory designed to capture the initial emissions of short gamma-ray bursts, electromagnetic signatures of neutron star mergers. This paper presents the final thermal design and analysis of the StarBurst Instrument, comprising the bus-to-instrument interface plate, control electronics, and twelve crystal detector units, which form the core of the mission’s science capability. The passive thermal control system design requires consideration of restrictive keep-out zones, unknown orbital parameters, and narrow temperature limits of the detectors. Also summarized is the instrument level thermal vacuum cycle test, correlated model refinements, and updated model results. Following successful completion of the instrument test campaign, the hardware was integrated with the spacecraft bus for spacecraft level testing, including additional thermal vacuum testing. The results from the spacecraft level thermal vacuum test will further inform the instrument thermal model, ensuring accurate flight temperature predictions. StarBurst launches as a secondary payload in 2027 and has a mission duration of at least one year.

StarBurst

Development and Validation of a High-Vacuum Thermal Conductivity Testbed for Aerospace Interface Materials

Spacecraft thermal margins depend on the temperature penalty of installed interfaces, yet catalog conductivity omits bondline thickness, mating surfaces, preload, and environment. The Testbed for Advanced Interface Materials in Vacuum (TAIMV) addresses this engineering-data gap by converting ambient/vacuum temperature fields from a four-coupon stack into quantities usable in spacecraft thermal models. Its 17-node steady-state model resolves axial conduction, parasitic fixture paths, grease-filled contact, radiation, and ambient convection. Four monolithic checks gave axial heat-rate ratios of 0.90–1.02, and Braycote calibration self-recovery gave 2.20% casewise mean absolute percentage error (MAPE) and 2.85% maximum difference. With transferable fixture terms frozen, Krytox LVP tested cross-material transfer against a manufacturer-derived external relation. The six Krytox cases gave 2.28% MAPE, 5.10% maximum difference, and full-field residual root-mean-square values of 0.553 °C ambient and 1.306 °C vacuum. Installed-joint resistance spanned 6.39–7.15 × 10⁻⁴ m²·K/W for Braycote and 4.09–4.49 × 10⁻⁴ m²·K/W for Krytox, corresponding to 2.33–6.66 K per modeled interface over the tested heat-flux range. TAIMV therefore supplies directly usable installed-joint resistance and conductance, plus apparent installed conductivity and explicitly model-conditioned grease conductivity/contact terms.

heat transfer

RTE: A Computer Code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Mohammad H N Naraghi

Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi: Simulations vs. Experiment

A recent microstructure-based FEM model that couples crystal-based plasticity, the B2<-> MB190 phase transformation and anisotropic elasticity at the grain scale is calibrated to recent data for polycrystalline NiTi (49.9 at.% Ni). Inputs include anisotropic elastic properties, texture and differential scanning calorimetry data, as well as a subset of recent isothermal deformation and load-biased thermal cycling data. The model is assessed against additional experimental data. Several experimental trends are captured - in particular, the transformation strain during thermal cycling monotonically increases and reaches a peak with increasing bias stress. This is achieved, in part, by modifying the martensite hardening matrix proposed by Patoor et al. [Patoor E, Eberhardt A, Berveiller M. J Phys IV 1996;6:277]. Some experimental trends are underestimated - in particular, the ratcheting of macrostrain during thermal cycling. This may reflect a model limitation that transformation-plasticity coupling is captured on a coarse (grain) scale but not on a fine (martensitic plate) scale.

Phase Transformations

Coalition for Community-Supported Affordable Geothermal Energy Systems (C2SAGES)

The C2SAGES project evaluated the feasibility of a community geothermal system for the planned Windy Ridge affordable housing development in Hinesburg, Vermont. Led by GTI Energy with Vermont Gas Systems, LN Consulting, NREL, and Frontier Energy, the work assessed technical design, energy performance, costs, business models, community engagement, maintenance, workforce development, and permitting. The proposed system was designed to serve 100% of the development’s heating, cooling, and domestic hot water loads. Compared with a baseline using air-source heat pumps and natural gas water heating, the geothermal system was estimated to reduce HVAC and domestic hot water energy use by about 45% to 48%, lower operating and maintenance costs, and reduce 30-year life-cycle costs by 37% for Phase 1 and 10% for Phase 2. Technical testing and modeling indicated that the Windy Ridge site is suitable for a community-scale geothermal system. The project also developed borehole field layouts, piping concepts, pump house designs, controls, maintenance plans, and supporting engineering drawings. The business model analysis found that first cost, ownership structure, and customer affordability remain major deployment challenges. Utility-led maintenance and operation were viewed favorably, but traditional utility cost-recovery models may require subsidy or revised financing structures to be practical for affordable housing. Community engagement highlighted the need for clear public education, transparent financing, reliable long-term maintenance, trained technicians, and the potential to pair geothermal systems with weatherization. Overall, the report concludes that community geothermal is technically feasible and offers meaningful energy, emissions, and life-cycle cost benefits, but broader deployment will depend on workable financing models and workforce readiness.

15 GEOTHERMAL ENERGY

Modeling Occupant Core Temperatures Across the Boston Building Stock to Advance Public Health

In recent history, extreme heat has been the cause of most deaths from a natural disaster. Exposure to extreme heat can aggravate preexisting conditions, increase hospitalization, and even cause death. Thus, modeling the thermal resilience of households across the United States will allow for a quantitative assessment of the health and safety risks posed by extreme heat. For the first time, we simulate occupant comfort in representative households across Boston by combining the granular results of the ResStock(TM) model with a two-node heat strain model. This model calculates occupants' core body temperature in each simulated household over a year. We compare the simulated core temperatures to two public health metrics: hyperthermia and heat stroke. Our results show for households in Boston that do not have or use aid conditioning, thousands potentially experience many dangerous heat events each summer and these events can last for more than a day at a time. These heat events peak during the late afternoon and evening, just as residents are coming home, cooking meals, and trying to go to sleep. We have found that multi-family buildings, renters, and low-income homes in more airtight and insulated homes are the most at risk for these events. These results demonstrate the scale and urgency of exposure to extreme heat.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

Mu2e Production Target Thermal Test Project

The Mu2e experiment requires a production target that is capable of operating under extreme thermal conditions caused by an 8 GeV proton beam. This project’s objective supports the development of the Mu2e Pro- duction Target by testing the Stickman model’s thermal behavior. Angel Flores Luviano has assisted Jonathan Williams in progressing this project by contributing to the development of a Radiative Cooling Test Fixture (RCTF) that will be used to evaluate the thermal behavior of the new production target model. Engineering calculations were performed to an- alyze thermal performance and pressure drop within the cold well and wa- ter cooling circuit. These calculations also determined the optimal sizing for key components of the water system. An engineering note was made to document these calculations. CAD models of the water circuit piping, thermocouple mounting bars, and radiator chimney were developed, and prototype test rig components were fabricated using 3D printing. Future work will focus on continued RCTF development which includes control system integration, heater hardware design, and interfaces that can be scaled up to increase thermal capacity.

Flores Luviano, Angel Augusto [Unlisted, US]

The thermal engine of Venus: Implications from the impact crater distribution

Questions concerning the state of Venus' volcanological activity are currently under investigation. These questions were raised by the imaging radar and altimetric data obtained by the Magellan spacecraft, which began mapping Venus in Sept. 1990. Clues about the answers to these questions are provided by the nature and distribution of impact craters on the surface of Venus. Statistical tests show that the distribution of craters on Venus cannot be distinguished from a completely spatially random population. Further, the majority of craters look relatively pristine; only about 5 percent of the observed craters appear to be embayed or flooded by lavas. The simplest interpretation of these observations is a model in which the impact craters lie on a surface that has been undisturbed by volcanological and tectonic processes since the time the surface was formed. Counting up the number of craters and estimating the rate at which meteoroids strike the surface then gives an age for the venusian surface of 500 million years; in this model the volcanic activity since that time has been essentially nil. This scenario is described as a 'production' model, or perhaps more appropriately, a 'catastrophic' model.

Phillips, Roger J.

Heat Flow and Thermal History of the Moon

An analysis was made of current heat flow data and thermal models of lunar evolution which satisfy the diverse information that has accumulated on internal processes.

Ye A Lyubimova

H 2 O and CO 2 Sorption in Ion-Exchange Sorbents: Distinct Interactions in Amine Versus Quaternary Ammonium Materials

This article examines how water (H₂O) and carbon dioxide (CO₂) interact with two classes of ion-exchange sorbents — a primary amine sorbent and a quaternary ammonium (QA⁺) sorbent — using calorimetry, thermal gravimetric analysis, gas analysis, and molecular modeling. Here, the QA⁺ sorbent exhibits stronger binding to both H₂O and CO₂ but also shows thermal stability limitations. Mixed-gas experiments reveal that humidity strongly influences CO₂ uptake and that moisture-driven sorbent regeneration enables cyclic moisture swing CO₂ capture, with implications for low-energy CO₂ separation from dilute gas streams.

36 MATERIALS SCIENCE

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

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence

WRF-Comfort: simulating microscale variability in outdoor heat stress at the city scale with a mesoscale model

Abstract. Urban overheating and its ongoing exacerbation due to global warming and urban development lead to increased exposure to urban heat and increased thermal discomfort and heat stress. To quantify thermal stress, specific indices have been proposed that depend on air temperature, mean radiant temperature (MRT), wind speed, and relative humidity. While temperature and humidity vary on scales of hundreds of meters, MRT and wind speed are strongly affected by individual buildings and trees and vary on the meter scale. Therefore, most numerical thermal comfort studies apply microscale models to limited spatial domains (commonly representing urban neighborhoods with building blocks) with resolutions on the order of 1 m and a few hours of simulation. This prevents the analysis of the impact of city-scale adaptation and/or mitigation strategies on thermal stress and comfort. To solve this problem, we develop a methodology to estimate thermal stress indicators and their subgrid variability in mesoscale models – here applied to the multilayer urban canopy parameterization BEP-BEM within the Weather Research and Forecasting (WRF) model. The new scheme (consisting of three main steps) can readily assess intra-neighborhood-scale heat stress distributions across whole cities and for timescales of minutes to years. The first key component of the approach is the estimation of MRT in several locations within streets for different street orientations. Second, mean wind speed and its subgrid variability are downscaled as a function of the local urban morphology based on relations derived from a set of microscale LES and RANS simulations across a wide range of realistic and idealized urban morphologies. Lastly, we compute the distributions of two thermal stress indices for each grid square, combining all the subgrid values of MRT, wind speed, air temperature, and absolute humidity. From these distributions, we quantify the high and low tails of the heat stress distribution in each grid square across the city, representing the thermal diversity experienced in street canyons. In this contribution, we present the core methodology as well as simulation results for Madrid (Spain), which illustrate strong differences between heat stress indices and common heat metrics like air or surface temperature both across the city and over the diurnal cycle.

Geology

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization

Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions

This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.

Chemistry