Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Projection-Based Reduced Order Modeling for Spacecraft Thermal Analysis

This paper presents a mathematically rigorous, subspace projection-based reduced order modeling (ROM) methodology and an integrated framework to automatically generate reduced order models for spacecraft thermal analysis. Two key steps in the reduced order modeling procedure are described: (1) the acquisition of a full-scale spacecraft model in the ordinary differential equation (ODE) and differential algebraic equation (DAE) form to resolve its dynamic thermal behavior; and (2) the ROM to markedly reduce the dimension of the full-scale model. Specifically, proper orthogonal decomposition (POD) in conjunction with discrete empirical interpolation method (DEIM) and trajectory piece-wise linear (TPWL) methods are developed to address the strong nonlinear thermal effects due to coupled conductive and radiative heat transfer in the spacecraft environment. Case studies using NASA-relevant satellite models are undertaken to verify the capability and to assess the computational performance of the ROM technique in terms of speed-up and error relative to the full-scale model. ROM exhibits excellent agreement in spatiotemporal thermal profiles (<0.5% relative error in pertinent time scales) along with salient computational acceleration (up to two orders of magnitude speed-up) over the full-scale analysis. These findings establish the feasibility of ROM to perform rational and computationally affordable thermal analysis, develop reliable thermal control strategies for spacecraft, and greatly reduce the development cycle times and costs.

0000↗
NASA NTRS1984

Transient behavior of a flare-associated solar wind. I - Gas dynamics in a radial open field region

A numerical investigation is conducted into the way in which a solar wind model initially satisfying both steady state and energy balance conditions is disturbed and deformed, under the assumption of heating that correspoonds to the energy release of solar flares of an importance value of approximately 1 which occur in radial open field regions. Flare-associated solar wind transient behavior is modeled for 1-8 solar radii. The coronal temperature around the heat source region rises, and a large thermal conductive flux flows inward to the chromosphere and outward to interplanetary space along field lines. The speed of the front of expanding chromospheric material generated by the impingement of the conduction front on the upper chromosphere exceeds the local sound velocity in a few minutes and eventually exceeds 100 million cm/sec.

Nagai, F.↗
NASA NTRS2014

Thermal Model Predictions of Advanced Stirling Radioisotope Generator Performance

This paper presents recent thermal model results of the Advanced Stirling Radioisotope Generator (ASRG). The three-dimensional (3D) ASRG thermal power model was built using the Thermal Desktop(trademark) thermal analyzer. The model was correlated with ASRG engineering unit test data and ASRG flight unit predictions from Lockheed Martin's (LM's) I-deas(trademark) TMG thermal model. The auxiliary cooling system (ACS) of the ASRG is also included in the ASRG thermal model. The ACS is designed to remove waste heat from the ASRG so that it can be used to heat spacecraft components. The performance of the ACS is reported under nominal conditions and during a Venus flyby scenario. The results for the nominal case are validated with data from Lockheed Martin. Transient thermal analysis results of ASRG for a Venus flyby with a representative trajectory are also presented. In addition, model results of an ASRG mounted on a Cassini-like spacecraft with a sunshade are presented to show a way to mitigate the high temperatures of a Venus flyby. It was predicted that the sunshade can lower the temperature of the ASRG alternator by 20 C for the representative Venus flyby trajectory. The 3D model also was modified to predict generator performance after a single Advanced Stirling Convertor failure. The geometry of the Microtherm HT insulation block on the outboard side was modified to match deformation and shrinkage observed during testing of a prototypic ASRG test fixture by LM. Test conditions and test data were used to correlate the model by adjusting the thermal conductivity of the deformed insulation to match the post-heat-dump steady state temperatures. Results for these conditions showed that the performance of the still-functioning inboard ACS was unaffected.

radioisotope power system↗
NASA NTRS2002

Wake Closure Characteristics and Afterbody Heating on a Mars Sample Return Orbiter

