Recent Advances in the Critical Heat Flux Amelioration of Pool Boiling Surfaces Using Metal Oxide Nanoparticle Deposition

Pool boiling is an effective heat transfer process in a wide range of applications related to energy conversion, including power generation, solar collectors, cooling systems, refrigeration and air conditioning. By considering the broad range of applications, any improvement in higher heat-removal yield can ameliorate the ultimate heat usage and delay or even avoid the occurrence of system failures, thus leading to remarkable economic, environmental and energy efficiency outcomes. A century of research on ameliorating critical heat flux (CHF) has focused on altering the boiling surface characteristics, such as its nucleation site density, wettability, wickability and heat transfer area, by many innovative techniques. Due to the remarkable interest of using nanoparticle deposition on boiling surfaces, this review is targeted towards investigating whether or not metal oxide nanoparticles can modify surface characteristics to enhance the CHF. The influence of nanoparticle material, thermo-physical properties, concentration, shape, and size are categorized, and the inconsistency or contradictions of the existing research results are recognized. In the following, nanoparticle deposition methods are presented to provide a worthwhile alternative to deposition rather than nanofluid boiling. Furthermore, possible mechanisms and models are identified to explain the amelioration results. Finally, the present status of nanoparticle deposition for CHF amelioration, along with their future challenges, amelioration potentials, limitations, and their possible industrial implementation, is discussed.


Pool Boiling Phenomenon
Pool boiling is proposed as a beneficial method for the transfer of high thermal energy in various applications of industry, such as refrigeration, cooling electronic chips as well as nuclear reactor cooling, the generation of thermal power and the conversion of energy. Regarding to the widespread use of this process in different fields, understanding the origin, process and existing developments in this field is of valuable importance [1,2]. In the pool boiling process, the heating surface is submerged in a pool with a great amount of stagnant liquid. To study the influence of diverse parameters, including the behavior of heater surfaces in addition to the use of various fluids, several investigations have been conducted on a laboratory scale [3][4][5].

CHF Importance in Energy Applications
CHF is a phenomenon relating to the point in which a persistent liquid contact cannot be preserved at the heated surface. Being thoroughly accurate, this special term shows the heat flux corresponding to the occurrence of the phenomenon. Among other terms that are often used are boiling crisis, departure from nucleate boiling (DNB), dryout and burnout. CHF is an interesting and substantial phenomenon from both fundamental and practical standpoints. From a fundamental viewpoint, CHF includes broad variations in heat transfer, pressure drop and flow regime. From the practical viewpoint, lots of application fields are in direct contact with CHF, including fusion applications, nuclear power plants, fossil power plants, steam generators and so on. In these fields, CHF notably affects the controllability, safety, efficiency, integrity and economic competitiveness of systems and subsystems [8].
CHF investigation is also precious in electronic cooling systems, fusion component heat removal, designing neutral particle beams to treat cancer patients, and removing heat from the

CHF Importance in Energy Applications
CHF is a phenomenon relating to the point in which a persistent liquid contact cannot be preserved at the heated surface. Being thoroughly accurate, this special term shows the heat flux corresponding to the occurrence of the phenomenon. Among other terms that are often used are boiling crisis, departure from nucleate boiling (DNB), dryout and burnout. CHF is an interesting and substantial phenomenon from both fundamental and practical standpoints. From a fundamental viewpoint, CHF includes broad variations in heat transfer, pressure drop and flow regime. From the practical viewpoint, lots of application fields are in direct contact with CHF, including fusion applications, nuclear power plants, fossil power plants, steam generators and so on. In these fields, CHF notably affects the controllability, safety, efficiency, integrity and economic competitiveness of systems and subsystems [8].
CHF investigation is also precious in electronic cooling systems, fusion component heat removal, designing neutral particle beams to treat cancer patients, and removing heat from the blades of high-temperature gas turbines [9]. Chio et al. [10] conducted CHF experiments in a three-rod bundle in a high-flux advanced neutron application reactor. So as to consider the geometric effect of the three-rod bundle, a three-pin correction factor was developed. An example of cooling electronic systems based on boiling mechanism utilizing PF-5060 dielectric liquid of high-power computer chips is exhibited, schematically, in Figure 2 in which the composite spreaders increased the total thermal power removed at relatively low chip temperature.
Energies 2020, 13, x FOR PEER REVIEW 3 of 50 blades of high-temperature gas turbines [9]. Chio et al. [10] conducted CHF experiments in a threerod bundle in a high-flux advanced neutron application reactor. So as to consider the geometric effect of the three-rod bundle, a three-pin correction factor was developed. An example of cooling electronic systems based on boiling mechanism utilizing PF-5060 dielectric liquid of high-power computer chips is exhibited, schematically, in Figure 2 in which the composite spreaders increased the total thermal power removed at relatively low chip temperature. Since vertical two-dimensional channels made of parallel plates are a frequently utilized element in thermal engineering facilities, Guo and Wu [11] systematically investigated thermal drag and CHF for natural convection of air in vertical channels. As an outcome, CHF may appear in the channel and the boiling process, because the possibility of heat transfer for the channel geometry is restricted. The following suggested correlation that can evaluate CHF presents a new criterion for thermal analysis and designing narrow-spaced equipment. , = 0.06 (1) where , , , , , , , and represent the critical heat flux, density of ambient fluid, specific heat capacity, gravitational acceleration, channel half-spacing, temperature of ambient fluid, kinematic viscosity of ambient fluid and channel height, respectively. Wu and Zhao [12] explored nanofluid heat transfer and CHF amelioration from the standpoint of a research gap in engineering applications. They expressed that the available empirical data of nanofluids' thermal characteristics are incomplete and invalid when engineering applications are considered to be important. To succeed in dealing with the research gaps for these applications, the suggested approaches include estimating nanofluid stability under quiescent and flow conditions, establishing a nanofluid database containing nanoparticle size, distribution and additives, which is helpful to advance the knowledge about the interaction of suspended nanoparticles and boundary layers. Hence, realizing convective heat transfer enhancement, and ultimately the empirical and numerical study of the bubble dynamics of boiling nanofluids, will be accessible.
Two-phase pressure drop, boiling heat transfer, and CHF to water were checked by Yu et al. [13] in small horizontal tubes installed in compact heat exchangers. It was pointed out that CHF occurred at relatively great qualities of water (between 0.5 and 1.0) in the horizontal channel being studied. These qualities are higher than those reported for large tubes (at higher pressures and mass fluxes). The CHF quality was concluded to decrease with mass flux reduction, while opposite behavior had been observed in large tubes.
Since enhancing CHF of industrial boilers through surface modification can result in vital energy savings and global reduction in greenhouse gas emissions, Dhillon et al. [14] scrutinized the influence of CHF maxima during boiling crisis on textured surfaces. They proposed a coupled Since vertical two-dimensional channels made of parallel plates are a frequently utilized element in thermal engineering facilities, Guo and Wu [11] systematically investigated thermal drag and CHF for natural convection of air in vertical channels. As an outcome, CHF may appear in the channel and the boiling process, because the possibility of heat transfer for the channel geometry is restricted. The following suggested correlation that can evaluate CHF presents a new criterion for thermal analysis and designing narrow-spaced equipment.
where q w,cri , ρ ∞ , C P , g, b, T ∞ , ν ∞ and l represent the critical heat flux, density of ambient fluid, specific heat capacity, gravitational acceleration, channel half-spacing, temperature of ambient fluid, kinematic viscosity of ambient fluid and channel height, respectively. Wu and Zhao [12] explored nanofluid heat transfer and CHF amelioration from the standpoint of a research gap in engineering applications. They expressed that the available empirical data of nanofluids' thermal characteristics are incomplete and invalid when engineering applications are considered to be important. To succeed in dealing with the research gaps for these applications, the suggested approaches include estimating nanofluid stability under quiescent and flow conditions, establishing a nanofluid database containing nanoparticle size, distribution and additives, which is helpful to advance the knowledge about the interaction of suspended nanoparticles and boundary layers. Hence, realizing convective heat transfer enhancement, and ultimately the empirical and numerical study of the bubble dynamics of boiling nanofluids, will be accessible. Two-phase pressure drop, boiling heat transfer, and CHF to water were checked by Yu et al. [13] in small horizontal tubes installed in compact heat exchangers. It was pointed out that CHF occurred at relatively great qualities of water (between 0.5 and 1.0) in the horizontal channel being studied. These qualities are higher than those reported for large tubes (at higher pressures and mass fluxes). The CHF quality was concluded to decrease with mass flux reduction, while opposite behavior had been observed in large tubes.

CHF Enhancement by Ameliorating Fluid Properties with Nanoparticles
Altering the working fluids can be considered as a primary way for CHF enhancement. Several research studies were conducted on the replacement of different fluids, including Freon-12 or Freon-134a with water in the early 1990s [15]. Improving the properties of the base fluid with nanoparticles is regarded as another useful solution of heat transfer enhancement. Nanofluids are a new category of heat transfer fluids, in which the stable suspending nanoparticles (with the approximate dimensions of 1-50 nm) are dispersed in common base fluids. Various fluids can be used as the base fluids, including water [16], refrigerant [17], or organic fluids [18]. among them, water, oil, ethylene, bio-fluids, and different polymer solutions can be mentioned. Moreover, nanoparticles can include diverse materials such as metals with chemical stability (copper, silver, gold), metal oxides (alumina, titania, silica, zirconia, bismuth oxide), carbon allotropes (diamond, fullerene, single-walled and multi-walled carbon nanotubes), as well as functionalized nanoparticles [19][20][21][22][23]. Even hybrid nanoparticles and nano-fluids can be used for heat transfer [24] and pool boiling enhancement [25,26]. Thus, it can be inferred that the CHF is enhanced by the modification of fluid properties through adding various particles to the base fluids [27,28].

CHF Enhancement by Modifying Surface Characteristics
Generally, surface characteristics can be improved by the deposition of the nanoparticles [29][30][31]. The heat transfer surface in pool boiling can be modified by uniform/modulated porous coatings and fabrication/installation of the diverse structures on the heated surface [32]. Different methods have been presented for the fabrication of surfaces with porous coating of uniform thickness, including welding, sintering, brazing of particles, flame spraying, electrolytic deposition, plating to bond of particle, galvanizing, plasma spraying of a polymer, metallic coating of a foam substrate, and the deposition of nanoparticles. In a different investigation, various processes have been implemented for producing a heater surface with enhanced CHF, including nano-wires [33,34], nano-porous structures [9,35], micro-porous surface layers [36,37], nano-tube arrays [28], carbon nano-tube coating [38,39] or composite porous surfaces [40].
Experimental studies on a modulated porous-layer coating revealed CHF enhancement about 2-3 times in comparison with a plain surface [41]. In some studies [42,43], innovative structures with separated paths of liquid and vapor were proposed as porous media. These separated paths may lead to a decline in flow resistance, which results in an enhancement in CHF. Further enhancement in CHF by micro-fabrication techniques can be achieved by the use of microstructures that are coated with zinc oxide nano-rods [44], silicon nanowires [45,46], sintered surface [47,48] and the tobacco mosaic virus [49].

CHF Enhancement by Changing Flow Channel Structures
Improving the structure of the flow passage is considered as an effective approach for CHF amelioration. The deposition of the nanoparticle affects the flow pattern of nanofluid in the vicinity of the structured surface. Through the machining of the micro-channels' surfaces, some bubbles may Energies 2020, 13, 4026 5 of 49 be formed, grow and detach to inhibit the instability during the transformation into film boiling. In a study by Kandlikar [50], pool boiling was performed in a substrate with fins. According to the obtained results, it was observed that CHF value was 3 MW/cm 2 at a superheating temperature equal to 4.8 • C. However, this value was 1.25 MW/cm 2 at a higher superheating temperature of 17 • C in a plain surface. In another study, Boziuk et al. [51] achieved 218% CHF amelioration in comparison with the plain surface through the use of micro-channels. According to [52], expanding the cross-section is another way to modify the micro-channels. By the expansion angle in the range of 0-2 • , a significant increase in CHF from 178 to 489 W/cm 2 can be achieved besides a low-pressure drop. In addition, Gheytaghi et al. [53] demonstrated the 65% increase in CHF with inclined micro-channels.

CHF Enhancement by Integration of Hybrid Approaches
Integration of the aforementioned approaches with the porous structures is regarded as an another approach for further CHF enhancement [54]. Bai et al. [55] explained the influence of a porous artery structure on the CHF enhancement. In their study, the environmental conditions for the buildup of capillary pressure was provided through the use of porous layers in the interface of liquid/vapor. This led to a movement of vapor venting through the artery to the side liquid rather than the upper bulk liquid. Thus, the possibility of phase transformation would be reduced during nucleation boiling. Three designs of porous coated surfaces were suggested by Joshi and Dede [56]. Moreover, jet orifices were used to evaluate the performance of two-phase flow. These designs included an open tunnel, a closed tunnel, and a pin-fin heat spreader on it. Among the three mentioned designs, the third design reached the best CHF of more than 218 W/cm 2 , while the reported CHF for other designs was about 200 W/cm 2 . Another example of a hybrid approach is the use of acoustic actuation, where interfacial forces would be induced, which motivates the detachment of the bubbles in a smaller scale and prevents the vapor buildup. It was reported that the integration of textured surface with acoustic actuation induced by a 1.7 MHz ultrasonic piezoelectric transducer [57] would lead to 65% and 30% CHF enhancement in plain and modified surfaces, respectively.

