Abstract

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. 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/m2 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.

1 Introduction

Post-critical heat flux (CHF) heat transfer is a critical phenomenon in the analysis of various industrial applications, such as cryogenic applications [1], metallurgical processing [2], conventional boilers [3], and water-cooled nuclear reactors [4]. The post-CHF regimes involve complicated flow and heat transfer processes. As an example, during the reflooding phase of a hypothetical loss of coolant accident in a light water-cooled nuclear reactor (LWR), post-CHF heat transfer may be encountered due to high wall temperatures of nuclear fuel rods. To prevent fuel rods from overheating, subcooled water from the reactor emergency core cooling system is injected into the reactor core. Inverted annular film boiling (IAFB) will usually occur just downstream of the quench front, where a thin vapor film develops on the fuel rods and a liquid core exists in the central region of the flow channel. As the void fraction increases and the liquid core breaks up into liquid slugs, the flow regime transitions to the inverted slug film boiling (ISFB) regime (for low mass flux conditions) or agitated inverted annular film boiling (AIAFB) regime (for high mass flux conditions) [5]. The IAFB and ISFB/AIAFB regimes are regarded as low-quality film boiling regimes. Further downstream, the liquid slugs eventually break up into small droplets and the flow enters a high-quality regime, i.e., dispersed flow film boiling (DFFB). These above-mentioned post-CHF processes are caused by the departure from nucleate boiling (DNB) mechanism, which usually occurs at subcooled and low-quality inlet conditions. Another major mechanism that leads to post-CHF is the dryout of a liquid film that contacts the heated wall at high-quality conditions, which directly leads to the DFFB regime without going through the IAFB or ISFB/AIAFB regimes. A common feature of these post-CHF regimes is that the heat transfer is significantly deteriorated, causing high wall temperatures that can potentially lead to damage or failure of the heated surface. Therefore, accurately predicting wall temperatures in the post-CHF regimes is of high importance in the design and safety analysis of light water-cooled nuclear reactors and a variety of other industrial applications.

In recent years, many experimental and theoretical studies have been performed on the heat transfer characteristics in the post-CHF flow regimes [611]. A commonly used method to calculate the heated wall temperatures is to use empirical heat transfer correlations, which are typically high-order nonlinear functions of multiple empirical parameters, including the flow channel geometry, fluid properties, and flow conditions. These correlations can be characterized into three groups: (1) modified Bromley's model [12], (2) modified Dittus-Boelter model [13], and (3) laminar flow model [10,14]. These heat transfer correlations are useful and convenient for engineering applications and, to some extent, reveal physical mechanisms of the film boiling phenomena. However, due to the complex physical nature and flow regime transitions in the post-CHF flow regimes, these correlations are usually complicated functions with multiple empirical parameters that are calibrated based on experimental measurements. As a result, most of these correlations are only applicable to relatively narrow ranges of flow conditions and cannot be readily extended to beyond their range of validity, limiting their applications in practice.

Considering the thermal nonequilibrium characteristics in film boiling regimes, one-dimensional (1-D) two-fluid model is widely applied to predict wall temperature and void fraction in post-CHF flow regimes by solving the conservation equations for the liquid and vapor phases [1518]. To solve the 1-D two-fluid model, closure models, including those for the wall heat transfer and friction, interfacial heat transfer and friction/shear/drag, etc., are needed as source terms to make the conservation equations solvable. However, due to the lack of local flow information in the post-CHF flow regimes, the empirical closure models currently used in the relevant computer codes, such as system-level and subchannel analysis codes, are usually extended from generally encountered two-phase flows. Such a mismatch in the closure models and flow regimes could cause large uncertainties in the simulations that involve post-CHF flow regimes. A robust and generalized model to predict the heat transfer and wall temperature in post-CHF regimes with good accuracy is, therefore, highly desired.

With the advances in algorithm and computing capabilities, data-driven modeling enabled by machine learning (ML) has received increasing attention and has been widely applied to classification and regression problems in many industrial applications in recent years [1922]. A variety of ML methods have been used for data-driven modeling in engineering applications, including support vector machines (SVMs), Gaussian processes (GPs), decision trees (DTs), and artificial neural networks (ANNs). These ML methods have unique advantages in complex and nonlinear engineering problems over the conventional empirical correlation method. Instead of first identifying a fixed function form and then calibrating a few empirical parameters/coefficients, ML methods can directly learn from large datasets and make predictions with good accuracy and generalizability. In the field of nuclear thermal hydraulics, ML methods have been widely used to analyze a large amount of data obtained from experiments and/or numerical simulations for heat transfer problems [2325]. ANNs and SVMs are two common ML methods used in this field [26,27]. However, both methods have some shortcomings. For example, it is difficult to determine the architecture and hyperparameters of ANNs, such as the number of hidden layers and the number of neurons in each hidden layer. In SVMs, the kernel function needs to be designed correctly to achieve good prediction results, however, selecting the right kernel could be tricky and largely depends on user's experience.

Compared to training a ML model with sophisticated architecture and many learnable parameters, such as deep neural networks, an alternative approach is to use the ensemble methods. The ensemble methods combine the predictions of multiple base estimators with relatively simple structure to improve generalizability and robustness over a single estimator. There are two groups of ensemble methods, i.e., averaging method (or bagging method) and boosting method. In the averaging method, the principle is to build many base estimators independently and form the final prediction by averaging their predictions with reduced variance, which is usually better than a single base estimator. In the boosting method, base estimators are built sequentially and fit on the residual of the former estimators to reduce the bias of the combined estimator. Random forest (RF) and gradient boosted decision trees (GBDT) are, respectively, representatives of the averaging method and boosting method, which combine a set of base estimators, i.e., DTs, to make more accurate predictions over a single DT. The model trained by the ensemble methods is robust and can give smoother predictions compared to a single complicated ML model, making it ideal for many engineering applications.

To the authors' best knowledge, although ML methods have been widely used for dryout point and CHF predictions [2426], limited ML-related research work has been performed on wall heat transfer problems in the post-CHF flow regimes, especially at subcooled and low-quality conditions. In this paper, RF and GBDT models are developed, as a first attempt, to predict wall temperatures in the post-CHF flow regimes. We collected and summarized 22,813 experimental data points obtained from steady-state subcooled and low-quality film boiling experiments with water as the working fluid in vertical tubular test sections in the literature and consolidated them into a Low-quality Water Film Boiling (LWFB) database, which is then compared with three conventional heat transfer correlations for the film boiling regimes and is also used to train and validate the ML models. In addition, a statistical analysis is performed to evaluate the performance of the ML models by analyzing the relative errors between the measured and predicted wall temperatures for the GBDT and RF models.

2 Data Description

2.1 Experimental Methods.