Aeroheating wind-tunnel tests were conducted on a 0.028 scale model of an orbiter concept considered for a possible Mars sample return mission. The primary experimental objectives were to characterize hypersonic near wake closure and determine if shear layer impingement would occur on the proposed orbiter afterbody at incidence angles necessary for a Martian aerocapture maneuver. Global heat transfer mappings, surface streamline patterns, and shock shapes were obtained in the NASA Langley 20-Inch Mach 6 Air and CF4 Tunnels for post-normal shock Reynolds numbers (based on forebody diameter) ranging from 1,400 to 415,000, angles of attack ranging from -5 to 10 degrees at 0, 3, and 6 degree sideslip, and normal-shock density ratios of 5 and 12. Laminar, transitional, and turbulent shear layer impingement on the cylindrical afterbody was inferred from the measurements and resulted in a localized heating maximum that ranged from 40 to 75 percent of the reference forebody stagnation point heating. Comparison of laminar heating prediction to experimental measurement along the orbiter afterbody highlight grid alignment challenges associated with numerical simulation of three- dimensional separated wake flows. Predicted values of a continuum breakdown parameter revealed significant regions of non-continuum flow downstream of the flow separation at the MSRO shoulder and in the region of the reattachment shock on the afterbody. The presence of these regions suggest that the Navier-Stokes predictions at the laminar wind-tunnel condition may encounter errors in the numerical calculation of the wake shear layer development and impingement due to non-continuum effects.

Horvath, Thomas J.↗
DOE OSTI2021

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 6061-T651 Aluminum Alloy from Plate

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The matters of interest here are the elastic-plastic response and ductile failure behavior of 6061-T651 aluminum alloy under the conditions described above. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of a modular $J_2$ plasticity model with isotropic hardening as well as a ductile failure model. Both included modules that accounted for temperature and strain rate dependence. The models were coupled with an adiabatic heating module to calculate the temperature rise due to the conversion of plastic work to heat. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s as well as four geometries of notched tension specimens and two tests on specimens with shear-dominated deformations. The test data collected allowed the calibration of both the plasticity and the ductile failure models. Most test specimens were extracted from a single piece of plate to maintain consistency. Notched tension tests came from a possibly different plate, but from the same lot. When using the model in structural finite element calculations, element formulations and sizes different from those used to model the test specimens in the calibration are likely to be used. A brief investigation demonstrated that the failure model can be particularly sensitive to the element selection and provided an initial guide to compensate in a specific example.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Thermal Transport Model for Heat Sink Design

A document discusses the development of a finite element model for describing thermal transport through microcalorimeter arrays in order to assist in heat-sinking design. A fabricated multi-absorber transition edge sensor (PoST) was designed in order to reduce device wiring density by a factor of four. The finite element model consists of breaking the microcalorimeter array into separate elements, including the transition edge sensor (TES) and the silicon substrate on which the sensor is deposited. Each element is then broken up into subelements, whose surface area subtends 10 10 microns. The heat capacity per unit temperature, thermal conductance, and thermal diffusivity of each subelement are the model inputs, as are the temperatures of each subelement. Numerical integration using the Finite in Time Centered in Space algorithm of the thermal diffusion equation is then performed in order to obtain a temporal evolution of the subelement temperature. Thermal transport across interfaces is modeled using a thermal boundary resistance obtained using the acoustic mismatch model. The document concludes with a discussion of the PoST fabrication. PoSTs are novel because they enable incident x-ray position sensitivity with good energy resolution and low wiring density.

Chervenak, James A.↗
NASA NTRS1977

Effect of haze absorption and scattering on radiative-convective-conductive heat flux divergence in Martian CO2 atmosphere and ground

The paper presents some guidelines of an improved numerical modeling effort developed to investigate the effect of an absorbing and scattering particulate phase on the temperature field of the Mars atmosphere and soil in its diurnal cycle and in response to a time-dependent convective heat transfer. Some guidelines are also formulated for the re-evaluation of Mariner 9 infrared radiometer or spectrometer inverted temperature measurements of the dust-laden atmosphere.

Pallmann, A. J.↗
DOE OSTI2025

Phenomena Identification and Ranking Table (PIRT) for Heat Pipes

This Phenomena Identification and Ranking Table (PIRT) report provides an evaluation of key phenomena affecting the performance and operational regimes of heat pipes, particularly in the context of heat pipe microreactors (HPMRs). Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in the heat pipe, including phase change, turbulent transition, and compressibility effects, among others, there is high uncertainty in identifying and ranking the important phenomena affecting the operation of heat pipes and the current knowledge for their modeling and simulation and experimental measurements and instrumentation. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. The report analyzes phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics, discussing challenges and future research directions for improving their modeling and simulation and experimental measurements. Additionally, the report addresses phenomena with low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geysering, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. This comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
NASA NTRS2026

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

