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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A Fast Algebraic Multigrid Solver and Accurate Discretization for Highly Anisotropic Heat Flux I: Open Field Lines

We present a novel solver technique for the anisotropic heat flux equation, aimed at the high level of anisotropy seen in magnetic confinement fusion plasmas. Such problems pose two major challenges: (i) discretization accuracy and (ii) efficient implicit linear solvers. We simultaneously address each of these challenges by constructing a new finite element discretization with excellent accuracy properties, tailored to a novel solver approach based on algebraic multigrid (AMG) methods designed for advective operators. We pose the problem in a mixed formulation, introducing the directional temperature gradient as an auxiliary variable. The temperature and auxiliary fields are discretized in a scalar discontinuous Galerkin space with upwinding principles used for discretizations of advection. We demonstrate the proposed discretization’s superior accuracy over other discretizations of anisotropic heat flux, achieving error 1000x smaller for anisotropy ratio of 10 9 , for closed field lines. The block matrix system is reordered and solved in an approach where the two advection operators are inverted using AMG solvers based on approximate ideal restriction, which is particularly efficient for upwind discontinuous Galerkin discretizations of advection. To ensure that the advection operators are nonsingular, in this paper we restrict ourselves to considering open (acyclic) magnetic field lines for the linear solvers. We demonstrate fast convergence of the proposed iterative solver in highly anisotropic regimes where other diffusion-based AMG methods fail.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

The effectiveness of D 2 pellet and gas injection in reducing intra-ELM tungsten erosion and heat flux in the DIII-D small angle slot divertor

Edge localized modes (ELMs) in H-mode plasmas can melt and erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes the use of D 2 pellet and gas injection for ELM mitigation on the DIII-D tokamak during the 2022 Small Angle Slot V-shaped Tungsten (W) (SAS-VW) divertor campaign, reducing W erosion and heat flux during ELMs. D α (656 nm) and WI (400.9 nm) filterscopes and Langmuir probes provide photon emission and electron density/temperature to estimate W atom erosion using the S/XB method, while surface eroding thermocouples measured ELM peak heat flux at the outer strike point (OSP). Thomson Scattering measurements of pedestal T e and n e provided input to predict W divertor erosion and heat flux during ELMs via the Free-Streaming plus Recycling Model (FSRM). While a greater D 2 mass injection rate decreased the ELM peak heat flux, the impact on W erosion was not monotonic. The average ‘large’ intra-ELM W erosion was lower for plasma shots with D 2 mass injection compared to the plasma shots without any D 2 mass injection in the SAS-VW divertor. However, plasma shots did not see significant changes in ‘large’ intra-ELM W divertor erosion as the D 2 mass injection rate increased. On average, the FSRM overestimated SAS-VW experimental intra-ELM W erosion by a factor of 6.9, but the overestimation is reduced to a factor of 3.6 with the implication of C deposition effects. Generally, experimentally measured and predicted quantities were worse for plasma shots with a lower D 2 mass injection rate and corresponding higher plasma carbon (C) impurity percentage (f C ). The discrepancies are postulated to be due to C/W material mixing, for which a simple analytic mixed-material model is presented. These results highlight the importance of incorporating and improving the robustness of a mixed-material layer model in the analysis of PFC erosion on present and future tokamak devices.

Carbon (C)↗

DOE OSTI2024

Development of a Hybrid Single/Two-Phase Capillary-Based Micro-Cooler using Copper Inverse Opals Wick with Silicon 3D Manifold for High-Heat-Flux Cooling Application

Previously, we reported two-phase capillary-based cooling using narrow (200 to 1000 ..mu..m) heater bridge copper inverse opal (CIO) wicks with heat flux level of 1400 Wcm -2 and superheat ~ 10 degrees C. Here, we demonstrate the area scaling of the proposed technology to a large-area micro-cooler for high-heat-flux cooling of microprocessors and power electronics. We developed a hybrid single/two-phase micro-cooler that relies on capillary-wicking in a 25-..mu..m-thick CIO with an open silicon microchannel 3D-manifold for liquid delivery and vapor extraction, achieving a high heat flux ~ 400 Wcm -2 over a heated area of 10 x 10 mm 2 . For a range of inlet water flowrates from 5 to 60 g(min) -1 , we achieved total thermal resistances and vapor qualities of 0.68 cm 2 degrees CW -1 to 0.2 cm 2 degrees CW -1 and 0.55 to 0.12, respectively. The flowrates are 10x smaller than those of conventional single- or two-phase microchannel cooling technology. The corresponding two-phase thermal resistances ranges from 0.05 to 0.02 cm 2 degrees CW -1 with temperature superheat of 8 to 6 degrees C, respectively. While the overall performance of the large-area (10 x 10 mm 2 ) capillary-based micro-cooler degraded compared to the previous demonstration of the technology for a heated area of 5 x 5 mm 2 , however, preliminary CFD modeling indicates that an improved manifold design will be able to achieve comparable performance.

