Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Flux”

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 901 records · Page 50

More Frequent and Persistent Heatwaves Due To Increased Temperature Skewness Projected by a High‐Resolution Earth System Model

Abstract Heatwaves are strongly associated with temperature distributions, but the mechanisms by which distributions are influenced by climate change remains unclear. Comparing simulations from a high‐spatial resolution Community Earth System Model (CESM‐HR) with those from low‐resolution models, we identify substantial improvements by CESM‐HR in reproducing observed Northern Hemisphere summer temperature skewness, as well as the frequency, intensity, persistence, and total heatwave days. Temperature skewness is strongly linked to land‐atmosphere interactions and atmospheric circulation. Under global warming projections, some regions exhibit enhanced temperature skewness, along with more frequent and persistent heatwaves of greater intensity. We find that in energy‐limited regimes, such as India, negative skewness in latent heat flux facilitates large positive skewness in sensible heat flux, which modulates near‐surface air temperatures. Skewness differences of latent and sensible heat fluxes are amplified under global warming, increasing the temperature skewness. We find that this contrasting flux mechanism is active in several heatwave‐prone regions.

Gao, Yang↗

Closed Gap Slug Calorimeter for Plasma Stream Characterization

Slug calorimeters are used in sheer and stagnation mode to characterize heat flux levels for high enthalpy streams. The traditional design features a gap between slug and holder, which can be of concern in these convective heat flux environments. The challenge is to develop a calorimeter that closes the gap to gas flow, but largely maintains thermal insulation of the slug. The work presented herein introduces two new slug calorimeter designs featuring a closed gap. This is done using either aerogel as a filler or press fitting the slug with a disk. The designs were verified and compared to the baseline calorimeter design under radiative heat flux. Building on this, the calorimeters were exposed to convective heat flux in the arc-jet facilities. Results from the new designs and conclusions on the impact of the gap in convective heat flux will be shown.

Nawaz, Anuscheh↗

Variations in lithium vapor cave performance predictions due to radial transport and recycling assumptions

The lithium vapor cave is a detached divertor design that uses a single private flux region baffle to contain a dense cloud of lithium vapor to dissipate heat flux. Plasma flows are created via fuel gas puffing in order to minimize lithium contamination of the main plasma. Significant modeling using the 2D edge code SOLPS-ITER has already been performed, predicting that sufficient target heat flux reductions (q$^{max}_{Target}$ MW m −2 ) with ( $n$ Li /$n$ e ) LCFS < 0.05 is possible in a case with 90 MW m −2 unmitigated heat flux in NSTX-U. Low heat flux and low upstream concentration was found with a variety of combinations of divertor geometries, target recycling coefficients, upstream plasma parameters, lithium evaporation locations and deuterium fueling locations, with variations in performance found for each design choice. However, the most universal uncertainty of SOLPS-ITER simulations has until now remained unaddressed systematically, namely the cross-field anomalous particle and heat diffusivities. This article aims to bound the uncertainty in lithium concentration prediction as a result of the assumed deuterium cross-field transport and recycling. For the simulations presented here, a factor of 2.1 increase in the upstream lithium density prediction across a factor of four decrease to the assumed deuterium radial particle diffusivity is found. This result is compared across different assumed deuterium recycling coefficients, known to reduce with lithium injection. Upstream lithium density is found to vary by a factor of 2.4 across a feasible range of recycling coefficients at similar $n$$^{OMP,sep}_{e}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Meteorological and Soil Data from Ecohydrology Sensor Towers at Pump House and Snodgrass Mountain in East River Watershed, Colorado, 2019-2025