To perform steady-state post-CHF experiments with water as the working fluid, especially at subcooled and low-quality conditions, Groeneveld [28] first applied an indirectly heated hot patch technique to achieve steady-state film boiling phenomena in a tubular test section in 1974. Then in 1984, Chen and Li [29] proposed a directly heated hot patch technique for their steady-state film boiling experimental investigation. Schematics of these two hot patch techniques are shown in Fig. 1. The test section is a Joule heated round tube. The indirectly heated hot patch is a copper block with high thermal inertia and a set of electrical cartridge heaters are included for the hot patch. In the directly heated hot patch technique, the wall thickness of the tubular test section is reduced over a short length near the tube ends to increase the local electrical resistance of the tube and therefore, increase the local heat flux at the reduced thickness section, which is called notch. In both techniques, high heat fluxes are realized for the hot patch, just upstream of the film boiling section, to form a dry patch in the hot patch and stabilize the quench front. The dry patch in the hot patch propagates to the downstream of the hot patch so that steady-state film boiling can be achieved in the test section. The function of the hot patch is to prevent the quench front from propagating along the test section and rewetting of the heated tube wall downstream of the quench front. With this technique, the average heat flux values required in the test section are relatively low and steady-state film boiling can be formed with the average wall heat flux considerably below the CHF values, which can avoid overheating the test section.

Fig. 1
Schematics of two hot patch techniques: (a) indirectly heated hot patch technique and (b) directly heated hot patch technique
Fig. 1
Schematics of two hot patch techniques: (a) indirectly heated hot patch technique and (b) directly heated hot patch technique
Close modal

2.2 Data Summary.

Both the indirectly heated and directly heated methods have been widely used in the literature for steady-state film boiling experiments with water as the working fluid in vertical tubular test sections. Relevant experimental data have been collected from steady-state post-CHF experiments that were performed in the 1980s and 1990s, available in the literature. All the experimental data include the wall temperature information but only part of the data have the void fraction information available. Table 1 summarizes some of these steady-state post-CHF experiments with water as the working fluid at subcooled and low-quality conditions [3040], including the geometrical parameters (inner diameter, ID; outer diameter, OD; and test section length), flow conditions (system pressure, p; mass flux G; inlet subcooling ΔTsub,in; wall temperature Tw; and wall heat flux qin), and the type of hot patch techniques. Tables 2 and 3 summarize the test matrices of these experiments. We collected a total of 22,813 experimental data points from the nine data sources in the accessible literature and consolidate them into our LWFB database, which covers a wide flow range of the system pressure from 0.1 to 9.0 MPa, mass flux from 25 to 2750 kg/m2 s, and inlet subcooling from 1 to 70 °C.

Table 1

Summary of some previous steady-state subcooled and low-quality water film boiling experiments

Refs.ID/OD (mm)Length (mm)Inlet subcooling (°C)Mass flux (kg/m2 s)System Pressure (MPa)Wall temperature (°C)Wall heat flux (kW/m2)Hot patch
Fung [30]11.81/12.70750–7801–7050–5000.1357 to 113924 to 258Indirect
11.94/13.08
Stewart [31]8.93/11.14197010–50120–27502.0–9.0365 to 79060 to 439
Costigan, et al. [32]9.85/14.09202–5025–1500.2, 0.4421 to 10874 to 139
Mosaad [33]9.0/12.02800–3050–5000.1485 to 115319 to 266Direct
Chen et al. [34,35]7.0/10.02,2000–2025–15000.1–6.058 to 223372 to 712
12.0/15.0
Swinnerton et al. [36,37]9.75/14.09205–2050–10000.2–2.05 to 380343 to 840Indirect
Savage et al. [38]9.75/13.979205–1550–20009 to 324351 to 862
Savage et al. [39]9.4/12.6892020–50200–10000.5–7.068 to 375520 to 927Direct
Savage et al [40]9.4/12.68920100–100052 to 345418 to 978
Refs.ID/OD (mm)Length (mm)Inlet subcooling (°C)Mass flux (kg/m2 s)System Pressure (MPa)Wall temperature (°C)Wall heat flux (kW/m2)Hot patch
Fung [30]11.81/12.70750–7801–7050–5000.1357 to 113924 to 258Indirect
11.94/13.08
Stewart [31]8.93/11.14197010–50120–27502.0–9.0365 to 79060 to 439
Costigan, et al. [32]9.85/14.09202–5025–1500.2, 0.4421 to 10874 to 139
Mosaad [33]9.0/12.02800–3050–5000.1485 to 115319 to 266Direct
Chen et al. [34,35]7.0/10.02,2000–2025–15000.1–6.058 to 223372 to 712
12.0/15.0
Swinnerton et al. [36,37]9.75/14.09205–2050–10000.2–2.05 to 380343 to 840Indirect
Savage et al. [38]9.75/13.979205–1550–20009 to 324351 to 862
Savage et al. [39]9.4/12.6892020–50200–10000.5–7.068 to 375520 to 927Direct
Savage et al [40]9.4/12.68920100–100052 to 345418 to 978
Table 2

Test matrices for the subcooled and low-quality film boiling experiments identified in Table 1 (0.1–0.7 MPa)

p (MPa)
ΔTsub,ina (°C)0.10.20.250.40.50.60.7
G (kg/m2 s)255,10,20,35,505,10,20,35,50,60
501,5,10,15,20,30,50,60,702,5,20,35,505,15,20,35,40,505
1001,2,5,10,15,20,30,40,50,605,20,355,15,20,35,405,10,15,20,25,30,40
1251,5,10,15
1505,10,20,30,405,205,15,20,30,4030,40
17015,25
2001,5,10,15,20,30,4055,10,20,25,35,45
2505,10,20
3001,5,10,20
3505,10,20
4001,5,10,20,252,25,30,35
45030-
5001,5,10,20,30510155,30,35
600
7005,10
75030
10005,105,305,10
150020
200010,155
p (MPa)
ΔTsub,ina (°C)0.10.20.250.40.50.60.7
G (kg/m2 s)255,10,20,35,505,10,20,35,50,60
501,5,10,15,20,30,50,60,702,5,20,35,505,15,20,35,40,505
1001,2,5,10,15,20,30,40,50,605,20,355,15,20,35,405,10,15,20,25,30,40
1251,5,10,15
1505,10,20,30,405,205,15,20,30,4030,40
17015,25
2001,5,10,15,20,30,4055,10,20,25,35,45
2505,10,20
3001,5,10,20
3505,10,20
4001,5,10,20,252,25,30,35
45030-
5001,5,10,20,30510155,30,35
600
7005,10
75030
10005,105,305,10
150020
200010,155
a

The numbers in the columns of the different pressures indicate the water inlet subcooling values.

Table 3

Test matrices for the subcooled and low-quality film boiling experiments identified in Table 1 (1.0 to 9.0 MPa)