capillary flow↗

DOE OSTI2021

Model-based real-time surface heat flux and temperature estimation for the DIII-D tokamak

A control-oriented model for monitoring of wall power flux densities on the DIII-D tokamak has been successfully implemented and validated experimentally. Future reactors will have to withstand severe steady state high heat flux loads on plasma-facing components (PFCs). Due to the difficulty of directly-measuring local heat fluxes on these components, monitoring and protection of PFCs during the plasma discharge can benefit from simplified physics-based real-time (RT) functional models to estimate and guide heat load control. As a first step into the development, a control-oriented model for monitoring of wall power flux densities and temperatures on DIII-D tokamak has been successfully implemented. The paper discusses the experimental demonstration and comparison of the 2-D model-based wall heat flux algorithm on the DIII-Dinner wall limiter (IWL) against infra-red camera heat flux measurements for limited plasma configurations. The paper also reports on the benchmarking of the field line tracing environment, SMITER, developed at ITER organization on DIII-D tokamak against experimental IR diagnostic data and the derivation of the component shaping weighting factors for the 2-D model-based approach. Here, the extension of the model-based approach for surface temperature estimation on the DIII-D IWL is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Capillary-Enhanced Two-Phase Micro-Cooler Using Copper-Inverse-Opal Wick with Silicon Microchannel Manifold for High-Heat-Flux Cooling Application

In this work, we demonstrate a two-phase capillary-fed boiling micro-cooler that consists of a ~ 25-..mu..m-thick copper inverse opal (CIO) porous wicking structure for high-heat-flux boiling and a silicon 3D-manifold for distributed liquid delivery and vapor extraction across a 0.5 cm x 0.5 cm heated area. At low inlet water mass flow rates of 1.5 to 1.9 g(min)-1, the micro-cooler displays nearly two-phase boiling with exit vapor quality ~ 1 and a high critical heat flux (CHF) of 253 to 320 W cm-2 with low superheat of ~ 10 degrees C resulting in a thermal resistance of boiling ~ 0.025 cm2 degrees C W-1 or heat transfer coefficient of 0.4 MW m-2 degrees C-1. For higher flow rates of 5, 10, and 15 g(min)-1, the micro-cooler exhibits a hybrid single-phase and two-phase cooling regime where the contribution of the sensible heat (single-phase) cooling is linearly added to that of the two-phase cooling. For the highest flow rate of 15 g(min)-1, the CHF is increased to ~ 500 W cm-2 resulting in an overall thermal resistance of ~ 0.18 cm2 degrees C W-1. However, the two-phase heat transfer effectiveness, which estimates the utilization level of the inlet mass flow rate for two-phase boiling, is reduced to ~ 0.11. To achieve the best cooling system performances, the micro-cooler must operate entirely within the two-phase boiling regime (exit vapor quality or two-phase heat transfer effectiveness ~ 1). Ideally, the "coolant" should be delivered near its saturation temperature (~ 100 degrees C for water), which provides significant advantages for the energy-efficient operation of data centers and power electronics. We present detailed analysis with Infrared and high speed camera images at various inlet flow rates and heat fluxes to understand complex heat transfer in the micro-cooler. Furthermore, a conjugate thermofluidic simulation model, which incorporates the physics of capillary-fed boiling in a porous copper layer, agrees well with the experimental data.

capillary-enhanced boiling↗

DOE OSTI2022

Development of coated heat flux gauges for fast responding measurements