This data package includes hourly meteorological and soil sensor data at eight ecohydrology monitoring sites in East River Watershed, Colorado as part of the Watershed Function Scientific Focus Area (WFSFA) research led by Lawrence Berkeley National Lab (LBNL). Four field sites were located on the hillslope of East River (ER) near Pump House (PH) at Mount Crested Butte (ER-PHS1 to 4), and the other four are in the Snodgrass Mountain (SG) area (SG-EHS5 to 8). In terms of vegetation cover, three sites are in montane grasslands (ER-PHS1, ER-PHS2, and SG-EHS5), three are below evergreen conifer canopy (ER-PHS3, SG-EHS6, and SG-EHS7), and two are below deciduous aspen canopy (ER-PHS4 and SG-EHS8). The monitoring period began in October 2019 at the East River sites, in October 2020 at SG-EHS5 and SG-EHS6, and in October 2021 at SG-EHS7 and SG-EHS8. In September 2024, all four East River sites were fully retired. The four Snodgrass Mountain sites remain active. Each site is equipped with a comprehensive suite of meteorological sensors on a tripod and soil sensors that measure weather, energy fluxes, and soil variables. This data package includes measurements from ten different types of sensors and up to thirteen individual sensors per site, including (1) a weather station (measurement height ranges from 2.8~3.8 meters (m) above ground), (2) a quantum sensor for photosynthetic active radiation (PAR) (2.4~3.3m), (3) a net radiometer (1.7~2.1m), (4) an infrared radiometer (1.6~2.2m), (5) a sonic distance sensor (1.5~1.9m), (6) a soil carbon dioxide (CO2) flux chamber (0m), (7) a soil heat flux plate (-0.05m below ground), (8) a soil oxygen sensor (-0.3m), (9) a soil water potential sensor (-0.3m), and (10) soil water content sensors at 3~4 depths (-1.15 ~ -0.1m). A total of twenty-three variables is reported in this data package, including (1) atmospheric variables: air temperature (TA), atmospheric pressure (PA), vapor pressure (VP), and vapor pressure deficit (VPD), (2) precipitation variables: rain precipitation (P) and snow depth (D_SNOW), (3) energy fluxes variables: four-component net radiation (NETRAD) (shortwave/longwave incoming/outgoing radiation, SW_IN, SW_OUT, LW_IN, LW_OUT), photosynthetic photon flux density (PPFD), and soil heat flux (G), (4) soil variables: soil water content (SWC), soil water potential (SWP), soil temperature (TS), soil bulk electrical conductivity (COND_SOIL), and soil gaseous oxygen concentration (O2_SOIL), (5) wind variables: two-dimensional wind speed (WS), gust speed (WS_MAX), and wind direction (WD), and (6) surface variables: surface infrared temperature (T_CANOPY) and soil CO2 flux (CO2_SOIL). Please see the Methods section for data processing and QA/QC steps taken to generate the hourly datasets. The following files are included in this data package (notes on version: v{x}-{y}, where x is the metadata version, and y is the data version, when applicable): (1) “metadata_site_v{x}-{y}.csv” - a site metadata file that summarizes location information of all sites, including site ID, description, coordinates, timeframe, elevation, and vegetation cover, (2) “metadata_instrument_v{x}-{y}.csv” - an instrument metadata file that summarizes sensor information of all sites, including sensor manufacturer and model, measurement height, and sampling and averaging interval of all variables, (3) "data_{SITE_ID}_v{x}-{y}.csv" - eight data files that contain hourly data of each site indicated by {SITE_ID} in the filename, (4) “/figure/data_{SITE_ID}_v{x}-{y}.png" - eight figures that help visualize data of each site indicated by {SITE_ID} in the filename, (5) “/photo/*” - photos of each site indicated by {SITE_ID} in the filename, and (6) four file level metadata (flmd.csv) and data dictionary (*_dd.csv) files that summarize file, header, column, and variable information of all files. Notes: (1) Measurement height: Each variable name is followed by conventional positional qualifiers “H_V_R”, where H indicates the relative horizontal positions of that specific variable, V the vertical positions, and R the replicates. In this data package, only the vertical qualifier V varies, and V increases from the highest vertical position (V=1) to the lowest. Variables with the same qualifier are not necessarily measured by the same sensor, and the same variable with the same qualifier across different sites are not necessarily measured at the same height. Please refer to “metadata_instrument.csv” for the sensor information and measurement heights, and whether a variable is measured below the canopy. (2) Variable availability: Snow depth is not available at ER-PHS3 and SG-EHS7. SWC, soil temperature, and soil bulk EC at the deepest depth (<-1m) are not available at SG-EHS6 and SG-EHS7. The missing value code for numeric variables is -9999, except for SWP. For SWP, the missing value code is +9999, because SWP values are negative. (3) Sampling frequency: Please refer to “metadata_instrument.csv” for the increase of sampling frequency of some variables from 30-min to 1-min at ER-PHS1 to 4 in July 2020. (4) Sensors: While the methods of each sensor are not detailed, all sensors are commercially available, and their methods can be found in their manuals. Please refer to “metadata_instrument.csv” for the sensor manufacturer and model information. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Ion thermal transport in the H-mode edge transport barrier on DIII-D

