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 325 records · Page 18

A novel 3-D pulsating heat pipe module for high heat-flux applications

This study proposed a novel pulsating heat pipes (PHP) for high-flux thermal module. According to the open literature, the maximum value of the heat flux is 23.7 W∙cm –2 handled by PHP, which is far behind industry requirements. Therefore, the present study proposes a novel 3D PHP suitable for high heat-flux applications. The PHP were experimentally tested for different working fluids, including methanol and water. Furthermore, the experiments were performed to analyze the influence of various operating parameters, such as filling ratio (25–60%), heat load (100–800 W), and air flow rate 100–250 cubic feet per minute (CFM). In this work, the results showed that this PHP can sustain a heat flux up to 70 W∙cm –2 with the maximum case temperature under 80 °C, which is approximately three times higher than the available literature. Furthermore, the performance of the PHP module was compared against the commercial module with wicked heat pipe. The PHP exhibited about 7% higher heat dissipation compared to the commercial module. Since the proposed module uses simple copper tubes in its construction over the wicked structure used in commercial module, the novel 3D PHP module provides additional economic and mass production capability.

3D pulsating heat pipe↗

A simple model for the interaction between vertical eddy heat fluxes and static stability

A numerical model for studying the interaction of vertical eddy heat fluxes and vertical temperature structure in midlatitude regions is described. The temperature profile for the model was derived from calculations of the equilibrium among heating rates in simplified representations of large-scale vertical eddy heat flux, moist convection and radiation. An eddy flux profile is calculated based on the quasi-geostrophic, liner baroclinic instability of a single wave. Model equilibrium states for summer and winter conditions are compared with observations, and the results are discussed in detail.

Gutowski, W. J., Jr.↗

Trends and Variations of Ocean Surface Latent Heat Flux: Results from GSSTF2c Data Set

Trends and variations of Goddard Satellite-based Surface Turbulent Fluxes (GSSTF) version 2c (GSSTF2c) latent heat flux (LHF) are examined. This version of LHF takes account of the correction in Earth incidence angle. The trend of global mean LHF for GSSTF2c is much reduced relative to GSSTF version 2b Set 1 and Set 2 for the same period 1988-2008. Temporal increase of GSSTF2c LHF in the two decades is 11.0%, in which 3.1%, 5.8%, and 2.1% are attributed to the increase in wind, the increase in sea surface saturated air humidity, and the decrease in near-surface air humidity, respectively. The first empirical orthogonal function of LHF is a conventional El Nino Southern Oscillation (ENSO) mode. However, the trends in LHF are independent of conventional ENSO phenomena. After removing ENSO signal, the pattern of LHF trends is primarily determined by the pattern of air-sea humidity difference trends.

global mean↗

A unique high heat flux facility for testing hypersonic engine components

A major concern in advancing the state-of-the-art technologies for hypersonic vehicles is the development of an aeropropulsion system capable of withstanding high thermal loads expected during hypersonic flights. Consequently, there is a need for experimental facilities capable of providing a high heat flux environment for testing compound concepts and verifying analyses. A hydrogen/oxygen rocket engine was developed to provide a high enthalpy/high heat flux environment for component evaluation. This Hot Gas Facility is capable of providing heat fluxes ranging from 200 (on flat surfaces) up to 8000 Btu per sq ft per sec (at a leading edge stagnation point). Gas temperatures up to 5500 R can be attained as well as Reynolds numbers up to 360,000 per ft. Test articles such as cowl leading edges, transpiration-cooled seals, fuel injectors, and cooled panel concepts can be evaluated with gaseous hydrogen as coolant. This facility and its configuration and test capabilities are discussed. Results from flow characterization experiments are also shown and their implications considered.

Melis, Matthew E.↗

Local Heat Flux Measurements with Single and Small Multi-element Coaxial Element-Injectors

