Results for “Heat Flux” — Page 48

NASA NTRS1999

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗

NASA NTRS1962

BOILING HEAT TRANSFER TO LIQUID HYDROGEN AND NITROGEN IN FORCED FLOW

Boiling heat transfer to liquid hydrogen and nitrogen was investigated experimentally. Results are presented from a study of bulk boiling inside a cylindrical tube under vertically upward forced-flow conditions. A 0.555-inch-inside-diameter and 16 (1/8)-inch-long electrically heated stainless-steel tube was used. The range of variables studied for hydrogen were mass velocity of 2850 to 17,000 pounds per hour per square foot, local heat flux of 3600 to 40,000 Btu per hour per square foot, inlet pressure of 30 to 74 pounds per square inch absolute, and inlet subcooling of 0° to 9° R. Nitrogen test conditions were mass velocity of 15,000 to 56,000 pounds per hour per square foot, local heat flux of 2300 to 40,000 Btu per hour per square foot, inlet pressure of 47 to 56 pounds per square inch absolute, and inlet sub cooling of 1° to 6° R. The axial distribution of the tube-wall temperatures is presented. A transition in the type of boiling heat transfer was obtained. The critical heat flux corresponding to this transition was determined over a range of flow and heating rates and local qualities. At specific combinations of flow and transition location, a range of critical-heat-flux values was obtained and maximum values were determined. The maximum critical heat flux increased with increasing fluid-flow rate and decreased with increasing length of tube before transition. Similar variations of the maximum critical heat flux have been reported for water. The tube-inner-wall temperatures upstream of transition were essentially uniform and were only slightly greater than the fluid saturation temperature. The Wall-temperature profiles downstream of transition generally resembled those obtained in film-boiling studies and appeared to be strongly dependent upon local quality at the point of transition. Maximum wall temperatures of 9000 and 18000 R were obtained with hydrogen and nitrogen, respectively. Fluctuations of pressure, flow rate, and temperature occurred during some of the boiling tests. Under some conditions, maximum critical-heat-flux values were attained during steady-state operation with fluctuations. In other cases the fluctuations became uncontrolled, and critical-flux values less than the maximum values were obtained upon restabilization of the test conditions. No measurable pressure drop across the test section was obtained at any condition.

LIQUID HYDROGEN↗

DOE OSTI2025

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Analysis of fast-ion losses measured in MAST-U via infrared thermography and a Fast Ion Loss Detector

Fast-ion losses need to be monitored to avoid damage to plasma facing components. In existing experimental devices, the scintillator-based fast-ion loss detector (FILD) is the most advanced diagnostic for measuring fast-ion losses. However, FILDs provide only local information about the losses. Infrared (IR) thermography can be used as a complementary tool for more global monitoring of the deposition of fast-ion losses on the wall, at the expense of no velocity-space resolution. IR cameras measure the temperature of the plasma facing components. This measurement, determined by a combined effect of the thermal plasma, radiation, neutrons and fast-ion losses, can be decomposed to infer the fast-ion load on the tokamak wall. In this manuscript, a workflow to estimate fast-ion losses via IR thermography is applied to the MAST-U spherical tokamak, using a 1D approximation to extract the experimental heat flux on the FILD front face from IR data. To numerically estimate the different contributions to this total heat flux, the field-line tracing environment SMITER is used to calculate the thermal plasma contribution, the orbit-following Monte-Carlo code ASCOT to estimate the fast-ion losses, and bolometry measurements for the radiation. To validate the workflow, two discharges, L-mode plasmas with low MHD activity, were executed using on and off-axis beams, respectively. The experimentally and numerically estimated heat flux are of the same order of magnitude for the on-axis heated scenario, with a strong dependence of the estimated fast-ion losses contribution on the fit to the kinetic profiles used as input. This is also true for the off-axis heated scenario, where the total numerically estimated heat flux is 2.1 or 1.3 times higher than the maximum experimentally estimated heat flux, depending on the ASCOT input used.

FILD↗