The power balance ion heat flux in the pedestal region on DIII-D increases and becomes increasingly anomalous (above conventional neoclassical) in experiments with higher temperature and lower density pedestals where the ion collisionality ($v^*_i$) is lowered toward values expected on ITER. Direct measurements of the main-ion temperature are shown to be essential on DIII-D when calculating the ion heat flux due to differences between the temperature of $D^+$ and the more commonly measured $C^{6+}$ impurity ions approaching the separatrix. Neoclassical transport calculations from NEO and non-linear gyrokinetic calculations using CGYRO are consistent with these observations and show that while neoclassical transport plays an important role, the turbulent ion heat flux due to ion scale electrostatic turbulence is significant and can contribute similar or larger ion heat fluxes at lower collisionality. Beam emission spectroscopy and Doppler backscattering measurements in the steep gradient region of the H-mode pedestal reveal increased broadband, long-wavelength ion scale fluctuations for the low $v^*_i$ discharges at the radius where the non-linear CGYRO simulations were run. Taken together, increased fluctuations, power balance calculations, and gyrokinetic simulations show that the above neoclassical ion heat fluxes, including the increases at lower $v^*_i$, are likely due to weakly suppressed ion scale electrostatic turbulence. These new results are based on world first inferred ion and electron heat fluxes in the pedestal region of deuterium plasmas using direct measurements of the deuterium temperature for power balance across ion collisionalities covering an order of magnitude from high $v^*_i$ values of 1.3 down to ITER relevant $v^*_i$ ~0.1.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Heating effects on jack pine pyrogenic organic matter properties from a pyrocosm study in 2022

This dataset contains data associated with the preprint “Fire removes preexisting pyrogenic organic matter from the ecosystem through the mechanisms of both direct combustion and increasing mineralizability” (Luo et al., 2025b), which is the complementary study to the published paper “Reburning pyrogenic organic matter: a laboratory method for dosing dynamic heat fluxes from above” (Luo et al., 2025a). We designed a full-factorial experiment with different burial depths of jack pine (Pinus banksiana Lamb) pyrogenic organic matter (PyOM) (Surface, 1 cm, and 5 cm) and different heat-flux profiles (High, Low, and Control) to examine how subsequent fires affect the properties of preexisting PyOM. We measured total carbon (C), pH, dissolved organic carbon (DOC), dissolved inorganic carbon (DIC), and mineralized C (as CO₂-C, from a 12-week incubation).We found that high heat flux and/or surface placement resulted in substantial direct C losses through combustion. Intermediate heat exposure produced both combustion losses and increases in DOC and mineralizability, which may have complex long-term implications: an increased dissolved fraction of PyOM may promote downward transport into mineral soils and potentially contribute to deeper, longer-term C storage, but it may also make PyOM more susceptible to microbial decomposition. Under the lowest heat flux and deepest burial, most PyOM was retained, and changes in DOC and C mineralization were minimal. Finally, PyOM pH, an important chemical property, decreased under low-temperature heating but increased under higher temperatures.We uploaded pH data for all samples (“pH_of_all_samples.csv”); pH and temperature-related data (peak temperature and degree hours) for samples in High and Low heat-flux treatments (“pH_vs_peakT_and_degree_hours_only_for_heated_samples.csv”); total C data (“CN_pct_C_stock_C_loss_in_samples.csv”); DOC and DIC data (“doc_dic.csv”); and mineralized C (CO₂-C) data (“CO2-C_all_original.csv”). Additional details can be found in the Methods & Sampling section.All datasets uploaded to ESS-DIVE are clearly labeled, cleaned, and include both raw and derived data, ready for reuse in other analyses. All analysis code and raw datasets are also available on GitHub: https://github.com/MengmengLuo/Fire-removes-preexisting-pyrogenic-organic-matter-from-the-ecosystem.