To support NASA's Vision for Space Exploration mission, the NASA Marshall Space Flight Center conducted a program in 2005 to improve the capability to predict local thermal compatibility and heat transfer in liquid propellant rocket engine combustion devices. The ultimate objective was to predict and hence reduce the local peak heat flux due to injector design, resulting in a significant improvement in overall engine reliability and durability. Such analyses are applicable to combustion devices in booster, upper stage, and in-space engines with regeneratively cooled chamber walls, as well as in small thrust chambers with few elements in the injector. In this program, single and three-element injectors were hot-fire tested with liquid oxygen and gaseous hydrogen propellants at The Pennsylvania State University Cryogenic Combustor Laboratory from May to August 2005. Local heat fluxes were measured in a 1-inch internal diameter heat sink combustion chamber using Medtherm coaxial thermocouples and Gardon heat flux gauges, Injector configurations were tested with both shear coaxial elements and swirl coaxial elements. Both a straight and a scarfed single element swirl injector were tested. This paper includes general descriptions of the experimental hardware, instrumentation, and results of the hot-fire testing for three coaxial shear and swirl elements. Detailed geometry and test results the for shear coax elements has already been published. Detailed test result for the remaining 6 swirl coax element for the will be published in a future JANNAF presentation to provide well-defined data sets for development and model validation.

Jones, Gregg↗

Tree-Based Ensemble Learning Models for Wall Temperature Predictions in Post-Critical Heat Flux Flow Regimes at Subcooled and Low-Quality Conditions

Accurately predicting post-critical heat flux (CHF) heat transfer is an important but challenging task in water-cooled reactor design and safety analysis. Although numerous heat transfer correlations have been developed to predict post-CHF heat transfer, these correlations are only applicable to relatively narrow ranges of flow conditions due to the complex physical nature of the post-CHF heat transfer regimes. In this paper, a large quantity of experimental data is collected and summarized from the literature for steady-state subcooled and low-quality film boiling regimes with water as the working fluid in vertical tubular test sections. In addition, a low-quality water film boiling (LWFB) database is consolidated with a total of 22,813 experimental data points, which cover a wide flow range of the system pressure from 0.1 to 9.0 MPa, mass flux from 25 to 2750 kg/m 2 s, and inlet subcooling from 1 to 70 °C. Two machine learning (ML) models, based on random forest (RF) and gradient boosted decision tree (GBDT), are trained and validated to predict wall temperatures in post-CHF flow regimes. The trained ML models demonstrate significantly improved accuracies compared to conventional empirical correlations. To further evaluate the performance of these two ML models from a statistical perspective, three criteria are investigated and three metrics are calculated to quantitatively assess the accuracy of these two ML models. For the full LWFB database, the root-mean-square errors between the measured and predicted wall temperatures by the GBDT and RF models are 5.7% and 6.2%, respectively, confirming the accuracy of the two ML models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Deployment of Traditional and Hybrid Machine Learning for Critical Heat Flux Prediction in the CTF Thermal-Hydraulics Code

Critical heat flux (CHF) marks the transition from nucleate to film boiling, where heat transfer to the working fluid can rapidly deteriorate. Accurate CHF prediction is essential for efficiency, safety, and preventing equipment damage, particularly in nuclear reactors. Although widely used, empirical correlations frequently exhibit discrepancies when compared to experimental data, limiting their reliability in diverse operational conditions. Traditional machine learning (ML) approaches have demonstrated potential for CHF prediction but often suffer from limited interpretability, data scarcity, and insufficient knowledge of physical principles. Hybrid model approaches, which combine data-driven ML with base models, mitigate these concerns by incorporating prior knowledge of the domain. This study integrates an externally trained purely data-driven ML model and two hybrid models (using the Biasi and Bowring CHF correlations) within the CTF subchannel code via a custom Fortran framework. Performance was evaluated using two validation cases: a subset of the Nuclear Regulatory Commission (NRC) CHF database and the Bennett dryout experiments. In both cases, the hybrid models demonstrated significantly lower error metrics compared to conventional empirical correlations, with the best models often reducing relative error by about 5 percentage points. The pure ML model achieved comparable accuracy, outperforming the hybrid Biasi model in the NRC test case (3.3% versus 5.5% relative error) but exhibiting slightly higher error against the hybrid Bowring model in the Bennett test case (7.7% versus 6.1%). Trend analysis of error parity indicated that ML-based models reduced the tendency for CHF overprediction, improving overall accuracy. These results demonstrate that ML-based CHF models can be effectively integrated into subchannel codes and could potentially increase performance compared to conventional methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Convective Heat Transfer Scaling of Ignition Delay and Burning Rate with Heat Flux and Stretch Rate in the Equivalent Low Stretch Apparatus