Thermal systems often exhibit transient behaviors that have important implications for the operation of the system and can be difficult to predict. For these reasons, experimental testing is often required to ensure system durability requirements are achieved. One important parameter governing the survivability of components in hot, high-stress environments is the heat flux into the part that dictates the temperature distribution for the component. However, sensors required to experimentally characterize heat fluxes in extreme environments must also be resilient. Here, this study presents the development of coated heat transfer gauges capable of robust, high-frequency measurements in turbine research facilities. The addition of a protective coating increases the durability of the gauge, but inherent of that coating is the attenuation of high-frequency temperature penetrations. As a result, this study first outlines the use of analytical solutions to define a gauge design for a specific frequency range and heat transfer, ensuring that subsurface signals can be rectified to surface conditions through inverse methods. Then, the fabrication of polyimide substrate sensors with a parylene-F coating is described. Micro surface heaters added to the custom sensors were used to determine important geometric and thermal properties necessary to calculate accurate surface heat flux. Ultimately, this work shows increased sensor robustness in a turbine test bed and experimentally validates that the frequency response of the fabricated sensors meet the design intent.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Predicting critical heat flux using localized sensing at invisible vapor-liquid interfaces

Predicting critical heat flux (CHF) in two-phase electronics cooling systems remains a significant challenge due to the sudden onset of boiling crisis and the difficulty in directly visualizing vapor-liquid interfaces. Existing sensing methods rely on lagging temperature measurements, optically accessible systems, or spatially averaged signals that cannot pinpoint CHF initiation at localized high-heat-flux regions. Here, we report a planar capacitive sensing approach that enables real-time, localized detection of vapor-liquid interface dynamics for CHF prediction in boiling heat transfer. The capacitive sensor exploits the dielectric constant difference between liquid and vapor phases to capture bubble nucleation, growth, and departure dynamics with a temporal resolution down to 2 ms. The capacitive sensing reveals distinct signals across boiling regimes: from high-frequency fluctuations during strong nucleate boiling to low-frequency fluctuations with increased amplitudes when approaching CHF. The multi-sensor array experiments demonstrate real-time localized sensing, where each sensor responds exclusively to boiling in its immediate vicinity without crosstalk from neighboring regions. This non-intrusive sensing approach provides predictive rather than lagging sensing signals of CHF occurrence, offering predictive diagnosis of two-phase liquid cooling for the thermal management of high-power-density electronics.

CHF↗

DOE OSTI2022

Unstable turbulent channel flow response to spanwise-heterogeneous heat fluxes: Prandtl's secondary flow of the third kind

Turbulent secondary flows are defined as Prandtl's secondary flow of the first or second kind, the former produced by stretching and/or tilting of vorticity, the latter produced via spatial heterogeneity of Reynolds stresses. Both mechanisms are instantaneously active within inertia-dominated wall turbulence; Reynolds stress spatial heterogeneity is required for Reynolds-averaged secondary flows. Spanwise-variable surface roughness can induce turbulent stress spatial heterogeneity in the spanwise–wall-normal plane and provide sustenance for streamwise-aligned mean secondary flows. Herein, we demonstrate that turbulent secondary flows can also be sustained by spanwise variability in the surface heat flux in unstably stratified turbulent channels, defined hereafter as Prandtl's secondary flow of the third kind. Support for this mechanism is established with scaling arguments, while large-eddy simulation is used to model inertia-dominated channel turbulence responding to a lower boundary with uniform aerodynamic/hydrodynamic roughness but spanwise-variable surface heat flux. Transport equations for streamwise vorticity and turbulent kinetic energy, k, outline the conditions needed for third-kind production: shear and buoyancy production over the elevated heat flux regions necessitates lateral entrainment of low- k fluid, inducing mean counter-rotating secondary cells aligned such that upwelling and downwelling occur over the high and low heat flux regions, respectively. Buoyancy-driven production of k alters aggregate flow response and thus is a distinctly different mechanism responsible for sustenance of secondary flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

A coupled ground heat flux–surface energy balance model of evaporation using thermal remote sensing observations

Abstract. One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Bhattacharya, Bimal K.↗

DOE OSTI2020

Impact of plasma density/collisionality on divertor heat flux width