54 ENVIRONMENTAL SCIENCES↗

Design and Analysis of Liquid Lithium Plasma Facing Components

In this study, liquid metal (LM) on plasma-facing components (PFCs) creates a renewable protective cover enhancing power exhaust and protecting the solid substrate. Additionally, energy confinement improvement through particle pumping can be achieved in the case of liquid lithium (LL). A new variant of PFCs was recently introduced at Princeton Plasma Physics Laboratory (PPPL) where a porous wall is used to stabilize the LM surface, while using magnetohydrodynamic (MHD) drive to push the LM flow inside the component. This arrangement allows efficient heat exhaust, and enhanced control of the LM surface temperature, leading to spatial control of evaporation and condensation of LL on the plasma interface. This feature is particularly attractive when vapor shielding is introduced to allow heat flux redistribution. This system has the advantage that as the heat flux increases, the evaporation rate will increase, while decreasing heat flux will decrease the evaporation rate of lithium, ideally creating a feedback effect which could self-regulate the amount of lithium evaporated. Analytical and numerical models for LL PFCs were developed in PPPL. To calculate the target temperature distribution for the case of evaporation from the divertor, we apply an iterative process allowing two-way coupling between the fluid-kinetic analysis of plasma using SOLPS-ITER code and the flow and heat transfer analysis of the PFC using an analytical model. In the last stage, the results are validated using computational fluid dynamics (CFD) analysis with customized version of the CFX code from ANSYS. CFX was modified at PPPL to allow MHD analysis at the high magnetic field typical for fusion applications. Results of the analysis for NSTX-U tokamak conditions will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theoretical study of the electron temperature in the high-latitude ionosphere for solar maximum and winter conditions

The T(e) variation in the high-latitude ionosphere at altitudes between 120 and 800 km has been modeled for solar maximum, winter solstice, and strong magnetic activity conditions. The calculated electron temperatures are consistent with the plasma densities and ion temperatures computed from a time-dependent ionospheric model. Heating rates for both solar EUV and auroral precipitation were included. In general, the predicted UT variation of the electron temperature that results from the displacement between the magnetic and geographic poles is only a few hundred degrees. However, in sunlit trough regions, T(e) hot spots develop, and these hot spots show a marked UT variation, by as much as 2500 K. The dominant parameter controlling the T(e) variation above 200 km is the magnetospheric heat flux into the ionosphere, which is essentially unknown. For realistic values of the magnetospheric heat flux, the maximum electron temperature ranges from 5000 to 10,000 K at 800 km. A magnetospheric heat flux is particularly effective in enhancing trough electron temperatures. In general, the electron heat flux at high altitudes is uniquely related to the electron temperature and gradient, except on auroral field lines where thermoelectric heat flow is important.

Schunk, R. W.↗

DESIGN AND TESTING OF DIMES CARBON ABLATION RODS IN THE DIII-D TOKAMAK

We present the design of ATJ graphite rods developed for ablation experiments under high heat flux (up to 50 MW/m2) in the lower divertor of the DIII-D tokamak, a magnetic plasma confinement device. This work is motivated by the need to test ablation models relevant to carbon-based thermal shields used in high-speed spacecraft atmospheric entries, where the heat fluxes encountered can be comparable to those achieved in the DIII-D divertor plasma. Several different designs for the flowfacing side of the rod are analyzed, including "sharp nose," "blunt," and "concave". The last shape is studied for its potential to lower heat fluxes at the rod surface by increased radiation from trapped neutrals and reduced parallel plasma pressure. We also analyze the possibility of applying a thin (approximately 30 microns) layer of silicon carbide (SiC) to the exposed part of several carbon ablation rods to benchmark its erosion calculations and lifetime predictions. Such calculations are of interest as SiC represents a promising material for both thermal protection systems (TPS) and a fusion plasma-facing material (PFM). Preliminary results from the DIII-D rod ablation experiments are also discussed.