Thermal Interface Materials (TIMs) are critical components in spacecraft thermal management systems, where thermal performance is strongly influenced by vacuum conditions, interface contact resistance, and layered metallic joint behavior. However, manufacturer-reported thermal conductivity values are often derived under idealized conditions and may not accurately represent performance within operational aerospace applications. To address this limitation, the Testbed for Advanced Interface Materials in Vacuum (TAIMV) was developed as a modular vacuum-compatible thermal conductivity characterization platform capable of evaluating aerospace-relevant TIM configurations under both ambient and high-vacuum environments. The testbed was derived from the ASTM C1044-16 guarded hot plate methodology and incorporates interchangeable layers of stainless steel coupon geometries, independently controlled main and guard heaters, embedded resistance temperature detectors (RTDs), thermocouples, multi-layer insulation (MLI), and a temperature-controlled cold plate to characterize through-thickness thermal gradients across layered interfaces. In the current configuration, interface compression is limited to the nominal contact pressure generated by the experimental stack assembly. Initial experimental campaigns were conducted at ambient pressure and below 1×10-5 torr for vacuum cases using multiple interface materials including Braycote 601EF and Krytox-based greases across a range of thermal operating conditions. In parallel, a coupled numerical Python thermal model was developed to predict temperature distribution throughout the stack while accounting for conduction, radiation, and parasitic heat transfer pathways and effective interface resistance effects. Experimental measurements and numerical predictions showed consistent thermal trends across multiple operating conditions and environmental states. Results also revealed measurable differences between ambient and vacuum thermal behavior, demonstrating the importance of interface resistance, parasitic heat transfer mechanisms, and stack geometry in determining effective thermal performance within layered thermal interfaces. The presented work establishes a foundation for future thermal model correlation efforts and expanded characterization of aerospace thermal interface materials under representative environmental conditions. Future work will focus on the integration of a load cell system to enable controlled pressure-dependent characterization of thermal interface materials under compressive loading. This capability will allow investigation of the influence of contact pressure on effective thermal conductivity, interface resistance, and thermal performance within layered aerospace thermal interfaces under representative operational conditions.

Thermal Development Testing↗
DOE OSTI2022

Experimental and numerical study on the impacts of various long carbon fiber inserts on the thermal performance of phase change materials

In thermal energy storage systems, phase change materials (PCMs) are widely used to provide efficient storing and extraction of thermal energy. However, most PCMs have low thermal conductivities, which limit the heat transfer of the latent energy inward/outward. In this study, the impacts of various fiber inserts on the solidification and melting processes with PCM were investigated by experiments and simulations. It was found that long anisotropic carbon fibers with high thermal conductivities can not only increase the heat transfer within PCM by enlarging hot-cold interface areas, but also act as crystallization nuclei and increase the solidification of PCM during the charging process. Recalescence of PCM during solidification was observed. COMSOL models were used to understand the thermal transfer mechanisms. COMSOL Multiphysics models of a single carbon fiber filament in PCM showed that the anisotropy in thermal conductivities controlled the radial and axial heat transfer ratios along its fiber axis and can enhance both the solidification and melting processes.

Zhou, Xinzhang↗
DOE OSTI2021

Engineering Graphene-Ceramic 3D Composite Foams by Freeze Drying

A 3D graphene foam produced by the chemical vapor deposition (CVD) technique is recognized as an effective nanofiller material, as it does not agglomerate in the matrix. Though CVD facilitates pristine graphene foam production, the method poses limitations in developing large-scale 3D graphene composite foams as advanced engineering materials. In this paper, a freeze-drying (FD) process is used to produce a composite 3D foam, which is a mixture of graphene nanoplatelet (GNP) and a low-temperature co-fired ceramic (LTCC). The freeze-dried GNP-LTCC reticulated 3D composite foam has an average pore size ranging from 70 to 100 μm. Pore size is varied by regulating the heat transfer rate during the freezing process using a thermally conductive aluminum (Al) mold and a thermally insulating acrylonitrile butadiene styrene (ABS) mold. Computational thermal modeling is used to visualize the heat transfer and its effect on foam pore size. Subsequently, a freeze-dried GNP-LTCC composite foam is embedded into the LTCC matrix to form a hierarchical assembly by a spark plasma technique without compromising the 3D structure of the FD foam. This study established that a simple, eco-friendly, and scalable processing methodology can produce advanced graphene-based 3D composite foams as the future tailorable nanofillers for designing multimatrix materials.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Results of phase change paint tests of 0.040 scale 50% forebody models (82-0) of the space shuttle orbiter in the AEDC VKF B hypersonic wind tunnel (OH75)