DOE OSTI2026

Concept of Dynamic Heat Insulation for Rotating Detonation Engines

This work introduces a new class of materials concept to dynamically reduce instantaneous heat fluxes in Rotating Detonation Engine (RDE) combustor chamber walls. The high-frequency and high-amplitude surface heat fluxes observed in RDEs arise from large instantaneous temperature differences between the detonation shockwave and chamber wall surface. These temperature gradients drive substantial energy loss and reduce the chamber gas pressure, ultimately limiting the cycle’s thermodynamic efficiency. This work introduces a concept for dynamically insulating the combustion chamber surfaces using surface layers or coatings with low thermal time scale. With such coatings, the surface temperature may follow the fluctuations of the cyclic detonation wave temperature, thus reducing the instantaneous heat flux therefore cycle-mean heat flux. To analyze these cyclic thermal phenomena, a one-dimensional analytical conduction solver was utilized with the capability to handle multilayered structures. Parametric modeling was performed using transient heat flux boundary conditions representative of a hydrogen–air RDE across a broad range of coating thermal properties and engine conditions. The coating effectiveness scaled with the product of thermal time constant and detonation wave frequency and the results were non-dimensionalized to guide future materials development. This strategy may offer benefits in increasing material survivability, reducing cooling requirements, and enhancing pressure gain.

heat transfer↗

DOE OSTI2023

Multi-scale Analysis of Surface Heterogeneity Induced Convection on Isentropic Coordinates

This research uses semi-idealized simulations to investigate multi-scale processes induced by the heterogeneity of soil moisture observed during the 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign. The semi-idealized simulations have realistic land heterogeneity, but large-scale winds are removed. Analysis on isentropic coordinates enables the tracking of circulation that transports energy vertically and facilitates the identification of the primary convective processes induced by realistic land heterogeneity. The isentropes associated with upward motion are found to connect the ground characterized by high latent heat flux to cloud bases directly over the ground with high sensible heat flux, while isentropes associated with downward motion connect precipitation to the ground characterized by high sensible heat fluxes. The mixing of dry, warm parcels ascending from the ground with high sensible heat fluxes and moist parcels from high latent heat regions leads to cloud formation. This new mechanism explains how soil moisture heterogeneity provides the key ingredients such as buoyancy and moisture for shallow cloud formation. We also found that the sub-mesoscale dominates upward energy transport in the boundary layer, while mesoscale circulations contribute to vertical energy transport above the boundary layer. Our novel method better illustrates and elucidates the nature of land atmospheric interactions under irregular and realistic soil moisture patterns.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2006

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/cm2 for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/cm2, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8- species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗

NASA NTRS2016

Preparation for Cassini Grand Finale - Validation of Saturn Ring Heating in TSS

Thermal analysis of spacecraft at Saturn have previously been restricted to three environmental heat fluxes: solar, Saturn planet IR, and Saturn planet albedo. Heat flux from Saturn’s rings were not available in any commercially available thermal analysis tools. While preparing for the Cassini Grand Finale, during which Cassini will travel extremely close to the rings of Saturn, it was determined the existing thermal model required updating to accurately predict heat fluxes from the rings. This is particularly important for Cassini’s Visual and Infrared Mapping Spectrometer infrared spectrometer (VIMS-IR), as it is cooled by a passive radiator which operates below the temperature of Saturn’s rings. This paper discusses the validation of the heat flux from solar, Saturn, and Saturn’s rings in Spacedesign’s Thermal Synthesizer System (TSS). The incorporation of heating from Saturn’s rings also led to the use of SPICE geometry to ensure accurate spacecraft positioning and pointing. Thermal model results are compared against historical heating events as well as predictive Grand Finale sequences.

Swanson, Kirsten S.↗

NASA NTRS2006

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/sq cm for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/sq cm, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8-species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗

NASA NTRSDate not supplied

Mars 2020 Reconstructed Aerothermal Environments and Design Margins