Orlov, Dmitri M.↗

Vaporization heat transfer of dielectric liquids on enhanced surfaces covered with screen wicks

Experiments were conducted to investigate the vaporization heat transfer characteristics for the dielectric liquid FC-72 on several wicking surfaces which may be used in flat-plate heat pipes. The wicking materials studied included microstructure enhanced surfaces and coarse surfaces covered with screen meshes. Experimental data for q versus deltaT curves and critical heat fluxes were obtained for the two different operating conditions of a heat pipe, evaporation, and shallow pool boiling. When the liquid level was above the heated surface, the height of the liquid level above the surface was varied from 0 to 10 mm. When the liquid level was below the heated surface, the distance from the liquid level to the edge of the surface was adjusted from 0 to 15 mm. Experimental results revealed that for shallow pool boiling when the heated surface was covered with a wire screen mesh, the heat transfer coefficient increased at lower heat fluxes but the critical heat flux (CHF) decreased for all the surfaces tested. In the case of evaporation, both CHF and the heat transfer coefficient increased as the microstructure surfaces were covered with screen meshes.

Gu, C. B.↗

Shock/shock interference on a transpiration cooled hemispherical model

Experimental results are presented which show the effectiveness of transpiration cooling in reducing the peak heat flux caused by an impinging shock on a bow shock of a hemispherical model. The 12-inch diameter hemispherical transpiration model with helium coolant was tested in the Calspan 48-inch Hypersonic Shock Tunnel at nominal Mach 12.1 and freestream unit Reynolds number of 0.33 x 10 to the 6th/ft. An incident shock wave, generated by a blunt flat-plate shock generator inclined at 10 deg to the freestream, intersected the bow shock of the model to produce shock/shock interference. The stagnation heat flux without coolant or shock/shock interference was about 1.6 times a smooth surface laminar prediction due to effective roughness of the coolant ejection slots. A coolant mass flux 31 percent of the freestream mass flux reduced the stagnation heat flux to zero without shock/shock interference. However, for the same coolant mass flux and with shock/shock interference the peak heat flux was only reduced 8.3 percent, even though the total integrated heat load was reduced.

Nowak, Robert J.↗

Insights on Simulating Summer Warming of the Great Lakes: Understanding the Behavior of a Newly Developed Coupled Lake-Atmosphere Modeling System

The Laurentian Great Lakes are the world's largest freshwater system and regulate the climate of the Great Lakes region, which has been increasingly experiencing climatic, hydrological, and ecological changes. An accurate mechanistic representation of the Great Lakes thermal structure in Regional Climate Models (RCMs) is paramount to studying the climate of this region. Currently, RCMs have primarily represented the Great Lakes through coupled one-dimensional (1D) column lake models; this approach works well for small inland lakes but is unable to resolve the realistic hydrodynamics of the Great Lakes and leads to inaccurate representations of lake surface temperature (LST) that influence regional climate and weather patterns. This work overcomes this limitation by developing a fully two-way coupled modeling system using the Weather Research and Forecasting model and a three-dimensional (3D) hydrodynamic model. The coupled model system resolves the interactive physical processes between the atmosphere, lake, and surrounding watersheds; and validated against a range of observational data. The model is then used to investigate the potential impacts of lake-atmosphere coupling on the simulated summer LST of Lake Superior. By evaluating the difference between our two-way coupled modeling system and our observation-driven modeling system, we find that coupled-lake atmosphere dynamics can lead to a higher LST during June-September through higher net surface heat flux entering the lake in June and July and a lower net surface heat flux entering the lake in August and September. The unstratified water in June distributes the entering surface heat flux throughout the water column leading to a minor LST increase, while the stratified waters of July create a conducive thermal structure for the water surface to warm rapidly under the higher incoming surface heat flux. This research provides insight into the coupled modeling system behavior, which is critical for enhancing our predictive understanding of the Great Lakes climate system.