Post-test information and data are presented from phase change paint, aerodynamic heating wind tunnel tests of a Rockwell International space shuttle orbiter forebody model. These tests were conducted in the Arnold Engineering and Development Center von Karman Facility Tunnel B Hypersonic Wind Tunnel. The purpose of these tests was to determine the effect of simulated orbiter protuberances and penetrations (including RCS nozzles) on aerodynamic heating rates during simulated entry conditions.

Dye, W. H.↗
DOE OSTI2020

City-Scale Building Anthropogenic Heating during Heat Waves

More frequent and longer duration heat waves have been observed worldwide and are recognized as a serious threat to human health and the stability of electrical grids. Past studies have identified a positive feedback between heat waves and urban heat island effects. Anthropogenic heat emissions from buildings have a crucial impact on the urban environment, and hence it is critical to understand the interactive effects of urban microclimate and building heat emissions in terms of the urban energy balance. Here we developed a coupled-simulation approach to quantify these effects, mapping urban environmental data generated by the mesoscale Weather Research and Forecasting (WRF) coupled to Urban Canopy Model (UCM) to urban building energy models (UBEM). We conducted a case study in the city of Los Angeles, California, during a five-day heat wave event in September 2009. We analyzed the surge in city-scale building heat emission and energy use during the extreme heat event. We first simulated the urban microclimate at a high resolution (500 m by 500 m) using WRF-UCM. We then generated grid-level building heat emission profiles and aggregated them using prototype building energy models informed by spatially disaggregated urban land use and urban building density data. The spatial patterns of anthropogenic heat discharge from the building sector were analyzed, and the quantitative relationship with weather conditions and urban land-use dynamics were assessed at the grid level. The simulation results indicate that the dispersion of anthropogenic heat from urban buildings to the urban environment increases by up to 20% on average and varies significantly, both in time and space, during the heat wave event. The heat dispersion from the air-conditioning heat rejection contributes most (86.5%) of the total waste heat from the buildings to the urban environment. We also found that the waste heat discharge in inland, dense urban districts is more sensitive to extreme events than it is in coastal or suburban areas. The generated anthropogenic heat profiles can be used in urban microclimate models to provide a more accurate estimation of urban air temperature rises during heat waves.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1972

An analysis of the booster plume impingement environment during the space shuttle nominal staging maneuver

An experimental study of the plume impingement heating on the space shuttle booster afterbody resulting from the space shuttle orbiter engine plumes was conducted. The 1/100-scale model tests consisted of one and two orbiter engine firings on a flat plate, a flat plate with a fin, and a cylinder model. The plume impingement heating rates on these surfaces were measured using thin film heat transfer gages. Results indicate the engine simulation is a reasonable approximation to the two engine configuration, but more tests are needed to verify the plume model of the main engine configuration. For impingment, results show models experienced laminar boundary layer convective heating. Therefore, tests at higher Reynolds numbers are needed to determine impingment heating.

Wojciechowski, C. J.↗
DOE OSTI2021

Multidimensional Numerical Modeling of Combustion Dynamics in a Non-Premixed Rotating Detonation Engine With Adaptive Mesh Refinement

In the present work, a novel computational fluid dynamics (CFD) methodology was developed to simulate full-scale non-premixed rotating detonation engines (RDEs). A unique feature of the modeling approach was the incorporation of adaptive mesh refinement (AMR) to achieve a good trade-off between model accuracy and computational expense. Here, unsteady Reynolds-averaged Navier–Stokes (RANS) simulations were performed for an Air Force Research Laboratory (AFRL) non-premixed RDE configuration with hydrogen as fuel and air as the oxidizer. The finite-rate chemistry model, along with a ten-species detailed kinetic mechanism, was employed to describe the H 2 -Air combustion chemistry. Three distinct operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel/air mass flowrates. For all conditions, the capability of the model to capture essential detonation wave dynamics was assessed. An exhaustive verification and validation study was performed against experimental data in terms of a number of waves, wave frequency, wave height, reactant fill height, oblique shock angle, axial pressure distribution in the channel, and fuel/air plenum pressure. The CFD model was demonstrated to accurately predict the sensitivity of these wave characteristics to the operating conditions, both qualitatively and quantitatively. A comprehensive heat release analysis was also conducted to quantify detonative versus deflagrative burning for the three simulated cases. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