p (MPa)
ΔTsub,ina (°C)1.02.03.03.54.06.07.08.09.0
G (kg/m2 s)505
1005,30,4530,45
12010,2010,2020
15030,4530,40,5030,45
2005,10,20,30,405,4030,35,40,45,5035,50
21510,2010,2020,35
2504030
30025,30,35,40
35010
35910,2010,20,3520,35,5020,35,50
4002,25,30,3525,30,35,4010305,1030
5005,30,35
52410,20,3510,2020,35,5020,35,5050
60025
70020
92310,202020,3520,35,5020,50
100010,1520,30
2000
275010,20
p (MPa)
ΔTsub,ina (°C)1.02.03.03.54.06.07.08.09.0
G (kg/m2 s)505
1005,30,4530,45
12010,2010,2020
15030,4530,40,5030,45
2005,10,20,30,405,4030,35,40,45,5035,50
21510,2010,2020,35
2504030
30025,30,35,40
35010
35910,2010,20,3520,35,5020,35,50
4002,25,30,3525,30,35,4010305,1030
5005,30,35
52410,20,3510,2020,35,5020,35,5050
60025
70020
92310,202020,3520,35,5020,50
100010,1520,30
2000
275010,20
a

The numbers in the columns of the different pressures indicate the water inlet subcooling values.

2.3 Data Reduction.

In the experiments summarized in Table 1, the outer surface temperature of the test section was measured by thermocouples. The inner surface temperature is directly related to the wall heat transfer characteristics in convective heat transfer problems. Based on the measured outer surface temperatures and wall heat fluxes under steady-state conditions, neglecting the axial heat conduction (except the quench front region) and heat losses to the ambient, the inner wall temperature can be derived by
Tw=Tw,outqindi(do2di2)kw[do22ln(dodi)do2di24]
(1)
where Tw and Tw,out are respectively the inner and outer surface temperatures of the test section (K); qin the heat flux from the test section to the fluid based on the inner surface area of the test section (W/m2); kw the thermal conductivity of the test section wall (W/m·K); do the outer diameter of the test section (m); and di the inner diameter of the test section (m). The wall convective heat transfer coefficient is defined based on the local liquid saturation temperature as
hc=qinqradTwTsat
(2)
where Tsat is the fluid saturation temperature (K). The radiative heat flux qrad is calculated as
qrad=σSB(Tw4Tsat4)1εl+1εw1
(3)

where σSB is Stefan-Boltzmann constant, 5.670367 × 10−8 (kg/s3-K4); and εl and εw are respectively the emissivity of the liquid water and heated wall inner surface.

3 Methodology

The prediction of the wall temperature in the post-CHF regimes can be characterized as a regression problem in the supervised learning. In other words, we aim to train a regression ML model with a set of input–output pairs. In this work, the input features to the ML model are the flow conditions of a certain post-CHF experiments, while the output is the predicted wall temperature associated with the conditions. To evaluate the overall generalization capability of the ML model, we divide the full LWFB database into two parts: a training dataset used to train the ML model and a testing dataset to evaluate the ML's performance on unseen observations. While a wide range of supervised learning algorithms can be used in regression problems, the ensemble methods have been proven to be simple but powerful approaches that combine the predictions from multiple base estimators to improve the predictive performance. RF and GBDT are two robust tree-based ensemble learning methods with DTs as the base estimators. In this work, these two methods are used to predict wall temperatures in the post-CHF flow regimes. Both methods are computationally efficient with good generalization capabilities and do not require fine tuning of many hyperparameters.

3.1 Random Forest and Gradient Boosted Decision Trees Methods.

The RF method was proposed by Breiman [41] and represents a prominent ensemble method that builds many individual trees in parallel. To grow multiple trees, RF uses the bagging or bootstrapping approach to generate samples with replacement from the training data. The sampled data are to grow DTs in parallel and all the DTs are independent and randomized. Since the data in each DT contain different portions of the original data, the corresponding trees differ across samples and therefore, form an ensemble of distinct DTs. After growing a RF, the final prediction is obtained by averaging the results over all the individual DTs for regression problems. Averaging usually achieves a better performance in prediction than a single DT because it decreases the variance and improves the accuracy of the model. Figure 2 shows a simplified RF structure for regression problems. The green nodes represent the optimal path. In this work, we train a RF model with 190 trees and the maximum depth of each tree is 16. The total number of the edges from the root node to the leaf node in the longest path is known as the depth of the tree. The tree depth controls the complexity of each tree and the computational cost.

Fig. 2
A simplified RF structure for regression problems
Fig. 2
A simplified RF structure for regression problems
Close modal

RF is an averaging method, in which each DT is grown separately and independent of each other. GBDT represents another family of ensemble methods based on sequential learning proposed by Friedman [42]. Figure 3 shows a simplified GBDT structure for regression problems. In the GBDT method, each new tree is grown in sequence based on the pseudo-residuals of the previous tree and thus each tree depends on the results of its predecessor and provides a more accurate prediction. By adding sequential trees to the ensemble model (combination of all the previous trees), the loss function is minimized, and the overall error of prediction reaches a global minimum. In this work, our GBDT model has 400 trees, and the maximum depth of each tree is 6.

Fig. 3
A simplified GBDT structure for regression problems
Fig. 3
A simplified GBDT structure for regression problems
Close modal

3.2 Model Training and Validation Procedure.

Both the RF and GBDT models are trained and validated using the scikit-learn [43], a Python-based open-source ML library. Figure 4 illustrates the overall workflow of training the two ML models for post-CHF wall temperature predictions. In this work, the input features are selected based on the understanding of the physical problem to construct the ML models. Based on the review of the water-based experimental data in Sec. 2 and understanding of the problem, five parameters, i.e., the system pressure, mass flux, inlet subcooling, hydraulic diameter (Dh), and local thermodynamic equilibrium quality (xe) are selected as the input variables and stored in a feature vector X. The thermodynamic equilibrium quality is obtained from each local measurement location, which represents the local thermodynamic state. Important two-phase flow terms, such as the fluid density, fluid enthalpy, local mixture Reynolds number, and so on, can be calculated from these five parameters. The inner wall surface temperature of the test section, i.e., Tw, is selected as the target output. The LWFB database described in Sec. 2.2 is formatted into a two-dimensional (2D) numerical matrix and randomly divided into two parts, i.e., 80% training data and 20% test data, as the input to the ML algorithms to train and validate the ML models, respectively. During the training process, the ML predicted target, i.e., Tw in this study, is compared with the results calculated from the training data through a loss function, ε, usually presented as the mean squared error (MSE) as
MSE=1Ni=1n(ymyp)2
(4)

where N is the number of data points, ym the observed value, and yp the predicted value. The objective of the training process is to minimize the loss function by gradually updating the trainable parameters of the ML model. The trained ML model is then validated by the testing dataset. If the trained ML model does not have good performance on either the training or testing dataset, the hyperparameters are tuned to repeat the training process. Through multiple iterations, the ML model with optimal performance can be obtained. The experimental data we collected and the two ML models we developed can be accessed by readers through the link provided in the Appendix.

Fig. 4
Overall workflow of training ML models
Fig. 4
Overall workflow of training ML models
Close modal

4 Results and Discussion

4.1 Comparison of Conventional Empirical Correlations and ML Models.