The Mars 2020 aeroshell's thermal protection system was nearly identical to the Mars Science Laboratory system that successfully completed its mission for the Curiosity rover's landing in 2012. It was predicted that, like Mars Science Laboratory, the Mars 2020 heatshield would experience boundary layer transition and that the thermal protection system would provide sufficient material thickness margins against the aerothermal environments, even after adding radiative heating that was not included for Mars Science Laboratory. The Mars 2020 flight instrumentation suite included heatshield sensors similar to Mars Science Laboratory and new backshell measurements. A full set of surface pressure and in-depth temperature data were collected during atmospheric entry. This paper provides an initial analysis of the flight data and focuses on the reconstructed total surface heat flux inferred from the measured temperatures. Turbulent boundary layer conditions again were observed on the Mars 2020 heatshield, and the temperatures at all eleven heatshield and six backshell thermocouple in-depth sensor locations were well within system capabilities due to favorable entry conditions and conservative uncertainties. New computational fluid dynamics results on the reconstructed entry trajectory are compared to the measured surface pressures and reconstructed total heat fluxes. The predicted heatshield surface pressures at six locations match the data qualitatively and quantitatively well, as expected. Smooth-wall laminar heating predictions prior to boundary layer transition fall above the reconstructed heating on the heatshield. After the observed boundary layer transition time, total heat flux based on algebraic turbulence model calculations generally match the heat flux trends relative to one another. The convective heat flux was predicted to constitute the majority of the total heating. On the backshell, smooth-wall laminar total heat flux predictions generally exceed the reconstructed total heating at six locations. At each of these locations, it is estimated that radiative heating provided almost all of the total heating. The Mars 2020 as-flown aerothermal environments were well below the design levels for all measurement locations due to conservative design assumptions and a stressing design trajectory.

Mars 2020↗

NASA NTRSDate not supplied

Mars 2020 Reconstructed Aerothermal Environments and Design Margins

The Mars 2020 aeroshell's thermal protection system was nearly identical to the Mars Science Laboratory system that successfully completed its mission for the Curiosity rover's landing in 2012. It was predicted that, like Mars Science Laboratory, the Mars 2020 heatshield would experience boundary layer transition and that the thermal protection system would provide sufficient material thickness margins against the aerothermal environments, even after adding radiative heating that was not included for Mars Science Laboratory. The Mars 2020 flight instrumentation suite included heatshield sensors similar to Mars Science Laboratory and new backshell measurements. A full set of surface pressure and in-depth temperature data were collected during atmospheric entry. This paper provides an initial analysis of the flight data and focuses on the reconstructed total surface heat flux inferred from the measured temperatures. Turbulent boundary layer conditions again were observed on the Mars 2020 heatshield, and the temperatures at all eleven heatshield and six backshell thermocouple in-depth sensor locations were well within system capabilities due to favorable entry conditions and conservative uncertainties. New computational fluid dynamics results on the reconstructed entry trajectory are compared to the measured surface pressures and reconstructed total heat fluxes. The predicted heatshield surface pressures at six locations match the data qualitatively and quantitatively well, as expected. Smooth-wall laminar heating predictions prior to boundary layer transition fall above the reconstructed heating on the heatshield. After the observed boundary layer transition time, total heat flux based on algebraic turbulence model calculations generally match the heat flux trends relative to one another. The convective heat flux was predicted to constitute the majority of the total heating. On the backshell, smooth-wall laminar total heat flux predictions generally exceed the reconstructed total heating at six locations. At each of these locations, it is estimated that radiative heating provided almost all of the total heating. The Mars 2020 as-flown aerothermal environments were well below the design levels for all measurement locations due to conservative design assumptions and a stressing design trajectory.

Mars 2020↗

NASA NTRS2023

Effects of Heating Configuration and Operating Parameters on Heat Transfer and Interfacial Physics of Microgravity Flow Boiling With Subcooled Inlet Conditions –Experiments Onboard the International Space Station