Previous Review Studies on CHF Enhancement
Lots of review papers have been issued in the field of enhancing heat transfer in pool boiling by surface modification. These papers can be categorized into two groups. In the first group, pool boiling was explored along with other heat transfer mechanisms. However, in the second group, pool boiling was studied exclusively. The investigations of Bergles et al. [58] and Webb et al. [59], in which 508 papers and reports as well as 59 U.S. patents were documented in the field of pool boiling and flow boiling prior to 1980 are classified in the first group. Furthermore, the studies of Bhavnani et al. [60], Shojaeian and Koşar [61], McCarthy et al. [62], Kim et al. [63], and Attinger et al. [64] are included in this group. Among the reviews focusing entirely on pool boiling enhancement, the paper of Honda and Wei [65] can be mentioned, which reported the studies for the improvement of cooling performance in electronic devices through the usage of surface microstructures prior to 2004. In another research, Lu and Kandlikar [66] reviewed the enhancement of pool boiling by using modification techniques for nanoscale surfaces. In addition, manufacturing methods for porous surfaces in pool boiling applications were reviewed by Patil and Kandlikar [67]. Moreover, Mori and Utaka [68] have recently published a review for the used techniques to CHF improvement in pool boiling. However, amelioration of HTC was not mentioned within the nucleate boiling region.
In a research study by Rioux et al. [69], nanostructure surfaces made by microscale acid etching through a sintering process were examined. Regarding their obtained results, the use of hierarchical multiscale modulated porous surfaces would lead to 200% augmentation in CHF in comparison with the polished plain surface. Silver-nanoparticles suspended in de-ionized water with plate heaters were scrutinized by Jo et al. [70]. They demonstrated the shift of boiling curves to the right side. Additionally, they inferred the CHF enhancement by the reduction in the nanoparticles size. Experimental studies on the surfaces of micro-channeled with micro-porous coatings were performed by Patil et al. [71] and Jaikumar et al. [72] via diverse methodologies (sintering, electrodeposition, etc.). In these investigations, remarkable enhancements in CHF were observed. Sarafraz et al. [73] compared the thermal performance of pool boiling in alumina nano-fluid for both plain and micro-structured surfaces. They pointed out the augmented rate of CHF enhancement by increasing nanofluid mass concentration, nanoparticle size and interspace between microstructures. Nucleate pool boiling of Al 2 O 3 -based aqueous nanofluid was prospected by Shahmoradi et al. [74] on a flat plate heater. They declared that by increasing CHF, HTC of nanofluid experiences a descending trend. Moreover, they reported an increase in CHF amelioration by a rising volume fraction of nanoparticles. Jaikumar and Kandlikar [47] revealed different nucleating and non-nucleating regions, which facilitates the separation of liquid-vapor pathways for the detachment of the vapor bubbles and returning liquid that would lead to CHF enhancement.
Heat transfer improvement with surfactants and polymeric additives was addressed in the comprehensive reviews of Cheng et al. [75] and Wasekar and Manglik [76]. It is worth mentioning that the use of nanofluids for the CHF amelioration of both pool boiling and forced convection was listed in the recent reviews. In a review by Ciloglu and Bolukbasi [77], the investigations associated with pool boiling were listed up to 2015. In the mentioned review, it was aimed to justify the conflicting findings for nanofluids pool boiling. Other reviews regarding the pool boiling of nanofluids, such as Yu et al. [ [90], and Devendiran and Amirtham [91], were concentrated on the preparation methods for nanofluid, instruments and methods for inspection, thermo-physical characteristics, as well as applications.

Nanoparticle Deposition Method
The deposition of nanoparticles is defined as the attachment of nanoparticles to the solid surfaces, which are called substrates, to form coatings of nanoparticles. The properties of different materials' surfaces will change if nanoparticles are deposed on them. This may be attributed to the improvement of roughness, porosity, wettability and topography of the surface by nanoparticle deposition.
Coatings can be classified into multilayer/monolayer and organized/unorganized structures considering the coating process. Nanoparticle deposition is usually encountered with some challenges according to the physical properties of nanoparticles. The available methods for coating with nanoparticle deposition can be distinguished through different parameters, including the ability of controlling particle density and layer thickness, the ability of using different particles, the complexity of the method as well as the required instrumentation. Various materials can be used for generating nanoparticles, such as metals, oxides of nanoparticles ceramics and polymers. The nanoparticles made of metal oxide are of great importance in chemistry, physics, mechanics and materials science. Among the applications of metal oxide nanoparticles in technology, production of microelectronic circuits, sensors, piezoelectric devices, fuel cells and anti-corrosion coatings can be mentioned. The unique physical and chemical properties of the metal oxide nanoparticles are related to the limited size and high density in corner or edge surface sites.
Stability is a significant issue in nanoparticles considering their tendency to lower their high surface energy attributed to their high surface-to-bulk ratio. Nanoparticles can reach stability through two processes: the absorption of molecules from the environments and lowering the surface area through coagulation and agglomeration [92,93]. The tendency of a nanoparticle to coagulate, which is unwanted in general, can be controlled by modifying the surface layer. In order to prevent coagulation in liquid medium, appropriate ligand molecules are generally attached to the nanoparticle surface, so that they can be solved in suitable solvents. The nanoparticle deposition on the layer leads to an improvement in the characteristics of the heated surface, including surface wettability, roughness and capillary wicking performance, which consequently results in the CHF amelioration.
Thermal conductivity of suspended nanoparticles is responsible for such amelioration in boiling heat transfer. Alumina, silica, titania and zirconia nanofluids have been widely investigated for CHF Energies 2020, 13, 4026 7 of 49 amelioration. The use of alumina and copper nanofluids in the experimental tests of pool boiling experiments resulted in the respective CHF amelioration of 171% and 176% compared to that of pure water [94]. Correspondingly, HTC was also augmented due to the nanofluid boiling. From the investigations concerning the nanoparticle deposition on the heater surface after boiling, it was inferred that heat flux improvement can be achieved at lower concentrations. However, for the higher concentrations of alumina nanoparticles, a decrease was observed in CHF enhancement. It can be concluded that increasing nanoparticle concentration would lead to an increase in the nanoparticle deposition over the surface with a thickness of a few micrometers. This increase is followed by a reduction in the active nucleation sites. From all of these studies, it can be inferred that the use of a hydrodynamic instability hypothesis is not sufficient for justifying CHF amelioration in nanofluids. An instability hypothesis is based on water contact angle measurement, and it is a challenge to quantify the amelioration of heat transfer over thin wire surfaces [95]. From the performed observation, the insufficiency of the existing theoretical models can be concluded for the prediction of the CHF with surface effects. The basic scheme of nanoparticle layer deposition is shown in Figure 3. for CHF amelioration. The use of alumina and copper nanofluids in the experimental tests of pool boiling experiments resulted in the respective CHF amelioration of 171% and 176% compared to that of pure water [94]. Correspondingly, HTC was also augmented due to the nanofluid boiling. From the investigations concerning the nanoparticle deposition on the heater surface after boiling, it was inferred that heat flux improvement can be achieved at lower concentrations. However, for the higher concentrations of alumina nanoparticles, a decrease was observed in CHF enhancement. It can be concluded that increasing nanoparticle concentration would lead to an increase in the nanoparticle deposition over the surface with a thickness of a few micrometers. This increase is followed by a reduction in the active nucleation sites. From all of these studies, it can be inferred that the use of a hydrodynamic instability hypothesis is not sufficient for justifying CHF amelioration in nanofluids. An instability hypothesis is based on water contact angle measurement, and it is a challenge to quantify the amelioration of heat transfer over thin wire surfaces [95]. From the performed observation, the insufficiency of the existing theoretical models can be concluded for the prediction of the CHF with surface effects. The basic scheme of nanoparticle layer deposition is shown in Figure 3.

Objectives of the Present Review
This review paper aims to highlight the essential concepts concerning recent advances in critical heat flux amelioration of pool boiling surfaces to access the vast literature available on research methodologies. The study is concentrated on the modification of surface characteristics through the use of metal oxide nanoparticle deposition process. By reviewing the related literature, the use of metal oxide nanoparticle for CHF enhancement, nanoparticle deposition methods and the mechanism of CHF enhancement in nanoparticle deposited surfaces are analyzed. Additionally, the existing gap in the research, future works, enhancement potentials, and limitations as well as their possible industrial implementation are debated in the present review article.

Nanoparticle Material
The nanoparticle dispersion in the base fluid plays a key role in the CHF amelioration. Diverse metals, metal oxides, CNTs and diamond have been employed for pool boiling heat transfer purposes. The CHF amelioration in some pool boiling investigations is summarized in Table 1. As it

Objectives of the Present Review
This review paper aims to highlight the essential concepts concerning recent advances in critical heat flux amelioration of pool boiling surfaces to access the vast literature available on research methodologies. The study is concentrated on the modification of surface characteristics through the use of metal oxide nanoparticle deposition process. By reviewing the related literature, the use of metal oxide nanoparticle for CHF enhancement, nanoparticle deposition methods and the mechanism of CHF enhancement in nanoparticle deposited surfaces are analyzed. Additionally, the existing gap in the research, future works, enhancement potentials, and limitations as well as their possible industrial implementation are debated in the present review article.

Nanoparticle Material
The nanoparticle dispersion in the base fluid plays a key role in the CHF amelioration. Diverse metals, metal oxides, CNTs and diamond have been employed for pool boiling heat transfer purposes. The CHF amelioration in some pool boiling investigations is summarized in Table 1. As it can be seen, the reported CHF enhancement was different in various studies, even for the same substance. For example, despite the extensive use of alumina nanoparticles in most of the investigations, the CHF amelioration was different. In fact, the nanoparticle material, size and concentration as well as the base fluid have different impacts on CHF. In spite of controversies regarding the effectuality of water-based oxide nanofluids, they have recently received a considerable amount of attention. In contrast, metal nanofluids have gained less attention, which could be attributed to difficulties due to their suspension within the base fluid. Kathiravan et al. [120,121] explored the nucleate pool boiling in a water-based Cu nanofluid (0.25-1.0 wt%) and reported the highest CHF amelioration (~50%) relative to the pure water. Incorporation of 9.0 wt% SDS surfactant, however, declined the CHF amelioration to one-third. Furthermore, the nucleate boiling HTC showed a descending trend by the enhancement of nanoparticle content. Based on Zhou [122], single-phase convection was increased; the nucleate boiling heat transfer, however, was compromised with acetone-based Cu nanofluids. Krishna et al. [123] reported a reduction in the HTC of water-based Cu nanofluids (0.01 vol.%) at the low-heat-flux nucleate boiling region in comparison with the pure water. The HTC was also enhanced by raising the concentration to 0.1 vol.%. Therefore, concentration enhancement can oppositely influence at the high-heat-flux region. Such a complex trend can be assigned to various parameters such as thermal conductivity of the microlayer, sorption layer formation and active nucleation site densities. Kole and Dey [124] observed a simultaneous increase in nucleate boiling HTC and CHF by raising the Cu content (0.005-0.5 wt%) in the water-based Cu nanofluids. Moreover, they recorded higher heating rise in the copper-made surfaces compared to the brass surfaces. The nucleate boiling HTC showed an amelioration by an increase in surface roughness (R a ) from 0.06 to 0.22 µm; however, a further increase in roughness by 0.7 µm, resulted in a decline, which can be explained by deactivation of nucleation sites by nanoparticle deposition.
Shi et al. [125] explored the impact of Fe nanoparticles and observed that the nucleate boiling heat transfer depended on nanoparticles thermal properties, size and concentration, in addition to heat flux. The rise in efficient thermal conductivity of the fluid as well as the decline in surface tension at the presence of nanoparticles are associated with improved nucleate boiling heat transfer. Besides, nanoparticles' entrapment in the surface voids will further smoothen the surface and hence reduce the nucleate boiling HTC. So, they concluded that any changes in heat transfer performance are the consequence of the two mentioned factors. The boiling performance of Al 2 O 3 and Cu water-based nanofluids was compared in the work of Cieslinski and Kaczmarczyk [126] on copper and stainless steel tubes. In the case of smooth copper tube, nucleate boiling HTC did not show dependence on the nanoparticle material; HTC was, however, reduced by nanoparticle concentration enhancement from 0.01 to 1.0 wt%. Besides, concerning the stainless-steel tube, nucleate boiling HTC exhibited improvement with both nanoparticle type and concentration. The substance used in heating surface seems to be ineffective on the boiling heat transfer as tested for 0.1 wt% Cu nanofluid.
The type of nanoparticles and the base fluids are demonstrated in Figure 4. Clearly, Al 2 O 3 , CuO and TiO 2 accounted for 33.1%, 14.5%, and 11.3%, respectively. They were followed by SiO 2 (8.2%), ZnO (6.6%), Fe 3 O 4 (3.1%) and ZrO 2 (3.3%), and that 85.1% of the applied base fluid was water, in which the EG (7.7%) and refrigerants (5.2%) had the subsequent ranks. As suggested Figure 4, Al 2 O 3 and CuO are the most generally utilized materials nanoparticles in water, which has the highest rate of application as the base fluid.

Nanoparticle Thermo-Physical Properties
Thermo-physical properties refer to properties with an impact on the heat transfer and storage. These properties may vary by the state variables, such as temperature, pressure and mixture components with no alternation in the chemical nature of the material. Density, thermal conductivity, heat capacity, thermal expansion, thermal radiative properties, vaporization heat, boiling point temperature, melting point, viscosity, mass and thermal diffusion coefficients, the speed of sound, surface and interfacial tension in fluids can be classified as the thermo-physical properties of a substance [127][128][129].
The type of nanoparticles and the base fluids are demonstrated in Figure 4. Clearly, Al2O3, CuO and TiO2 accounted for 33.1%, 14.5%, and 11.3%, respectively. They were followed by SiO2 (8.2%), ZnO (6.6%), Fe3O4 (3.1%) and ZrO2 (3.3%), and that 85.1% of the applied base fluid was water, in which the EG (7.7%) and refrigerants (5.2%) had the subsequent ranks. As suggested Figure 4, Al2O3 and CuO are the most generally utilized materials nanoparticles in water, which has the highest rate of application as the base fluid.

Nanoparticle Thermo-Physical Properties
Thermo-physical properties refer to properties with an impact on the heat transfer and storage. These properties may vary by the state variables, such as temperature, pressure and mixture components with no alternation in the chemical nature of the material. Density, thermal conductivity, heat capacity, thermal expansion, thermal radiative properties, vaporization heat, boiling point temperature, melting point, viscosity, mass and thermal diffusion coefficients, the speed of sound, surface and interfacial tension in fluids can be classified as the thermo-physical properties of a substance [127][128][129].
Numerous studies have addressed the impact of nanofluids' thermo-physical properties on their boiling performance. Nanofluid properties such as heat capacity, density and viscosity showed no sensitivity to the nanoparticles' presence; the thermal conductivity and the surface tension, however, exhibited a significant dependence on the nanoparticles. Studies have revealed that the thermal conductivity of nanofluids is greater than the conventional fluids [106,130], which can be attributed to the large specific surface area of small-sized particles, as well as an enhancement in the fluid heat capacity, Brownian motion and interfacial liquid layering [120]. As Numerous studies have addressed the impact of nanofluids' thermo-physical properties on their boiling performance. Nanofluid properties such as heat capacity, density and viscosity showed no sensitivity to the nanoparticles' presence; the thermal conductivity and the surface tension, however, exhibited a significant dependence on the nanoparticles. Studies have revealed that the thermal conductivity of nanofluids is greater than the conventional fluids [106,130], which can be attributed to the large specific surface area of small-sized particles, as well as an enhancement in the fluid heat capacity, Brownian motion and interfacial liquid layering [120]. As stated by Das et al. [130] and Wen and Ding [131], an augmentation in the thermal properties of a stable nanofluid can elevate its boiling heat transfer. In fact, the conduction heat transfer plays a key role in the thin fluid layer above the heating surface [132].
Lee et al. [133], Kwak and Kim [134] and Suganthi and Rajan [135] expressed the resistance of metal oxides toward oxidation which can result in their chemical stability. Moreover, some metal oxides have densities inferior to their corresponding metals which can be helpful in resolving the problem of particle settling during nanofluid formulation. This can explain the popularity of metal oxides in nanofluid formulation in spite of their inferior thermal conductivity (compared to metals). Owing to its superior thermal conductivity and low density, alumina has been widely explored for nanofluid applications.
Liu et al. [103] empirically scrutinized the nucleate boiling heat transfer in the CuO/water nanofluid at different pressures and nanoparticle contents. Their results revealed a significant impact of operation pressure on nucleate boiling profile. Furthermore, they reported an increase in the HTC and CHF of nanofluids by pressure reduction when compared to the pure water. A slow enhancement of HTC and CHF by nanoparticle content elevation (until its optimal value, 1.0 wt%) was also observed for the entire tested pressures. Beyond the optimal concentration, HTC showed a deteriorating trend, whereas CHF maintained a fixed value. This phenomenon was explored by investigating the impact of nanofluid thermo-physical properties and variations in the surface features. In the case of former parameter, nanofluid thermal conductivity, viscosity and the surface tension were considered, which showed 102%, 101% and 88% variation in comparison with the pure water, respectively. These data were then employed in the Kutateladze correlation [136]: where σ shows the liquid vapor interface tension and ρ l and ρ v represent the liquid and vapor densities, respectively; ν l denotes the liquid kinematic viscosity and P stands for the pressure. Furthermore, g signifies the gravitational acceleration, λ represents the fluid thermal conductivity, h f g is the latent evaporation heat and, finally, Pr l denotes the saturated liquid Prandtl number. Variations in the fluid thermo-physical features were found responsible for the HTC degradation in nanofluid.
Concerning the latter case, it was found that as the nanoparticles' size is one or two orders of magnitude smaller than the surface roughness, their deposition will decline the surface roughness, which in turn may decline the contact leading to a reduction in boiling HTC as well as an enhancement in CHF. Besides, it was assumed that a steady porous layer is only present in the sub-atmospheric pressure. Moreover, the bubble formation under sub-atmospheric pressures is largely different from the case at atmospheric pressure. Yang and Liu [22] empirically explored the pool boiling heat transfer on a copper bar under atmospheric and sub-atmospheric pressures for the functionalized and traditional SiO 2 /water nanofluids and reported a substantial difference in their boiling properties. According to them, traditional nanofluid showed an increasing trend in CHF and deteriorating pattern in HTC, whereas a slight increase was observed in the HTC of the functionalized nanofluid, while its CHF remained unchanged. In contrast to the traditional nanofluids, no nanoparticles' layer covered the heated surface during the boiling in case of a functionalized nanofluid. Such dramatic differences could be assigned to modifications in the thermo-physical properties of the functionalized nanofluid, for instance, a slight increase in its thermal conductivity, accompanied with substantial enhancement in its viscosity and a decline in its surface tension.
Fe nanoparticles were employed to formulate an EG-based nanofluid by Hong et al. [137]. According to them, Fe resulted in a higher increment in the thermal conductivity of the nanofluid when compared to Cu nanofluids.
Quan et al. [138] recently published an article on the impact of nanoparticle wettability on nucleate boiling performance. According to Figure 5, nanoparticles with moderate hydrophilicity are adsorbed to the vapor-liquid interface, preventing from liquid drainage among the adjacent bubbles. Thus, the bubble coalescence will be resisted; hence, a decline will be observed in the bubble departure diameter as well as an amelioration in the nucleate boiling HTC and CHF. Highly hydrophilic nanoparticles, however, were not adsorbed to the interface. Therefore, the bubble coalescence will remain unchanged.
Energies 2020, 13, x FOR PEER REVIEW 12 of 50 are adsorbed to the vapor-liquid interface, preventing from liquid drainage among the adjacent bubbles. Thus, the bubble coalescence will be resisted; hence, a decline will be observed in the bubble departure diameter as well as an amelioration in the nucleate boiling HTC and CHF. Highly hydrophilic nanoparticles, however, were not adsorbed to the interface. Therefore, the bubble coalescence will remain unchanged. As demonstrated in Figure 5, nanoparticle wettability affected the deposited layer morphology in which the high hydrophilic nanoparticle layers were relatively smooth, with uniform particle dispersion while the layers out of nanoparticles with medium hydrophilicity possessed higher roughness and irregularity. As demonstrated in Figure 5, nanoparticle wettability affected the deposited layer morphology in which the high hydrophilic nanoparticle layers were relatively smooth, with uniform particle dispersion while the layers out of nanoparticles with medium hydrophilicity possessed higher roughness and irregularity.
According to Bourdon et al. [139][140][141], surface hydrophobicity will decline the incipience boiling superheat of water. Jo et al. [142] also tested the identical nano/microstructures (micro-patterns surrounded by nano-grass structure) with two different sorts of surface wettability (i.e., silicon-oxidized, with a contact angle of~0 • , and Teflon-coated, with contact angle of 162 • ), where both wettability conditions managed to significantly improve the CHF. Although the Teflon-coated sample exhibited earlier boiling incipience and larger HTC, its CHF was less than the silicon-oxidized specimen. In the work of Kim et al. [44,143], microstructural hybrid surfaces modified by nano-rods were applied. The application of these nano-rods substantially enhanced the surface wettability and liquid spreadability in such a way that the CHF improvement of the hybrid surface was equal to the sum of enhancements achievable by use of microstructures and nano-rods.
According to Sun et al. [144], oxygen vacancies in the crystal structure of TiO 2 can enhance the wettability as these vacancies tend to absorb water molecules. Such phenomenon can be observed under ultraviolet (UV) exposure, heat treatment and Ar + sputtering. Similarly, Takata et al. [145,146] reported a 200% improvement in the CHF of water and nucleate boiling heat transfer in the Cu-sputtered TiO 2 surface in comparison with the bare surface. Furthermore, the Leidenfrost temperature showed a significant enhancement by the decline of the contact angle. CHF improvement in water was also detected by Liaw and Dhir [147] upon wettability increment in a copper surface. They controlled the contact angle through variation in the surface oxidation degree. In contrast, Maeng et al. [148] stated that plasma sputtering induced enhanced wettability in TiO 2 on a nickel surface, which did not lead to the creation of a porous nanostructure; it also did not affect the CHF.
Kim et al. [149] recently indicated that temperature-induced wetting transitions in TiO 2 coatings can substantially influence the CHF. Their developed surface exhibited hydrophobic behavior at low wall temperatures, which altered to hydrophilic characteristics by the temperature rise. For instance, temperature increment from 100 to 200 • C caused a decline in the water contact angle from 83.1 to 32.7 • . In addition to these temperature effects, a time dependency was observed in the wetting transitions [150]; thus, the boiling performance could be further enhanced by boiling time augmentation at high heat flux zones.
The thermal transport characteristics (i.e., thermal conductivity, density, specific heat, viscosity and surface tension) of the nanofluids are of crucial significance in the single-phase heat transfer purposes. Interestingly, numerous studies (with a few exceptions) have stated that these properties are not highly effective on CHF enhancement. This feature, however, should be first explored prior to examinations for the pool boiling heat transfer purposes.
The major underlying reason for incorporation of nanoparticles in the conventional fluid is the thermal conductivity enhancement. The thermal conductivity values of some common metallic and non-metallic solids as well as liquids are presented in Figure 6. Experimental or analytical methods can be applied to calculate the thermal conductivity of nanofluids. Experimental approaches such as transient hot wire [151], steady state parallel plate [152] and temperature oscillation methods [153] have been widely employed to quantify the nanofluid thermal conductivity. Among the analytical techniques for the prediction of thermal conductivity [154][155][156][157], Hamilton-Crosser [158] and Yu and Choi [159] are almost comparable while Maxwell model [160] underestimates the values with difference increase by nanofluid concentration increment [161]. Interestingly, the experimental results were far higher than the analytical outcomes reflecting the presence of a neglected heat transport mechanism in nanofluid.
Energies 2020, 13, 4026 13 of 49 [153] have been widely employed to quantify the nanofluid thermal conductivity. Among the analytical techniques for the prediction of thermal conductivity [154][155][156][157], Hamilton-Crosser [158] and Yu and Choi [159] are almost comparable while Maxwell model [160] underestimates the values with difference increase by nanofluid concentration increment [161]. Interestingly, the experimental results were far higher than the analytical outcomes reflecting the presence of a neglected heat transport mechanism in nanofluid. For a deeper understanding of the nanofluid thermal conductivity and CHF improvement, You et al. [97] incorporated a very small nanoparticle content to enhance the CHF. Their employed nanoparticle concentration was substantially smaller than the practical nanoparticle content used to enhance the nanofluid thermal conductivity. Based on their observations, CHF improvement cannot be assigned to the nanofluids thermal conductivity, which was in line with the previous studies on pool boiling heat transfer [163,164]. On the other hand, Ahmed et al. [165] stated the more profound impact of nanofluids thermal conductivity compared to the rate of nanoparticles deposition at lower concentrations, which may result in CHF increment. Their observation, however, requires deeper exploration regarding the particles' size and shape, nanofluid preparation For a deeper understanding of the nanofluid thermal conductivity and CHF improvement, You et al. [97] incorporated a very small nanoparticle content to enhance the CHF. Their employed nanoparticle concentration was substantially smaller than the practical nanoparticle content used to enhance the nanofluid thermal conductivity. Based on their observations, CHF improvement cannot be assigned to the nanofluids thermal conductivity, which was in line with the previous studies on pool boiling heat transfer [163,164]. On the other hand, Ahmed et al. [165] stated the more profound impact of nanofluids thermal conductivity compared to the rate of nanoparticles deposition at lower concentrations, which may result in CHF increment. Their observation, however, requires deeper exploration regarding the particles' size and shape, nanofluid preparation procedure as well as its stability. A recent investigation by Ahn et al. [166] indicated that a small enhancement in the thermal conductivity of reduced graphene oxide (RGO) colloids combined to surface characteristics can significantly affect the CHF delay. Yang and Liu [22] study also reported a drastic thermal conductivity increment in the functionalized silica nanofluids. Such an increment did not affect the CHF improvement, as the resulting CHF did not differ from that of pure water. Therefore, it can be inferred that the nanofluid thermal conductivity can significantly affect the CHF enhancement at lower nanoparticle contents. A comprehensive investigation on a wide spectrum of nanoparticles with various sizes is inevitable for a definite statement on CHF.
The difference in the density and specific heat values was negligible or even identical to those of pure water at the entire tested nanoparticle contents [163,167]. A comparison was made between the analytical and empirical specific heat and density of nanofluids by Pantzali et al. [168], which revealed 5% and 2% deviations in density and specific heat capacity, respectively. Therefore, nanofluids' density and specific heat can be determined by the following equations [169] in the majority of the pool boiling heat transfer studies.
where subscripts of n f and b f are respective expressions for nanofluid and base fluid, while ρ, φ and C P denote density, concentration and specific heat capacity, respectively. In spite of a negligible variation in the analytically predicted density and specific heat values, a limited number of investigations reported a decline in specific heat by elevation of the nanoparticles content [164,170]. Regarding the lack of a systematic study addressing the impact of density and specific heat on the CHF augmentation, a detailed investigation sounds essential to deeply explore the effect of these parameters at high nanoparticles' concentration.
The viscosity of nanofluids should be also addressed to elucidate its role in CHF improvement. The nanofluid viscosity can be determined by Brinkman equation [171].
in which subscripts of n f and b f are respective expressions for nanofluid and base fluid, while µ and φ denote dynamic viscosity and concentration, respectively. The mentioned equation holds for charge-free particles; it, however, underestimates the viscosity. This equation was modified by the inclusion of electro-viscous force emerged during electro stabilization [172].
in which subscripts of n f and b f are respective expressions for nanofluid and base fluid, while µ and φ denote dynamic viscosity and concentration, respectively. The nanofluid viscosity generally grows by nanoparticle concentration and declines with temperature rise, as suggested by various studies [131,173,174]. The underlying mechanisms are, however, unknown. Some studies examined the relationship between nanofluid viscosity and nanoparticle size. The viscosity can be well estimated by classical models and other formulae at lower volume fractions. The estimation, however, fails at different temperatures [175].
Concerning the impact of viscosity on CHF increment, Kwark et al. [106] evaluated the viscosity and observed a minimum difference between the viscosity of nanofluid and pure water, which could be attributed to low nanoparticle content. An incremental trend was found in the kinematic viscosity by concentration enhancement [164], which could be assigned to the electro stabilization-induced and electro viscous forces. Das et al. [130] assessed the viscosity of alumina-water nanofluid as a function of the shear rate. They recorded an ascending trend in viscosity by the growth of nanoparticles' content, while it reduced by temperature. The viscosity showed an incremental behavior by shear rate. They concluded that such patterns cannot be exploited for deteriorating the boiling heat transfer features. Currently, the viscosity values of the pool boiling studies are similar to those of pure water in all the evaluated concentrations [176,177]. Further studies are required in the following areas: • Empirical and analytical studies to establish a database for the nanofluid viscosity.

•
Application of more accurate models and/or correlations for the simulation or analysis of the nanofluid heat transfer.
The nanofluid surface tension is another key factor in CHF increase that has been extensively addressed in recent decade. Vafaei et al. [178] studied the impact of bismuth telluride (Bi 2 Te 3 ) nanoparticle size and content on the efficient gas-liquid surface tension. Surprisingly, the surface tension declined by the enhancement of nanoparticles' content until a critical value. A further increase of nanoparticle concentration resulted in an inverse behavior. Furthermore, at similar mass concentrations, the nanofluid encompassing 10.4 nm nanoparticles exhibited higher surface tension compared to the one including 2.5 nm nanoparticles, reflecting the impact of nanoparticles size and concentration on the liquid-gas surface tension.
An empirical research on CNT nanofluids by Kumar and Milanova [179] indicated the surface tension relaxation (the difference between surface tension of nanofluid and base fluid) as the main underlying mechanism in the CHF delay. The optimization of CNT and surfactant concentrations for the heat flux maximization requires more studies. According to Kathiravan et al. [121], the incorporation of surfactant in Cu-water nanofluid will suppress the CHF increment as the result of a declined surface tension. Jeong et al. [180] also reported a decline in surface tension by Tri-sodium phosphate (TSP) and nanofluid contents' increase. Indeed, a reduction in surface tension will lower the CHF values. The impact of nanoparticle deposition on the surface (which reduced the contact angle and hence improved the wettability) could explain the augmentation in CHF of TSP and nanofluids. Contrarily, the nanofluid surface tension remained nearly unchanged relative to pure water at all evaluated concentrations [176,177].
The surface tension generally shows a decreasing trend with a rise in nanoparticles' concentrations and temperature. The following fields should be more investigated to gain a deeper insight into the CHF enhancement: (1) synergetic influence of surface tension and surface wettability; (2) synergetic impact of nanoparticle concentrations and surfactant type.
To conclude, variations in the thermal transport characteristics do not significantly affect the CHF delay. Its role is, however, more pronounced at lower concentrations. In this regard, thermal transport parameters should be further investigated at lower volume fractions and finer particle sizes for a more reliable conclusion on the CHF enhancement.

Nanoparticle Shape and Size
The morphology of the individual particles has a significant impact on the boiling CHF and HTC. This is related to the fact that different parameters, including post-sintering pore shape, permeability, surface roughness, coating's effective conductivity and diffusivity are affected by particles morphology. The effect of particles' morphology on pool boiling has not been addressed in the literature. However, in some other applications for heat transfer, the influence of this parameter on geometrical and thermo-physical properties has been studied [181][182][183]. Deng et al. [184] investigated the impact of particle parameters such as shape and size on flow boiling performance of microchannels. The studied microchannel was a porous copper heat sinks being made of particles with spherical and irregular shapes (50-150 µm effective diameter). According to the obtained results, the highest HTC was observed for the irregular-shaped particles with smaller sizes (50-75 µm). However, for the spherical-shaped particles with smaller sizes, the lowest HTC was achieved.
It can be inferred that the nanoparticle size plays a significant role in the amelioration of pool boiling CHF [185]. As the size of nanoparticles increase, better boiling heat transfer performance would be achieved for the nanofluids [186,187].
The effect of nanoparticle concentration in the base fluid was investigated by Kole et al. [110]. Two diverse concentrations of ZnO nanoparticles with a size of 30-40 nm in ethylene glycol base fluid were evaluated. It was observed that, after pool boiling, nanoparticles were deposited over the heater surface, which prevented the active nucleation sites, and, therefore, HTC decreased. By measuring the CHF values through the use of thin Constantan wire, they observed the significant increase in CHF by increasing the concentration of ZnO. They reported a maximum CHF enhancement of 117% for ZnO nanoparticle volume fractions of 2.6%. Moreover, they inferred that the CHF amelioration is associated with the increased surface roughness of the heating wire as a result of nanoparticle deposition.
Peng et al. [188] performed an empirical study on the impact of Cu nanoparticle dimensions on nucleate pool boiling heat transfer of R113/oil mixture. According to the obtained results, the maximum increase of 23.8% was achieved for HTC by reducing the nanoparticle dimension from 80 to 20 nm. Additionally, smaller Cu nanoparticles would lead to a higher nucleate pool boiling HTC.
An empirical analysis of nanoparticle size impact on pool boiling was conducted in some studies [189,190]. However, satisfactory evidence and validation were not provided for the obtained results. The influence of size and concentration of the nanoparticle on the boiling performance was investigated by Hu et al. [191] for SiO 2 nanofluid. They reported the increasing trend of HTC by decreasing the nanoparticle size from 120 to 84 nm.

Nanoparticle Concentration
The impact of diverse concentrations of Al 2 O 3 /water nanofluids on the pool boiling of a smooth horizontal flat surface was studied by Bang and Chang [99]. According to the obtained results, as the concentration of the nanoparticles increased, worse heat transfer performance was observed compared to pure water. However, the CHF enhancements for horizontal and vertical pool boiling were reported to be 32% and 13%, respectively.
In an experimental investigation by Golubovic et al. [105], the effect of the concentration, size and type of nanoparticles on CHF was evaluated under saturated conditions. The Al 2 O 3 -BiO 2 /water nanofluids were used on a Ni-Cr wire, which was submerged in a chamber with dimensions of 70 × 70 × 100 mm. Based on Figure 7 CHF enhances by the particle concentration increment. However, after reaching the maximum heat flux, a further increase in the particle concentration does not show any significant impact on CHF. The maximum CHF enhancements were 50% and 33% for alumina and bismuth oxide nanofluids, respectively. investigated by Hu et al. [191] for SiO2 nanofluid. They reported the increasing trend of HTC by decreasing the nanoparticle size from 120 to 84 nm.

Nanoparticle Concentration
The impact of diverse concentrations of Al2O3/water nanofluids on the pool boiling of a smooth horizontal flat surface was studied by Bang and Chang [99]. According to the obtained results, as the concentration of the nanoparticles increased, worse heat transfer performance was observed compared to pure water. However, the CHF enhancements for horizontal and vertical pool boiling were reported to be 32% and 13%, respectively.
In an experimental investigation by Golubovic et al. [105], the effect of the concentration, size and type of nanoparticles on CHF was evaluated under saturated conditions. The Al2O3-BiO2/water nanofluids were used on a Ni-Cr wire, which was submerged in a chamber with dimensions of 70 × 70 × 100 mm. Based on Figure 7 CHF enhances by the particle concentration increment.
However, after reaching the maximum heat flux, a further increase in the particle concentration does not show any significant impact on CHF. The maximum CHF enhancements were 50% and 33% for alumina and bismuth oxide nanofluids, respectively. The CHF performance of TiO2-Al2O3-SiO2/water nanofluids was investigated by Kim and Kim [192] by the use of Ni-Cr wire with diameter of 0.2 mm. As can be seen in Figure 8, they reported The CHF performance of TiO 2 -Al 2 O 3 -SiO 2 /water nanofluids was investigated by Kim and Kim [192] by the use of Ni-Cr wire with diameter of 0.2 mm. As can be seen in Figure 8, they reported CHF enhancement by the use of nanofluids. Besides, they declared the strong dependency of CHF on the nanoparticle type and concentration.
The effect of different nanofluids, including alumina, titania and magnetite nanofluids on improving the CHF was explored by Lee et al. [193]. The experiments were conducted in saturated conditions and under atmospheric pressure by the use of Ni-Cr wire, which was embedded in a rectangular stainless-steel vessel with the dimensions of 250 × 100 × 230 mm. Regarding the results, they reported the maximum CHF by concentration increment for the magnetite-water nanofluid. They justified the CHF increase by the increased surface wettability as a result of a higher frequency of departing bubbles in magnetite-water nanofluid (MWNF).
In an experimental study, Jung et al. [194] evaluated the effect of a stabilizer on CHF and boiling HTC. The considered nanofluid was Al 2 O 3 /water nanofluids, which was studied upon a plate Cu heater in the rectangular test chamber (1650 mL). In addition, polyvinyl alcohol (PVA) was utilized as a steric stabilizer. They stated the lower boiling HTC despite the higher CHF of nanofluids compared to the base fluid by increasing the nanoparticle concentration. Moreover, they declared that up to a critical concentration of nanoparticle (10 −3 vol.% in this investigation), a higher effective boiling surface area was formed by the use of a stabilizer. However, further increase in the concentration would lead to a reduction in the effective boiling surface area. They found that the effective boiling surface area was increased by the deposition of the nanoparticles, which would result in the CHF amelioration, consequently. CHF enhancement by the use of nanofluids. Besides, they declared the strong dependency of CHF on the nanoparticle type and concentration. The effect of different nanofluids, including alumina, titania and magnetite nanofluids on improving the CHF was explored by Lee et al. [193]. The experiments were conducted in saturated conditions and under atmospheric pressure by the use of Ni-Cr wire, which was embedded in a rectangular stainless-steel vessel with the dimensions of 250 × 100 × 230 mm . Regarding the results, they reported the maximum CHF by concentration increment for the magnetite-water nanofluid. They justified the CHF increase by the increased surface wettability as a result of a higher frequency of departing bubbles in magnetite-water nanofluid (MWNF).
In an experimental study, Jung et al. [194] evaluated the effect of a stabilizer on CHF and boiling HTC. The considered nanofluid was Al2O3/water nanofluids, which was studied upon a plate Cu heater in the rectangular test chamber (1650 mL). In addition, polyvinyl alcohol (PVA) was utilized as a steric stabilizer. They stated the lower boiling HTC despite the higher CHF of nanofluids compared to the base fluid by increasing the nanoparticle concentration. Moreover, they declared that up to a critical concentration of nanoparticle (10 −3 vol.% in this investigation), a higher effective boiling surface area was formed by the use of a stabilizer. However, further increase in the concentration would lead to a reduction in the effective boiling surface area. They found that the effective boiling surface area was increased by the deposition of the nanoparticles, which would result in the CHF amelioration, consequently.
In another study by Shahmoradi et al. [74], an increase in the surface roughness was reported by the increase in nanoparticle concentration as a result of higher particle size compared to the surface roughness (as can be seen in atomic force microscopy (AFM) images in Figure 9). Besides, the CHF enhancement was obtained according to the improved wettability. In another study by Shahmoradi et al. [74], an increase in the surface roughness was reported by the increase in nanoparticle concentration as a result of higher particle size compared to the surface roughness (as can be seen in atomic force microscopy (AFM) images in Figure 9). Besides, the CHF enhancement was obtained according to the improved wettability.  As was obvious from Figure 10, a significant impact of nanofluid concentration on CHF augmentation was reported by You et al. [97]. In this figure, it can be observed that the CHF values for nanofluids are almost three times the corresponding value for water at a lower saturation temperature of T = 60 ℃ (Al O ) . In addition, intensive CHF amelioration was observed by increasing the nanoparticle concentration at low concentrations, which becomes almost constant at concentrations higher than 0.01 g/l. A similar trend was observed in several investigations [101,103,192,195]. In the study by Yang and Liu [22], no amelioration in CHF was reported by increasing the concentration of functionalized silica nanoparticles [21] from 0.5 to 2.5 wt%. This was justified by the grafting silanes to the surface of silica nanoparticles. As was obvious from Figure 10, a significant impact of nanofluid concentration on CHF augmentation was reported by You et al. [97]. In this figure, it can be observed that the CHF values for nanofluids are almost three times the corresponding value for water at a lower saturation temperature of T sat = 60 • C (Al 2 O 3 ). In addition, intensive CHF amelioration was observed by increasing the nanoparticle concentration at low concentrations, which becomes almost constant at concentrations higher than 0.01 g/L. A similar trend was observed in several investigations [101,103,192,195]. In the study by Yang and Liu [22], no amelioration in CHF was reported by increasing the concentration of functionalized silica nanoparticles [21] from 0.5 to 2.5 wt%. This was justified by the grafting silanes to the surface of silica nanoparticles. augmentation was reported by You et al. [97]. In this figure, it can be observed that the CHF values for nanofluids are almost three times the corresponding value for water at a lower saturation temperature of T = 60 ℃ (Al O ) . In addition, intensive CHF amelioration was observed by increasing the nanoparticle concentration at low concentrations, which becomes almost constant at concentrations higher than 0.01 g/l. A similar trend was observed in several investigations [101,103,192,195]. In the study by Yang and Liu [22], no amelioration in CHF was reported by increasing the concentration of functionalized silica nanoparticles [21] from 0.5 to 2.5 wt%. This was justified by the grafting silanes to the surface of silica nanoparticles.
The impact of utilizing Al O as nanofluid for improving the pool boiling features of a copper surface with the average roughness of 177.5 nm-292.8 nm was studied by Ham et al. [198]. The alumina nanoparticles with the concentrations of 0, 10 −3 , 10 −2 , 0.05, and 0.1 vol.% were added to the water (as the base fluid). It was observed that CHF was increased up to 224.8% and 138.5%, respectively. When the increase in the concentration of nanoparticles reached 0.1 vol.%, a reduction was observed in the CHF behavior. This behavior is related to the variation in surface roughness. Regarding these findings, the significant impact of surface roughness, as well as nanoparticle concentration in heat transfer properties, can be inferred. The effect of diverse nanoparticle concentrations in CHF behavior for [95,199,200] is depicted in Figure 11.  [106] compared to [196,197]. Adapted from Mori and Utaka. Reproduced from [68], Elsevier: 2017.
The impact of utilizing Al 2 O 3 as nanofluid for improving the pool boiling features of a copper surface with the average roughness of 177.5 nm-292.8 nm was studied by Ham et al. [198]. The alumina nanoparticles with the concentrations of 0, 10 −3 , 10 −2 , 0.05, and 0.1 vol.% were added to the water (as the base fluid). It was observed that CHF was increased up to 224.8% and 138.5%, respectively. When the increase in the concentration of nanoparticles reached 0.1 vol.%, a reduction was observed in the CHF behavior. This behavior is related to the variation in surface roughness. Regarding these findings, the significant impact of surface roughness, as well as nanoparticle concentration in heat transfer properties, can be inferred. The effect of diverse nanoparticle concentrations in CHF behavior for [95,199,200] is depicted in Figure 11.
The impact of utilizing Al O as nanofluid for improving the pool boiling features of a copper surface with the average roughness of 177.5 nm-292.8 nm was studied by Ham et al. [198]. The alumina nanoparticles with the concentrations of 0, 10 −3 , 10 −2 , 0.05, and 0.1 vol.% were added to the water (as the base fluid). It was observed that CHF was increased up to 224.8% and 138.5%, respectively. When the increase in the concentration of nanoparticles reached 0.1 vol.%, a reduction was observed in the CHF behavior. This behavior is related to the variation in surface roughness. Regarding these findings, the significant impact of surface roughness, as well as nanoparticle concentration in heat transfer properties, can be inferred. The effect of diverse nanoparticle concentrations in CHF behavior for [95,199,200] is depicted in Figure 11.  Sarafraz et al. [116] studied the impact of different concentrations of Zirconium oxide nanofluids in a water-ethylene glycol (50:50) mixture on heat transfer behavior. Four volumetric concentrations of 0.025%, 0.05%, 0.075% and 0.1% were considered for the mentioned nanofluid. It was observed that the CHF was increased up to 29% according to the particle deposition over the heater surface, which tincreases the roughness and wettability. Moreover, the contact angle was changed from 86.2 • to 58.2 • by increasing the volumetric concentration of nanoparticles. The ameliorated heat transfer was attributed to the capillary wicking into the porous layer, which affects the dry-spot rewetting considerably.
The heat transfer features of water-based nanofluids were scrutinized by Sulaiman et al. [202] through performing pool boiling experiments. Both material and concentration were considered in this study. The working fluid was obtained by the addition of TiO 2 , Al 2 O 3 and SiO 2 with the respective mass concentrations of 0.04, 0.4 and 1 kg/m 3 in water. It was observed that both parameters of material and concentration play important roles in increasing the CHF. In higher concentrations, almost 2.5-3 times CHF enhancement was observed compared to the CHF of water.
CHF enhancement by the utilization of TiO 2 nanoparticles was studied in many studies [94,200,[203][204][205]. In Figure 12, the TEM image of TiO 2 nanoparticle dispersed in water can be seen. For the nanoparticle concentration equal to 0.001 vol. %, the maximum CHF amelioration equal to 30% was observed. A further increase in concentration was followed by a slow decrease in CHF [206]. This trend can be justified by the fact that increasing the nanoparticle concentration suppresses the distribution of active nucleation sites, which results in lesser dry spot rewetting.
Hiswankar and Kshirsagar [199] performed a study for evaluating the effect of nanoparticle concentration on CHF behavior. It was found that, for smaller nanoparticle concentrations, CHF was improved up to 170%. This amelioration was because of the increase in the surface roughness as a result of nanoparticle deposition on the heater surface. Therefore, it can be concluded that a trade-off should be performed between different parameters, such as nanoparticle concentration, surface roughness and a reduction in nucleation sites. For this purpose, an optimization can be conducted on the size of nanoparticle. CHF enhancement by the utilization of TiO2 nanoparticles was studied in many studies [94,200,[203][204][205]. In Figure 12, the TEM image of TiO2 nanoparticle dispersed in water can be seen. For the nanoparticle concentration equal to 0.001 vol. %, the maximum CHF amelioration equal to 30% was observed. A further increase in concentration was followed by a slow decrease in CHF [206]. This trend can be justified by the fact that increasing the nanoparticle concentration suppresses the distribution of active nucleation sites, which results in lesser dry spot rewetting. Hiswankar and Kshirsagar [199] performed a study for evaluating the effect of nanoparticle concentration on CHF behavior. It was found that, for smaller nanoparticle concentrations, CHF was improved up to 170%. This amelioration was because of the increase in the surface roughness as a result of nanoparticle deposition on the heater surface. Therefore, it can be concluded that a trade-off should be performed between different parameters, such as nanoparticle concentration, surface roughness and a reduction in nucleation sites. For this purpose, an optimization can be conducted on the size of nanoparticle.
According to the results of [207], 80% CHF amelioration can be achieved by the dispersion of 0.001 wt% ZnO nanoparticle in the water. However, the CHF performance deteriorates in higher concentrations due to the deposition of a thick layer of nanoparticle deposition on the heater surface. This investigation introduced a preliminary idea for the selection of the type and concentration of nanoparticles in the working fluid.
An empirical study on the impacts of different parameters of nanofluids was performed by Ciloglu and Bolukbasi [77] for pool boiling CHF. They declared the significant effect of nanoparticle concentration on CHF behavior. In addition, they found that the maximum CHF amelioration According to the results of [207], 80% CHF amelioration can be achieved by the dispersion of 0.001 wt% ZnO nanoparticle in the water. However, the CHF performance deteriorates in higher concentrations due to the deposition of a thick layer of nanoparticle deposition on the heater surface. This investigation introduced a preliminary idea for the selection of the type and concentration of nanoparticles in the working fluid.
An empirical study on the impacts of different parameters of nanofluids was performed by Ciloglu and Bolukbasi [77] for pool boiling CHF. They declared the significant effect of nanoparticle concentration on CHF behavior. In addition, they found that the maximum CHF amelioration occurred in an optimum value of nanoparticle concentrations. Before this optimum value, increasing the concentration would lead to an increase in CHF, while, after this point, CHF remains constant by a further increase in concentration. In the investigation of Kole and Dey [124] for the water-based Cu nanofluids, it was found that the CHF increases by raising the concentration of nanoparticle in the range of 0.005-0.5 wt%.
Sarafraz et al. [116] scrutinized the ethylene glycol/water-based with 0.025-0.1 vol.% ZrO 2 nanofluids and 50 vol.% glycol, experimentally. They also confirmed the increasing trend of CHF with increasing the nanoparticle concentration. The boiling performance of H 2 O/LiBr-based binary Al 2 O 3 nanofluids was studied by Jung et al. [208]. They used polyvinyl alcohol as a stabilizer and reported the CHF amelioration by increasing Al 2 O 3 concentration from 0.01 to 0.1 vol.%. Moreover, Lotfi and Shafii [209] reported CHF reduction for 0.5-4 wt% Ag and 0.125-1.0 wt% TiO 2 nanofluids as the nanoparticle concentration decreases. The boiling properties of a honeycomb porous plate, which covers a 30-mm diameter heating surface that was placed in water-based TiO 2 nanofluids was studied by Mori et al. [210,211]. They achieved considerable CHF amelioration equal to 320 W/cm 2 by increasing the nanoparticle concentration to 0.1 vol.%.

Nanoparticle Dispersion Method
The method for the dispersion of a nanoparticle into the base fluid is of great significance as it affects the characteristic as well as the agglomeration of nanoparticles [192,212,213]. By reviewing the available studies, it was found that the existing experimental data were different even for the same nanoparticle material, which may be attributed to the different dispersion methods of the nanoparticles in the base fluids [77,137].
Nanofluids can be prepared using either two-step or one-step processes [161,214]. In the prior method, nanoparticles, nanotubes, or nanofibers are formed as dry powder through the physical or chemical treatments such as condensation of inert gas and deposition of chemical vapor in the first step. In the subsequent step, the powder is dispersed into the base liquid. Aggregation of nanoparticles in the fluid is considered as an important challenge in two-step processes. Moreover, inefficiency of the two-step process in dealing with heavier metal nanoparticles despite the fairly good performance for oxide nanoparticles can be considered as the other challenge. In the one-step process, the synthesis and dispersion of nanoparticles into the base fluid are conducted at the same time. However, even though this process provides superior dispersibility and stability in comparison with the two-step processes, the complexity of the process due to the stringent preparation requirements is regarded as the challenge. It should be noted that most of the conducted studies concerning nanofluid pool boiling were performed using two-step methods.
Nanofluid preparation is an important factor. Particle deposition and high thermal resistance can be generated as a result of improper preparation. Stability and appropriate distribution are considered as the effective parameters for preparing nanofluids. The appropriate preparation can be obtained through the use of the following processes [215]: • Changing the pH of the solution via the addition of acid to keep the nanoparticles away from their isoelectric point; • The addition of surfactants and/or dispersants; • Stabilization through the use of ultrasonic vibration or electrostatic stabilization.
Based on the performed studies, the rheological behavior of the fluid and the nucleate boiling heat transfer can be affected by the first and second processes [131,216,217]. Pool boiling heat transfer of the aqueous-based alumina nanofluids with adjustable pH solution was conducted by Wen and Ding [131], experimentally. In order to improve the CHF performance, solvable additives such as TSP, sodium lauryl sulfate (SLS), and sodium lauryl benzene sulphonate (SLBS) were utilized in most of the studies in the field of pool boiling heat transfer [218,219]. In an empirical investigation by Jeong and Kwon [220], the impact of ultrasonic vibration on CHF was examined under natural convection condition. The CHF measurements were conducted in the presence and absence of ultrasonic vibration applied to the distilled water as the working fluid. Regarding the obtained results, it was found that applying the ultrasonic wave to water results in the CHF improvement. As visual observation showed, it can be inferred that there is a close relation between the CHF amelioration and dynamic behavior of bubble generation and departure in the acoustic field.

Key Surface Factors of Deposited Nanoparticle for CHF Enhancement
The processes by which the nanoparticles are attached to the solid surfaces that are called substrates and the resultant coatings of nanoparticles are formed are referred to as nanoparticle deposition. Different structures of coatings, including monolayer/multilayer and organized/unorganized structures, would be achieved due to the used coating method. Several coating methods have been proposed for nanoparticle deposition. The difference of these methods is in their capability of controlling the density and layer thickness of particle packing, their ability to use different particles, the complexity of the method as well as the required instrumentation [221][222][223]. As mentioned in the previous sections, the use of nanofluid in pool boiling heat transfer would lead to CHF amelioration as a result of nanoparticle deposition on the heating surface. Recently, different heating surfaces were utilized in the boiling heat transfer. The empirical results of flat plates and thin wires were reviewed by Kim [224]. It was inferred that flat plates are more reliable; however, thin wires would lead to the better amelioration of CHF.
It is worth understanding the mechanism of nanoparticle deposition on the heating surface. Microlayer evaporation was considered as the reason behind nanoparticles' deposition on the heater surface by Kim et al. [102]. They declared that the formation of a microlayer at the base of the bubbles is generally because of the nanoparticles left behind the evaporating liquid. Therefore, nanoparticles are bonded to the heater surface as a result of microlayer evaporation.
In another study by Kim et al. [100], it was reported that the growth of the deposited layer was attributed to the natural convection or gravitational sedimentation. Additionally, they stated that the CHF amelioration of nanofluids was due to the deposited layer. Moreover, Bang and Chang [99] introduced the particles deposition (not increased concentration) as the responsible factor for CHF amelioration. For a better understanding of the deposited layer effect, two sets of experiments were performed by Kim and Kim [203]. The first test was conducted with nanofluid and bare heater wire, while the second one was performed with nanoparticle-deposited wire and pure water. According to the obtained empirical data, they concluded that the CHF value in the second experiment was 1.35 times greater than the corresponding value in the first experiments. These results were consistent with the result of Kwark et al. [106]. In the latter study, they confirmed the good bonding strength of the deposited layer even after 16 consecutive experiments, which emphasizes the great role of nanoparticle deposition on CHF enhancement. Moreover, it was found that not only the nanoparticles were deposited but also solvable salts may be deposited during the boiling of nanofluids, which may be responsible for the CHF amelioration [225].
Surprisingly, in the boiling of functionalized silica nanofluids, no deposited layer was formed in all concentrations and pressure ranges. This was followed by an insignificant CHF variation [22]. In a recent research by Kim et al. [100], two key parameters were stated to be influential in the physical mechanism of CHF amelioration for deposited wire. The thickness of porous structure and capillarity phenomena; the former is responsible for holding the additional liquid microlayer, while the latter supplies the liquid toward dry area underneath growing bubbles. In the following section, the surface roughness, surface wettability and capillary wicking structure that are considered to play key roles for the CHF enhancement are discussed.

Surface Roughness
As a result of increasing the concentration of nanofluids, the deposited layer grows on the heater surface. On the other hand, a direct relationship is established between the surface roughness and deposited layer, which has been confirmed in several studies [113,226]. For example, for alumina-water nanofluid, the surface roughness of 0.5 vol.% was lower than that of 0.01 vol.%. This could be due to the lesser uniformity of the deposited layer on the surface, which would lead to the higher heat transfer rates compared to the boiling pure water with nanoparticle deposited surfaces [165].
For a better understanding of the surface roughness effects, rough and smooth surfaces were studied experimentally by Wen et al. [227]. According to the obtained results, particle deposition on the heating surface would result in the increased roughness on the smooth surface. However, no significant variation was illustrated for the rough surface. Moreover, they claimed that the surface modification was performed by the inherent feature nanoparticles every time after boiling. In contrast, for rough surfaces, the arithmetic surface roughness was increased by increasing wall superheat. The inverse trend was observed for the smooth surfaces [164]. In fact, surface roughness is a dominant factor over the surface wettability. Furthermore, the boiling heat transfer of nanofluids was about 2 times that of the base fluid for the rough surface. For the smooth surface, approximately a 30% reduction was observed when boiling with nanofluid in comparison with base fluid.
Although the average surface roughness reduced from 0.167 to 0.099 µm after boiling experiments, the CHF amelioration was accomplished. This can be attributed to the changes in surface microstructure and improved topography [121]. Utilizing ionic additive (nitric acid) is an alternative way, which affects the formation of self-assembled nanoparticle structures on the heater surfaces. In this way, a more uniform and smoother structure would be formed, which reduces the CHF enhancement [228].
It can be stated that, based on the particle deposition, the surface roughness of the heating surface was not constant. Furthermore, CHF was reduced as the result of thermal insulation caused by the continuous growth of the deposited layer. In order to obtain the optimum value for the thickness of deposited layer, which corresponds on the maximum CHF, supplementary investigations are required.
It was also inferred that the boiling performance of nanofluids was affected considerably by the transient nature of nanofluid.

Surface Wettability
The surface wettability amelioration is one of the prominent parameters in the pool boiling heat transfer of nanofluids [229]. Based on Kolev's model [230], variations of contact angle significantly affect the CHF. Accordingly, Figure 13 presents a decrease in contact angle from π/2 to π/8 giving rise to a 2-fold CHF increment. In the mentioned figure, θ = π/2 was regarded as a reference on which all other cases were normalized. Not only static contact angle, but dynamic contact angle also have an effect on bubble dynamic and pool boiling performance [231]. The empirical findings of Kim et al. [102] were consistent with the impact of contact angle for all the evaluated nanofluids, including alumina, zirconia and SiO 2 . Furthermore, surface wettability and contact angle are inversely correlated [146]. The influence of the surface tension and contact angle of TSP and nanofluids has been addressed by Jeong et al. [180]. Accordingly, a decline in the contact angle managed to ameliorate the surface wettability, and hence led to CHF amelioration in both fluids. These results are in line with other investigations in the field of pool boiling [100,104,108,192].  Figure 14 illustrates the contact angle variations in distilled water and alumina-water nanofluid. The decrease in contact angle can be attributed to the nanoparticles residing on the heater surface. The wettability can also alter the nanoparticle content of the base fluid in a way that the concentration of alumina-water nanofluid experienced a 2-fold amelioration by the contact angle decline from 46.5° to 33° [105,232]. The influence of surface wettability was also addressed by Park et al. [233] in alumina-water and graphene/graphene oxide nanofluids, where no definite surface wettability-CHF amelioration relationship was observed. Moreover, they found that the variation of CHF ameliorationcorresponding wavelength was more profound than the impact of contact angle in all studied nanofluids. A recent study by Ahn et al. [234] indicated that the wettability alone fails to completely justify the CHF amelioration upon the use of graphene oxide (RGO) colloids. Additionally, they concluded the CHF amelioration can be assigned to the synergetic impacts of surface wettability and the contact angle dynamics during water absorption as well as the liquid-saturated porous environment as the consequence of water absorption onto the base graphene multilayer and the self-assembled inter-connected graphene layer.   Figure 14 illustrates the contact angle variations in distilled water and alumina-water nanofluid. The decrease in contact angle can be attributed to the nanoparticles residing on the heater surface. The wettability can also alter the nanoparticle content of the base fluid in a way that the concentration of alumina-water nanofluid experienced a 2-fold amelioration by the contact angle decline from 46.5° to 33° [105,232]. The influence of surface wettability was also addressed by Park et al. [233] in alumina-water and graphene/graphene oxide nanofluids, where no definite surface wettability-CHF amelioration relationship was observed. Moreover, they found that the variation of CHF ameliorationcorresponding wavelength was more profound than the impact of contact angle in all studied nanofluids. A recent study by Ahn et al. [234] indicated that the wettability alone fails to completely justify the CHF amelioration upon the use of graphene oxide (RGO) colloids. Additionally, they concluded the CHF amelioration can be assigned to the synergetic impacts of surface wettability The influence of surface wettability was also addressed by Park et al. [233] in alumina-water and graphene/graphene oxide nanofluids, where no definite surface wettability-CHF amelioration relationship was observed. Moreover, they found that the variation of CHF amelioration-corresponding wavelength was more profound than the impact of contact angle in all studied nanofluids. A recent study by Ahn et al. [234] indicated that the wettability alone fails to completely justify the CHF amelioration upon the use of graphene oxide (RGO) colloids. Additionally, they concluded the CHF amelioration can be assigned to the synergetic impacts of surface wettability and the contact angle dynamics during water absorption as well as the liquid-saturated porous environment as the consequence of water absorption onto the base graphene multilayer and the self-assembled inter-connected graphene layer.

Capillary Wicking
The impact of surface wettability cannot sufficiently describe the further enhancement of CHF beyond the attainable values in cases where the contact angle is close to zero (reported by Kim et al. [100]). Moreover, the liquid rise on the capillary-induced deposited nanoparticles can be determined by: in which σ, θ, R c , ρ and g are liquid vapor interface tension, contact angle, density and gravitational acceleration, respectively. The phenomenon of capillary wicking on the nano-porous layers will delay the dry-out and cool the dried spots giving rise to an ameliorated CHF [204]. The retarded CHF is not only attributed to the ameliorated surface wettability but also the capillary impacts originated from the nanoparticle deposition. Consideration of capillarity cannot completely elucidate the CHF enhancement for grapheme/grapheme-oxide nanofluids [233] as other factors, such as wavelength modulation during the nanofluid boiling are also involved. The influence of water absorption phenomenon was addressed by Ahn et al. for the saturated porous layers during the RGO nucleate boiling [235]. The substantial CHF amelioration (320%) cannot be fully justified considering the capillary wicking solely. In a recent investigation by Kim et al. [100], the capillary rise of all nanoparticle deposited wires was explored, where the CHF value exhibited a strong dependence on the water rise on the fouled heater wires, although all the deposited wires possessed the same contact angles. Capillarity trend showed an initial increment by the thickness of the deposited nanoparticle layer, which can be attributed to the larger number of micro-scaled structures. Further augmentation of the nanoparticles layer thickness resulted in no capillarity elevation since the valleys were already filled by the nanoparticles. According to the available literature, the capillary wicking during the boiling of nanofluids can significantly affect the CHF enhancement. The capillary rise alone, however, does not suffice to completely elucidate the underlying mechanism of CHF improvement in the majority of the studies. The general census of the phase change heat transfer community recognizes the nanoparticle deposition as the most influential factor in the CHF amelioration. The impact of nanoparticle-induced surface modifications should be deeper studied by investigating the combined influences of surface roughness and wettability (contact angle) along with the capillary wicking. These results can be helpful in clarifying the CHF enhancement mechanism for a wide range of nanofluid contents with various particle sizes.

Mechanisms of CHF Enhancement
Possible mechanisms of CHF amelioration upon nanoparticle deposition on the surfaces are discussed in this section. Despite the many types of research during the last century on CHF phenomena, the establishment of a robust general correlation between empirical results and the physical process underlying the CHF amelioration has remained elusive. As discussed in the previous section, nanoparticle deposition can change the surface properties and it has an effect on CHF. Five different mechanisms have been introduced to connect the surface properties on CHF prediction: hydrodynamic instability [236], macrolayer dry-out [237], bubble interference [238], hot/dry spot [239] Energies 2020, 13, 4026 25 of 49 and interfacial lift-off [240], among which, the theory of hydrodynamic instability by Zuber [236] has gained the highest popularity as the first theoretical and mechanistic modeling for CHF. Most of subsequent pool boiling CHF studies tried to enhance the predictive ability of Zuber's model by considering the parametric impacts neglected in the original model.

Hydrodynamic Instability Model
Based on Zuber [241], the hydrodynamic instability model is based on the assumption that the latency in CHF is highly dependent on the hydrodynamics of the vapor and liquid countercurrent flows at a relatively large distance from the top of the heated surface. This theory expresses that an ameliorated superficial vapor velocity will result in an upward movement of the vapor, which will inhibit the downward flow of liquid to the heater surface for increased heat fluxes. Such a condition will lead to an instability, which finally will give rise to CHF. The superficial vapor velocity ( j v ) causing the instability above the flat upward-facing heated surface was determined by Zuber as: where the subscripts of l and v are respective expressions for liquid and vapor, while σ, ρ and g denote liquid vapor interface tension, density and gravitational acceleration, respectively. Thus, the CHF of a saturated boiling can be presented by: where h f g represents the latent heat of evaporation. Many studies have confirmed the validity of the hydrodynamic instability model to elaborate the CHF amelioration mechanism in a nanofluid. For example, the study by Kim et al. [102] can be mentioned, which used the hydrodynamic model as an interpretative tool for CHF amelioration in nanofluids since its stability could consider the surface impacts including wettability and roughness. Moreover, an earlier study by You et al. [97] declared that some key but unknown factors missing in Zuber's theory were responsible for CHF amelioration in nanofluids. For CHF enhancement interpretation in nanofluids, Zuber's theory was further modified by Lienhard and Dhir [242], who included the size and geometrical impacts. They also theorized that the variations in the surface contact angle can alter the size and spacing of the vapor jets located above the heater surface. The hydrodynamic theory is recently widely questioned, even in the case of pure fluids [243].

Macrolayer Dry-Out Model
Developed by Haramura and Katto, the macrolayer dry-out model [237] considers large mushroom-shaped bubbles above the heated surface for a sufficiently long time prior to their departure. Based on this model, CHF occurs at an adequately high heat flux for the macrolayer evaporation before the mushroom bubble departure. Kim et al. [102] showed that the macro-layer dry-out theory indicates a strong correlation between CHF amelioration and the surface wettability augmentation due to nanoparticle deposition. Calculating the macrolayer thickness by the Sadasivan model [244], they detected a decline in contact angle from 70 • to 20 • , which could be attributed to nanoparticle deposition giving rise to a four-fold elevation of the liquid layer thickness and hence a 400% increment in CHF. The dent-like structures of the porous layers of TiO 2 on wires were also addressed [100]. This layer is capable of holding liquid and increasing the efficient liquid layer thickness in the macrolayer model for CHF. A limit was, however, indicated in CHF regardless of the liquid rise and/or TiO 2 -deposited wires' thickness; it requires deeper investigations to elucidate the underlying mechanism.

Bubble Interaction Model
Rohsenow and Griffith [238] suggested bubble interaction theory for the heated surface, which postulates the occurrence of CHF, particularly at high wall superheats. The number of bubbles and the departure frequency become very high, meaning that the bubbles coalesce radially. This prevents the liquid flow towards the surface. The impact of shear force created by the mutual interactions of the growing and departing bubbles is contained in most of the developed models. The bubble growth cycle was cut prematurely by shear force, resulting in a decrease in bubble departure diameter, leading to lower latent heat removal upon departure. Based on Kolev model [230], the heat flux in the nucleate boiling region and the wall superheat can be stated as: where subscripts of w, sat and d are respective expressions for wall, saturation and departure, while n and T stand for active nucleation site density and temperature, respectively. Furthermore, the nucleation site density will reduce by contact angle, as observed by Wang and Dhir [245]. Therefore, the shear stress intensity due to the mutual interaction of bubbles will have a slower growth on a low-contact angle heated surface compared to the one having a larger contact angle. According to Kim et al. [102], the bubble interaction theory accounts for surface wettability amelioration as a probable mechanism of CHF amelioration.

Yagov Model
The hot/dry spot procedure associates CHF to abundant small dry spots located on the heater surface during nucleate boiling. The mentioned model was first developed by Yagov [246] in 1988 and then amended by Yagov [239] in 2014. Some studies tried to extrapolate CHF from nucleate boiling models [247][248][249][250]. According to Van Ouwerkerk [251], the initiation mechanism of CHF based on the hot/dry spot model can be regarded as an irreversible growth of the dry spot area on a surface (see Figure 15).
where subscripts of , and are respective expressions for wall, saturation and departure, while and stand for active nucleation site density and temperature, respectively. Furthermore, the nucleation site density will reduce by contact angle, as observed by Wang and Dhir [245]. Therefore, the shear stress intensity due to the mutual interaction of bubbles will have a slower growth on a low-contact angle heated surface compared to the one having a larger contact angle. According to Kim et al. [102], the bubble interaction theory accounts for surface wettability amelioration as a probable mechanism of CHF amelioration.

Yagov Model
The hot/dry spot procedure associates CHF to abundant small dry spots located on the heater surface during nucleate boiling. The mentioned model was first developed by Yagov [246] in 1988 and then amended by Yagov [239] in 2014. Some studies tried to extrapolate CHF from nucleate boiling models [247][248][249][250]. According to Van Ouwerkerk [251], the initiation mechanism of CHF based on the hot/dry spot model can be regarded as an irreversible growth of the dry spot area on a surface (see Figure 15).
where the subscripts of g, f and sat are respective expressions for gas, fluid and saturation, while h, σ, ρ, k, g, v, C P and T stand for latent heat, liquid vapor interface tension, density, thermal conductivity, gravitational acceleration, kinematic viscosity, specific heat capacity and temperature, respectively. Moreover, R i shows the individual gas constant, and the Prandtl number function f Pr f can be determined by: where the subscript of f stands for fluid and Pr denotes Prandtl number. For high reduced pressures P P c > 0.03 , the following equation can be stated: where the subscripts of f and g are respective expressions for fluid and gas, while h, ρ, σ, g and µ stand for latent heat, density, liquid vapor interface tension, gravitational acceleration and dynamic viscosity, respectively. The following weighting function was suggested by Yagov [239] to obtain the CHF for intermediately diminished pressures, 0.001 ≤ P P c ≤ 0.03 : where the subscripts of hp and lp denote high pressure and low pressure, respectively. Notably, in contrary to the hydrodynamic instability model, the dry spot model is dependent on liquid viscosity. The impact of viscosity will be pronounced at lower pressures (Equations (11) and (12)) while getting much weaker at high pressures (Equation (13)). Yagov [239] concluded that Equation (13) predicts CHF values closer to the results of hydrodynamic instability models since the majority of experiments were conducted at moderate to high reduced pressures in which the influence of viscosity was rather weak. Kam et al. [253] developed the observation-based model by dividing the dry patch into two kinds of reversible (quenchable) or irreversible (unquenchable).

Theofanous and Dinh Formulation
With the help of a high-speed, high resolution infrared imaging, Theofanous et al. [240,254] reported a temperature pattern on the surface during nucleate boiling. Lower temperatures were observed at cold spots (known as active nucleation sites), which might be due to effective heat removal by the bubble growth. By elevation of the heat flux, the cold spot pattern got more regular, despite an augmentation in the pool's two-phase motion chaos. Moreover, enhancement of the heat flux will result in the development of the hot spots inside the originally cold spots; at this time, CHF will be initiated by irreversible growth of hot spots, as illustrated in Figure 16. Moreover, Jung et al. [255] expressed that CHF occurrence is upon the faster growth of dry spot size compared to the heat transfer increment through the wetted area; thus, higher CHF is achievable by either the enhancement of the wetted fraction of the surface or improvement of heat transfer within the wetted parts. Presuming that CHF is governed by the liquid microlayer dynamics and instability, a followup investigation by Theofanous and Dinh [256] resulted in a CHF relation similar to that of Kutateladze-Zuber: In the above relation, k denotes a function of surface conditions that descend by an increase in the surface wettability. No approach was, however, proposed to calculate this parameter. An analytical equation was then developed for by Kim et al. [102] based on Lord Rayleigh's formula for a volume of static liquid meniscus [257]: The above relation can be applied for hydrophilic surfaces, whose contact angles are less than 90°.

Other Observations
On the basis of the total reflection technique, Chu et al. [258] lately proved that, in contrast to the general belief in the presence of a thin, stable liquid film beneath large coalescent hovering bubbles, these bubbles reside on an almost dry surface even prior to CHF onset. The hydraulic and thermal criteria were developed by Choi et al. [259] for the irreversible dry spot initiation. Based on the hydraulic criterion, the irreversible condition occurs in cases where the number of dry spots grows to a critical level, at which point the spots are in contact with each other. According to the thermal criterion, the irreversible condition initiates when the temperature of dry spot environment reaches the Leidenfrost temperature preventing the surface rewetting even after the bubble departure. The CHF is achievable upon the fulfillment of both criteria. Kim et al. [260] observed that Presuming that CHF is governed by the liquid microlayer dynamics and instability, a follow-up investigation by Theofanous and Dinh [256] resulted in a CHF relation similar to that of Kutateladze-Zuber: In the above relation, k denotes a function of surface conditions that descend by an increase in the surface wettability. No approach was, however, proposed to calculate this parameter. An analytical equation was then developed for k by Kim et al. [102] based on Lord Rayleigh's formula for a volume of static liquid meniscus [257]: The above relation can be applied for hydrophilic surfaces, whose contact angles are less than 90 • .

Other Observations
On the basis of the total reflection technique, Chu et al. [258] lately proved that, in contrast to the general belief in the presence of a thin, stable liquid film beneath large coalescent hovering bubbles, these bubbles reside on an almost dry surface even prior to CHF onset. The hydraulic and thermal criteria were developed by Choi et al. [259] for the irreversible dry spot initiation. Based on the hydraulic criterion, the irreversible condition occurs in cases where the number of dry spots grows to a critical level, at which point the spots are in contact with each other. According to the thermal criterion, the irreversible condition initiates when the temperature of dry spot environment reaches the Leidenfrost temperature preventing the surface rewetting even after the bubble departure. The CHF is achievable upon the fulfillment of both criteria. Kim et al. [260] observed that the irreversible dry spots will be formed in water at 134 • C, which is far lower than the Leidenfrost and maximum liquid-contact temperatures. Kim et al. [260] and Chu et al. [261] agreed that an irreversible dry spot will be formed as the consequence of vigorous bubble nucleation close to the triple phase contact line, which will prevent from the dry spot wetting.
Advances to stronger flow visualization and measurement methods will drastically improve readers' insight into interfacial mechanisms near the dry spot. This will further enhance the dry-spot-based CHF model's prediction capabilities.

Vertical and Near-Vertical Orientations
Galloway and Mudawar [262] conducted a high-speed video motion investigation and hypothesized that CHF within flow boiling occurs upon sufficient growth of vapor momentum in a way that it can lift the bulk liquid away from the surface. They applied this mechanism to develop a phenomenological interfacial lift-off model [263]. Subsequently, Mudawar et al. [264,265] reported the similarity of the CHF interfacial behavior in pool boiling along the vertical surface to the flow boiling CHF. Based on pool boiling illustrations, a series of incidents will cause the CHF. In cases where heat flux is slightly less than CHF, intensive vaporization may lead to the formation of a vapor layer propagating along the surface. Considerable waves may be formed in this layer as the result of Helmholtz instability which allow liquid contact with the surface only in the 'wetting fronts' corresponding to the wave troughs. The mentioned wetting fronts sweep along the surface and provide the last opportunity for the surface cool down. CHF can be regarded as a consequence of wetting front loss in cases where the severe surface-perpendicular vapor momentum rises beyond the opposing interfacial curvature-induced pressure force.
where q l denotes the localized heat flux within the wetting front. Additionally, subscripts of g, f and sub are respective expressions for gas, fluid and subcooled, while ρ, h, C P and T stand for density, latent heat, specific heat capacity and temperature, respectively. The average interfacial curvature-induced pressure difference can be determined by: where δ represents the mean vapor layer thickness, whereas σ stands for liquid vapor interface tension; γ c is the critical wavelength corresponding to the initiation of the Helmholtz instability. By the substitution of Equation (18) in Equation (17) and considering the wetting front span as 1 4 of the critical wavelength (i.e., q CHF = q l /4), one can present: Different flow models were used to derive the analytical equations of δ and γ c : where In the above equation, f i denotes the interfacial friction factor (considered as 0.5). By the substitution of Equations (20) and (21) in Equation (19), the following equation will be obtained: Under saturated condition and for noncritical pressure values, the mentioned equation could be further simplified as: Equation (23) is similar to the Kutateladze-Zuber equation. Notably, the gravity term in Equation (23), however, originates from the buoyancy-driven vapor flow rather than the Taylor instability in Zuber's model. Subsequent studies by Howard and Mudawar [265] showed that the CHF interfacial lift-off trigger process is not restricted to the vertical orientation, as it encompasses all the near-vertical orientations varying in θ = 60-165 • (see Figures 17b and 18). The mentioned model can effectively predict the near-vertical-downward-facing transition angle of 165 • , as illustrated in Figure 19. By the application of various reported data, they also indicated that, at 0 < θ < 90 • , the impact of orientation on CHF is very weak, implying that Equation (23) is equally valid for h = 0 • and 90 • . Recent studies by Kim et al. [266] and Zhong et al. [267] also confirmed the interfacial lift-off model in pool boiling CHF with various orientations.
The interfacial lift-off model is highly versatile and precise in the prediction of flow boiling CHF at diverse conditions, including microgravity [268,269]. In particular, it confirms the gravity dependence of CHF. Overall, in comparison with pool boiling, it is more successful in the prediction of decreased gravity CHF for flow boiling [270].
In the above equation, denotes the interfacial friction factor (considered as 0.5). By the substitution of Equations (20) and (21) in Equation (19), the following equation will be obtained: Under saturated condition and for noncritical pressure values, the mentioned equation could be further simplified as: Equation (23) is similar to the Kutateladze-Zuber equation. Notably, the gravity term in Equation (23), however, originates from the buoyancy-driven vapor flow rather than the Taylor instability in Zuber's model. Subsequent studies by Howard and Mudawar [265] showed that the CHF interfacial lift-off trigger process is not restricted to the vertical orientation, as it encompasses all the near-vertical orientations varying in = 60 − 165° (see Figure 17b and Figure 18). The mentioned model can effectively predict the near-vertical-downward-facing transition angle of 165°, as illustrated in Figure 19. By the application of various reported data, they also indicated that, at 0 < < 90°, the impact of orientation on CHF is very weak, implying that Equation (23) is equally valid for ℎ = 0° and 90°. Recent studies by Kim et al. [266] and Zhong et al. [267] also confirmed the interfacial lift-off model in pool boiling CHF with various orientations.
The interfacial lift-off model is highly versatile and precise in the prediction of flow boiling CHF at diverse conditions, including microgravity [268,269]. In particular, it confirms the gravity dependence of CHF. Overall, in comparison with pool boiling, it is more successful in the prediction of decreased gravity CHF for flow boiling [270].

Horizontal, Upward-Facing Orientation
The interfacial lift-off process in pool boiling CHF was also addressed along with the upward-facing orientation by Guan et al. [271]. The lift-off process for pentane is illustrated in Figure 20 through a series of images. At t = 0 s, a small vapor patch can be seen at the right edge of the surface which propagated from right to left until t = 2.8 s when the whole surface was covered by vapor. Accordingly, they proposed a CHF model for pool boiling on the basis of a liquid macrolayer lift-off process, which relies on the Galloway's and Mudawar's interfacial lift-off model [263]. The major difference is that the Guan et al. [271] model ( Figure 21) does not require a separate flow model for the determination of the vapor layer parameters. It rather employs the momentum conservation and Laplace condition on the vapor-liquid interface to obtain an equation for CHF-corresponding upward vapor velocity.
In the above equation, r represents the local interface curvature radius. Maximum vapor velocity can be achieved at minimal local radius: where δ * shows the liquid macrolayer thickness predictable by the theoretical formulation proposed by Rajvanshi et al. [272]. Using this approach, the following equation can be obtained for pool boiling CHF: The above equation is relatively similar to the Kutateladze-Zuber equation. The interfacial lift-off process in pool boiling CHF was also addressed along with the upwardfacing orientation by Guan et al. [271]. The lift-off process for pentane is illustrated in Figure 20 through a series of images. At = 0 s, a small vapor patch can be seen at the right edge of the surface which propagated from right to left until = 2.8 s when the whole surface was covered by vapor. Accordingly, they proposed a CHF model for pool boiling on the basis of a liquid macrolayer lift-off process, which relies on the Galloway's and Mudawar's interfacial lift-off model [263]. The major difference is that the Guan et al.   Altogether, the interfacial lift-off mechanism for pool boiling CHF relies on microscopic observation with phenomenological appeals that may not be found in the experimental CHF correlations. Moreover, owing to its dependence on hydrodynamic instability theory, it will result in a dimensionless formulation similar to the Kutateladze-Zuber equation. Altogether, the interfacial lift-off mechanism for pool boiling CHF relies on microscopic observation with phenomenological appeals that may not be found in the experimental CHF correlations. Moreover, owing to its dependence on hydrodynamic instability theory, it will result in a dimensionless formulation similar to the Kutateladze-Zuber equation.

Challenges and Suggestions
Based on the review of available studies, enhanced CHF not only increases heat transfer efficiency but also augments the safety of heat management. A wide diversity of pool boiling experiments have been carried out at various concentrations, preparation methods and nanofluids stabilities to significantly enhance the CHF. Furthermore, different surface properties, including porosity, capillarity, wettability and roughness, can increase CHF through a complicated mechanism. CHF can be also elevated through modifying flow channel by either altering its size or incorporating surface modifiers. The application of hybrid or hierarchical structures can substantially improve the CHF for some particular purposes, among which the hydrophilic membrane with micro-channels and a combination of high porosity with the super-hydrophilic nature of the surface can be mentioned. Moreover, CHF is under the substantial influence of surface wettability and capillary wicking. The synergetic impact of the mentioned parameters should be extensively examined for diverse nanoparticle sizes and concentrations to elucidate the mechanism of CHF amelioration. Despite a wide diversity in boiling heat transfer theories to explain the variations in nanofluids' properties and surface wettability, no comprehensive theory has been introduced to mechanistically justify the CHF amelioration for a wide spectrum of nanoparticle size and contents. A lack of such an inclusive framework has limited the nanofluids commercialization in various fields including nuclear reactor, fossil fuel boiler, spray cooling and high-power electronics. Resolving the key issues can pave the way for experiments on thermal engineering systems. Researchers have recognized the nanoparticle deposition-induced enhancement in surface wettability and capillarity as the main cause of CHF augmentation. Table 2 summarizes the impacts on nanofluids' CHF.

Challenges and Suggestions
Based on the review of available studies, enhanced CHF not only increases heat transfer efficiency but also augments the safety of heat management. A wide diversity of pool boiling experiments have been carried out at various concentrations, preparation methods and nanofluids stabilities to significantly enhance the CHF. Furthermore, different surface properties, including porosity, capillarity, wettability and roughness, can increase CHF through a complicated mechanism. CHF can be also elevated through modifying flow channel by either altering its size or incorporating surface modifiers. The application of hybrid or hierarchical structures can substantially improve the CHF for some particular purposes, among which the hydrophilic membrane with micro-channels and a combination of high porosity with the super-hydrophilic nature of the surface can be mentioned. Moreover, CHF is under the substantial influence of surface wettability and capillary wicking. The synergetic impact of the mentioned parameters should be extensively examined for diverse nanoparticle sizes and concentrations to elucidate the mechanism of CHF amelioration. Despite a wide diversity in boiling heat transfer theories to explain the variations in nanofluids' properties and surface wettability, no comprehensive theory has been introduced to mechanistically justify the CHF amelioration for a wide spectrum of nanoparticle size and contents. A lack of such an inclusive framework has limited the nanofluids commercialization in various fields including nuclear reactor, fossil fuel boiler, spray cooling and high-power electronics. Resolving the key issues can pave the way for experiments on thermal engineering systems. Researchers have recognized the nanoparticle deposition-induced enhancement in surface wettability and capillarity as the main cause of CHF augmentation. Table 2 summarizes the impacts on nanofluids' CHF. Table 2. Summary of the parameter variation impacts on CHF.

Parameter Variations CHF Reason
Increase in the surface roughness Increase The CHF causes a higher surface superheat. The efficient contact area increases.
Increase in the surface wettability Increase The liquid traps in the porous structure and the vapour layer formation is prevented. The number of active nucleation sites decreases because of the low contact angle and the vapour layer formation is prevented.
Increase in the capillarity Increase The dry regions under the vapour bubbles are wetted by fresh liquid. So far, several questions have remained unanswered. Thus, highly creative and practical investigations are required to establish the advanced approaches for CHF amelioration. The following recommendations can be considered for further understanding the fundamental mechanisms of CHF amelioration: (1) Despite the relative maturity of nanofluid-related studies, issues such as long-period stability, probable erosion and required maintenance measures as well as cleaning procedures have remained obstacles to nanofluid commercialization. In this context, studies to extend the available predictive correlations or numerical simulation tools for CHF prediction are highly encouraged. Furthermore, the next generation of nanofluids, including shell-core nanofluids, microfluidics and hybrids of more than two nanofluids, are interesting topics owing to their phase-change characteristics and optimized thermo-physical properties, which may result in higher heat transfer efficiency and outstanding stability. (2) Regarding the depositing tendency of nanofluids during boiling, texture fabrications and coating particles on the surface have been intensively explored to improve CHF. The probable micro/nanostructure detachment or failure during the boiling process needs to be further addressed experimentally. Moreover, in the case of applying nanotubes or nanowires as the coating substance, the choking issue must be avoided by the systematic optimization of their diameter and length. Besides, the complexity of the surface properties' relationships should be addressed to further clarify the fundamental mechanisms of CHF amelioration. (3) In spite of huge improvements to surface modifications for flow boiling, a considerable gap still remains between the studies and engineering applications, particularly considering the size of system. Large engineering scales with different shapes were explored by micro-channels, making the coating process an obstacle in practical engineering scales. (4) The literature review showed that the integration of hybrid and hierarchical structures can significantly enhance the CHF. It is, however, accompanied by difficulties for separately investigating the impact of surface features on CHF amelioration. Putting the safe CHF enhancement as the major objective, the integration of advanced approaches, such as acoustic methods, magnetic fields and other creative techniques through novel preparation methods and materials are highly recommended. (5) Nanofluid stability has remained one of the major challenges. Instead of experimental studies with varying mass fractions, various approaches have been tested (i.e., pH control, surfactants and surface functionalization) to prevent from nanoparticles' agglomeration and sedimentation, which are, albeit, at a research level. In this regard, the industrial application of nanofluids will be realized only when their long-term physical and chemical stabilities are guaranteed at the mass production level. (6) Nanoparticle size is another important factor in the heat transfer field. The use of smaller nanoparticles is recommended. Thus, cost-effective synthesis methods should be established to prepare relatively small nanoparticles. (7) The relative size effect of nanoparticles in the liquid phase requires further examination to prevent from particle clustering. (8) To the best of our knowledge, a limited number of studies have attributed the CHF amelioration to alternations in nanofluids' thermal transport properties. Therefore, a database including the thermal transport characteristics along with the detailed specification of nanoparticle sizes and dispersion stability with/without surfactant is recommended, in which the promising nanofluids are prioritized. (9) Regarding the key role of the deposited layer on the heating surface in CHF amelioration, the thickness should be optimized to induce the maximum latency in the CHF occurrence. Moreover, the stability of the deposited layer should be experimentally checked in several replications.  (10) The synergetic impact of surface wettability and capillary wicking should be addressed to mechanistically elucidate the CHF enhancement for various particle sizes and concentrations. (11) Application of nano-coats on the heating surface by the physical/chemical vapor deposition method is a promising approach compared to nanofluids. This technique, however, demands extensive investigations regarding its stability and optimal thickness to examine the delay in the CHF occurrence. (12) The bubble dynamics should be empirically and numerically examined to determine the exact share of deposited layer and nanoparticle in the CHF amelioration. (13) Regarding the significance of pressure, the irreversible growth of dry patches should be explored at various pressure levels to elucidate the CHF enhancement in nanoparticles deposited layer.

Concluding Remarks
The present study is an overview on the CHF amelioration of pool boiling surfaces through the use of metal oxide nanoparticles deposition. The impact of metal oxide nanoparticles and their thermo-physical properties, concentration, shape and size were addressed and the inconsistencies or conflicts in the available literature were highlighted. In addition, nanoparticle deposition techniques were reviewed to find an efficient substitution for deposition rather than nanofluid boiling. Moreover, possible mechanisms and models were discussed to explore the enhancement findings. This review study had the following observations: (1) Besides the innovative dimensional analysis-based Kutateladze's CHF formulation, five different CHF mechanisms have been widely employed: bubble interference, hydrodynamic instability, macrolayer dry-out, hot/dry spot, and interfacial lift-off, among which the Zuber's hydrodynamic instability theory has gained the highest popularity owing to its mechanistic formulation and theoretical attractiveness. Lately, the theoretical-based interfacial lift-off mechanism has been widely confirmed by experimental results that can deal with various surface orientations. (2) The impacts of thermo-physical properties, concentration, shape and size on CHF have been extensively addressed in heat transfer field. Accordingly, CHF rises with wall thickness enhancement but reaches an asymptotic value beyond a thickness threshold. Thus, the data corresponding to wall thicknesses beyond this threshold are highly essential. Altogether, a limited amount of data covers the entire relevant parameters of the fluids with various thermal features. A severe shortage is also felt regarding the horizontal, downward-facing surface orientations reflecting the need for more advanced strategies to plan the future studies including the micro-photographic analysis of near-wall interfacial phenomena to confirm or reject the introduced CHF mechanism. (3) Considering the complexity of the CHF phenomenon and its dependence on a diverse range of factors, the available data should be combined in a comprehensive database to assess diverse models and correlations. Such database will be expanded by the inclusion of new data to fill the vacancies regarding the relevant parameters. (4) Nanofluids can undoubtedly improve the CHF, which can be attributed to the enhanced surface wettability upon nanoparticles deposition. The data regarding the impact of nanofluids on nucleate boiling HTC are, however, contradictory, which may be due to involvement of numerous complex factors, such as the type of liquid, initial surface roughness, and heat flux in addition to nanoparticle type, size, concentration, preparation and functionalization procedure, which can substantially affect the thermo-physical characteristics of the nanofluid and surface characteristics (i.e., surface finish, active nucleation site density, wettability and changes in the triple line). Such complexities can definitely restrict the theoretical attempts to model the nanofluid boiling. (5) Bath quenching the metal portion of a nanofluid is related to the cooling rates similar to, or even weaker than, those of the base liquid. The cooling rate can be accelerated by successive quenching as it will thicken the nanoparticle layer on the surface. At elevated surface temperatures, this layer will further destabilize the vapor film, leading to the premature disruption of the film boiling pattern. (6) The nanofluid-induced nucleate boiling enhancement can be attributed to nanoparticle deposition on the surface, resulting in capillary wicking in the porous layer and enhancement in surface wettability and bubble dynamics. These impacts are competing with each other depending on the size of nanoparticles relative to the surface roughness. (7) Although the nanofluids have shown high potentials in improving the boiling performance, some practical issues should be closely considered prior to the utilization of nanofluids in practical cooling purposes, among which the clustering, sedimentation and precipitation of nanoparticles, the clogging of intricate features, the heating surface erosion of the temporal variation of cooling performance, and the cost of quality can be mentioned. (8) Despite the improving impacts of nanofluids on the thermal conductivity of the boiling fluid, the majority of their usefulness is rooted in their surface modification. Other approaches of surface modifications (such as micro/nano studs, nanotube/nanowire arrays, microporous structure, and nanoparticle pre-deposition) are capable of offering comparable or even better heat transfer enhancements with no practical problems relevant to use of nanofluid boiling.