To better evaluate the buoyant contributions to the convective cooling (or heating) inherent in normal-gravity material flammability test methods, we derive a convective heat transfer correlation that can be used to account for the forced convective stretch effects on the net radiant heat flux for both ignition delay time and burning rate. The Equivalent Low Stretch Apparatus (ELSA) uses an inverted cone heater to minimize buoyant effects while at the same time providing a forced stagnation flow on the sample, which ignites and burns as a ceiling fire. Ignition delay and burning rate data is correlated with incident heat flux and convective heat transfer and compared to results from other test methods and fuel geometries using similarity to determine the equivalent stretch rates and thus convective cooling (or heating) rates for those geometries. With this correlation methodology, buoyant effects inherent in normal gravity material flammability test methods can be estimated, to better apply the test results to low stretch environments relevant to spacecraft material selection.

Olson, Sandra↗

Data Efficiency Assessment of Generative Adversarial Networks for Critical Heat Flux Synthetic Data Generation

This study investigates the application of generative artificial intelligence techniques, particularly conditional generative adversarial networks (cGAN), in real-world engineering contexts, with a specific focus on synthetic data generation for critical heat flux (CHF). Utilizing a dataset comprising more than 20,000 real experimental CHF measurements, we conduct a series of experiments to examine cGAN’s behavior. These experiments encompass varying sizes of the training dataset, training cGAN on data from diverse experimental sources to generate new data on unseen experimental setups, and assessing the impact of excluding various input features on cGAN’s data generation accuracy. Our findings underscore the pronounced data dependency of cGAN for reliable performance, with decreased efficacy observed with smaller training dataset sizes. Notably, cGAN exhibits varying performance when trained on data from different experiments, with superior predictive capabilities observed for certain experiment sources compared to others. For instance, when cGAN was trained on data from Smolin et al.’s experiments or Zenkevich et al., it exhibited relatively good performance in generating the data from Becker et al., Kirillov et al., and Alekseev et al. experiments. In contrast, when trained with Alekseev et al.’s data and tasked with generating other experimental setups, cGAN showed notably poor performance. In both scenarios, cGAN’s performance was inferior compared to training on samples from all experiments concurrently. A feature importance analysis highlights the significant influence of parameters such as mass flux and heated length on accurate CHF generation, while other parameters like diameter and pressure have less impact. Inlet temperature is identified as a moderating factor by cGAN.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Estimation of the average surface heat flux over an inhomogeneous terrain from the vertical velocity variance

An indirect method of estimating the surface heat flux from observations of vertical velocity variance at the lower mid-levels of the convective atmospheric boundary layer is described. Comparison of surface heat flux estimates with those from boundary-layer heating rates is good, and this method seems to be especially suitable for inhomogeneous terrain for which the surface-layer profile method cannot be used.

Eilts, M. D.↗

Unified descriptor for enhanced critical heat flux during pool boiling of hemi-wicking surfaces

Boiling heat transfer is dictated by interfacial phenomena at the three-phase contact line where vapor bubbles form on the surface. Structured surfaces have shown significant enhancement in critical heat flux (CHF) during pool boiling by tailoring interfacial phenomena. This CHF enhancement has been primarily explained by two structural effects: roughness, which extends the contact line length at the bubble base, and wickability, the ability to imbibe liquid through surface structures by capillary pumping. Here, in this work, we show that CHF enhancement on structured surfaces cannot be described by roughness or wickability alone. This result was confirmed using systematically designed micropillar surfaces with controlled roughness and wickability. Further, we performed a scaling analysis and derived a unified descriptor, which represents the combined effects of thin film density and volumetric wicking rate. This unified descriptor shows a reasonable correlation with CHF values with our experiments and literature data. This work provides important insights in understanding the role of surface structures on CHF enhancement, thereby providing guidelines for the systematic design of surface structures for enhanced pool boiling heat transfer.

