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 379 records · Page 21

Time and Space Resolved Heat Transfer Measurements Under Nucleate Bubbles with Constant Heat Flux Boundary Conditions

Investigations into single bubble pool boiling phenomena are often complicated by the difficulties in obtaining time and space resolved information in the bubble region. This usually occurs because the heaters and diagnostics used to measure heat transfer data are often on the order of, or larger than, the bubble characteristic length or region of influence. This has contributed to the development of many different and sometimes contradictory models of pool boiling phenomena and dominant heat transfer mechanisms. Recent investigations by Yaddanapyddi and Kim and Demiray and Kim have obtained time and space resolved heat transfer information at the bubble/heater interface under constant temperature conditions using a novel micro-heater array (10x10 array, each heater 100 microns on a side) that is semi-transparent and doubles as a measurement sensor. By using active feedback to maintain a state of constant temperature at the heater surface, they showed that the area of influence of bubbles generated in FC-72 was much smaller than predicted by standard models and that micro-conduction/micro-convection due to re-wetting dominated heat transfer effects. This study seeks to expand on the previous work by making time and space resolved measurements under bubbles nucleating on a micro-heater array operated under constant heat flux conditions. In the planned investigation, wall temperature measurements made under a single bubble nucleation site will be synchronized with high-speed video to allow analysis of the bubble energy removal from the wall.

Myers, Jerry G.↗

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗

An empirical determination of the polytropic index for the free-streaming solar wind using Helios 1 data

Observations of solar wind proton temperatures indicate that the solar wind is heated as it moves outward toward the orbit of Earth. This heating, which may be the results of electron heat conduction and perhaps MHD waves, has proven difficult to quantify and hence is often neglected in MHD models of the solar wind. An alternate approach to finding explicit heating terms for the MHD energy equation is to use a polytropic approximation. This paper discusses the properties of the polytropic approximation and its application to the solar wind plasma. By using data from the Helios 1 spacecraft, an empirical value for the polytropic index of the free-streaming solar wind is determined. Various corrections to the data are made to account for velocity, nonuniformity in radial sampling, and stream interaction regions. The polytropic index, as derived from proton data, is found to indepedent of speed state, within statistical error, and has an average value of 1.46. If magnetic pressure is included, the polytropic index has an average value of 1.58.

Totten, T. L.↗

Liquid Hydrogen Tank Chill and No-Vent Fill Prediction using Computational Fluid Dynamics

Cryogenic tank chill and fill is an important cryogenic fluid management technology that supports and enables many of NASA’s long-duration space missions. For no-vent fill, the receiver tank pressure remains below the supply tank pressure during the entire duration of the fill so that the tank does not require venting. This is especially advantageous for tank fill operations in low gravity where the position of the liquid is not always known and venting the tank may cause loss of propellant by venting liquid. In lieu of expensive tests conducted on-orbit, accurate computational models capable of predicting receiver tank pressure during cryogenic propellant tank fill may be used to reduce system and propellant mass as well as mission risk. However, these numerical models must be validated or anchored to test data. This study presents computational fluid dynamics (CFD) models with conjugate heat transfer that are used to predict a liquid hydrogen tank chill and no-vent fill ground test conducted at Lewis Research Center (now Glenn Research Center) in 1991. The specific test case chosen for model validation implemented an upward-facing jet near the bottom of a cylindrical 34-liter tank. CFD predictions are compared to experimental measurements of tank pressure, fill level, and wall temperature. The CFD results show reasonable agreement to the test data but overpredict the pressure collapse near the end of the fill despite the liquid jet penetrating the liquid-vapor interface.

Computational Fluid Dynamics↗

A model for analysis of the temperature field downstream of a heated jet injected into an isothermal crossflow at an angle of 90 deg

The temperature distribution downstream of a heated jet entering an isothermal crossflow at an angle of 90 deg is predicted using two conduction models with energy sources above the point of injection, in one case a point source and in the second a line source. The models use effective turbulent diffusivities that are determined empirically from previous measurements. Temperatures predicted by the models are compared to experimental results.

Eriksen, V. L.↗

Assessment of CFD Hypersonic Turbulent Heating Rates for Space Shuttle Orbiter