54 ENVIRONMENTAL SCIENCES↗

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↗

Reconstruction of Thermal Protection System Aeroheating using a Green’s Function Approach

Inverse heat transfer (IHT) techniques are often used to reconstruct the surface heating conditions on spacecraft thermal protection systems (TPS) during atmospheric entry. Current IHT techniques for entry spacecraft applications, however, demand substantial computational resources, and are impractical for analyses such as uncertainty quantification and real-time health monitoring. In this paper, a Green’s function sensor fusion approach is used to reconstruct the TPS surface aeroheating conditions on experimental spaceflight and ground test systems from collocated temperature and heat flux sensors embedded in the TPS. The algorithm leverages Green’s functions to model the heat conduction within the spacecraft TPS and stabilizes the recovery of the surface heating condition using the direct heat flux sensor measurement. The algorithm is validated using arc-jet ground test data and applied to the reconstruction of the Mars 2020 backshell heating during Martian atmospheric entry. The performance of the algorithm is benchmarked against a current state-of-the-art IHT framework, FIAT_Opt. The Green’s function-based reconstruction algorithm recovers the net hot-wall heat flux absorbed by the TPS and the incident heat flux from the atmospheric entry environment in close agreement with FIAT_Opt. Notably, computation of the surface heating condition is completed in three orders of magnitude less time with the Green’s function sensor fusion approach using a consumer-grade PC, versus with FIAT_Opt running on a high performance computer cluster. The efficiency of the algorithm is leveraged to compute the uncertainty contributions of input parameters to the total uncertainty in reconstructed Mars 2020 backshell heating for the full atmospheric entry heat pulse. The sensitivity analysis uncovers that, at different times throughout the entry heat pulse, uncertainties in the TPS specific heat, thermal conductivity, and emissivity are all dominant drivers of the reconstruction uncertainty. These results demonstrate Green’s functions and sensor-fusion techniques as promising IHT approaches to reconstruct atmospheric entry environments from TPS-embedded measurements, and highlight how these techniques may give access to post-flight analyses previously hindered by the prohibitive cost of current methods.

Kenneth McAfee↗

Reconstruction of Thermal Protection System Aeroheating using a Green’s Function Approach

Inverse heat transfer (IHT) techniques are often used to reconstruct the surface heating conditions on spacecraft thermal protection systems (TPS) during atmospheric entry. Current IHT techniques for entry spacecraft applications, however, demand substantial computational resources, and are impractical for analyses such as uncertainty quantification and real-time health monitoring. In this paper, a Green’s function sensor fusion approach is used to reconstruct the TPS surface aeroheating conditions on experimental spaceflight and ground test systems from collocated temperature and heat flux sensors embedded in the TPS. The algorithm leverages Green’s functions to model the heat conduction within the spacecraft TPS and stabilizes the recovery of the surface heating condition using the direct heat flux sensor measurement. The algorithm is validated using arc-jet ground test data and applied to the reconstruction of the Mars 2020 backshell heating during Martian atmospheric entry. The performance of the algorithm is benchmarked against a current state-of-the-art IHT framework, FIAT_Opt. The Green’s function-based reconstruction algorithm recovers the net hot-wall heat flux absorbed by the TPS and the incident heat flux from the atmospheric entry environment in close agreement with FIAT_Opt. Notably, computation of the surface heating condition is completed in three orders of magnitude less time with the Green’s function sensor fusion approach using a consumer-grade PC, versus with FIAT_Opt running on a high performance computer cluster. The efficiency of the algorithm is leveraged to compute the uncertainty contributions of input parameters to the total uncertainty in reconstructed Mars 2020 backshell heating for the full atmospheric entry heat pulse. The sensitivity analysis uncovers that, at different times throughout the entry heat pulse, uncertainties in the TPS specific heat, thermal conductivity, and emissivity are all dominant drivers of the reconstruction uncertainty. These results demonstrate Green’s functions and sensor-fusion techniques as promising IHT approaches to reconstruct atmospheric entry environments from TPS-embedded measurements, and highlight how these techniques may give access to post-flight analyses previously hindered by the prohibitive cost of current methods.