This study is part of the Flow Boiling and Condensation Experiment (FBCE), a collaborative effort between the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center. The FBCE fitted with the Flow Boiling Module (FBM) was launched to the International Space Station (ISS) in August 2021 and experiments were successfully performed from February to July 2022 to amass a large microgravity-flow-boiling database. This study is focused on heat transfer and flow visualization of microgravity flow boiling of n-Perfluorohexane in a rectangular channel of 5.0 mm height, 2.5 mm width (heated), and 114.6 mm length, with subcooled inlet conditions. High-speed-video photography is utilized to present flow patterns and temporal interfacial behavior. Heat transfer results are presented in the form of flow boiling curves and both parametric curves and streamwise profiles of wall temperature and heat transfer coefficient. Firstly, the parametric effects of mass velocity (199.4 – 3200.0 kg/m 2 s), inlet subcooling (0.2 – 46.0°C), and inlet pressure (124.2 – 176.7 kPa), on the aforementioned aspects are assessed for double-sided heating to establish them for a microgravity environment. Of these three parameters, mass velocity and inlet subcooling mostly determine the microgravity flow boiling behavior, while inlet pressure plays an insignificant role. Flow patterns for double-sided heating are more complex than those for single-sided heating due to interaction between the two vapor layers. Vapor interaction is minimized at high subcoolings and high mass velocities due to strong condensation offered by the subcooled bulk liquid layer separating them. Despite the different flow patterns, both single- and double-sided heating generally result in similar parametric trends and local heat transfer coefficients for similar operating conditions. Flow instabilities manifest as temporal flow anomalies and temperature oscillations, and their severity increases with increasing boiling number. Secondly, the effects of heating configuration are analyzed by comparing and contrasting several aspects of single- and double-sided heating data. The heat fluxes at which onset of nucleate boiling degradation (ONBD) and critical heat flux (CHF) occur are distinctly different for single- and double-sided heating. There exists a threshold inlet subcooling demarcating the dominance of flow acceleration and condensation effects in vapor removal from the near-wall region and replenishment of fresh liquid for boiling. Above the threshold, condensation from the near-wall region is dominant and single-sided heating yields higher heat fluxes, and below it, acceleration is dominant and double-sided yields higher heat fluxes. At mass velocity in the range of 200 – 2400 kg/m 2 s, the threshold inlet subcooling lies in the approximate range of 20 – 30°C (corresponding inlet quality of roughly -0.40 – -0.20).

Microgravity↗

NASA NTRS1993

Turbulent heat transfer from a sparsely vegetated surface - Two-component representation

The conventional calculation of heat fluxes from a vegetated surface involving the coefficient of turbulent heat transfer which increases logarithmically with surface roughness, is inappropriate such highly structured surfaces as desert scrub or open forest. An approach is developed here for computing sensible heat flux from sparsely vegetated surfaces, where the absorption of insolation and the transfer of absorbed heat to the atmosphere are calculated separately for the plants and for the soil. This approach is applied to a desert-scrub surface in the northern Sinai, for which the turbulent transfer coefficient of sensible heat flux from the plants is much larger than that from the soil below, as shown by an analysis of plant, soil, and air temperatures. The plant density is expressed as the sum of products (plant-height) x (plant-diameter) of plants per unit horizontal surface area. The solar heat absorbed by the plants is assumed to be transferred immediately to the airflow. The effective turbulent transfer coefficient k(g-eff) for sensible heat from the desert-scrub/soil surface computed under this assumption increases sharply with increasing solar zenith angle, as the plants absorb a greater fraction of the incoming irradiation. The surface absorptivity (the coalbedo) also increases sharply with increasing solar zenith angle, and thus the sensible heat flux from such complex surfaces is a much broader function of time of day than when computed under constant k(g-eff) and constant albedo assumptions.

Otterman, J.↗

NASA NTRS2017

Mars 2020 Entry, Descent, and Landing Instrumentation 2 (MEDLI2) Sensor Suite