42 ENGINEERING↗

Intercomparison of Latent Heat Fluxes Over Global Oceans

Turbulent fluxes of momentum, moisture, and heat at the air-sea interface are essential for climate studies. Version 2 Goddard Satellite-based Surface Turbulent Fluxes (GSSTF2) has been derived from the Special Sensor Microwave/Imager (SSM/I) radiance measurements. This dataset, covering the period July 1987-December 2000 over global oceans, has a spatial resolution of 1 deg x 1 deg lat-long and a temporal resolution of 1 day. Turbulent fluxes are derived from the SSM/I surface winds and surface air humidity, as well as the 2-m air and sea surface temperatures (SST) of the NCEP/NCAR reanalysis, using a bulk aerodynamic algorithm based on the surface layer similarity theory. The GSSTF2 bulk flux model, and retrieved daily wind stress, latent heat flux, wind speed, and surface air humidity validate well with ship observations of ten field experiments over the tropical and midlatitude oceans during 1991-99. The global distributions of 1988-2000 annual- and seasonal-mean turbulent fluxes show reasonable patterns related to the atmospheric general circulation and seasonal variations. Latent heat fluxes and related input parameters over global oceans during 1992-93 have been compared among GSSTF1 (version 1), GSSTF2, HOAPS (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data), NCEP/NCAR reanalysis, and one based on COADS (Comprehensive Ocean-Atmosphere Data Set). Our analyses suggest that the GSSTF2 latent heat flux, surface air humidity, surface wind, and SST are quite realistic compared to the other four flux datasets examined. However, significant differences are found among these five flux datasets. The GSSTF2, available at http://daac.gsfc.nasa.gov/CAMPAIGN_DOCS/hydrology/hd_gsstf2.O.html, is useful for climate studies.

Chou, Shu-Hsien↗

Heat flux in a non-Maxwellian plasma

A hybrid numerical scheme is applied to solve the Landau equation for the electron distribution function over all velocity space. Evidence is presented for the first time of the degree and character of the failure of the classical Spitzer-Haerm heat flux approximation in a realistic solar coronal loop structure. In the loop model used, the failure is so severe at some points that the role of the heat flux in the plasma's energy balance is completely misinterpreted. In the lower corona the Spitzer-Haerm approximation predicts that the heat flux should act as an energy source, whereas the more accurate distribution functions calculated here show this to be an energy sink.

Ljepojevic, N. N.↗

Three-Dimensional Thermal Boundary Layer Corrections for Circular Heat Flux Gauges Mounted in a Flat Plate with a Surface Temperature Discontinuity

Three-dimensional Navier-Stokes computational fluid dynamics (CFD) analysis has been performed in an effort to determine thermal boundary layer correction factors for circular convective heat flux gauges (such as Schmidt-Boelter and plug type)mounted flush in a flat plate subjected to a stepwise surface temperature discontinuity. Turbulent flow solutions with temperature-dependent properties are obtained for a free stream Reynolds number of 1E6, and freestream Mach numbers of 2 and 4. The effect of gauge diameter and the plate surface temperature have been investigated. The 3-D CFD results for the heat flux correction factors are compared to quasi-21) results deduced from constant property integral solutions and also 2-D CFD analysis with both constant and variable properties. The role of three-dimensionality and of property variations on the heat flux correction factors has been demonstrated.

Kandula, M.↗

Viscous hydrodynamic instability theory of the peak and minimum pool boiling heat fluxes