42 ENGINEERING↗
NASA NTRS1999

Integrated Modeling Tools for Thermal Analysis and Applications

Integrated modeling of spacecraft systems is a rapidly evolving area in which multidisciplinary models are developed to design and analyze spacecraft configurations. These models are especially important in the early design stages where rapid trades between subsystems can substantially impact design decisions. Integrated modeling is one of the cornerstones of two of NASA's planned missions in the Origins Program -- the Next Generation Space Telescope (NGST) and the Space Interferometry Mission (SIM). Common modeling tools for control design and opto-mechanical analysis have recently emerged and are becoming increasingly widely used. A discipline that has been somewhat less integrated, but is nevertheless of critical concern for high precision optical instruments, is thermal analysis and design. A major factor contributing to this mild estrangement is that the modeling philosophies and objectives for structural and thermal systems typically do not coincide. Consequently the tools that are used in these discplines suffer a degree of incompatibility, each having developed along their own evolutionary path. Although standard thermal tools have worked relatively well in the past. integration with other disciplines requires revisiting modeling assumptions and solution methods. Over the past several years we have been developing a MATLAB based integrated modeling tool called IMOS (Integrated Modeling of Optical Systems) which integrates many aspects of structural, optical, control and dynamical analysis disciplines. Recent efforts have included developing a thermal modeling and analysis capability, which is the subject of this article. Currently, the IMOS thermal suite contains steady state and transient heat equation solvers, and the ability to set up the linear conduction network from an IMOS finite element model. The IMOS code generates linear conduction elements associated with plates and beams/rods of the thermal network directly from the finite element structural model. Conductances for temperature varying materials are accommodated. This capability both streamlines the process of developing the thermal model from the finite element model, and also makes the structural and thermal models compatible in the sense that each structural node is associated with a thermal node. This is particularly useful when the purpose of the analysis is to predict structural deformations due to thermal loads. The steady state solver uses a restricted step size Newton method, and the transient solver is an adaptive step size implicit method applicable to general differential algebraic systems. Temperature dependent conductances and capacitances are accommodated by the solvers. In addition to discussing the modeling and solution methods. applications where the thermal modeling is "in the loop" with sensitivity analysis, optimization and optical performance drawn from our experiences with the Space Interferometry Mission (SIM), and the Next Generation Space Telescope (NGST) are presented.

Milman, Mark H.↗
NASA NTRS1982

Ionospheric hot spot at high latitudes

Schunk and Raitt (1980) and Sojka et al. (1981) have developed a model of the convecting high-latitude ionosphere in order to determine the extent to which various chemical and transport processes affect the ion composition and electron density at F-region altitudes. The numerical model produces time-dependent, three-dimensional ion density distributions for the ions NO(+), O2(+), N2(+), O(+), N(+), and He(+). Recently, the high-latitude ionospheric model has been improved by including thermal conduction and diffusion-thermal heat flow terms. Schunk and Sojka (1982) have studied the ion temperature variations in the daytime high-latitude F-region. In the present study, a time-dependent three-dimensional ion temperature distribution is obtained for the high-latitude ionosphere for an asymmetric convection electric field pattern with enhanced flow in the dusk sector of the polar region. It is shown that such a convection pattern produces a hot spot in the ion temperature distribution which coincides with the location of the strong convection cell.

Schunk, R. W.↗
NASA NTRS1995

Helioseismology: A probe of the solar interior, atmosphere, and activity cycle

Helioseismology began in earnest in the mid 1970's. In the two decades which have elapsed since that time this branch of solar physics has become a mature field of research. Helioseismology has demonstrated that the solar convection zone is about twice as deep as was generally thought to be the case before 1977. Helioseismology has also provided measurements of the solar internal angular velocity over much of the sun's interior. Helioseismology has also ruled out models which would solve the solar neutrino problem by a lowering of the temperature of the core. Recently, some of the seismic properties of the sun have been demonstrated to vary with changing levels of solar activity. Also, helioseismology has recently provided evidence for helical flow patterns in the shallow, sub-photosphere layers. The techniques of helioseismology are also expanding to include seismic probes of solar active regions. Some work is also being conducted into the possible contributions of the solar acoustic models to the heating of the solar atmosphere. In this talk I will highlight a few of the above results and concentrate on current areas of research in the field.

Rhodes, E. J., Jr.↗