The Mars 2020 Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite seeks to address the aerodynamic, aerothermodynamic, and thermal protection system (TPS) performance issues during atmospheric entry, descent, and landing of the Mars 2020 mission. Based on the highly successful instrumentation suite that flew on Mars Science Laboratory (MEDLI), the new sensor suite expands on the types of measurements and also seeks to answer questions not fully addressed by the previous mission. Sensor Package: MEDLI2 consists of 7 pressure transducers, 17 thermal plugs, 2 heat flux sensors, and one radiometer. The sensors are distributed across both the heatshield and backshell, unlike MEDLI (the first sensor suite), which was located solely on the heat-shield. The sensors will measure supersonic pressure on the forebody, a pressure measurement on the aftbody, near-surface and in-depth temperatures in the heatshield and backshell TPS materials, direct total heat flux on the aftbody, and direct radiative heating on the aftbody. Instrument Development: The supersonic pressure transducers, the direct heat flux sensors, and the radiometer all were tested during the development phase. The status of these sensors, including the piezo-resistive pressure sensors, will be presented. The current plans for qualification and calibration for all of the sensors will also be discussed. Post-Flight Data Analysis: Similar to MEDLI, the estimated flight trajectory will be reconstructed from the data. The aerodynamic parameters that will be reconstructed will be the axial force coefficient, freestream Mach number, base pressure, atmospheric density, and winds. The aerothermal quantities that will be determined are the heatshield and backshell aero-heating, turbulence transition across the heatshield, and TPS in-depth performance of PICA. By directly measuring the radiative and total heat fluxes on the back-shell, the convective portion of the heat flux will be estimated. The status of the current tools to perform the post-flight data analysis will be presented, along with plans for model improvements.

Hwang, Helen H.↗

NASA NTRS2017

Mars 2020 Entry, Descent, and Landing Instrumentation 2 (MEDLI2) Sensor Suite

The Mars 2020 Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite seeks to address the aerodynamic, aerothermodynamic, and thermal protection system (TPS) performance issues during atmospheric entry, descent, and landing of the Mars 2020 mission. Based on the highly successful instrumentation suite that flew on Mars Science Laboratory (MEDLI), the new sensor suite expands on the types of measurements and also seeks to answer questions not fully addressed by the previous mission. Sensor Package: MEDLI2 consists of 7 pressure transducers, 17 thermal plugs, 2 heat flux sensors, and one radiometer. The sensors are distributed across both the heatshield and backshell, unlike MEDLI (the first sensor suite), which was located solely on the heat-shield. The sensors will measure supersonic pressure on the forebody, a pressure measurement on the aftbody, near-surface and in-depth temperatures in the heatshield and backshell TPS materials, direct total heat flux on the aftbody, and direct radiative heating on the aftbody. Instrument Development: The supersonic pressure transducers, the direct heat flux sensors, and the radiometer all were tested during the development phase. The status of these sensors, including the piezo-resistive pressure sensors, will be presented. The current plans for qualification and calibration for all of the sensors will also be discussed. Post-Flight Data Analysis: Similar to MEDLI, the estimated flight trajectory will be reconstructed from the data. The aerodynamic parameters that will be reconstructed will be the axial force coefficient, freestream Mach number, base pressure, atmospheric density, and winds. The aerothermal quantities that will be determined are the heatshield and backshell aero-heating, turbulence transition across the heatshield, and TPS in-depth performance of PICA. By directly measuring the radiative and total heat fluxes on the back-shell, the convective portion of the heat flux will be estimated. The status of the current tools to perform the post-flight data analysis will be presented, along with plans for model improvements.

Hwang, Helen↗

DOE OSTI2022

Wettability and CHF limits of Accident-Tolerant nuclear fuel cladding materials in light water reactor conditions