Liquid viscosity was included in the Bellman-Pennington theory of the Taylor wave in a liquid vapor interface. Predictions of the most susceptible wavelength, and of the wave frequency, were made as a function of a liquid viscosity parameter and the Bond number. The stability of a gas jet in a viscous liquid was studied and the result is used to predict the peak heat flux on large horizontal heaters. Experimental measurements of the dominant Taylor wave and its growth rate were made during the film boiling of cyclohexanol on cylindrical heaters. The results bear out the predictions quite well. The thickness of the vapor blanket surrounding a cylindrical heater was measured and a correlation suggested. The effect of large fluxes of vapor volume on the dominant wavelength was also noted. Theoretical results of the peak heat flux are compared with the experimental data, and the effect of finite geometry of flat plate heaters on the peak heat flux is also discussed.

Dhir, V. K.↗

Peak pool boiling heat flux in viscous liquids

The stability of a gas jet in a surrounding viscous liquid is studied. An expression is developed for the critical velocity at which the jet becomes unstable in a returning viscous liquid. The stability analysis is made to correspond with the geometrical configuration of gas jets and liquid columns similar to those observed near the peak pool boiling heat flux. The critical velocity of the gas jet is then used to obtain the functional form of the peak heat flux on flat plates and cylindrical heaters. The expressions are compared with original observations of the peak heat flux in very viscous liquids on flat plate, and cylindrical, heaters at both earth-normal, and elevated, gravities.

Dhir, V. K.↗

Three-Dimensional Heat Flux and Thermal Analysis of Angled Tungsten Samples on DIII-D

ITER-grade tungsten and dispersoid-strengthened tungsten samples with the top surface angled at ~15° towards the incident plasma flux were exposed to 9 H-mode discharges with edge-localized modes (ELMs) in the lower divertor of DIII-D tokamak using the Divertor Material Evaluation System (DiMES). Surface damage included cracking and flaking of material on the two samples farthest away from the plasma strike point, and significant melting of the two samples closest to the strike point. Heat flux and thermal analysis tools new to DIII-D have been applied to better understand this material response and to help optimize the exposure conditions for future experiments. SMITER field-line tracing simulations based on IRTV data and EFIT equilibria estimate an average inter-ELM perpendicular heat flux, 𝑞⊥,𝑖nter−𝐸LM , on the angled surfaces of 10.1 – 19.6 MW/m² for a majority of the 9 discharges, increasing to 15.6 – 24.5 MW/m² for the single, higher-power shot where samples melted. Fast camera data showed shallow intra-ELM melting and re-solidification, which transitioned to bulk inter-ELM melting with melt motion in the 𝐽⃗ 𝑥 𝐵⃗ direction. About 50% of the protruding volume of the most affected sample was displaced via melt-motion. SIERRA thermal modeling software was able to reproduce an onset time of melting consistent with fast camera data and final sample conditions, within < 200 ms. Maximum surface temperatures of 3122 K and 2787 K are estimated for the samples farthest away from the strike point, while the closest samples achieve melting at 4067 ms and 4750 ms into the ~5000 ms plasma exposure. A +10% increase in both the SMITER 𝑞⊥,𝑖nter−𝐸LM calculations and the estimated ELM heat loads 𝑞⊥, 𝐸LM was required to achieve this result, which is within the uncertainty of the diagnostic data but likely accounts for non-ideal geometry effects plus other physics uncertainties not included in this first iteration of modeling. This work provided valuable estimates of the 3D temperature evolution to help better understand the observed surface morphology and internal recrystallization of samples, which are discussed in detail in a complementary manuscript [1]. Benchmarking efforts with more diagnosed DIII-D experiments are underway to further refine the SMITER and SIERRA models for DiMES. Future use of these tools will enable researchers to precisely target heat flux exposure conditions in DIII-D to test, but not exceed, the thermomechanical limitations of novel plasma-facing materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Heat flux observations and the location of the transition region boundary of Giacobini-Zinner

Electron heat flux observations and associated plasma phenomena upstream from Comet Giacobini-Zinner show many similarities to observations upstream from the earth's magnetosheath. Two similarities are discussed. Heat flux events are used to compute the transition region boundary location. The computed subsolar standoff distance of 4,000 km is commensurate with a comet gas production rate, G, of 3 x 10 to the 28th molecules/s in accord with ground-based determination of G.

Fuselier, S. A.↗