Turbulent CFD codes are assessed for the prediction of convective heat transfer rates at turbulent, hypersonic conditions. Algebraic turbulence models are used within the DPLR and LAURA CFD codes. The benchmark heat transfer rates are derived from thermocouple measurements of the Space Shuttle orbiter Discovery windward tiles during the STS-119 and STS-128 entries. The thermocouples were located underneath the reaction-cured glass coating on the thermal protection tiles. Boundary layer transition flight experiments conducted during both of those entries promoted turbulent flow at unusually high Mach numbers, with the present analysis considering Mach 10{15. Similar prior comparisons of CFD predictions directly to the flight temperature measurements were unsatisfactory, showing diverging trends between prediction and measurement for Mach numbers greater than 11. In the prior work, surface temperatures and convective heat transfer rates had been assumed to be in radiative equilibrium. The present work employs a one-dimensional time-accurate conduction analysis to relate measured temperatures to surface heat transfer rates, removing heat soak lag from the flight data, in order to better assess the predictive accuracy of the numerical models. The turbulent CFD shows good agreement for turbulent fuselage flow up to Mach 13. But on the wing in the wake of the boundary layer trip, the inclusion of tile conduction effects does not explain the prior observed discrepancy in trends between simulation and experiment; the flight heat transfer measurements are roughly constant over Mach 11-15, versus an increasing trend with Mach number from the CFD.

Wood, William A.↗

Model for energy transfer in the solar wind: Model results

A description is given of the results of solar wind flow in which the heating is due to (1) propagation and dissipation of hydromagnetic waves generated near the base of the wind, and (2) thermal conduction. A series of models is generated for fixed values of density, electron and proton temperature, and magnetic field at the base by varying the wave intensity at the base of the model. This series of models predicts the observed correlation between flow speed and proton temperature for a large range of velocities. The wave heating takes place in a shell about the sun greater than or approximately equal to 10 R thick. We conclude that large-scale variations observed in the solar wind are probably due mainly to variation in the hydromagnetic wave flux near the sun.

Barnes, A. A., Jr.↗

Liquid Hydrogen Tank Chill and No-Vent Fill Prediction using Computational Fluid Dynamics

Cryogenic tank chill and fill is an important cryogenic fluid management (CFM) technology that supports and enables many of NASA’s long-duration space missions. For no-vent fill, the receiver tank pressure remains below the supply tank pressure during the entire duration of the fill so that the tank does not require venting. This is especially advantageous for tank fill operations in low gravity where the position of the liquid is not always known and venting the tank may cause loss of propellant by venting liquid. In lieu of expensive tests conducted on-orbit, accurate computational models capable of predicting receiver tank pressure during cryogenic propellant tank fill may be used to reduce system and propellant mass as well as mission risk. However, these numerical models must be validated or anchored to test data. This study presents a computational fluid dynamics (CFD) model with conjugate heat transfer that is used to predict a liquid hydrogen tank chill and no-vent fill ground test conducted at Lewis Research Center (now Glenn Research Center) in 1991. The specific test case chosen for model validation implemented an upward-facing jet near the bottom of a cylindrical 34 liter tank. CFD predictions are compared to experimental measurements of tank pressure, fill level, fluid temperatures, and wall temperatures. The CFD results show reasonable agreement to the test data but overpredict the pressure collapse near the end of the fill despite the liquid jet penetrating the liquid-vapor interface. Several sensitivity studies are considered due to notable uncertainties in the experiment.

Computational Fluid Dynamics↗

Liquid Hydrogen Tank Chill and No-Vent Fill Prediction Using Computational Fluid Dynamics

Cryogenic tank chill and fill is an important cryogenic fluid management (CFM) technology that supports and enables many of NASA’s long-duration space missions. For no-vent fill, the receiver tank pressure remains below the supply tank pressure during the entire duration of the fill so that the tank does not require venting. This is especially advantageous for tank fill operations in low gravity where the position of the liquid is not always known and venting the tank may cause loss of propellant by venting liquid. In lieu of expensive tests conducted on-orbit, accurate computational models capable of predicting receiver tank pressure during cryogenic propellant tank fill may be used to reduce system and propellant mass as well as mission risk. However, these numerical models must be validated or anchored to test data. This study presents a computational fluid dynamics (CFD) model with conjugate heat transfer that is used to predict a liquid hydrogen tank chill and no-vent fill ground test conducted at Lewis Research Center (now Glenn Research Center) in 1991. The specific test case chosen for model validation implemented an upward-facing jet near the bottom of a cylindrical 34 liter tank. CFD predictions are compared to experimental measurements of tank pressure, fill level, fluid temperatures, and wall temperatures. The CFD results show reasonable agreement to the test data but overpredict the pressure collapse near the end of the fill despite the liquid jet penetrating the liquid-vapor interface. Several sensitivity studies are considered due to notable uncertainties in the experiment.

Computational Fluid Dynamics↗

Nuclear winter - Three-dimensional simulations including interactive transport, scavenging, and solar heating of smoke

A reexamination is conducted of the 'nuclear winter' hypothesis with a three-dimensional global model modified to allow for localized injection of smoke, its transport by the simulated winds, its absorption of sunlight, and its removal by model-simulated precipitation. Smoke injected into the troposphere is driven upward by solar heating. The tropopause, initially above the smoke, reforms below the heat smoke layer and separates it from precipitation below. Although much smoke is scavenged while the thermal structure is being altered, the residence time of the remaining smoke is greatly increased. Particularly for July conditions, a longer-lasting 'nuclear winter' effect is observed than was found in earlier modeling studies in which normal tropospheric residence times were assumed. In January the smaller solar flux in the northern hemisphere allows faster removal of smoke than in July. Significant cooling of the northern hemisphere continents is predicted; its dependence on season and injected smoke mass is described.

Malone, R. C.↗

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↗

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.↗

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↗

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.↗

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.↗

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.↗

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↗

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.↗