Kenneth McAfee↗

Stagnation Point Convective Heating Correlations for Entry Into H2/He Atmospheres

New convective heating correlations have been developed over a broad range of entry conditions into H2/He atmospheres. The correlations are based on calculations of stagnation point convective heat fluxes over a large range of velocities, densities, and effective nose radii. The heat fluxes and corresponding flowfields are computed using NASA’s hypersonic computational fluid dynamics code, DPLR. Cases in this study have been computed using both a temperaturelimited, radiative equilibrium wall and a cold wall temperature boundary condition. The results for the scaled heat flux, as a function of the difference in freestream and wall enthalpies, show a clear deviation from the linear trend predicted by the method of Sutton and Graves. Therefore, updated non-linear correlations are proposed for the prediction of stagnation point convective heat fluxes. Generally, these updated correlations predict lower heat loads and therefore smaller required TPS thicknesses when compared to previous correlations for a range of test cases in the Neptune, Saturn, and Uranus atmospheres.

Alex T. Carroll↗

Uncertainty Analysis of Slug Calorimeters in the HyMETS Arc-Jet Facility

The objective of this work is to perform an uncertainty analysis of the deduced stagnation heat flux environment on a slug calorimeter for conditions that span the performance envelope of the Hypersonic Materials Environmental Test System arc-jet facility located at NASA Langley Research Center. Analytical solutions are developed for boundary-value problems on the slug element accounting for non-ideal effects, including spatial variation in the slug heat flux, multi-dimensional thermal conduction, and back-face losses, which departs from the state-of-the-art method derived from the American Society of Testing and Materials. Boundary-value problem definitions are informed by preliminary finite element thermal analysis of the slug calorimeter assembly (including both slug and housing) and just the slug element. The analytical solutions are presented in a general sense and in a truncated form from error analysis. Results are shown in optimizing and validating the analytical models against available slug back-face thermal data. The optimization results indicate that the appropriate epistemic uncertainty of the deduced stagnation heat flux on the slug calorimeter is at most±2.5% for both a high-and low-enthalpy test condition. In addition, a numerical approach is used to determine the aleatory (probabilistic) uncertainty component in the slug stagnation heat flux by applying a marching least-squares slope routine through the steady-state portion of the slug back-face thermal response. Results indicate a compromise between the number of samples and the filter frequency of slug back-face thermal data points when evaluating the standard deviation of the deduced stagnation heat flux statistics. When combining the mixed uncertainty, both aleatory and epistemic, the interval of uncertainty in the deduced stagnation heat flux is determined to be up to ±4%, which is at least a 60% reduction from the standard uncertainty used in the state-of-the-art method.

uncertainty↗

Stagnation Point Convective Heating Correlations for Entry Into H2/He Atmospheres

New convective heating correlations have been developed over a broad range of entry conditions into H2/He atmospheres. The correlations are based on calculations of stagnation point convective heat fluxes over a large range of velocities, densities, and effective nose radii. The heat fluxes and corresponding flow fields are computed using NASA’s hypersonic computational fluid dynamics code, DPLR. Cases in this study have been computed using both a temperature limited, radiative equilibrium wall and a cold wall temperature boundary condition. The results for the scaled heat flux, as a function of the difference in freestream and wall enthalpies, show a clear deviation from the linear trend predicted by the method of Sutton and Graves. Therefore, updated non-linear correlations are proposed for the prediction of stagnation point convective heat fluxes. Generally, these updated correlations predict lower heat loads and therefore smaller required TPS thicknesses when compared to previous correlations for a range of test cases in the Neptune, Saturn, and Uranus atmospheres.

Alex T. Carroll↗