In the past several decades, many empirical heat transfer correlations have been developed to calculate wall temperatures in the post-CHF flow regimes. These empirical correlations usually have limited application ranges of flow conditions. In this section, three widely-used empirical correlations, i.e., Bromley correlation [12], Berenson correlation [44], and the correlation adopted by the nuclear reactor thermal hydraulics analysis code COBRA-TF or CTF [45], are compared with the LWFB database to assess their prediction capabilities for a wide range of flow conditions. The Bromley correlation for film boiling is given as
hBromley=0.62[ρvg(ρlρv)hfgkv3μvdi(TwTsat)]1/4
(5)
where hBromely is the convective heat transfer coefficient calculated using the Bromley correlation (W/m2·K); ρv the vapor density (kg/m3); ρl the liquid water density (kg/m3); g the gravitational acceleration (m/s2); hfg the enthalpy difference between the vapor (evaluated at the arithmetic mean temperature of the wall temperature and saturation temperature) and the liquid at the saturation temperature, kJ/kg; kv the vapor thermal conductivity (W/m-K); and μv the vapor dynamic viscosity (kg/m s). The Berenson correlation for film boiling is given as
hBerenson=0.425[ρvg(ρlρv)hfgkv3μvσ/g(ρlρv)(TwTsat)]1/4
(6)
where hBerenson is the convective heat transfer coefficient calculated using the Berenson correlation (W/m2·K) and σ is the surface tension between the liquid and vapor (kg/s2). The major difference between these two correlations is the choice of the characteristic length scale, i.e.,, using the Laplace length scale σ/g(ρlρv) in the Berenson correlation to replace the tube/pipe inner diameter di in Eq. (5). In CTF, the critical wavelength of the Rayleigh-Taylor instability, λt, was introduced into the Bromley correlation as
hCTF=0.62(diλt)0.172[ρvg(ρlρv)hfgkv3μvdi(TwTsat)]1/4
(7)
where hCTF is the convective heat transfer coefficient calculated using the correlation in CTF for film boiling (W/m2·K) and the critical wavelength of the Rayleigh-Taylor instability is given as
λt=2πσg(ρlρv)
(8)

Based on the measured heat flux in the experiments, the wall temperature can be calculated from Eqs. (2) and (3) using these three empirical correlations.

Figure 5 shows comparisons of the wall temperatures calculated by these three correlations against those in the LWFB database, as well as the relative errors. The predictions of these correlations are similar and large discrepancies are observed between the experimental results and the three empirical correlations. The percentage of the data points that are within the ±50% error band is 68.9%, 69.8%, and 71.3% for the Bromley correlation, Berenson correlation, and CTF correlation, respectively. The relatively poor performance of these three correlations suggests that the conventional empirical correlations do not have good generalization capability to cover the wide range of flow conditions in the LWFB database. Moreover, these conventional heat transfer coefficient correlations were developed based on the fluid saturation temperature instead of the vapor bulk temperature, which may be considerably higher than the saturation temperature in the post-CHF flow regimes. However, the vapor temperature was not measured in the experiments included in Table 1, primarily due to the difficulties in accurately measuring the (high) vapor temperature. In summary, it is difficult for the conventional empirical correlations to accurately predict the wall temperatures in the post-CHF flow regimes over a wide range of flow conditions and therefore, a more robust model with improved generalization capabilities is needed.

Fig. 5
Comparisons of the experimental results with the Bromley correlation [12] (a) and (b), Berenson correlation [44] (c) and (d), and the correlation in CTF [45] (e) and (f)
Fig. 5
Comparisons of the experimental results with the Bromley correlation [12] (a) and (b), Berenson correlation [44] (c) and (d), and the correlation in CTF [45] (e) and (f)
Close modal

In the present work, two data-driven tree-based ML models, based on GBDT and RF, are developed for wall temperature predictions in the post-CHF flow regimes using the LWFB database. The full LWFB database is split 80% and 20%, respectively, into a training dataset and a test dataset for model training and validation, as illustrated in Fig. 4. Figure 6 shows the predictions of the GBDT and RF models against the training data and test data, respectively. In general, the wall temperatures calculated by the RF and GBDT models match well with the experimental data. 98.7% of the training data points and 96.8% of the test data points are within the ±20% error band for the GBDT model and 98.4% and 96.1%, respectively, for the RF model. We observe that the ML models perform slightly better for the training data than the test data. This could be caused by the overfitting of the models, but such an overfitting is controlled well through the training process and can be deem as acceptable. Figures 6(b) and 6(d) show that the GBDT model performs slightly better than the RF model. The relative error slightly decreases as the wall temperature increases. Comparing Fig. 5 with Fig. 6, the predictive capabilities of these two ML models demonstrate significant improvement over the three conventional empirical correlations.

Fig. 6
Comparisons of the wall temperature results by the two ML models with the corresponding experimental data: (a) wall temperature by the GBDT model, (b) relative error of the GBDT model, (c) wall temperature by the RF model, and (d) relative error of the RF model
Fig. 6
Comparisons of the wall temperature results by the two ML models with the corresponding experimental data: (a) wall temperature by the GBDT model, (b) relative error of the GBDT model, (c) wall temperature by the RF model, and (d) relative error of the RF model
Close modal

The ML models can be implemented into computational thermal-hydraulic codes of interest to replace the empirical model(s) for improved modeling capability. As a result, the wall temperature, vapor temperature, and heat transfer coefficient can be evaluated through an iterative process. A detailed procedure for such an implementation will not be discussed in this work, but the implementation can be similar to the subcooled flow boiling heat flux partitioning algorithm [46]. It should be noted that although the current ML models are trained based on a Python script and the scikit-learn library, no issue is expected to convert the models to another computer language, such as Fortran or C++, so the ML models can be compiled with the thermal-hydraulic codes as a data-driven closure relation. In this way, the ML models will be naturally fit into the physics-based model and its solution process. One such an example is the SAM-ML work [47], in which a Python-based neural network model is converted into C++ and is then integrated into the System Analysis Module (SAM) code [48] for sodium-cooled fast reactor safety analysis.

4.2 Statistical Analysis of the Machine Learning Models.

To further evaluate the performance of the GBDT and RF models, a statistical analysis is performed in this section. In general, any ML model will inevitably involve uncertainties that can be decomposed into two parts: the epistemic uncertainty and aleatoric uncertainty. The epistemic uncertainty stems from the lack of knowledge or reliable data, which we aim to minimize through a robust training process and a comprehensive training dataset; while the aleatoric uncertainty stems from the stochastic nature of the physical process and is irreducible. To investigate the accuracy of the trained ML models, the aleatoric part of the prediction, represented by the relative error of the wall superheat, is calculated as
δi=Tw,pTw,mTw,mTsat
(9)

where Tsat is the saturation temperature (°C). The relative error of a good regression model should meet three criteria as follows:

  1. The relative errors follow a normal distribution;

  2. The relative errors are independent of the input features X (p, G, ΔTsub,in, xe, and Dh in this study). In other words, the relative errors should have no pattern relative to the input features X;

  3. Most of the relative errors are close to zero.