Both ASDEX-Upgrade (AUG) data and the generalized HD (GHD) model showed that the scrape-off width broadens as the density/collisionality increases [1, 2]. A series of BOUT++ transport simulations are performed to study the physics of the scaling characteristics of the divertor heat flux width vs density/collisionality via a plasma density scan with either fixed pressure profile or fixed temperature profile inside separatrix. Additionally, the simulations show that even in the drift dominated regime, the divertor heat flux width can be broadened due to the transition of the SOL residence time from the parallel particle flow time to the enhanced parallel conduction time as the collisionality/density increases as posited in the GHD model. In addition, the heat flux width is found to be proportional to the square root of ion mass for low collisionality while it has a weakly dependence on ion mass for high collisionality. Furthermore, our simulations show that as the density increases, the radial electric field (E r ) well shallows, which potentially weakens E r × B flow shear stabilization of turbulence at high density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code

Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun.228 69–82). Here the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χ c and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λ q . From transport simulations, we obtain the critical thermal diffusivity χ c = 0.5 m 2 s –1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current I p = 5 MA, and χ c = 0.3 m 2 s –1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current I p = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity ${\chi }_{\text{c}}\propto {A}^{1/2}/(Z{\left(1+Z\right)}^{\frac{1}{2}}{B}_{\text{p}}^{2})$ using ITER scenarios with fixed safety factor q 95 , major radius R, aspect ratio R/a, and the separatrix temperature T sep , and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q 95 , R, R/a, and T sep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χ c , possibly yielding a transition from turbulence- to drift-dominant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Two-color pyrometry based flame to fuel surface radiative heat flux diagnostic using flamelets

In this study, a non-intrusive two-color pyrometry based radiative heat flux diagnostic is presented that accounts for the self-absorption effects. This diagnostic relies on looking up the flame temperature and soot volume fraction using the ratio of color intensities obtained from digital cameras and numerical solutions of 1D steady-state diffusion flames. Virtual two color pyrometry of computed flame show a unique one-to-one mapping of strained flames with intensity ratio; thereby allowing flames to be indexed and creation of a two-color flamelet manifold (TCFM). The TCFM is a function of two variables; the ratio of red to green intensities (flame index) and physical distance to the stoichiometric surface defined using a level-set function. The developed diagnostic is used in upward flame spread experiments where an approximate flame hull reconstruction method is applied to provide the 3D level-set function and intensity ratio for use of TCFM to provide radiative absorption properties for 3D ray tracing. Accounting for self- absorption effects provide improved soot volume fraction estimates. In conclusion, the radiative heat flux estimates are found to be in very good agreement with literature and improved results are obtained near the flame base.

42 ENGINEERING↗

DOE OSTI2020

Subcloud and Cloud-Base Latent Heat Fluxes during Shallow Cumulus Convection

Doppler and Raman lidar observations of vertical velocity and water vapor mixing ratio are used to probe the physics and statistics of subcloud and cloud-base latent heat fluxes during cumulus convection at the ARM Southern Great Plains (SGP) site in Oklahoma, United States. The statistical results show that latent heat fluxes increase with height from the surface up to ~0.8 Z i (where Z i is the convective boundary layer depth) and then decrease to ~0 at Z i . Peak fluxes aloft exceeding 500 W m −2 are associated with periods of increased cumulus cloud cover and stronger jumps in the mean humidity profile. These entrainment fluxes are much larger than the surface fluxes, indicating substantial drying over the 0–0.8 Z i layer accompanied by moistening aloft as the CBL deepens over the diurnal cycle. We also show that the boundary layer humidity budget is approximately closed by computing the flux divergence across the 0–0.8 Z i layer. Composite subcloud velocity and water vapor anomalies show that clouds are linked to coherent updraft and moisture plumes. The moisture anomaly is Gaussian, most pronounced above 0.8 Z i and systematically wider than the velocity anomaly, which has a narrow central updraft flanked by downdrafts. This size and shape disparity results in downdrafts characterized by a high water vapor mixing ratio and thus a broad joint probability density function (JPDF) of velocity and mixing ratio in the upper CBL. We also show that cloud-base latent heat fluxes can be both positive and negative and that the instantaneous positive fluxes can be very large (~10 000 W m −2 ). However, since cloud fraction tends to be small, the net impact of these fluxes remains modest.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

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↗

DOE OSTI2023

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↗

DOE OSTI2026

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↗

DOE OSTI2024

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↗

DOE OSTI2021

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↗