Here, we present the results of experimental investigations aimed at evaluating the thermal–hydraulic performance of chromium-coated zircaloy, i.e., one of the most promising accident tolerant fuel (ATF) cladding material for light water nuclear reactors. Precisely, we investigate the wettability and critical heat flux (CHF) limits of chromium-coated and conventional zircaloy surfaces in prototypical reactor conditions. For both surface types, we measure the contact angle in a vapor-saturated atmosphere from ambient pressure to the operating pressure of pressurized water reactors (PWRs), i.e., ~15 MPa. We measure the ambient-pressure steady-state flow boiling CHF with a spatially uniform heat flux. We measure the high-pressure steady-state flow boiling CHF with a cosine shape heat flux (up to 20 MPa) and with a uniform heat flux (up to 15 MPa), also exploring the effect of CRUD deposits on the chromium-coated surface. Our results reveal that the chromium surface and the bare zircaloy surface have similar wettability and both become super-hydrophilic in PWR conditions, and that there is practically no difference in the steady-state CHF limits, both at low-pressure and high-pressure conditions, also when the chromium-coated surface is covered by a CRUD deposit. However, while the chromium-coating does not improve the CHF compared to the bare zircaloy surface, it improves the post-CHF behavior. The chromium coating prevents the reaction between zircaloy and steam, which results in the formation of a brittle zirconium oxide through the surface of the cladding. We also measure the transient CHF under exponentially escalating heat flux inputs of a nano-smooth and a rough surface mimicking a chromium-coated zircaloy cladding. Interestingly, the results of the transient heat flux tests suggest that the CHF limit for very short periods (i.e., fast transients) is independent of the surface finish, being the same for a rough chromium surface or a nano-smooth surface.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Two-layer transient heat transfer using impulse response methods

Solutions to the inverse heat conduction problem (IHCP) are methods that can be used to quantify surface heat flux in multi-layer materials for components in which there are limited subsurface (internal) temperature measurements, such as coated components. A critical consideration is to capture high frequency fluctuations using a practical heat flux sensor. To that end, this paper highlights key parameters for calculating accurate surface heat transfer. Specifically, this research extends the available solutions to the IHCP for multi-substrate structures through an impulse response methodology. The sensitivity of the impulse method was quantified with respect to practical measurements. Further, when compared to the inverse case, the impulse method resulted in lower errors when calculating surface heat flux over a range of conditions. Overall, this work provides a foundation for deducing heat flux from a subsurface heat flux sensor while maintaining a high-frequency response.

42 ENGINEERING↗

NASA NTRS1990

Heat Pipe Vapor Dynamics

The dynamic behavior of the vapor flow in heat pipes is investigated at startup and during operational transients. The vapor is modeled as two-dimensional, compressible viscous flow in an enclosure with inflow and outflow boundary conditions. For steady-state and operating transients, the SIMPLER method is used. In this method a control volume approach is employed on a staggered grid which makes the scheme very stable. It is shown that for relatively low input heat fluxes the compressibility of the vapor flow is low and the SIMPLER scheme is suitable for the study of transient vapor dynamics. When the input heat flux is high or the process under a startup operation starts at very low pressures and temperatures, the vapor is highly compressible and a shock wave is created in the evaporator. It is shown that for a wide range of input heat fluxes, the standard methods, including the SIMPLER scheme, are not suitable. A nonlinear filtering technique, along with the centered difference scheme, are then used for shock capturing as well as for the solution of the cell Reynolds-number problem. For high heat flux, the startup transient phase involves multiple shock reflections in the evaporator region. Each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe. Furthermore, shock reflections cause flow reversal in the evaporation region and flow circulations in the adiabatic region. The maximum and maximum-averaged pressure drops in different sections of the heat pipe oscillate periodically with time because of multiple shock reflections. The pressure drop converges to a constant value at steady state. However, it is significantly higher than its steady-state value at the initiation of the startup transient. The time for the vapor core to reach steady-state condition depends on the input heat flux, the heat pipe geometry, the working fluid, and the condenser conditions. However, the vapor transient time, for an Na-filled heat pipe is on the order of seconds. Depending on the time constant for the overall system, the vapor transient time may be very short. Therefore, the vapor core may be assumed to be quasi-steady in the transient analysis of a heat pipe operation.

Issacci, Farrokh↗