To study the characteristics of the relative error of the ML models, the above three criteria are investigated by plotting the relative error in various ways. To study the first criterion, the histograms of the relative errors of the wall temperatures from the training and test datasets for the GBDT and RF models are shown in Fig. 7. The probability density function in the histogram is calculated by each bin's raw count divided by the total number of counts and the bin width. The normal distributions in Fig. 7 are obtained based on the mean and variance of the relative errors in the training and test datasets, respectively. The relative errors from these two ML models, in general, follow a normal distribution for both the training and test data, but the relative errors for the GBDT model are closer to the normal distribution than those for the RF model, indicating the GBDT model performs marginally better than the RF model.

Fig. 7
Distributions of the relative errors of the predicted wall temperatures (a) GBDT model (b) RF model
Fig. 7
Distributions of the relative errors of the predicted wall temperatures (a) GBDT model (b) RF model
Close modal

We further use the quantile–quantile (QQ) plot, a commonly used visual method, to evaluate the relative error distribution of the predictions by the two ML models by plotting their quantiles against the ideal normal distributions, as depicted in Fig. 8. The values on the x axis are based on the normal distribution and those on the y axis are relative errors of the wall superheats. If the two distributions are similar, the points in the QQ plot will approximately lie on the line y = x. It can be seen from Fig. 8 that each of these two plots for the GBDT and RF models have an “S” shape. In the relative error range of [−0.1 0.2], the plots approximately lie on the line y = x for both the GBDT and RF models, which indicates that the distribution of the relative error of the predicted wall superheats follows a normal distribution in this range. Near the two ends of the plots, where the relative error is larger than 0.2 or smaller than −0.1, the trend of the QQ plot is steeper than the line y = x, indicating that distributions of the relative errors of the wall superheats have heavier tails than a normal distribution, which can also be seen in Fig. 7. Over 91% of the relative error data points are within [−0.1, 0.2] for the training data and over 86% for the test data for both the GBDT and RF models. Therefore, the relative errors can be regarded as, in general, following normal distributions based on the QQ plots for both the GBDT and RF models.

Fig. 8
The Q–Q plots for the relative errors of the wall temperatures: (a) the GBDT model (b) the RF model
Fig. 8
The Q–Q plots for the relative errors of the wall temperatures: (a) the GBDT model (b) the RF model
Close modal

To investigate the second criterion, Fig. 9 shows the relative error of the wall superheats relative to five variables, i.e., p, G, ΔTsub,in, xe, and Dh, for the GBDT and RF models. The relative errors are randomly distributed along the x axis and do not follow any patterns relative to the x axis, which indicates the relative errors meet the second criterion. For the thermodynamic equilibrium quality in Fig. 9(d), the relative errors have slightly higher peaks near thermodynamic equilibrium quality value of zero. From Figs. 9(a) and 9(b), the absolute relative errors have large values, slightly over 60%, for low-pressure (<0.7 MPa) and low-mass flux (<550 kg/m2 s) conditions while for high pressures and high mass fluxes, most of the absolute relative errors are within 20%. This is because most of the data in the LWFB database were obtained from low-pressure and low-mass flux conditions, i.e., 57.6% of the data points fall in the pressure range of [0.1, 0.7] MPa and 86.3% of the data in the mass flux range of [25, 550] kg/m2 s. As described in Sec. 2, the LWFB database is consolidated based on nine different data sources, and this introduce uncertainties resulting from the use of different test facilities and the data collection by different researchers. Therefore, it appears reasonable that the more data we have, the more scattering the data will be, which will result in relatively larger prediction errors. From Fig. 9(d), the absolute relative errors have large values in the thermodynamic equilibrium quality range of [0, 0.1]. One of the reasons for this large relative errors is that similar to the above hypothesis, 54.6% of the data points fall in this range. In addition, the flow regime transition from the IAFB to the ISFB regimes may occur in this thermodynamic equilibrium quality range [49], which may also contribute to the observed large prediction uncertainties.

Fig. 9
Relative errors of the predicted wall superheat versus the feature variables for all the collected data: (a) pressure, (b) mass flux, (c) inlet subcooling, (d) thermodynamic equilibrium quality, and (e) Hydraulic diameter
Fig. 9
Relative errors of the predicted wall superheat versus the feature variables for all the collected data: (a) pressure, (b) mass flux, (c) inlet subcooling, (d) thermodynamic equilibrium quality, and (e) Hydraulic diameter
Close modal

To assess the two models with respect to the third criterion, the cumulative probability distribution function (CPDF) of the absolute relative error for the two ML models are calculated and shown in Fig. 10. The CPDF curves for the training data rise slightly steeper than those for the test data. The curves rise rapidly when the absolute relative error is smaller than 10% and then tend to become flat. For both the GBDT and RF models, over 89% of the relative errors are within [−0.1, 0.1] for the training data and over 80% for the test data. Therefore, it is believed that the relative errors of the two models meet the third criterion.

Fig. 10
Cumulative probability distribution of the absolute relative error
Fig. 10
Cumulative probability distribution of the absolute relative error
Close modal
It should be noted that the above discussions are based on the data split of 80% training data and 20% test data. To study the sensitivity of the data split ratio (DSR) and prevent overfitting, the LWFB database is split based on another two different DSRs, i.e., 90% training data and 10% test data, and 70% training data and 30% test data. In addition to the above three criteria, three metrics are calculated to assess the accuracy of the GBDT and RF models, i.e., the root-mean-square error (RMSE), the correlation coefficient (R), and the 95% confidence interval (CI) of the relative error based on three different DSRs. The RMSE is calculated as
RMSE=1Ni=1Nδi2
(10)
The correlation coefficient, i.e., R, is defined by
R=±1i=1N(ym,iyp,i)2i=1N(ym,iy¯m)2
(11)

The range of R is [−1:1].

The results for these three metrics are summarized in Table 4. Minor differences of these three metrics are observed for the three different DSRs. The two trained ML models showed good predictive capabilities for all different DSRs. Taking the (80%, 20%) split as an example, the RMSEs of both models are very small, i.e., 5.7% for the GBDT model and 6.2% for the RF model. The correlation coefficient is 97.1% for the GBDT model and 96.6% for the RF model. The 95% CI of the relative error is [−17.5%, 15.5%] for the GBDT model and [−18.8%, 16.7%] for the RF model. All these three metrics for the test data indicate that the GBDT and RF models have good predictive capabilities for the whole LWFB dataset while the GBDT model performs slightly better than the RF model, which is consistent with the discussions in Sec. 4.1. The applicable range of the developed GBDT and RD models is recommended to ensure the validity of the two ML models as: 0.1 p 9.0 MPa, 25 G 2750 kg/m2 s, and 1 ΔTsub, in 70 °C.

Table 4

Metrics to evaluate the performance of two ML models

ModelGBDT modelRF model
DSRa (%)90, 1080, 2070, 3090, 1080, 2070, 30
RMSEa (%)4.5, 5.54.4, 5.74.3, 5.84.4, 6.04.3, 6.24.3, 6.3
Ra (%)98.2, 97.398.2, 97.198.3, 97.098.4, 96.998.4, 96.698.4, 96.5
95% CIb (%)[−14.1, 12.5][−13.8, 12.3][−13.4, 12.0][−13.7, 12.3][−13.4, 12.0][−13.6, 12.1]
[−16.8, 14.9][−17.5, 15.5][−17.8, 15.9][−18.2, 16.1][−18.8, 16.7][−19.2, 17.2]
ModelGBDT modelRF model
DSRa (%)90, 1080, 2070, 3090, 1080, 2070, 30
RMSEa (%)4.5, 5.54.4, 5.74.3, 5.84.4, 6.04.3, 6.24.3, 6.3
Ra (%)98.2, 97.398.2, 97.198.3, 97.098.4, 96.998.4, 96.698.4, 96.5
95% CIb (%)[−14.1, 12.5][−13.8, 12.3][−13.4, 12.0][−13.7, 12.3][−13.4, 12.0][−13.6, 12.1]
[−16.8, 14.9][−17.5, 15.5][−17.8, 15.9][−18.2, 16.1][−18.8, 16.7][−19.2, 17.2]
a

The first number is for the training data and the second number is for the test data.

b

The first row is for the training data and the second row is for the test data.

5 Conclusions

In this study, we propose a data-driven approach to develop machine learning models for wall temperature prediction in the post-CHF flow regimes. A general low-quality water film boiling database with 22,813 data points is consolidated from nine different data sources in the literature. Two tree-based ensemble models, i.e., GBDT and RF, are developed for the wall temperature prediction in the post-CHF flow regimes based on the LWFB database. The ML models take five flow features as the inputs, namely, the system pressure, mass flux, inlet subcooling, hydraulic diameter of the flow channel, and thermodynamic equilibrium quality, and predict the wall temperature as the output. The LWFB database is split into training data and test data to train and validate the GBDT and RF models, respectively. The performance of the two trained ML models is thoroughly evaluated and compared with three widely used empirical correlations with the full LWFB dataset. The conventional correlations show relatively poor performance in predicting the LWFB database with large errors, while the ML models show good accuracies for the entire database. A statistical analysis is further performed to evaluate the prediction uncertainty of the ML models. We found out that the relative errors between the measured and predicted wall superheats by the ML models approximately follow normal distributions. The relative errors are randomly distributed and have no clear patterns relative to the input variables. The cumulative probability distribution function within 10% error band covers over 80% of the test data points for both the GBDT and RF models. The RMSEs are 5.7% and 6.2% for the GBDT and RF models for the test data, respectively. The correlation coefficient is 97.1% for the GBDT model and 96.6% for the RF model. The 95% CI of the relative error is [-17.5%, 15.5%] for the GBDT model and [−18.8%, 16.7%] for the RF model. All these three metrics for the test data indicate that the GBDT and RF models have good predictive capabilities over a broad range of flow conditions for wall temperature in the post-CHF flow regimes, while the GBDT model has slightly better performance than the RF model. These two trained ML models have good predictive capabilities covering a broad range of conditions for vertical upward flow, i.e., the system pressure from 0.1 to 9.0 MPa, mass flux from 25 to 2750 kg/m2 s, and inlet subcooling from 1 to 70 °C, and can obtain significantly more accurate wall temperature predictions than the currently used empirical correlations in the post-CHF flow regimes at subcooled and low-quality conditions.

Acknowledgment

The authors at the University of Michigan appreciate the financial support as well as the technical and administrative support provided by Mr. Andrew Ireland, Dr. Chris Hoxie, and Jennifer Dudek of the U.S. NRC.

Funding Data

  • U.S. Nuclear Regulatory Commission (Award ID: 31310020C0009; Funder ID: 10.13039/100005187).

Data Availability Statement

The data and information that support the findings of this article are freely available online.2

Appendix: The Collected Post-CHF Data and the Tree-Based ML Models

Part of the LWFB database used for training and the two tree-based ensemble models can be accessed through the following link given in footnote2.

Footnotes

References

1.
Sindhuja
,
R.
,
Balakrishnan
,
A. R.
, and
Murthy
,
S. S.
,
2010
, “
Post-CHF Heat Transfer During Two-Phase Upflow Boiling of R-407C in a Vertical Pipe
,”
Appl. Therm. Eng.
,
30
(
2–3
), pp.
167
173
.10.1016/j.applthermaleng.2009.07.016
2.
Jha
,
J. M.
,
Ravikumar
,
S. V.
,
Tiara
,
A. M.
,
Sarkar
,
I.
,
Pal
,
S. K.
, and
Chakraborty
,
S.
,
2015
, “
Ultrafast Cooling of a Hot Moving Steel Plate by Using Alumina Nanofluid Based Air Atomized Spray Impingement
,”
Appl. Therm. Eng.
,
75
, pp.
738
747
.10.1016/j.applthermaleng.2014.10.005
3.
Kefer
,
V.
,
Käler
,
W.
, and
Kastner
,
W.
,
1989
, “
Critical Heat Flux (CHF) and Post-CHF Heat Transfer in Horizontal and Inclined Evaporator Tubes
,”
Int. J. Multiphase Flow
,
15
(
3
), pp.
385
392
.10.1016/0301-9322(89)90008-6
4.
Nelson
,
R. A.
,
1982
, “
Forced-Convective Post-CHF Heat Transfer and Quenching
,”
ASME J. Heat Mass Transfer-Trans. ASME
,
104
(
1
), pp.
48
54
.10.1115/1.3245067
5.
Hammouda
,
N.
,
Groeneveld
,
D. C.
, and
Cheng
,
S. C.
,
1996
, “
An Experimental Study of Subcooled Film Boiling of Refrigerants in Vertical Up-Flow
,”
Int. J. Heat Mass Transfer
,
39
(
18
), pp.
3799
3812
.10.1016/0017-9310(96)00062-2
6.
Yagov
,
V. V.
,
Zabirov
,
A. R.
,
Kanin
,
P. K.
, and
Denisov
,
M. A.
,
2017
, “
Heat Transfer in Film Boiling of Subcooled Liquids: New Experimental Results and Computational Equations
,”
J. Eng. Phys. Thermophys.
,
90
(
2
), pp.
266
275
.10.1007/s10891-017-1564-3
7.
Ebrahim
,
S. A.
,
Chang
,
S.
,
Cheung
,
F. B.
, and
Bajorek
,
S. M.
,
2018
, “
Parametric Investigation of Film Boiling Heat Transfer on the Quenching of Vertical Rods in Water Pool
,”
Appl. Therm. Eng.
,
140
, pp.
139
146
.10.1016/j.applthermaleng.2018.05.021
8.
Oliveira
,
A. V. S.
,
Carrillo
,
J. P.
,
Labergue
,
A.
,
Glantz
,
T.
, and
Gradeck
,
M.
,
2020
, “
Experimental Study of Dispersed Flow Film Boiling at Sub-Channel Scale in LOCA Conditions: Influence of the Steam Flow Rate and Residual Power
,”
Appl. Therm. Eng.
,
172
, p.
115143
.10.1016/j.applthermaleng.2020.115143
9.
Ni
,
P.
,
Wen
,
Z.
,
Su
,
F.
,
Huang
,
J.
,
Liu
,
X.
,
Lou
,
G.
, and
Dou
,
R.
,
2020
, “
Film Boiling Collapse in a Solid Hot Sphere Immersed in Subcooled Forced Convection
,”
Appl. Therm. Eng.
,
166
, p.
114630
.10.1016/j.applthermaleng.2019.114630
10.
Liu
,
Q.
, and
Sun
,
X.
,
2020
, “
Wall Heat Transfer in the Inverted Annular Film Boiling Regime
,”
Nucl. Eng. Des.
,
363
, p.
110660
.10.1016/j.nucengdes.2020.110660
11.
Liu
,
Q.
,
Sun
,
H.
,
Liu
,
Y.
,
Kelly
,
J.
, and
Sun
,
X.
,
2021
, “
Experimental Study of Post-CHF Heat Transfer in a Vertical Tubular Test Section
,”
Int. J. Heat Mass Transfer
,
166
, p.
120697
.10.1016/j.ijheatmasstransfer.2020.120697
12.
Bromley
,
L. A.
,
1950
, “
Heat Transfer in Stable Film Boiling
,”
Chem. Eng. Prog.
,
46
(
5
), pp.
221
227
.https://scholar.google.com/scholar_lookup?title=Heat%20transfer%20in%20stable%20film%20boiling&publication_year=1950&author=L.A.%20Bromley
13.
Dittus
,
F. W.
, and
Boelter
,
L. M. K.
,
1985
, “
Heat Transfer in Automobile Radiators of the Tubular Type
,”
Int. Commun. Heat Mass Transfer
,
12
(
1
), pp.
3
22
.10.1016/0735-1933(85)90003-X
14.
U.S. NRC
,
2010
, “
TRACE V5.0 Theory Manual-Field Equations, Solution Methods, and Physical Models
,”
U.S. Nuclear Regulatory Commission
,
Washington, DC
.
15.
Analytis
,
G. T.
, and
Yadigaroglu
,
G.
,
1987
, “
Analytical Modeling of Inverted Annular Film Boiling
,”
Nucl. Eng. Des.
,
99
, pp.
201
212
.10.1016/0029-5493(87)90121-X
16.
Hammouda
,
N.
,
Groeneveld
,
D. C.
, and
Cheng
,
S. C.
,
1997
, “
Two-Fluid Modelling of Inverted Annular Film Boiling
,”
Int. J. Heat Mass Transfer
,
40
(
11
), pp.
2655
2670
.10.1016/S0017-9310(96)00278-5
17.
El Nakla
,
M.
,
Groeneveld
,
D. C.
, and
Cheng
,
S. C.
,
2010
, “
Modeling Subcooled Flow Film Boiling in a Vertical Tube
,”
ASME J. Therm. Sci. Eng. Appl.
,
2
(
2
), p.
021002
.10.1115/1.4002526
18.
Liu
,
Q.
,
Shi
,
S.
,
Sun
,
X.
, and
Kelly
,
J.
,
2018
, “
Thermal Hydraulic Performance Analysis of a Post-CHF Heat Transfer Test Facility
,”
Nucl. Eng. Des.
,
339
, pp.
53
64
.10.1016/j.nucengdes.2018.08.020
19.
Liu
,
Y.
,
Wang
,
D.
,
Sun
,
X.
,
Liu
,
Y.
,
Dinh
,
N.
, and
Hu
,
R.
,
2021
, “
Uncertainty Quantification for Multiphase-CFD Simulations of Bubbly Flows: A Machine Learning-Based Bayesian Approach Supported by High-Resolution Experiments
,”
Reliab. Eng. Syst. Saf.
,
212
, p.
107636
.10.1016/j.ress.2021.107636
20.
Liu
,
Y.
,
Hu
,
R.
, and
Balaprakash
,
P.
,
2021
, “
Uncertainty Quantification of Deep Neural Network-Based Turbulence Model for Reactor Transient Analysis
,”
ASME
Paper No. VVS2021-65045. 10.1115/VVS2021-65045
21.
Dawahdeh
,
A.
,
Oh
,
J.
,
Zhai
,
T.
, and
Palazzolo
,
A.
,
2021
, “
Computational Fluid Dynamics—Machine Learning Prediction of Machinery Coupling Windage Heating and Power Loss
,”
ASME J. Heat Mass Transfer-Trans. ASME
,
143
(
8
), p.
082201
.10.1115/1.4051351
22.
Levy
,
Y.
,
Fan
,
H. Y.
, and
Sherbaum
,
V.
,
2005
, “
A Numerical Investigation of Mixing Processes in a Novel Combustor Application
,”
ASME J. Heat Mass Transfer-Trans. ASME
,
127
(
12
), pp.
1334
1343
.10.1115/1.2103090
23.
Liu
,
Y.
,
Dinh
,
N.
,
Sato
,
Y.
, and
Niceno
,
B.
,
2018
, “
Data-Driven Modeling for Boiling Heat Transfer: Using Deep Neural Networks and High-Fidelity Simulation Results
,”
Appl. Therm. Eng.
,
144
, pp.
305
320
.10.1016/j.applthermaleng.2018.08.041
24.
Wei
,
H. M.
,
Su
,
G. H.
,
Tian
,
W. X.
,
Qiu
,
S. Z.
,
Ni
,
W. F.
, and
Yang
,
X. B.
,
2010
, “
Study on Dryout Point by Wavelet and GNN
,”
Appl. Therm. Eng.
,
30
(
6–7
), pp.
664
672
.10.1016/j.applthermaleng.2009.11.014
25.
Zhao
,
X.
,
Shirvan
,
K.
,
Salko
,
R. K.
, and
Guo
,
F.
,
2020
, “
On the Prediction of Critical Heat Flux Using a Physics-Informed Machine Learning-Aided Framework
,”
Appl. Therm. Eng.
,
164
, p.
114540
.10.1016/j.applthermaleng.2019.114540
26.
Jiang
,
B. T.
, and
Zhao
,
F. Y.
,
2013
, “
Combination of Support Vector Regression and Artificial Neural Networks for Prediction of Critical Heat Flux
,”
Int. J. Heat Mass Transfer
,
62
, pp.
481
494
.10.1016/j.ijheatmasstransfer.2013.03.025
27.
Yarahmadi
,
M.
,
Robert Mahan
,
J.
, and
McFall
,
K.
,
2020
, “
Artificial Neural Networks in Radiation Heat Transfer Analysis
,”
ASME J. Heat Mass Transfer-Trans. ASME
,
142
(
9
), p.
092801
.10.1115/1.4047052
28.
Groeneveld
,
D. C.
,
1974
, “
Effect of a Heat Flux Spike on the Downstream Dryout Behavior
,”
ASME J. Heat Mass Transfer-Trans. ASME
,
96
(
2
), pp.
121
125
.10.1115/1.3450150
29.
Chen
,
Y.
, and
Li
,
J.
,
1984
, “
Subcooled Flow Film Boiling of Water at Atmospheric Pressure
,”
X. J.
Chen
, and
T. N.
Veziroglu
,
Two-Phase Flow and Heat Transfer
,
Hemisphere Pub. Co
., London, UK, pp.
141
150
.
30.
Fung
,
K. K.
,
1981
, “
Subcooled and Low-Quality Film Boiling of Water in Vertical Flow at Atmospheric Pressure
,” Ph.D. dissertation,
University of Ottawa
,
Ottawa, ON, Canada
.
31.
Stewart
,
J. C.
,
1981
, “
Low Quality Film Boiling at Intermediate and Elevated Pressures
,” M.S. thesis,
University of Ottawa
,
Ottawa, ON, Canada
.
32.
Costigan
,
G.
,
Ralph
,
J. C.
, and
Wade
,
C. D.
,
1987
, “
Steady State Heat Transfer Measurements Relevant to Reactor Reflood Conditions, Part 1: Upflow
,” UKAEA Report No. AERE-R12561.
33.
Mosaad
,
M.
,
1988
, “
Subcooled Film Boiling Heat Transfer to Flowing Water in a Vertical Tube
,” Ph.D. dissertation,
Technische Universitat Berlin
,
Berlin, Germany
.
34.
Chen
,
Y.
,
Fu
,
X.
, and
Chen
,
S.
,
1988
, “
Experimental and Analytical Study of Inverted Annular Film Boiling of Water
,”
Experimental Heat Transfer Fluid Mechanics and Thermodynamics
,
M. D.
Kelleher
,
K. R.
Sreenivasan
,
R. K.
Shah
, and
Y.
Joshi
eds.,
Elsevier Sci. Publ
., Amsterdam, The Netherlands, pp.
1438
1443
.
35.
Chen
,
Y.
,
Cheng
,
P.
,
Wang
,
J.
, and
Yang
,
M.
,
1989
, “
Experimental Results of Subcooled and Low Quality of Film Boiling Heat Transfer of Water in Vertical Tubes at Moderate Pressure
,”
Proceedings of the 4th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics
,
Karlsruhe, Germany
, Oct. 10–13, Vol.
2
, pp.
1106
1110
.
36.
Swinnerton
,
D.
,
Hood
,
M. L.
, and
Pearson
,
K. G.
,
1988
, “
Steady State Post-Dryout Experiments at Low Quality and Medium Pressure Data Report
,” UKAEA Report No. AEEW-R2267.
37.
Swinnerton
,
D.
,
Pearson
,
K. G.
, and
Hood
,
M. L.
,
1988
, “
Steady State Post-Dryout Experiments at Low Quality and Medium Pressure
,” UKAEA Report No. AEEW-R2192.
38.
Savage
,
R. A.
,
Swinnerton
,
D.
, and
Pearson
,
K. G.
,
1989
, “
Heat Transfer and Voidage Measurements in Steady State Post-Dryout at Low Quality and Medium Pressure
,” UKAEA Report No. AEEW-R2486.
39.
Savage
,
R. A.
,
Archer
,
D.
, and
Swinnerton
,
D.
,
1992
, “
Heat Transfer and Voidage Measurements in Steady State Post-Dryout at Low Quality and High Pressure
,”
IChemE 3rd UK National Heat Transfer Conference
,
Birmingham
, Vol.
1
, Sep. 16–18, pp.
607
615
.
40.
Savage
,
R. A.
,
Archer
,
D.
, and
Swinnerton
,
D.
,
1993
, “
Flow Visualization, Heat Transfer, and Voidage Measurements in Steady-State Post-Dryout at High Pressure
,”
Proceedings of the 6th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics
,
Grenoble, France
, Vol.
2
, Oct. 5–8, pp.
1311
1318
.
41.
Breiman
,
L.
,
2001
, “
Random Forests
,”
Mach. Learn
,
45
(
1
), pp.
5
32
.10.1023/A:1010933404324
42.
Friedman
,
J. H.
,
2001
, “
Greedy Function Approximation: A Gradient Boosting Machine
,”
Ann. Stat.
,
29
(
5
), pp.
1189
1232
.
43.
Pedregosa
,
F.
,
Varoquaux
,
G.
,
Gramfort
,
A.
,
Michel
,
V.
,
Thirion
,
B.
,
Grisel
,
O.
,
Blondel
,
M.
,
Prettenhofer
,
P.
, et al.,
2011
, “
Scikit-Learn: Machine Learning in Python
,”
J. Mach. Learn Res.
,
12
, pp.
2825
2830
.https://www.jmlr.org/papers/volume12/pedregosa11a/pedregosa11a.pdf?ref=https:/
44.
Berenson
,
P. J.
,
1961
, “
Film Boiling Heat Transfer From a Horizontal Surface
,”
ASME J. Heat Mass Transfer-Trans. ASME
,
83
(
3
), pp.
351
356
.10.1115/1.3682280
45.
Salko
,
R.
,
Wysocki
,
A.
,
Avramova
,
M.
,
Toptan
,
A.
,
Porter
,
N.
,
Blyth
,
T.
,
Dances
,
C.
,
Gomez
,
A.
,
Jernigan
,
C.
, and
Kelly
,
J.
,
2017
,
CTF Theory Manual
,
The North Carolina State University
,
Raleigh, NC
.
46.
Kurul
,
N.
, and
Podowski
,
M.
,
1990
, “
Multi-Dimensional Effects in Forced Convection Subcooled Boiling
,”
Proceedings of the 9th International Heat Transfer Conference
,
Jerusalem, Israel
, Aug. 19–24, pp.
21
26
.10.1615/IHTC9.40
47.
Liu
,
Y.
,
Hu
,
R.
,
Zou
,
L.
, and
Nunez
,
D.
,
2022
, “
SAM-ML: Integrating Data-Driven Closure With Nuclear System Code SAM for Improved Modeling Capability
,”
Nucl. Eng. Des.
,
400
, p.
112059
.10.1016/j.nucengdes.2022.112059
48.
Hu
,
R.
,
Zou
,
L.
,
Hu
,
G.
,
Nunez
,
D.
,
Mui
,
T.
, and
Fei
,
T.
,
2021
, “
SAM Theory Manual
,”
Argonne National Lab. (ANL)
,
Argonne, IL
, Report No. ANL/NSE-17/4 Rev.1.
49.
Liu
,
Q.
,
Kelly
,
J.
, and
Sun
,
X.
,
2021
, “
Flow Regime Transition in the Post-CHF Flow Regimes Under Subcooled and Low-Quality Conditions
,”
Int. J. Multiphase Flow
,
136
, p.
103543
.10.1016/j.ijmultiphaseflow.2020.103543