1 s2.0 S030645492100178X Main
1 s2.0 S030645492100178X Main
1 s2.0 S030645492100178X Main
Review
a r t i c l e i n f o a b s t r a c t
Article history: All over the world, the nuclear reactor technology has been in the critical stage of transformation and
Received 24 January 2021 upgrading. Also, there have been the higher requirements for the key technical indicators, such as the
Received in revised form 26 February 2021 power level, the miniaturization and the security. According to analyzing the engineering requirements
Accepted 25 March 2021
of typical nuclear reactor devices, the authors demonstrated the scientific problems and the physical
Available online 19 April 2021
mechanisms behind the limitation of thermal-hydraulic parameters, such as the fuel assembly, the steam
generator, and the special safe system. Thus, the technical route of heat transfer enhancement was put
Keywords:
forward in this paper. The technology of heat transfer enhancement referred to the application of various
Nuclear reactor
Interdisciplinary
engineering means to improve the heat transfer performance of component and reduce the power con-
Heat transfer enhancement technology sumption of coolant, the component temperature, as well as the device size. In this paper, the heat trans-
Passive technology fer enhancement technology was divided into the passive technology and the active technology.
Active technology Meanwhile, the passive technology could be subdivided into the structure innovation technology and
the surface modification technology. The active technology could be subdivided into the magnetic field
technology and the electric filed technology. Based on the analyses of this paper, it was considered that
the effects of different heat transfer enhancement technologies on the single-phase and two-phase con-
ditions were quite different. Moreover, the corresponding technical development maturity and the man-
ufacturing process were also much different. First of all, the structure innovation technology mainly
included two categories: the micro-channel technology and the longitudinal vortex technology. By
increasing the specific surface area of thermal component, enhancing the turbulent mixing degree of flow
field, and disturbing the near-wall thermal boundary layer, the heat transfer performance of single-phase
and two-phase conditions could be significantly improved. The above technology had been widely
applied to many heat transfer fields, and the manufacturing technology was mature as well as diversified,
which had the strong application prospect in the field of reactor core and heat exchanger without nuclear
conditions. It should be noted that the technical difficulty of structure innovation technology lied in the
optimization of structure design to obtain the best enhanced heat transfer performance. As for the surface
modification technology, it mainly included the surface micro-nano structure technology and the surface
coating technology. Moreover, it was suggested that the structural dimensions of above two technologies
were far smaller than that of the near-wall thermal boundary layer, and the corresponding effects on the
sing-phase heat transfer condition was not obvious. However, the above two technologies could signifi-
cantly affect the bubble dynamics and the two-phase interface evolution under boiling condition, which
could obviously improve the heat transfer performance under two-phase condition. In addition, the mag-
netic field technology and the electric field technology belonged to the active heat transfer enhancement
technology, which relied on the coolant physical properties and the intervention technology of flow field.
It should be noted that the above two technologies could effectively improve the heat transfer perfor-
mance under the single-phase and two-phase conditions. However, this kind of technology was still in
the stage of scientific research, and there was no practical engineering application. Based on the current
technology development, the application of these two technologies was weak in the reactor core, and the
application in non-nuclear environment was also limited.
Ó 2021 Elsevier Ltd. All rights reserved.
⇑ Corresponding authors at: Science and Technology on Reactor System Design Technology Laboratory, Chengdu 610041, China.
E-mail addresses: [email protected] (Q. Lu), [email protected] (J. Deng).
https://doi.org/10.1016/j.anucene.2021.108302
0306-4549/Ó 2021 Elsevier Ltd. All rights reserved.
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Nomenclature
Re Reynolds number
General symbol
g Gravitational acceleration (m/s2) Greek letters
h Heat transfer coefficient (W/m2) e Dielectric constant (F/m)
k Thermal conductivity (W/mK) q Density (kg/m3)
qe Heat flux (W/m2) r Surface tension (N/m)
De Equivalent diameter (m)
E Electric field intensity (V/m) Subscripts
Fe Electric field force (N) l Liquid-phase
Nconf Confinement coefficient g Gas-phase
Nu Nusselt number
Bond Bond number
Contents
1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.1. Fuel assembly . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2. Steam generation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3. Special safe system . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2. Theory of heat transfer enhancement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.1. Field synergy theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.2. Configuration theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.3. Entropy minimization theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.4. Entransy dissipation theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.5. Process enhancement theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3. Passive heat transfer enhancement technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.1. Structural innovation technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.1.1. Micro-channel technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.1.2. Longitudinal vortex technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.2. Surface modification technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.2.1. Micro-nano structure technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.2.2. Surface coating technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
4. Actively heat transfer enhancement technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
4.1. Magnetic field heat transfer enhancement technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
4.2. Electric field enhanced heat transfer technology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
5. Comprehensive evaluation. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
6. Conclusions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
Declaration of Competing Interest . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
Acknowledgements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
1. Introduction reactor system. Moreover, the fuel assembly plays as the key heat
release and exchange role, and its importance is self-evident. As is
With the continuous development of reactor technology around know to all, there are three key technical indicators for the fuel
the world, the higher requirements have been put forward for the assembly.
key thermal-hydraulic indicators, such as the power level, the
miniaturization and the security. Therefore, the need to improve (1) Power-flow ratio
the heat transfer performance of various components becomes
stronger. In recent decades, the heat transfer enhancement tech- The power-flow ratio is one of the important indicators to eval-
nology has been paid much attention in the nuclear reactor field. uate the thermal-hydraulic performance of reactor system, which
It should be noted that the heat transfer enhancement technology directly determines the average temperature of core coolant and
refers to one which can increase the heat transfer coefficient under the quality of steam generator. Also, it is closely related to the com-
a certain temperature difference and heat transfer area. Mean- pactness and the miniaturization of reactor system. One of the
while, it aims to improve the thermal efficiency, reduce the power most direct ways to increase the core power-flow ratio is to
consumption, limit the component temperature, and reduce the increase the contacting area between the fuel element and the
device size. coolant. Therefore, how to further increase the specific surface area
(surface area/volume) of fuel assembly is much important. In addi-
1.1. Fuel assembly tion, the reactor outlet is generally in a sub-cooled state under the
steady-state operating conditions. Thus, the mechanism of single-
The thermal-hydraulic performance of reactor core is directly phase convection heat transfer is generally dominant. Also, it is
related to the power level, the safety, and the miniaturization of believed that the heat transfer coefficient under the single-phase
2
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
condition can be improved by increasing the continuous distur- two categories, one is DO (Dry Out) type, and the other is DNB
bance of coolant thermal boundary layer near the heating surface, (Departure from Nucleate Boiling) type. In this case, the DO type
so as to achieve the purpose of improving the power-flow ratio. boiling crisis generally occurs under the large quality condition,
and the corresponding flow pattern is generally the annular flow.
(2) Limit quality Also, the continuous liquid film near the heating surface is contin-
uously evaporated or transfers to liquid droplets with vapor phase
The two-phase flow instability under boiling conditions has an
disturbance. Finally, the liquid film is completely evaporated at a
important impact on the stability and safety of reactor system.
specific position of heating surface. DNB type boiling crisis gener-
Meanwhile, it has always been an important safety criterion in
ally occurs under the low quality conditions. The generation rate
the thermal-hydraulic design process. At this time, the two-phase
of bubbles is significantly greater than the departure rate, which
flow changes from the previous stable state changes to another
results in the vapor phase aggregation and prevents the surround-
stable state, or into an unstable state with significant parameter
ing liquid from wetting the heating surface. By now, although lots
oscillation characteristics. As the thermal-hydraulic parameters
of experimental studies and theoretical analyses have been carried
(flow, pressure, temperature, etc.) significantly oscillate, the
out, there is still no theoretical model that can fully explain the
mechanical components of system can present fluctuation, which
occurrence mechanism of boiling crisis, as shown in Fig. 3. As for
may further lead to the heat transfer deterioration of reactor core.
the fuel assembly, the current method of predicting boiling crisis
The authors conducted a systematic study on flow instability and
is generally dependent on the empirical relationship (Liu et al.,
formed a relatively complete classification of two-phase flow
2020). According to the research foundation of authors (Lu et al.,
instability (Lu et al., 2019), as shown in Fig. 1. It can be seen that
2019a, 2019b), it is believed that CHF can be effectively increased
the more attention should be paid to the dynamic two-phase flow
by suppressing the vapor accumulation on the heating surface
instability in the actual engineering, which reflects the dynamic
before boiling crisis occurring, and promoting the rewetting effect
interaction process among parameters such as flow, pressure,
of heating surface, as shown in Fig. 4.
and temperature. In addition, it is believed that the pressure drop
proportion of two-phase and single-phase sections can signifi-
1.2. Steam generation
cantly affect the occurrence of flow instability (Lu et al., 2020;
Hu, 2017), as shown in Fig. 2. Therefore, the authors believes that
The steam generator is a key energy transfer device connecting
the two-phase flow instability can be suppressed or eliminated by
the primary and secondary loop systems of nuclear reactor.
optimizing the phase distribution to weak the sudden variation of
According to the fluid flow mode of secondary loop in the steam
pressure drop and reduce the resistance proportion of two-phase
generator, the steam generators of pressurized water reactor can
section, or by achieving the horizontal mixing of coolant between
be divided into the natural circulation type and the once-through
different sub-channels.
type (Yan and Cao, 2010)。
(3) Critical heat flux
(1) Natural circulation steam generator
The boiling crisis phenomenon refers to the variation of boiling
Natural circulation steam generators can be divided into two
heat transfer mechanism, and the temperature of heating surface
types: the vertical type and the horizontal type. Firstly, the vertical
rapidly rises, which may leads to the burning of fuel element. At
U-type tube natural circulation steam generator is commonly used.
this time, the heat flux is called as the critical heat flux (CHF). As
The typical vertical U-type tube natural circulation steam genera-
is know to all, the boiling crisis phenomenon is usually divided into
tion is generally composed of a lower head, a tube sheet, a U-
Rayleigh-Taylor Ledinegg
instability Fundamental Pressure drop Compound instability
dynamic instability oscillations dynamic instability
Flow distribution
Bubble collapse instability
Thermal oscillations
Flow pattern
Density-wave Thermo-acoustic transition
oscillations oscillations
Oscillations in
parallel channels Boiling crisis
3
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Fig. 2. The Two-phase flow instability and the two-phase interface evolution characteristics (Hu, 2017).
Droplet
After CHF sedimentation After CHF Bubbl e
boundary
Liquid Film
Droplet
Before CHF entrainment Before CHF
V a p o r
V a p o r
After CHF
Micro-liquid
layer
D r y w a l l
p h a s e
p h a s e
Before CHF
B e f o r e C H FA f t e r C H FI n t e r f a
s e p e r a t
type tube bundle, a steam-water separator and a cylinder. The heat and boiling heat exchange, the steam-water mixture flows upward
flux of primary loop is driven by the circulating water pump to through the steam-water separation device, and the saturated
flow in the U-type tube bundle to transfer the heat to the sec- steam flows into the steam turbine. As the steam generator is run-
ondary working fluid. After two stages of convection heat exchange ning, the secondary fluid is driven by buoyancy to establish the
4
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Fig. 4. Evolution of the liquid film at the bottom of the bubble (Lu et al., 2019).
natural circulation flow due to the density difference between the ent, and the pressure of secondary side steam is sensitive to the
single-phase fluid in descending channel and the steam-water external load, which may present the flow instability and the over-
mixture in ascending channel. As for the natural circulation steam all pulsation.
generator, there is no high requirement for the automatic control Considering the current development trend of once-through
of feed water and steam due to the buffering effect of large water steam generator, the design of smaller size and higher heat transfer
storage volume on the secondary side. However, the natural circu- efficiency is one of the most important directions in the future.
lation steam generator can just produce saturated steam. To ensure Meanwhile, the focus is to improve the heat transfer efficiency
the steam quality, the steam-water separator is necessary, which by changing the structure of heat transfer tube. The common struc-
leads to the large size. As for the miniaturization requirement of tures of heat transfer tube for once-through steam generator
nuclear reactor, the application of natural circulation steam gener- include the straight tube, the sleeve, the spiral tube, and the spiral
ator is much limited. sleeve, as shown in Fig. 5. Compared with the straight tube design,
the spiral tube not only has the significant heat exchange area
(2) Once-through steam generator between the primary and secondary working fluids, and the cen-
trifugal force generated by spiral structure is also conducive to
Compared with the natural circulation steam generator, the the uniform two-phase distribution of secondary side, which is
once-through steam generator has obvious advantages, such as conducive to the substantial improvement of boiling heat transfer
the small size, the light weight and the high heat exchange perfor- performance. At present, KLT-40S reactor of Russia (V. J., B. R. M., S.
mance. Meanwhile, the flow of secondary working fluid is driven R, 2009) and the IRIS reactor of USA (S. G. D., M. D. Carelli, 2008)
by the water pump head, and the secondary working fluid once use the once-through steam generator with spiral tube. As for
flows through the heat transfer tube, which then flows through the spiral sleeve structure, it is composed of the spiral inner tube
the preheating section, the evaporation section and the super- and the outer tube, which can form an obvious arc-shaped flow
heated section to generate the superheated steam. According to channel structure. This structure has the following advantages:
the secondary working fluid flowing outside or inside the heat 1) the inner tube and the outer tube can form the fulcrum, thus
transfer tube, the once-through steam generator can be divided the thermal expansion difference between them can be compen-
into the out-tube type and the in-tube type. Also, the in-tube sated; 2) The spiral flow of working fluid can increase the heat
once-through steam generator is the most common in the nuclear transfer coefficient. 3) The spiral inner tube can significantly
reactor. increase the tube length, thus the heat exchange area can be fur-
As for the once-through steam generator, it can directly gener- ther increased. 4) As for the arc-shaped structure, the sharp corner
ate the superheated steam with the superheat of about 20 ℃~30 ℃ can form an artificial vaporization core, which is easy to generate
and the stable pressure. Also, there is no need for dehumidification, the vapor bubbles. In the bubble form process with working fluid
which greatly reduces the requirement on the steam turbine and flowing upwards, the bubbles can continue to coalescence with
effectively improves the thermal efficiency of nuclear reactor. each other. Under the special arc-shaped structure, the coalesced
Since there is no complicated steam-water separation device, the bubbles trend to the areas with less resistance and larger cross-
structure of once-through steam generator is simple and compact. section. As a result, the working fluid in these areas correspond-
Compared with the natural circulation steam generator, the once- ingly flows to the narrow section, and fills the space left by bub-
through steam generator is more available for the limited space. In bles. In addition to flowing upward with the main flow, the
addition, the once-through steam generator has small heat storage bubbles can also produce the strong circumferential flow, this dis-
capacity and water storage capacity, thus the heating and cooling turbance in the circumferential direction is much stronger than the
processes of heating surface are easy to achieve the uniformity. circumferential flow of working fluid, which can enhance the boil-
Therefore, the once-through steam generator can quickly start ing heat transfer performance. In addition to the above advantages,
and stop, which can well meet the mobility requirements of the large arc-shaped channel has the supplement capacity of work-
nuclear reactor. However, the heat transfer situation of once- ing liquid, which can effectively prevent the heating surface from
through steam generator is more complicated than the natural cir- drying out.
culation steam generator. The once-through steam generator has
six heat transfer conditions in the secondary side, including the 1.3. Special safe system
single-phase convection, the subcooled boiling, the nucleated boil-
ing, the forced convection with evaporation, the film boiling and The special safety system is a general term for various safety
the superheated steam. Meanwhile, each heat transfer mechanism systems, which are specially set up for nuclear reactor under vari-
between the working fluid and the heating surface is quite differ- ous accident conditions and perform the safety functions to control
5
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Outer tube
Arc-shaped
structure
Fig. 5. The schematic diagram of heat transfer tube structure for steam generator (V. J., B. R. M., S. R, 2009; S. G. D., M. D. Carelli, 2008).
the accident consequences. The special safety system mainly matic diagram of typical passive waste heat removal system is
includes: the waste heat removal system, the safety injection sys- shown in Fig. 6 (Fei and Gu, 2016).
tem, and the molten material retention system. By means of summarizing the existing passive waste heat
removal system, it can be found that the existing design ideas of
(1) Waste heat removal system passive waste heat removal system mainly include two types
(Yan and Cao, 2010): 1) The direct cooling of primary working fluid
The waste heat removal system is a system used to remove the is used to export the waste heat, and the heat of primary working
waste heat of reactor core as the reactor shuts down. It is also fluid is discharged to the final heat sink according to the heat
called the shutdown cooling system, which mainly includes: the exchange device. 2) The waste heat is exported according to the
active waste heat removal system and the passive waste heat secondary side feed water cooling, namely the two-phase flow cir-
removal system. culation in steam generator. Also, it can be seen that no matter
The main function of waste heat removal system is shown as what design idea is adopted, the heat removal capacity of final heat
follows: the remaining heat of nuclear core and the sensible heat sink directly determines the performance of passive waste heat
of working fluid as well as devices can be released as for the emer- removal system. As shown in Fig. 6, there is the natural circulation
gency shutdown under normal and accident conditions. Mean- flow in heat exchange tube, and the water tank outside the heat
while, the temperature of reactor working fluid can decrease to exchange tube includes two heat transfer mechanisms, namely
the target temperature with a certain rate. The active waste heat the pool boiling heat transfer and the natural convection heat
removal system needs to rely on the supply of external energy. transfer in large space (Men, 2015).
Under the condition of the whole plant power failure, all the power
supply systems fail, and the active waste heat can not be available. (2) Safe injection system
Compared with the active waste heat removal system, the passive
waste heat removal system uses the natural principles such as the The safety injection system is also known as the emergency
gravity and convection of working fluid to make the nuclear reac- core cooling system, and its main function is to provide the cooling
tor safety, which is less dependent on the active systems such as to the nuclear reactor core under the abnormal conditions, so as to
pump and diesel engine (Zhang, 2012). Therefore, the passive maintain the integrity of fuel assembly. The safety injection system
waste heat removal system has become an important part of generally consists of the safety injection box, the high-pressure
advanced special safety system for nuclear reactor. Also, the sche- safety injection subsystem and the low-pressure safety subsystem.
6
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
a) The schematic diagram of passive waste heat removal system on the secondary side
Moreover, the action of safety injection system is shown as follows As for LOCA (Lack of Coolant Accident), the coolant is continu-
(Kong, 2014), as shown in Fig. 7: the high-pressure safety injection ously discharged from the primary loop into the containment, thus
system is firstly put into operation. The high-pressure safety injec- the primary loop is continuously depressurized, and the coolant
tion pump takes water from the refueling water tank and inject the volume is continuously reduced, which correspondingly causes
boron-containing water into the nuclear reactor core. As the the water level of pressure vessel to decrease below the lower
nuclear reactor pressure decreases below the pressure storage end of nuclear reactor core. As the pressure of primary loop
tank, the pressure storage tank safety injection system starts to decreases below the nitrogen pressure of safety injection tank,
feed the water injection. As the nuclear reactor pressure decreases the safety injection system starts to inject the water into the pri-
below the injection pressure of low-pressure safety injection tank, mary loop. The emergency cooling water injected by safety injec-
the pressure storage tank is emptied. Then, the pressure storage tion system reaches the bottom of nuclear reactor core. Initially,
tank is withdrawn from operation, and the low-pressure safety the nuclear reactor core is basically exposed, and a small amount
injection system is put into operation, while the low-pressure of heat is taken away by thermal radiation and natural convection.
safety injection pump draws water from the refueling water tank. Due to the release of decay heat, the temperature of nuclear reactor
It should be noted that the flow rate of low-pressure safety injec- core correspondingly rises. As for the pressurized water reactor,
tion tank is much greater than that of high-pressure safety injec- the heating rate at this stage is about 8–12 °C/s, the temperature
tion pump, namely the injection flow of high-pressure safety of fuel cladding may reach above 1100 °C. As the emergency cool-
injection system is small. As the water level of refueling water tank ing water reaches the bottom of nuclear reactor core, the reflood
reaches the low water level, the water source of low-pressure stage appears. Due to the high temperature of fuel cladding, the
safety injection pump is switched to the containment pit, thus water has been vaporized before contacting the cladding surface,
the system enters the recirculation cooling stage. which provides the initial cooling. As the temperature of fuel clad-
7
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
away, which can effectively cool the lower head of pressure vessel.
Under the severe accident conditions, there are three key condi-
tions for adopting the external cooling scheme of pressure vessel
and successfully implementing the retention strategy of molten
materials (Asmolov et al., 2001): 1) the reactor coolant system
has been depressurized to prevent the high-temperature creep fail-
ure and the integrity loss. 2) The molten materials of reactor core is
repositioned in the lower chamber of pressure vessel, and the
water level of reactor has reached a certain height. 3) It is neces-
sary for the reasonable design of external cooling channel, thus
ERVC circulating flow can ensure that the outer wall of vessel
can be sufficiently cooled. Therefore, the decay heat of molten
materials can be reliably discharged, and the pressure vessel can
be prevented from failing. It should be noted that this is also the
most direct link for the successful implementation of ERVC.
Since the 1950s, the international and domestic heat transfer
fields have carried out lots of researches on the heat transfer
enhancement technology. In general, the heat transfer enhance-
ment technology can be divided into the passive technology and
the active technology. The passive technology refers to the heat
transfer enhancement by just changing the physical properties of
fluid, the shape of channel and the heat structure design without
Fig. 7. The schematic diagram of safety injection system (Kong, 2014). additional power. Moreover, the active technology requires the
external power to change the fluid flow characteristics, so as to
intensityen the heat transfer performance. The passive heat trans-
ding decreases to a certain level, the water can contact the cladding
fer enhancement technology is subdivided into the structural inno-
surface, and the heat transfer coefficient is greatly increased, which
vation technology and the surface modification technology. Also,
can better reduce the surface temperature of fuel cladding. In this
the structural innovation technology is mainly divided into the
process, the surface characteristics of fuel cladding can signifi-
micro-channel technology and the longitudinal vortex technology.
cantly affect the contacting behavior of water, and the thermo-
The surface modification technology mainly includes the micro-
physical properties of water can also affect the heat taken away
nano structure technology and the surface coating technology. In
from the nuclear reactor core. In a word, these two factors deter-
addition, the authors subdivide the active heat transfer enhance-
mine the emergency cooling performance of safety injection
ment technology into the magnetic field heat transfer enhance-
system.
ment technology and the electric field heat transfer enhancement
technology. The technology route of this study is shown as Fig. 9.
(3) Molten material retention system
As the severe accidents occurs in the nuclear reactor, the 2. Theory of heat transfer enhancement
nuclear reactor core is exposed due to the loss of cooling water
and begins to heat up. The fuel assemblies are melted due to the The basic theory of heat transfer enhancement is an important
insufficient cooling, and the molten core falls into the lower cham- foundation for the application of heat transfer enhancement tech-
ber of pressure vessel. Once the pressure vessel melts and the mol- nology, which has received continuous attention in recent decades.
ten materials flow into the reactor chamber, the steam explosion As is know to all, eliminating the dead zone of flow backflow,
may happen, which can cause the temperature and pressure of enhancing the fluid disturbance and mixing, destroying the devel-
containment to increase. Finally, the containment integrity can opment of laminar boundary layer, promoting the departure of
be threatened. To mitigate the consequences of severe accidents,
the various strategies have been proposed to deal with severe acci-
dents. In-Vessel Retention (IVR) is one of the most important mit-
igation schemes of serious accident. The current IVR measures are
mainly divided into Internal Reactor Vessel Cooling (IRVC) and
External Reactor Vessel Cooling (ERVC). IRVC is to keep the integ-
rity of pressure vessel by directly injecting water into the pressure
vessel after severe accidents happening. ERVC aims to inject water
to the outside pressure vessel after severe accidents happening, so
as to fully cool the lower head and keep the molten materials in
pressure vessel.
At present, the external cooling implementation of the lower
head of pressure vessel has become the most widely used IVR mea-
sure based on the natural circulation to achieve the retention of
molten materials, which is generally called IVR-ERVC measure,
and the schematic diagram is shown in Fig. 8 (Guo, 2015). After
the severe accidents occurring, the water injection program is
started to flood the pressure vessel, and the cooling water enters
the flow channel between the lower head of pressure vessel and
the insulation layer. The cooling water is vaporized on the surface
of lower head to form the two-phase flow. In addition, the natural
circulation is established, and the heat is correspondingly taken Fig. 8. The schematic diagram of IVR-ERVC (Guo, 2015).
8
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
bubbles from heating surface, and inhibiting the vapor film cover- Then, Tao (Tao et al., 2002) from Xi’an Jiaotong University con-
ing heating surface can achieve the heat transfer enhancement. ducted lots of numerical calculations and experiments to verify
Additionally, increasing the heat exchange area, increasing the the field synergy theory. It was pointed out that the field synergy
average heat exchange temperature difference, and changing the theory could be used to guide the design of heat transfer devices.
heating surface characteristics can also improve the heat transfer According to the study of Liu et al. (Liu et al., 2009), the field syn-
coefficient. Although the heat transfer enhancement technology ergy theory of more physical quantities in convective heat trans-
can enhance the heat exchange between the heating surface and fer was analyzed. Therefore, the multi-field synergy theory was
the fluid, it can also increase the pressure consumption of flow pro- further proposed: the gradient coordination angle a of velocity
cess. Moreover, the increase of pressure consumption can lead to and mainstream direction velocity, the coordination angle b of
an increase of required pump power, which is not conductive to velocity field and temperature gradient field, the coordination
the energy-saving benefits. Some scholars have conducted theoret- angle h of velocity and pressure gradient could reflect the distur-
ical researches to reveal the scientific mechanisms of heat transfer bance degree of fluid, the synergy of heat consumption and the
enhancement technology, and tried to provide the theoretical guid- synergy of flow power consumption. Meanwhile, the synergy
ance for heat transfer enhancement theories. The current heat angles of convective heat transfer process could restrict each
transfer enhancement theories mainly include: the field synergy other, which explained the laws of heat transfer, flow and power
theory, the configuration theory, the entropy production mini- consumption. Recently, Liu et al. (Liu et al., 2017) constructed the
mization theory, the entransy dissipation theory, and the process constitutive relationship between the fluid flow and the pressure
enhancement theory based on flow field optimization. gradient, and established the conservation equation of mechanical
energy expressed in the pressure form. Thus, a physical quantity
2.1. Field synergy theory describing the functional force of fluid was proposed, namely the
available potential. Tao et al. (Tao et al., 2002) investigated the
Since the 1990s, Guo (Guo and Zhuang, 1992) from Tsinghua applicability of field synergy theory in the process of elliptical
University established the synergy equation based on the energy flow heat transfer, which was successfully applied to optimize
conservation, which revealed the relationship between the dot the arrangement of flat plates and the structure of bellows. More-
product of velocity vector as well as temperature gradient and over, Tao et al. (Tao et al., 2002) suggested the three major mech-
the Nusselt number (Nu). The above method was called the field anisms in the heat transfer enhancement, namely thinning the
synergy theory, which revealed the heat transfer enhancement thermal boundary layer, increasing the fluid disturbance, and
from a new perspective. Under the same Reynolds number (Re) increasing the wall velocity gradient, which could improve the
and Prandtl number (Pr), the effect of convection heat transfer synergy of the velocity field and the temperature gradient field.
was dependent on the synergy degree between the velocity vector Obviously, it could be seen that the field synergy theory was con-
and the temperature gradient, and the synergy degree between sistent with the experience of heat transfer enhancement. The
them was related to the following three aspects (Guo et al., experience is consistent. As proposed by Zhu et al. (Xiao and
2005): 1) the synergy degree is inversely proportional to the angle Jing, 2016), the local synergy angle was used to analyze the con-
value between the speed vector and the temperature gradient, vective heat transfer, which could overcome the limitations of
namely proportional to the direction consistence of them. 2) As average synergy angle evaluation method, and he also proposed
for the local position, the speed value, the temperature gradient the application criteria of incompressible flow at the conventional
value and the cosine value of them were proportional to the filed scale. In addition, there were also good application of field syn-
synergy degree. 3) The velocity and temperature distribution on ergy theory in the electric vehicle cabin (Xu et al., 2017), the liq-
each section should have good uniformity. After the field synergy uefied petroleum gas vehicle (Ou et al., 2014), and the diesel
theory proposed, a large number of scholars have conducted in- particulate filter (Deng et al., 2007).
depth researches.
9
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
2.2. Configuration theory temperature heat source with infinite volume, the entransy of
object would be fully derived. Meanwhile, the entransy repre-
In the 1990s, Bejan (B. A., K. J., 1988) proposed the configuration sented the heat transfer performance of object, and simultaneously
theory based on the biological evolutionary law in nature. This the- had the characteristics of energy and transfer ability. Then, Guo
ory points out that the flow system trends to provide the easier (Guo et al., 2007) proposed the entransy dissipation theory, which
path for the ‘‘flow” flowing through it, which undergoes the struc- could also be considered as the theory of minimum thermal resis-
tural evolution along the direction of time arrow. The flow system tance. As the heat transfer temperature difference was given, the
of configuration theory includes the various devices in the field of optimal result was that the entransy dissipation was the largest,
engineering technology. In addition, the ‘‘flow” of configuration which referred to maximum the heat flux. As the heat flux was
theory includes the fluid flow, the heat flow, the mass flow, the given, the optimal result was that the entransy dissipation was
current, the vehicle flow, and the information flow. According to the smallest, and the corresponding temperature difference was
this view, the configuration theory can unify lots of problems in the smallest. As for the above two cases, the equivalent thermal
different fields into the structural or path optimization problem. resistance was the smallest.
The configuration theory was firstly applied to the cooling prob- Furthermore, Liu et al (Liu et al., 2011) derived the entransy
lem of electronic component (B. A., K. J., 1988). The maximum tem- equation of transmission process, and the flow field of square cav-
perature of heating body was taken as the objective function, and ity was optimized by minimizing the entransy dissipation. Wu
the basic structural unit was optimized. Then, the several basic et al. (W. J, L. X, 2008) applied this theory to the field of radiation
structural units were assembled to obtain the first-level structural heat transfer, and it was proposed that there was partial entransy
unit, which was further optimized until the thermal conductive dissipation in the process of radiation heat transfer. Also, the
structure could completely cover the entire heating region. This entransy dissipation theory could guide the optimization of radia-
gradual assembly optimization method could achieves the goal of tion heat transfer process. Wang et al. (W. S, J. G, W. J, 2018) stud-
minimizing temperature. Moreover, the configuration theory was ied the geometric parameter effect of spiral wound heat exchanger
also applied to the transportation economy, the urban water sup- on the equivalent thermal resistance. The results showed that the
ply network and the power network optimization. four geometric parameters of spiral angle, layer spacing, tube spac-
ing, and outer diameter were inversely proportional to the equiva-
2.3. Entropy minimization theory lent thermal resistance caused by flow, namely the entransy. Also,
Yuan et al. (Yuan and Chen, 2011) optimized the energy utilization
As is known to all, the second law of thermodynamics is widely efficiency of indirect evaporative cooling system by means of
applied to the various thermal systems for the irreversible analy- entransy dissipation theory.
ses. Thus, Bejan et al. (B. A., 1979) applied it to the analyses of heat
transfer process, the expression of volume entropy production rate 2.5. Process enhancement theory
was proposed. The irreversibility of convective heat transfer pro-
cess included the temperature difference heat transfer and the vis- In recent years, the process enhancement method based on flow
cous flow resistance. These two sources were different and could filed optimization has been applied to the convection transfer pro-
affect each other. Meanwhile, the reduction of one irreversibility cess. With the specific constraints and the objective function as the
would lead to the increase of another irreversibility. Therefore, theory core, the flow field is adjusted and the transfer process is
the purpose of convection heat transfer optimization was to reduce artificially made. By now, this optimization method has been
the sum of these two irreversibility. applied to the fluid flow resistance reduction and the mass transfer
After the principle of volume entropy generation rate in the enhancement (Jia et al., 2012; C. Q., W. M., P. N., 2009).
field of heat transfer, the irreversible analyses of heat transfer pro- As the core of flow field structure optimization, the objective
cess have been widely studied. Zhu et al. (Zhu and Yu, 2017) opti- function directly determines the optimization direction of flow
mized the heat sink structure of thermoelectric cooler by the structure. The convective heat transfer process is taken as an
principle of minimum entropy generation, and the cooling liquid example, the typical thermodynamic optimization methods
with suitable heat capacity was selected. Also, Zhao et al. (Zhao include the entropy minimization theory and the entransy dissipa-
et al., 2017) optimized the geometric structure of rectangular tion theory. The flow field optimization needs to consider the
finned elliptical tube by the entropy generation analysis. The fluid mechanical energy loss (namely the pressure drop) caused by the
flow and the heat transfer characteristics of finned elliptical tube in flow. Therefore, the mechanical energy loss is generally considered
different arrangements were studied, and a structure was obtained as a constraint condition in the optimization process. From the
based on the minimization entropy generation. Vargas et al. above two views, it can be seen that the optimization of flow field
(Vargas and Bejan, 2000) studied the boiling heat transfer with dif- can be transformed into a target extremum problem with con-
ferent temperatures, and the fluid mass flow ratio was optimized straints. By solving the target extremum problem, the correspond-
by this principle. Sahiti et al. (Sahiti et al., 2008) applied the ing optimal flow field can be obtained. The specific flow field
entropy minimization theory to the structural optimization of obtained by this theory could provide a direct basis for the optimal
pin–fin heat exchanger, including the length of pin–fin and the design of device structure (Li et al., 2014; M. J. A., L. X. G., L. Z. X.,
number of subchannels. Cui et al. (Cui et al., 2010) analyzed the 2005). This process enhancement method based on optimizing
heat exchange performance with multiple fluids according to the flow field can provide a new optimization method. Firstly, the opti-
entropy minimization theory, and the influence of channel mal flow field structure is constructed from the objective function,
arrangement on the entropy generation was studied. and then the best equipment structure is searched from the per-
spective of optimizing flow field, so as to avoid the empirical com-
2.4. Entransy dissipation theory paring method in the past.
From the above analyses, it can be seen that the reasonable flow
Guo (Guo et al., 2007) also defined a heat transfer potential field structure can promote the heat exchange between the fluid
capacity, namely entransy, whose physical meaning was explained and the wall. In addition, the second law of thermodynamics points
as follows: an object could be regarded as a container for storing out that the irreversibility is the decline result of heat conversion
heat energy, and its potential energy could be analogous to the efficiency. Although the convective heat transfer process is not a
electric potential energy. As the object contacted with the low- typical heat conversion process, the entropy production (namely
10
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
the irreversibility degree) of system or process is closely related to the channels with diameters of 1.10 mm and 2.01 mm presented
the convective heat transfer. Thus, the process can be obtained the typical characteristics of micro-channel, while the channels
indirectly from the entropy production information of system. with diameters of 2.88 mm and 4.26 mm presented the obvious
Guo (Guo, 2015) suggested that the entropy production of convec- characteristics of conventional channel. Nevertheless, Ribatski
tive heat transfer process consisted of two parts: the flow entropy et al. (Ribatski et al., 2006) believed that the above criteria could
production and the heat transfer entropy production. The flow not reflect the scale influence of channel on the physical mecha-
entropy production was related to the flow field structure and nism, and the difference of two-phase flow and heat transfer char-
was caused by fluid overcoming shear stress. Also, the heat transfer acteristics between them should be taken into account. In addition,
entropy production was related to the temperature field, which Thome (Thome, 2004) suggested that the criterion for the micro-
was caused by the irreversibility of heat transfer. Therefore, Guo channel should consider the separation diameter of bubble and
(Guo, 2015) proposed that under a given viscosity dissipation, the restriction effect of channel wall. Furthermore, Kew and Corn-
the continuity equation and the energy balance equation were well (Kew and Cornwell, 1997) selected the restriction coefficient
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
simultaneously satisfied, and the heat transfer entropy production (N conf ¼ gðq rq ÞD2 ) as the criterion for these two channels. Under
l g e
could reach the extreme value by adjusting the flow field structure.
The heat transfer tube structure and the temperature field distribu- the atmospheric pressure, the channel could be defined as the con-
tion characteristics obtained by this optimization theory were ventional channel for Nconf < 0.5 (namely De > 5 mm). Also, the
shown in Fig. 10. channel could be defined as the micro-channel for Nconf > 0.5
(namely De < 5 mm) . Then, Cheng (Cheng and Wu, 2006) et al. pro-
posed the boundary between these two channels based on Bond
3. Passive heat transfer enhancement technology number. As for Bond > 3 (Dh > 4.34 mm) under the atmospheric
pressure, the channel could be defined as the conventional runner,
3.1. Structural innovation technology which could be conversely considered as the micro-channel.
Although there were many methods to divide the micro-channel
3.1.1. Micro-channel technology and the conventional flow channel, no matter which kind of defini-
It is believed that the micro-channel can significantly enhance tion method, the two-phase characteristics in micro-channel were
the heat transfer performance compared with the conventional different from the conventional channel. For example, since the
channel. On the one hand, since the volume is the cube of charac- surface tension in the micro-channel was more significant than
teristic length, and the surface area is the square of characteristic other forces, the integrity and the fluctuation of two-phase inter-
length, it can be seen that the ratio of surface area to volume is face in micro-channel were better than the conventional channel
inversely proportional to the characteristic length. Therefore, the (Lu et al., 2015). In addition, the bubble could quickly grow up in
micro-channel has the large ratio of surface area to volume. On the micro-channel to contact the channel wall, which further pro-
the other hand, for the single-phase laminar flow in geometrically moted the coalescence of surrounding bubbles. Moreover, the coa-
similar channels, Nu = hDe/k = constant, the heat transfer coeffi- lesced vapor slugs could cover the heating wall and block the flow
cient can be obtained, which is inversely proportional to the char- channel.
acteristic length, so the micro-channel has the large convective As is well known, there are two main mechanisms for flow boil-
heat transfer coefficient. In addition to the above-mentioned ben- ing heat transfer in micro-channels: the nucleate boiling heat
efits under the single-phase flow condition, the micro-channel transfer and the forced convection heat transfer with evaporation.
technology exhibits the significant heat transfer enhancement per- Charnay et al. (Charnay et al., 2015) suggested that the nucleate
formance under the boiling conditions (Zan, 2013). boiling heat transfer was closely related to the bubble formation
However, there is still no consensus on the scale definition of on the heating wall, which was significantly affected by the heat
micro-channel (Kandlikar, 2002). It was believed that the scale flux and the system pressure. Also, the forced convection heat
limit between the conventional channel and the micro-channel transfer with evaporation mainly depended on the heat transfer
was 3 mm. However, Mehendal et al. (Mehendale et al., 2000) con- of micro liquid film and the evaporation of two-phase interface.
sidered the size of 6 mm as the boundary between them. According Thus, the heat transfer performance was mainly dependent on
to the research of Chen et al. (Chen et al., 2006), it was found that the flow rate and the quality. Meanwhile, Tran et al. (Tran et al.,
a) The heat transfer tube structure b) The temperature field distribution characteristics
Fig. 10. The heat transfer tube structure and the temperature field distribution characteristics (Guo, 2015).
11
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
1996) proposed that the transition point of nucleate boiling heat influence on the heat transfer coefficient. By means of experimen-
transfer and forced convection heat transfer with evaporation tal investigation in micro-channel, Pan et al. (L., M., C. He, 2002)
was mainly determined by the wall superheat. Moreover, the concluded that the boiling heat transfer enhancement in micro-
nucleate boiling heat transfer mainly occurred under the condition channel was mainly due to the disturbance of bubble on the sur-
of high wall superheat, while the forced convection heat transfer rounding fluid. Moreover, Guo and Zhu (Guo et al., 2007) proposed
with evaporation mainly appeared with low wall superheat. How- that the channel scale reduction could restrict the bubble growth,
ever, Vlasie et al. (Vlasie et al., 2004) and Collier-Thome (Collier which could promote the coalescence with the neighboring bub-
and Thome, 1996) found that the these two boiling heat transfer bles. As the bubbles took away lots of heat in the sliding process,
mechanisms could coexist under the specific working conditions, they also took away the bubbles attached on the heating wall,
and the heat transfer coefficient was simultaneously determined which further promote the bubble renewal. However, a large
by parameters, such as flow rate, system pressure, and quality. amount of bubbles would block the flow channel and worsen the
Additionally, some researchers (Tran et al. (Tran et al., 1996), Bao heat transfer performance under the condition of higher heat flux.
et al (Bao et al., 2000) and Ali et al. (Ali et al., 2011) concluded that This conclusion was inconsistent with the results of Passos et al. (J.
the heat transfer coefficient in micro-channel just depended on the C. P., F. R. H., L. F. B. P, M. B., M. M, 2004) and Misale et al. (M. M, P.
heat flux, which was not significantly related to the flow rate of A. G. G, 2011). In addition, Xia et al. (Xia et al., 1996) proposed that
working fluid and the quality. That is to say, the nucleate boiling there was an optimal scale for the flow boiling heat transfer in the
heat transfer in micro-channel could play a major role. However, micro-channel, which was not conducive to the heat transfer per-
Lin et al (Lin et al., 2001) suggested that both of the nucleate boil- formance as the channel scale was too large or too small.
ing heat transfer and the forced convection heat transfer with
evaporation would occur in the micro-channel. Similar to the con- 3.1.2. Longitudinal vortex technology
clusions of Lin et al. (Lin et al., 2001), Huo et al. (Huo et al., 2004) The longitudinal vortex technology refers to the improvement
used R134a as the experimental fluid to study the boiling heat of the flow channel structure or the thermal component form
transfer in micro-channel. It was found that the nucleate boiling according to series of processing and manufacturing processes, so
heat transfer could play a major role with the channel diameter as to achieve the forced mixing of coolant and the significant dis-
of 2.01 mm and the quality less than 20% to 30%. As for the turbance of thermal boundary layer. As the fluid flows across the
micro-channel with diameter of 4.26 mm and the quality less than obstacle, the vortexes can be generated in the downstream space,
40%~50%, the nucleate boiling heat transfer could be dominant. which includes the transverse vortex whose rotation axis is per-
With the further in-depth research, Kew and Cornwell (Kew and pendicular to the main flow direction, and the longitudinal vortex
Cornwell, 1997) found that the nucleate boiling heat transfer whose rotation axis is consistent with the main flow direction. The
played a major role in the isolated bubbly flow area. As the bubbles above phenomena can promote the momentum and energy mixing
occupied the entire channel, namely the flow pattern transforming between the mainstream and the fluid near wall, which can further
into the slug flow, the relationship between the heat transfer coef- strengthen the turbulence and the disturbance of thermal bound-
ficient and the heat flux was no longer obvious. As the flow pattern ary layer. Therefore, the heat transfer enhancement can be
gradually transferred from the slug flow to the annular flow, the achieved. Moreover, the studies have shown that if the lateral size
boiling heat transfer mechanism mainly presented the characteris- of obstacle is limited and is suitable for the incoming fluid angle,
tics of forced convection heat transfer with evaporation. Karayian- the generated vortex will follow the mainstream movement to
nis et al. (Karayiannis et al., 2012) also found that the boiling heat form a series of ordered longitudinal vortexes, which will greatly
transfer mechanism gradually evolved from the nucleate boiling enhance the heat transfer performance. Moreover, the longitudinal
heat transfer to the forced convection heat transfer with evapora- vortex technology has been widely used in the field of nuclear
tion as the heating section length increased. Also, Sumith et al. reactor.
(Sumith et al., 2003) suggested that as the flow pattern transferred
to the annular flow in micro-channel, the boiling heat mechanism (1) Fuel assembly field
transferred from the nucleate boiling heat transfer to the forced
convection heat transfer with evaporation. As the liquid film of At present, most of the fuel assemblies in commercial nuclear
annular flow was thick enough, there would be the nucleate boil- power plants are designed with the cylindrical fuel bundle struc-
ing inside the liquid film. Kasza et al. (K. K. E, D. T, W. M. W, ture, as shown in Fig. 11a. The sub-channels of this open bundle
1997) found there was the obvious thin liquid film between the type fuel assembly are connected to each other, which can promote
vapor phase and the heating wall in the micro-channel under the the lateral fluid mixing and the homogenization of fluid tempera-
conditions of slug flow as well as annular flow, and the existence ture distribution. It should be noted that this type of fuel assembly
of nucleate boiling phenomenon could be observed. Moreover, structure can also eliminate the two-phase flow instability from
Balasubramanian and Kandlikar (B. P, K.S. G, 2005) concluded that the perspective of physical mechanism. Since increasing the heat
the nucleate boiling could occur in both of the mainstream and the transfer area per unit volume of fuel assembly is an effective
micro liquid layer in the micro-channel, and the bubble generation way, the annular fuel assembly with double-sided cooling charac-
frequency in the micro liquid layer was higher than the indepen- teristics has been proposed (as shown in Fig. 11b). It can be seen
dent bubble generation frequency in the mainstream. that the annular fuel rod can effectively reduce the radial thermal
In addition, a large number of studies had shown that the resistance and the maximum temperature of pellet. Meanwhile, it
micro-channel could significantly improve the flow boiling heat can also increase the convection heat transfer area between the
transfer performance. Owhaib et al. (Owhaib et al., 2004) found fuel assembly and the fluid. Thus, the power density and the safety
that the reduction would lead to the increase of heat transfer coef- margin of the fuel assembly can be improved. By means of further
ficient. Meanwhile, Yan and Lin (Yan and Lin, 1998) pointed out reducing the fuel rod spacing, the tight-lattice fuel assembly (as
that the heat transfer coefficient in micro-channel was 30% to shown in Figu. 11c) can be formed, which thereby achieves the lar-
80% higher than the conventional runner. Also, it was suggested ger heat transfer area compared with the traditional rod bundle
that the heat transfer coefficient in micro-channel was about 15% fuel assembly. According to the previous discussion, another way
higher than that of conventional channel (Shen et al., 2001). Mean- to effectively improve the thermal-hydraulic performance of fuel
while, Sobierska et al. (Sobierska et al., 2007) proposed that the assembly is to enhance mixing effect of fluid in the sub-channels
reduction of channel scale would significantly enhance the quality of fuel assembly. Thus, most of the rod-bundle fuel assemblies
12
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Dw
Dr
e) HCF assembly
Fig. 11. The application of longitudinal vortex technology in the field of fuel assembly (McKrell et al., 2010).
are designed with the positioning grids with mixing vanes. Also, basis of tight-lattice fuel assembly, the spatial self-positioning fuel
these special structure can position the fuel rod and can effectively rod can be realized by means of the design of wire structure, which
enhance the mixing degree of fluid, as shown in Fig. 11a. On the further realizes the continuous mixing along the flowing direction,
13
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
as shown in Fig. 11d. Combined with the above advantages of fuel tics and CHF, a large number of experiments and numerical studies
assembly structures, MIT (Massachusetts Institute of Technology) had been carried out (Koroush, 2016). The centrifugal effect gener-
proposed an innovative design of HCF (Helical Cruciform Fuel) ated by swirling flow could promote the bubble separation from
(McKrell et al., 2010), and the schematic diagram is shown in the heating wall, which was beneficial to enhance CHF of HCF
Fig. 11e. Obviously, it can be seen that during the spiral period of assembly. As the flow pattern in HCF assembly was the annular
fuel rod (b = 0° ~ 90°), the tips of adjacent fuel rods are contact with flow, the swirling effect would cause the liquid film to be destroyed
each other to achieve the spatial self-positioning without addi- in advance, which could worsen the heat transfer performance.
tional positioning grid. Moreover, HCF assembly has the open Based on VOF (Volume Of Fluid) two-phase flow model, Shirvan
structure characteristics of traditional rod-bundle fuel assembly, (K. S, 2016) compared the motion characteristics of a single bubble
which can fully eliminate the two-phase flow instability. Mean- in the conventional rod-bundle fuel assembly and HCF assembly.
while, the helical cladding structure of fuel rod can produce the The research showed that the bubble dynamics on the surface of
continuous mixing effect on the fluid in sub-channels in the whole HCF fuel rod were affected by the swirling flow field and the helical
flow process. In addition, the HCF fuel rod has the larger heat trans- heating surface. Moreover, the most typical feature was that the
fer area compared with the traditional rod-bundle fuel assembly. bubble on the surface of HCF fuel rod would slip over a long dis-
The relevant studies had shown that if HCF assemblies were used tance, and this bubble behavior would introduce a new boiling
to replace the those of Westinghouse AP1000 PWR (Pressurized heat transfer mechanism and enhance the nucleation core. Accord-
Water Reactor), the core temperature at maximum linear power ing to the existing boiling heat transfer experiments for HCF
density could be reduced by about 600 °C, which could greatly assembly, the swirling flow field generated by helical structure
increase the thermal safety margin of fuel assembly. Conversely, could effectively increase CHF (Yan, 2019). Also, Bolshakov et al.
the HCF assemblies could increase the core power density by about (B. V. V., B. S. M., K. L. L., 2007) used the 3 3 rod bundle of HCF
20%-45% with the same safety system and safety limits, which assembly to conduct CHF experiment under the working condition
could significantly improve the nuclear reactor economy (C. T. of 15.7 MPa. The study showed that CHF of HCF assembly with high
M., M. T. J., K. M. S., 2014). In view of the unique thermal perfor- quality was about 10% larger than that of the cylindrical rod-
mance advantages of HCF assembly, Lighting Bridge company bundle assembly, but CHF of HCF assembly with low quality was
and Framaton company had conducted in-depth research on it, about 10%-25% smaller than that of the cylindrical rod-bundle
and they planned to use the HCF assembly for the large PWR of assembly. The reason for the above phenomenon might be that
Westinghouse and the small nuclear reactor of NuScale (M. J, T. the temperature measurement point could note accurately capture
A, S. N, 2012). Also, some companies in USA and France had com- the dry point, or it might be attributed to the uneven heat flux due
pleted the design and the manufacture of HCF assembly, which to the heating rod circumference.
was expected for engineering application in recent years.
At present, a large number of numerical and experimental stud- (2) Non-nuclear heat exchange device
ies had been carried out on the swirl flow of HCF assembly and its
influence on the mixing characteristics between sub-channels There are also many applications of structural innovation tech-
(Agbodemegbe et al., 2016). According to the numerical research nology in the field of non-nuclear heat exchange device. Ma et al.
conducted by Krauss and Meyer (Krauss and Meyer, 1996), it was (Ma et al., 2011) applied the longitudinal vortex in the laminar flow
found that the obvious fluctuations of fluid velocity and temperature zone to increase the heat transfer coefficient by up to 100.9%, while
in the side channels were mainly attributed to the heat and momen- the frictional resistance coefficient increased by only 11.4%. In the
tum exchange with the central channels. By means of PIV (Particle turbulent flow zone, the longitudinal vortex could achieve the heat
Image Velocimetry) technology, Xiong et al. (Xiong et al., 2020) transfer enhancement by more than 87.1%, while the friction resis-
investigated the development process of swirl vortex induced by tance obviously increased by 100.3%. Huang et al. (Huang et al.,
helical wire. Also, McClusky et al. (Holloway, 2002) suggested that 2011, 2013) suggested that the longitudinal vortex could increase
the attenuation rate of vortex in the rod bundle channel was signif- the boiling heat transfer coefficient by 1.1%~25.8% compared with
icantly greater than that in HCF assembly. Then, Shen et al. (Shen the smooth channel. Also, the corresponding bubble distribution
et al., 1991) observed the significant cross-flow phenomenon caused was more uniform under the action of longitudinal vortex, and
by helical wire based on LDV (Laser Doppler Velocimetry) technol- the boundary was is obviously thinned or destroyed. These reasons
ogy, and the generated swirl vortex would significantly enhance were the main factors leading to the heat transfer enhancement.
the lateral mixing between sub-channels. Moreover, Ikeda and Meanwhile, Wang et al. (Wang et al., 2007) proposed that under
Hoshi (I. K, M. H, 2006) installed LDV in a experimental rod coated the same mass flow, pressure drop, and pump power, the longitu-
with optical fiber to observe the swirling vortex and the lateral mix- dinal vortex structure of alternate arrangement had better heat
ing caused in HCF assembly. Conboy et al. (C. T. M., 2007) measured transfer performance than the central arrangement. Also, Meng
the single-phase flow resistance and the mixing characteristics test et al. (M. J. A., L. X. G., L. Z. X., 2005) designed the heat transfer tube
of HCF assembly, and it was found that the corresponding drag coef- with double diagonal rib, as shown in Fig. 12a. The results showed
ficient was about 10% lower than that of rod-bundle assembly. How- that the fluid was disturbed by the double diagonal rib to produce
ever, the above phenomena were affected by the swirling flow field, the longitudinal swirling flow, and the convective heat transfer
and the pressure distribution on the cross section of HCF assembly performance was significantly improved. Moreover, the tapered
was relatively uneven. As for the sub-channels of tight-lattice fuel tube (as shown in Fig. 12b) had been investigated, and its heat
assembly, the experimental and numerical studies had confirmed transfer performance was nearly 4 times that of smooth round
that there was a large-scale periodic vortex structure in the gap tube (R. L. W., 2009). Xu et al. (Xu et al., 2016) found that there
regions between fuel rods, which could also effectively enhance was a vortex structure in the rear area of helical rib tube, which
the inter-rod mixing intensity (Lee et al., 2014). As for the single- could periodically disturb the fluid, which thereby improved the
phase flow pressure drop of 5 5 rod-bundle channel for HCF heat transfer performance, as shown in Fig. 12c. Wang et al. (W.
assembly, the flow resistance coefficient was obtained by Conboy W., Z. Y., L. B., 2018) researched the heat transfer enhancement
(McKrell et al., 2013), and the mixing coefficient of HCF assembly performance of laterally convex tube, as shown in Fig. 12d. It
was calculated by the thermal mixing method. was found that the fluid would present a jet effect while flowing
Since the two-phase distribution characteristics in the fuel through the ripples, and these ripples would regularly disturb
assembly could directly affect the boiling heat transfer characteris- the flow boundary layer. In addition, the formed swirl flow could
14
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
suppress the pulsation and the secondary flow, which correspond- showed that the cross-sectional variation of this tube could induce
ingly reduced the flow resistance. Sun et al. (S. M, Z. M, 2018) stud- the multiple longitudinal vortexes, which was in favor of the heat
ied the flow and heat transfer performance in the transversely transfer enhancement. Compared with the above heat transfer
concave tube were studied, and its structure was shown in tubes, the cross-zoomed elliptical tube could enhance the heat
Fig. 12e. It was found that the heat transfer and the flow synergy transfer while the flow resistance less increasing. Also, Sajadi
of it was better than the smooth tube. et al. (S. A. R, Y. D. S, A. D, 2014) compared the heat transfer perfor-
In addition to the above mentioned ribbed heat transfer tube, mance and the flow resistance of cross-zoomed elliptical tube, flat
the surface groove is also an effective structural innovation tube, and round tube. It was found that the cross-zoomed elliptical
method. In general, the flow resistance of surface groove tubes is tube presented the best heat transfer enhancement performance.
relatively small. Pirbastami et al. (P. S, M. S. F, M. S. G, 2016) Moreover, Tan et al. (T. X. H, Z. D. S, Z. G. Y, 2013) studied the flow
numerically studied the heat transfer enhancement performance and heat transfer performances of tube side and shell side for the
of helically outer grooved tube, as shown in Fig. 13a, and the influ- twisted tube, as shown in Fig. 14b. The results showed that there
ence of groove spacing on the heat transfer as well as the friction was a secondary flow in the twisted elliptical tube, which changed
factor was paid much attention. Also, Bilen et al. (B. K, C. M, G. H, the speed and temperature of fluid in this tube, and reduced the
2009) investigated the effects of circular groove, rectangular synergy angle between the temperature gradient and the velocity
groove, and trapezoidal groove on the heat transfer performance, vector. Thus, the heat transfer performance could be obviously
as shown in Fig. 13b. It was found that the groove structure could improved. Next, Jin et al. (J. Z. J, C. F. Q, G. Z. X, 2017) proposed a
effectively improve the heat transfer performance of heat transfer new type of six-head spiral corrugated tube, as shown in
tube, which had been increased by 63%, 58% and 47% respectively. Fig. 14c. Also, the effects of geometric parameters, Reynolds num-
Moreover, Aroonrat et al. (A. K, J. C, 2013) researched the flow and ber and fluid physical properties on the heat transfer performance
the heat transfer performance of fluid in the helically inner were studied. It was found that the heat transfer coefficient
grooved tube, as shown in Fig. 13c. Compared with the smooth decreased as the pitch increased. Meanwhile, the secondary flow
tube and the straight grooved tube, the heat transfer performance velocity and the longitudinal vortex intensity increased as the spi-
of helically grooved tube greatly increased. Also, the heat transfer ral depth increased.
performance enhancement was much more significant compared Moreover, inserting a spoiler element in the tube is regarded as
with the resistance increase. a common method to enhance the heat transfer performance. Since
Additionally, the special-shaped tube has begun to be widely the spoiler element is easy to install and remove, this method has
used in the various heat exchange devices in recent years. Meng been widely used in the field of heat transfer enhancement. The
et al. (M. J. A., L. X. G., C. Z. J., 2005) designed a new type of inserted spoiler element has a disturbing effect on the flow fluid
cross-zoomed elliptical tube, and the corresponding structure and promotes the fluid mixing. However, the turbulence can also
was shown in Fig. 14a. Meanwhile, the experimental results increase the flow resistance. At present, there are many types of
15
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
16
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Fig. 15. The heat transfer tube with inserted twisted tape (E. E, A. H, N. A, 2014; E.-A. S, W. K, E.-A. P, 2010; E.-A. S, T. C, E.-A. P, 2010; E.-A. S, Y. K, N. K, 2012; S. P, C. W, T. C,
2018; B.M.M. K, C.M.S. U, S. M, 2013).
hensive performance of interpolated double helical coil tube, as immersion time included 0–600s to form different zirconium oxide
shown in Fig. 16c. It was found that the interpolated double helical microscopic surfaces. Then, the samples were baked at the temper-
coil tube had the increased heat transfer performance of 74%-178% ature of 300 °C for 6 h to remove the fluoride. Also, the scanning
and the increased resistance of 283–879% compared to the smooth electron microscope presented the significant micro-structure vari-
tube. According to the above discussions, it can be seen that ation, and the surface with small contact angle had the structural
although the above-mentioned heat transfer tubes apply the characteristics similar to the porous material. After that, the pool
macro-structural innovation technology on both sides of the tube boiling was carried out for CHF experiment. The results showed
side and the shell side, most of they increase the heat transfer area, that: as for the low contact angle (0°-10°), the capillary action could
enhance the enhance turbulence intensity, or disturb the near-wall lead to the significant improvement of droplet diffusion perfor-
boundary layer to enhance the heat transfer performance. mance, and could effectively delay the occurrence of CHF. In
2018, Haas et al. (Haas et al., 2018) treated the surface of Zr-4 alloy
3.2. Surface modification technology tube by different manufacturing processes, as shown in Fig. 17. The
morphological characteristics were obtained by using the physical
The surface modification refers to the use of chemical treat- vapor deposition, the electroplating, the mechanical rolling, the
ment, mechanical processing, coating treatment, and micro-nano wet powder sintering, and the oxidation treatment. After that, the
electroplating processing to process the surface of heat transfer low-pressure forced convection boiling heat transfer was con-
element, which can presents different wettability, roughness, or ducted for each micro-nano structure zirconium alloy tube. It was
nucleation core density. Also, the existing studies show that the shown that the increase of CHF was not only related to the surface
surface modification technology can significantly improve the heat micro-structure, but also affected by the thermal parameters such
transfer coefficient and CHF. Meanwhile, the micro-nano structure as the system pressure and the flow rate. Compared with the
technology and the surface coating technology are the hot research smooth surface before treatment, the variation range of CHF for dif-
focuses in the field of surface modification technology. ferent surface micro-nano structures was about 13% to 29%. More-
over, it was believed that the bubble dynamic behaviors led to the
3.2.1. Micro-nano structure technology above phenomena (including slipping, departure, and coalescence),
Micro-nano structure is a general term for micro-scale and which were mainly affected by material properties and surface
nano-scale structures, which mainly includes: the micro array hole, modification processes. In addition, Fong et al. (Fong et al., 2001)
the micro array pillar, the carbon nano-tube, the metal oxide nano- used the glass particles (60–180 lm) to impact the surface of zirco-
wire, the silicon nano-wire, and the nano-porous material. Ahn nium alloy tube and then applied the oxidation treatment to it.
et al. (Ahn, 2010) used the anodic oxidation treatment method to Thus, the micro-nano structure surfaces with different roughness
change the surface characteristics of Zr-4 alloy. Firstly, the sample were obtained. The pool boiling experimental results showed that:
17
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Porous Oxidized
Fig. 17. The Micro-nano structure surfaces with different morphology characteristics (Haas et al., 2018).
ing to prepare the aluminum-based micro-channel with super- transfer performance. Firstly, it can effectively ensure the struc-
hydrophilic and super-hydrophobic surfaces, as shown in tural integrity of original heating surface, which mainly include
Fig. 23a. As shown in Fig. 23b and Fig. 23c, the super- as follows: the metal oxide, the nano-particle, SiO2, carbon nano-
hydrophobic surface required less energy to generate the bubbles, tube, and the graphene. Also, the coating methods generally
and the nucleation core was easier to be activated. Therefore, the include: the vapor deposition, the atomic layer deposition, the
wall superheat required for ONB (Onset of Nucleate Boiling) was splash coating, the nano-fluid boiling, the calcination, and the arc
the smallest among these three surfaces, which was only 2.86 °C. ion plating.
As the heat flux was lower than 16 kW/m2, the boiling heat transfer Kim et al. (Kim et al., 2017) achieved the hydrophobic proper-
performance of super-hydrophobic surface was the best, and the ties of heating surface according to the PTFE (Polytetrafluo-
average heat transfer coefficient was 1.15 ~ 1.32 times that of roethylene) coating. Takata et al. (Takata et al., 2006) obtained
the smooth surface. As the heat flux was higher than 25 kW/m2, the super-hydrophobicity of heating surface by nickel spraying
the generated bubbles were easy to coalesce, which prevented technology and PTFE particle. Phan et al. (Phan et al., 2009) pre-
the boiling heat transfer process, so as to resulting in the signifi- pared the hydrophilic properties by adding the silicon oxide coat-
cant reduction of thermal coefficient. As for the super- ing to the heating surface. Kim et al. (kim, 0000) made the super-
hydrophilic surface with the contact angle close to 0°, the bubble hydrophilicity of heating surface by the titanium dioxide coating.
diameter decreased, and the departure frequency increased. Thus, According to the research of Takata et al. (Takata et al., 2006), it
it could be concluded that the super-hydrophilic surface could pro- was suggested that as the surface roughness of coating reached
vide the great capillary suction, so that the heating wall could be 10 lm or less, the effect of surface roughness on boiling heat trans-
replenished in time by the liquid after the bubble departure hap- fer was much weak. Also, Hanley et al. (OHanley et al., 2013)
pening, finally the boiling heat transfer process was correspond- applied the ITO (Indium Tin Oxide) to make the porous and hydro-
ingly strengthened. philic heating surface. The results showed that the porous nature of
Dong et al. (Dong, 2014) fabricated the micro-nano structure ITO coating provided more nucleation cores. Meanwhile, the fluid
surface of micro-pillar, micro-hole, nano-pillar, and nano-pore could be transported between the nucleation core under the capil-
according to MEMS technology, as shown in Fig. 24. The results lary action, thereby enhancing the liquid rewetting effect after the
of boiling heat transfer experiment showed that the micro- bubbles departure. Wang and Dhir (Wang and Dhir, 1993) found
structure surface had the higher nucleation core compared with that the heating surface coated with TiO2 would exhibit the
the smooth surface at low heat flux. Also, the diameter decrease super-hydrophilic properties under the ultraviolet light irradiation.
of micro-cavity would increase the nucleation core, while the sur- Then, Takata et al. (Takata et al., 2006) studied the droplet evapo-
face of nano-scale structure had unobvious variation. In addition, ration on the super-hydrophilic surface. It was shown that the
the surface of nano-scale structure could increase the bubble super-hydrophilic surface coated with TiO2 had much higher CHF
departure frequency and reduce the bubble departure diameter, than the heating surface without coating anything. Moreover, the
which was mainly attributed to the capillary effect of nano- Leidenfrost temperature (namely the transition temperature corre-
structure. Under the high heat flux, the nano-structure could pre- sponding to the transitional boiling region and the film boiling
vent the formation of vapor film by delaying the bubble coales- region) could also be greatly increased. Later, Takata et al.
cence, which thereby greatly increased CHF. (Takata et al., 2006) electroplated the coating composed of nickel
and PTFE particles on the copper surface, and achieved the
3.2.2. Surface coating technology super-hydrophobic heating surface. The experimental results show
In the past few decades, many studies have used different mate- that the pool boiling characteristics of super-hydrophobic surface
rials to coat the heat transfer surface for improving the boiling heat were much different from the conventional surface. Firstly, the
19
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
was less than 5 nm, and the bubble nucleation always occurred
at the junction of hydrophilic and hydrophobic regions, as shown
in Fig. 25. Moreover, Betz et al. (Betz et al., 2013) conducted the
pool boiling experiment based on the hydrophilic silicon oxide
wafer and the hydrophobic silicon oxide wafer. Obviously, the
mixed hydrophilic-hydrophobic surface presented the best pool
boiling heat transfer performance, including the enhancement of
pool boiling heat transfer coefficient and CHF. Suroto et al.
(Suroto et al., 2013) also researched the pool boiling heat transfer
performance of the mixed hydrophilic-hydrophobic surface, and
obtained the similar research results as Betz et al. (Betz et al.,
2013). Meanwhile, Jo et al. (Jo et al., 2011) further investigated
the pool boiling heat transfer on the mixed hydrophilic-
hydrophobic silicon surface, and it was proposed that the
hydrophobic points on the hydrophilic surface were in favor of
the boiling heat transfer enhancement with CHF varying little.
According to the above discussions, it can be seen that whether
a) Surface structure characteristics
it is the micro-nano structure technology or the surface coating
technology, the contribution to boiling heat transfer enhancement
is mainly attributed to the hydrophilicity and hydrophobicity vari-
ation of heating surface, which can lead to the variation of bubble
dynamic behaviors. As the bubble is attached to the heating sur-
face, the bubble boundary attached to the heating wall is called
the vapor–liquid-solid three-phase contact line. Also, the wettabil-
ity of heating surface, namely the hydrophilicity or hydrophobicity,
can be measured by the contact angle between the bubble and the
heating surface (Gong, 2015). As the contact angle h is less than
90°, the heating surface is called the hydrophilic surface; on the
contrary (h > 90°), it is called the hydrophobic surface, as shown
in Fig. 26. As the contact angle is close to 0°, the heating surface
is called the super-hydrophilic surface; as the contact angle is
greater than 150°, it is called the super-hydrophobic surface
(Betz et al., 2013). Based on the previous analyses, the following
conclusions can be drawn: the hydrophilic surface is conducive
to the increase of CHF, while the hydrophobic surface can promote
the bubble nucleation and has the higher boiling heat transfer per-
formance in the low superheat region. Moreover, Hanley et al.
(OHanley et al., 2013) proposed that since the hydrophilic surface
could promote the contact area rewetting, thus CHF was signifi-
cantly improved. Kousalya et al. (Kousalya et al., 2015) found that
the nucleation core of super-hydrophobic surface had the stronger
vapor phase capture ability, so the density of nucleation point
could be greatly increased.
However, the long-lasting stability and corrosion resistance of
micro-nano surface structure under high temperature, high pres-
sure, and high flow rate coolant scouring need to be further verified
combined with the actual engineering background, especially the
corrosion product deposition and the hydrochemical effect in the
nuclear reactor. Since the thickness of surface coating is generally
only micrometers, it needs to focus on the structure stability in the
harsh fluid environment such as the thermal stress, the chemical
corrosion, and the high pressure. According to the research founda-
tion of authors, it is believed that how to improve the bonding
force of surface coating and reduce the thermal-mechanical stress
b) Boiling heat transfer curve between different materials should be paid much attention in
Fig. 18. The boiling heat transfer characteristics of different surface micro-nano
future for the engineering application.
structures (Zheng, 2017).
Fig. 19. The boiling heat transfer characteristics of different surface structures (Zheng, 2017).
fluid (E. E, C. R. G, F. S, R. S. A, M. N, 2016) (K. E, R. M, A. N, 2018; B. the convection effect, which in turn increased the convective heat
M, H. M, 2015). transfer coefficient. Moreover, Asfer et al. (A. M, M. B, K. A, K. S, P.
PK, 2016) researched the heat transfer characteristics of magnetic
(1) Single phase flow heat transfer fluid under the condition of laminar flow in the microchannel
(with the diameter of 2 mm) . It was suggested that the convective
Sundar et al. (S. L. S, S. M. K, S. A. C. M, 2013) conducted the sys- heat transfer coefficient of magnetic fluid under the action of mag-
tematic thermal conductivity study based on the ethylene glycol netic field was likely to increase or decrease, which was dependent
aqueous solution + Fe3O4 nano-particle magnetic fluid. The work- on many factors, mainly including the ratio of magnetic force to
ing conditions included different mass ratios of ethylene glycol inertial force, the interaction between the magnetic fluid and the
and water (20%:80%, 40%:60%, 60%:40%), different temperatures magnetic particle, and the effect of magnetic particle chain-like
(20 ~ 60 ℃), different nano-particle concentrations (0.2 ~ 2.0 vol cluster, as shown in Fig. 27.
%). The results showed that the thermal conductivity of magnetic Then, Hejazian and Nguyen (H. M, N.-T. N, 2017) proposed a
fluid increased with the concentration and the temperature. As miniature magnetic fluid device. This device contained a circular
the concentration was 2.0 vol%, the ratio of ethylene glycol to cavity (the channel depth was 500 lm, and the channel width
water was 20%:80%, the temperature was 60 °C, the thermal con- was 2 mm) with a heat source on its top, which could be used
ductivity of magnetic fluid could be increased by about 46%. In for the magnetic field heat transfer enhancement according to
addition, Newcastle University and MIT had cooperated for the the external permanent magnet. As for three cases of the ionized
effects of magnetic field intensity and uniformity on the heat trans- water, the magnetic fluid (magnetic NdFeB nano-particles + deio
fer performance under the condition of laminar flow (A. R, D. E, M. nized water, without magnetic field) and the magnetic fluid (with
T, B. J, H. LW, M. B, 2014). Also, it was shown that the local heat magnetic field), the introduction of magnetic nano-particles could
transfer coefficient increased with the increase of magnetic field improve the heat transfer performance without the magnetic field,
intensity and gradient. Meanwhile, the enhancement degree of and the enhancement degree was proportional to Re number.
convective heat transfer coefficient increased with the increase of Moreover, the addition of magnetic field was not good for the heat
Re number. Zhang et al (Z. X, H. H, W. H, 2014) numerically studied transfer at low Re number because of the magnetic secondary flow.
the magnetic turbulence with the low Re number of 600. The However, the adverse effect gradually disappeared with the
results showed that the interaction in wake area could lead to increase of Re number. Meanwhile, Kaneda et al. (K. M, T. A, O.
the variety of thermo-hydrodynamic forces under the condition H, S. K, 2015) investigated the heat transfer performance of laminar
of low separation ratio (the ratio of circumfluence distance to flow magnetic fluid (Re = 10) in the horizontal channel by the
width). However, this phenomenon did not occur under the condi- numerical simulation. It was found that the position of electromag-
tion of high separation ratio. Also, the magnetic field could make netic coil could significantly affect the flow and heat transfer effect,
21
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Fig. 20. Micro-nano surface structure characteristics with femtosecond laser technology (Liu et al., 2019a, 2019b).
which could lead to the strengthening or weakening of local heat the local convective heat transfer coefficient could significantly
transfer performance. It should be noted that the strengthening enhanced as the upper edge of magnet was located at the lower
effect could only be achieved as the coil was located at the edge edge of heater. Also, the heat transfer enhancement degree would
of heating region. Also, the heat transfer enhancement was mainly weaken as the magnet position moved up, which was mainly
attributed to that the magnetic-thermal effect could make the cold attributed to the combined effect of buoyancy, gravity and mag-
fluid close to the electromagnetic coil, thereby thinning the bound- netic force. As for two magnets, the increase of local convective
ary layer. Under the condition of no gravity field, the heat transfer heat transfer coefficient was similar to the single magnet working
coefficient could increases by 10%, and the addition of gravity field condition
would increase the thermal stratification effect and further thin the
boundary layer. (2) Boiling heat transfer
Zonouzi et al. (Z. S. A, K. R, S. H, A. H, T. Y, M. M, 2018)
researched the forced convection heat transfer performance of Fe3- Naphon et al. (N. P., 2015) used the magnetic fluid of TiO2 +-
O4 + magnetic fluid in the vertical channel, and it was shown that R141b to conduct the pool boiling experiment, and the main focus
the local heat transfer coefficient increased by 23.4–48.9% while was on the factors such as the magnetic field intensity, the particle
the pressure drop just increased by about 1%. Since the magnetic concentration, the boiling pressure, and the bubble characteristics.
field could accelerate the magnetic fluid entering the region near The results showed that the boiling heat transfer coefficient was
heating wall. Also, Wang studied (Wang et al., 2019) the effect of proportional to the particle concentration as there was no mag-
magnetic field distribution on the heat transfer performance with netic field. As the magnetic field existed, the boiling heat transfer
low Re number. It was suggested that the number, the spacing, coefficient was proportional to the particle concentration and the
and the distribution symmetry of multiple magnets could affect magnetic field intensity. Moreover, Li et al. (L. S-Y., J. W-T., Z. C-
the heat transfer performance. As the magnets were adjacently Y., Z. H., T. W-Q., 2019) used the magnetic fluid of a-Fe2O33c-
arranged, the heat transfer performance was weakened. As the Fe2O3 + water (with the concentration of 0.005 ~ 0.1 g/L) to study
spacing was increased and the symmetry was maintained, the heat the heat transfer enhancement. It was shown that the magnetic
transfer performance could be enhanced. Since the the spacing field could change the concentration of local magnetic fluid and
increase could reduce the degree of particle aggregation and was conducive to the bubble generation, so as to enhance the boil-
improve the convective heat transfer performance, as shown in ing heat transfer. Also, the strengthening effect was not only
Fig. 28. Kaneda et al. (K. M, N. K, F. H, W. K, S. K, 2019) investigated dependent on the intensity and the distribution of magnetic field,
the natural convection heat transfer characteristics of magnetic but also affected by the magnetic particle magnetism and the con-
fluid in the partially vertical channels. The results showed that centration of magnetic fluid.
22
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Fig. 21. The micro-nano surface structure characteristics (Zhang et al., 2018).
Fig. 22. The micro-nano surface structure with etching process (Zhou et al., 2019).
23
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
c) The variation characteristics of boiling heat transfer coefficient with heat flux
Fig. 23. The flow boiling heat transfer characteristics of different micro-nano surface structures (Wu, 2017).
Then, Aminfar et al. (A. H, M. M, M. R, 2014) investigated CHF of 2019) developed a new type of pulsating heat pipe based on the
Fe3O4 + water magnetic fluid in the vertical channel, and it was magnetic fluid (CoFe2O4 + water), whose heat transfer performance
suggested that the magnetic fluid had a positive effect on CHF. could be increased by 71%.
Since the particles were deposited on the wall surface to improve
the surface wettability, and the magnetic field generated by quad- 4.2. Electric field enhanced heat transfer technology
rupole magnet had the improvement effect on CHF. In general, the
magnetic effects included the nature variation of working fluid, the As another one active heat transfer enhancement technology,
unidirectional convection, the nucleate boiling suppression, the the electric field heat transfer enhancement was a method which
boiling state stabilization, as shown in Fig. 29. Subsequently, Amin- used the coupling phenomena between the electric field, the flow
far et al. (A. H, M. M, M. R, 2014) further studied the influence of field and the temperature field to enhance the heat transfer perfor-
magnetic field on the bubble departure diameter of magnetic field. mance. As early as 1916, Chubb (Chen, 2016) found that the elec-
As there was the magnetic field, the bubble departure diameter tric field could enhance the boiling heat transfer performance
was reduced by 5–10% compared with the absence of magnetic and increase the steam generation amount. In the 1950s, Kronig
field. Monroe et al. (M. J. G, K. S, F. J. D, W. K. B, B. M. J, T. S. M, et al. (Lin et al., 2007) suggested that the electric field could
24
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
a) Hydrophilic surface
for the working fluid to form the convection, and the electric field
force could be established to form the convection heat transfer. 5)
The low energy consumption. As for the electric field enhancement
heat transfer, the generated current was relatively weak, just at the
level of lA. It could be seen that the electric field heat transfer
enhancement technology had much large engineering application
prospect, such as the thermal control application of spacecraft.
The numerous theoretical analyses and experimental studies
suggested that although the electric field could increase the heat
transfer between molecules, the movement of single-phase dielec-
tric liquid under the electric field should be attributed to the elec-
tric convection effect (Chen, 2016) (B. E., V. J., W. L., 1962). As early
as the 1960s, Bonjour (B. E., V. J., W. L., 1962) proposed that the
electric conductivity could strongly affect the natural convection
heat transfer of organic dielectric. Later, Panofsky (Chen, 2016)
derived the equation of electric field force
h i
(F e ¼ qe E 2 E re þ 2 r E @ q q ), which laid the theoretical foun-
1 2 1 2 @e
dation for the electric field heat transfer enhancement. As for the
above equation, the first term on the right was the electrophoretic
force or the coulomb force; the second term was the dielec-
trophoretic force; the third term was the electro-induced contrac-
tion force. Moreover, these three electric field forces could affect
Fig. 25. The bubble behavior characteristics of mixed hydrophilic- hydrophobic
surface (Betz et al., 2010). the flow phenomena for the natural convection heat transfer of
single-phase fluid. If there was a temperature gradient for the
working fluid, the conductivity and dielectric constant of working
increase the convective heat transfer coefficient of insulating fluid would change with temperature, and the working fluid would
dielectric liquid. By now, the electric field heat transfer enhance- be subjected to the electrophoretic force as well as the dielec-
ment technology has attracted the much attention of scholars. trophoretic force. If the unevenness of electric field distribution
The literatures shown that the electric field heat transfer enhance- was significant, the working liquid element would move from
ment technology had the following advantages: 1) the heat trans- the location with low electric field intensity to that with high elec-
fer effect enhancement was significant. The electric field jet tric field intensity. Also, this movement speed was much larger
method could increase the heat transfer coefficient by about 100 than the movement speed of liquid during the natural convection
times, and the heat transfer enhancement effect of small tempera- process and was not affected by the temperature. It should be
ture difference was more obvious. 2) The electric field generator noted that the above phenomenon was called liquid jet under
was simple. There was just a high-voltage conversion device and the electric filed. Fig. 30 shown the liquid jet as the highly uneven
the additional electrodes required. 3) The heat flux and the wall electric field formed by the ring electrode. Up to now, the work of
temperature were easy to control. By means of controlling the volt- electric field heat transfer enhancement mainly includes two direc-
age, the local temperature and heat flux could be controlled, which tions: 1) one is the natural convection of of pure working fluid,
could also realize the automatic and remote control. 4) The electric including the heat transfer experiment and the dimensionless cri-
field could be applied to some extreme situations where the grav- terion equation; 2) the other is the phenomenon of single-phase
ity did not exist or was much weak. In this situation, it was difficult liquid jet under the action of high electric field.
25
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
Fig. 27. The phenomenon of magnetic particle chain cluster (A. M, M. B, K. A, K. S, P. PK, 2016).
According to the relevant references (Y. K., Yabe, 1997; was greater than that in the low electric field area. Therefore, the
SunadaK., Y. A., T. T., 1991), it could be seen that the electric field direction of resultant force pointed to the high electric field area,
could enhance the heat transfer performance of many medium, and the single-phase liquid would move to the high electric field
including R-123, R-113 and other organic substances. Didkovsky area. It should be noted that this formed force due to the non-
and Bologa (Didkovsky, 1981) investigated three types of working uniform electric field was called electrostrictive force. Meanwhile,
fluids (including ether, ethane and Freon), and it was found that Yabe (Y. K., Yabe, 1997) solved the Laplace equation, and obtained
the electric field had little effect on ethane and Freon, while the the potential distribution between the ring electrode and the plate.
electric field had a greater increase for the heat transfer of ether. It was concluded that the electric field intensity outside the ring
In addition, Yamashita et al. (Y. K., Yabe, 1997) studied the effect electrode was greater than the electric field intensity between
of electric field on the heat transfer of R123, R114 and C6F144. It them. From the perspective of charge relaxation time, since the
was suggested that the heat transfer performance of C6F144 under charge relaxation time was relatively small, the Coulomb force
the electric field was increased by nearly 4 times, and the heat could be considered as approximately zero. In the absence of tem-
transfer performance of R114 and R123 was nearly 6 times com- perature gradient, the electric field force of liquid dielectric was
pared with the heat transfer without electric field. Later, Paschke- just the electrostrictive force. Finally, Yabe (Y. K., Yabe, 1997) pro-
witz and Pratt (P. J. S, P. D. M, 2000) researched the influence of posed that the liquid jet was mainly caused by the electrostrictive
fluid properties on the electric flied heat transfer enhancement. force.
Considering the heat transfer, the pressure drop, the electric power Based on the electric field heat transfer enhancement technol-
demand, and the transition point where viscosity and electric field ogy discussed above, Chen (Chen, 2016) combined it with the
dominated, three cooling oils with different physical properties nano-fluid technology, thus a new technical route of heat transfer
were studied. The results showed that under the given electric enhancement was proposed. It was proposed that the nanoparti-
power input, the working fluid with low viscosity and low conduc- cles with electric charge would produce the directional migration
tivity had the greatest heat transfer enhancement performance. As under the action of electric field, which would move to the liquid
for the working fluid with small charge relaxation time (defined as phase, and this phenomenon was called electrophoresis. Therefore,
the ratio of dielectric constant to conductivity), it was necessary to the nano-fluid particles would be subjected to electric field force
obtain a given heat transfer coefficient with more electric power. and resulted in the migration phenomenon. Since the Brownian
Moreover, Wang et al. (Wang et al., 2005) believed that under motion of nanoparticles could not be ignored, the irregular motion
the high-voltage electric field, the heat transfer characteristics of caused by Brownian force could increase the collision probability
working fluid could be enhanced, and the intensityening effect between particles, which could also make the micro-convection
could gradually increase with the increase of voltage until the between the nanoparticle and the fluid, so as to achieve the heat
maximum value reached. After that, even if the voltage continued transfer enhancement. According to the above analyses, the
to increase, the intensityening effect could no longer change. nanoparticles would be subjected to the electrophoresis force
Meanwhile, Liu et al. (L. Z. H, Y. J, 2004) found that electric field and the dielectrophoresis force under the high electric field. More-
heat transfer enhancement had a very excellent heat transfer over, these forces would strengthen the movement of nanoparti-
enhancement effect on the weakly polar dielectric, such as the cles and increase the collision as well as heat exchange between
insulating oil. As for the electric field between the ring electrode particles. Jung et al. (M. J. Y., J. S. P., C. S., U. S., 2005) studied the
and the plate, Yabe (Y. K., Yabe, 1997) solved the Navier-Stokes dynamic characteristics of nanoparticles (Al2O3) under the electric
equation containing the electric field force, and determined the field, which was generated by the alternating current. In the force
velocity distribution and the influence of liquid jet on the natural model, the electric field force considered the electrostatic force and
convection heat transfer. First of all, the single-phase liquid work- the dielectrophoretic force, which also included the viscous force
ing fluid under the electric field would be polarized, but the total and the Brownian force of particles. The results showed that the
polarization charge of working fluid was zero. However, the elec- average speed of nanoparticles decreased with the increase of par-
tric field force of polarized charge in the high electric field area ticle size, and increased with the increase of electric potential. In
26
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
a) Experiment process
Fig. 29. The boiling heat transfer characteristics of magnetic fluid (A. H, M. M, M. R, 2014).
27
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
5. Comprehensive evaluation
Table 1
The comprehensive evaluation of different heat transfer enhancement technologies.
Evaluating Micro-channel Longitudinal Micro-nano Surface coating Magnetic field Electric field
indicator technology vortex technology structure technology technology technology
technology
Single-phase heat www www w w ww ww
transfer
performance
Two-phase heat www Affected by wwwAffected by www www ww Dependent on ww
transfer structural structural Magnetofluid &
performance magnetic field
Dependent on
electric field
Application in ww ww w w
reactor core
Application in www www ww ww w w
non-nuclear
environment
Manufacturing www Mature & www Mature & ww Mature but ww
maturity Diversified Diversifie limited size
Mature but w
limited size
Scientific w
research stage
Scientific
research stage
Technology Microelectronics Various fields New energy New energy Scientific research
maturity stage
【Space field】 Scientific research
stage
【Space field】
Technology Impurities block Heat transfer Corrosion product Corrosion product Dependent on dependent on fluid
bottleneck channel, resulting in performance is deposition and deposition and magnetic fluid and properties and electric
heat transfer significantly affected hydrochemical hydrochemical magnetic field control field control
deterioration by structure effect effect technology technology
Comprehensive Reactor core www www ww ww w w
evaluation
Non-nuclear www www www www ww ww
environment
28
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
29
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
(2) The future development trend of steam generator is smaller Lu, Q., Zhou, L., Shen, C., Liu, W., Liu, L., He, H., 2020. Experimental study on the
critical trigger mechanism of boiling caused by flow instability. J. Eng.
volume and higher heat transfer efficiency. Also, the heat
Thermophys. 41, 966–975.
transfer efficiency of waste heat removal system mainly L. Hu, Research on instability mechanism of boiling two-phase flow in parallel
depends on the heat carrying capacity of final heat sink. rectangular narrow channels coupled with cavitation evolution characteristics,
Meanwhile, the surface characteristics of fuel cladding and Chongqing University, 2017.
Liu, W., Peng, S., Jiang, G., Liu, Y., Du, S., Zhang, Y., Deng, J., 2020. Development and
the physical properties of coolant jointly determine the assessment of a new rod-bundle CHF correlation for China fuel assemblies. Ann.
emergency cooling capacity of safety injection system. The Nucl. Energy 138, 171–175.
reasonable design of external cooling channel can ensure Lu, Q., Liu, Y., Huang, H., Shen, C., He, H., Zhou, L., Chen, D., 2019a. Visual
investigation on the boiling crisis post the flow instability in the channel with
the sufficient cooling of pressure vessel, which is the most bypass. Ann. Nucl. Energy 131, 171–184.
direct link for the successful implementation of ERVC. Lu, Q., Chen, D., Zhou, L., Li, Z., He, X., 2019b. Experimental investigation on the
(3) The heat transfer enhancement theory can effectively guide distribution characteristics of bubble contact areas in the narrow channel. Int. J.
Heat Mass Transf. 137, 876–881.
the design and application of heat transfer enhancement Yan, C., Cao, X., 2010. Safe heat transfer in nuclear peactor, Harbin Engineering.
technology, including: the field synergy theory, the configu- University Journal.
ration theory, the entropy minimization theory, the entransy V. J., B. R. M., S. R, 2009. Safety design features of the KLT-40S. International Atomic
Energy Agency 6, 94–118.
dissipation theory, and the process enhancement S. G. D., M. D. Carelli, 2008. The IRIS (International Reactor Innovative and Secure)
technology. and its unique advancement of the state-of the art. Power and Energy Society
(4) The micro-channel technology and longitudinal vortex tech- General 75, 1–5.
Zhang, X., 2012. Research on the characteristics of the passive waste heat removal
nology are the structural innovation technologies, both of
system of advanced reactors. Harbin Engineering University.
which belong to the passive heat transfer enhancement Fei, Y., Gu, C., 2016. Hualong-1 active and passive safety system introduction.
technology. This type of technology can effectively improve Equipment application and research 12, 56–57.
the heat transfer performance under the single-phase and Men, Q., 2015. Research on heat transfer mechanism and calculation method of
passive waste heat removal heat exchanger. East China University of Science
two-phase conditions, but which depends on the optimiza- and Technology.
tion of structural parameters. In addition, the manufacturing Kong, F., 2014. Dynamic reliability analysis of AP1000 safety injection system based
of this technology is mature, which has the broad applica- on Go-FLow. Harbin Engineering University.
Guo, R., 2015. Experiments on the Critical Heat Flux of the Outer Wall of the Lower
tion prospects in both of reactor core and non-nuclear fields. Head of a Pressure Vessel and Research on the Mechanism Prediction Model
(5) The surface modification technology includes the micro- Shanghai. Jiaotong University.
nano structure technology and the surface coating technol- Asmolov, V., Ponomarev-Stepnoy, N., Strizhov, V., 2001. Challenges left in the area
of in-vessel melt retention. Nucl. Eng. Des. 209, 87–91.
ogy, both of which belong to the passive enhancement heat Guo, Z., Zhuang, W., 1992. Physical mechanism analysis of convective heat transfer
transfer technology. This type of technology can effectively and application. Acta thermophysics Sinica 13, 52–56.
improve the heat transfer performance under the two- Guo, Z.Y., Wen, W.Q., Shah, R.K., 2005. The field synergy (coordination) principle
and its applications in enhancing single phase convective heat transfer. Int. J.
phase condition. The current manufacturing technology Heat Mass Transf. 48, 1797–1807.
faces the size limitation. Moreover, this type of technology Tao, W.Q., He, Y.L., Wang, Q.W., Qu, Z.G., Song, F.Q., 2002. A unified analysis on
has the broad application prospect in the non-nuclear field, enhancing single phase convective heat transfer with field synergy principle.
Int. J. Heat Mass Transf. 45, 4871–4879.
but it is not suitable for the reactor core.
Liu, W., Liu, Z., Ming, T., Guo, Z., 2009. Physical quantity synergy in laminar flow
(6) Both of the magnetic field heat transfer enhancement tech- field and its application in heat transfer enhancement. Int. J. Heat Mass Transf.
nology and the electric field heat transfer enhancement 52, 4669–4672.
technology belong to the active enhanced heat transfer tech- Liu, W., Wang, J.B., Liu, Z.C., 2017. A method of fluid dynamic analysis based on
Navier-Stokes equation and conservation equation on fluid mechanical energy.
nology, which are dependent on the physical properties of Int. J. Heat Mass Transf. 109, 393–396.
fluid and the field control technology. The above two tech- Tao, W.Q., Guo, Z.Y., Wang, B.X., 2002. Field synergy principle for enhancing
nologies can effectively improve the heat transfer perfor- convective heat transfer––its extension and numerical verifications. Int. J. Heat
Mass Transf. 45, 3849–3856.
mance under the single-phase and two-phase conditions. Xiao, W., Jing, Q., 2016. Improvement in field synergy principle: more rigorous
However, these two technologies are still in the scientific application, better results. Int. J. Heat Mass Transf. 100, 347–354.
research stage, whose applicability in the reactor core is Xu, X.M., Li, R.Z., Fu, J.Q., Jiang, H.B., 2017. Research on the heat flow field synergy of
electric vehicle power cabin at different charge and discharge rates. Appl.
weak, and practicality in the non-nuclear environment is Therm. Eng. 117, 397–408.
also limited. Ou, J., Li, L., Tao, C., Chen, Z., 2014. Application of field synergy principle to analysis
of flow field in underhood of LPG bus. Comput. Fluids 103, 186–192.
Deng, Y., Cui, J., J. E, B. Zhang, X. Zhao, Z. Zhang, D. Han, 2007. Investigations on the
temperature distribution of the diesel particulate filter in the thermal
Declaration of Competing Interest regeneration process and its field synergy analysis. Appl. Therm. Eng. 123,
92–102.
B. A., K. J., 1988. Advanced Engineering Thermodynamics, 1st Entropy Generation
The authors declare that they have no known competing finan- through heat and fluid flow. Wiley, New York.
cial interests or personal relationships that could have appeared B. A., 1979. A study of entropy generation in fundamental convective heat transfer.
J. Heat Transfer 101, 718–725.
to influence the work reported in this paper.
Zhu, L., Yu, J., 2017. Optimization of heat sink of thermoelectric cooler using entropy
generation analysis. Int. J. Therm. Sci. 118, 168–175.
Acknowledgements Zhao, L., Gu, X., Wang, R., 2017. Entropy generation analysis and structure
optimization of rectangular finned elliptical tube heat exchanger. Therm. Sci.
Technol. 16, 337–344.
The authors are grateful for the support of the National Natural Vargas, J.V.C., Bejan, A., 2000. Thermodynamic optimization of the match between
Science Foundation of China (No. 11905212), the Applied Basic two streams with phase change. Energy 25, 15–33.
Sahiti, N., Kransniqi, F., Fejzullahu, Xh., Bunjaku, J., Muriqi, A., 2008. Entropy
Research Program of Sichuan Province (No. 2021YJ0514), and the generation minimization of a double-pipe pin fin heat exchanger. Appl. Therm.
China National Nuclear Corporation Science Fund for Talented Eng. 28, 2337–2344.
Young Scholars. Cui, G., Zhou, X., Zhao, Y., 2010. Channel arrangement analysis and design of multi
stream heat exchanger based on minimum entropy generation. Progress in
chemical industry 29, 642–645.
References Guo, Z.Y., Zhu, H.Y., Liang, X.G., 2007. Entransy, A physical quantity describing heat
transfer ability. Int. J. Heat Mass Transf. 50, 2545–2556.
Liu, W., Liu, Z.C., Jia, H., Fan, A.W., Nakayama, A., 2011. Entransy expression of the
Lu, Q., Zhang, Y., Liu, Y., Zhou, L., Shen, C., Chen, D., 2019. An experimental
second law of thermodynamics and its application to optimization in heat
investigation on the characteristics of flow instability with the evolution of two-
transfer process. Int. J. Heat Mass Transf. 54, 3049–3059.
phase interface morphology. Int. J. Heat Mass Transf. 138, 468–482.
30
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
W., J, L. X, 2008. Application of entransy dissipation extremum principle in radiative C. T. M., M. T. J., K. M. S., 2014. Evaluation of helical-cruciform fuel rod assemblies
heat transfer optimization. Sci. China Ser. E: Technol. Sci. 51, 1306–1314. for high-power-density LWRs. Nucl. Technol. 188, 139–153.
W. S, J. G, W. J, 2018. Application ofentransy-dissipation-based thermal resistance M. J, T. A, S. N, 2012. Lightbridge corporation’s advanced metallic fuel for light water
for performance optimization of spiral-wound heat exchanger. Int. J. Heat Mass reactors. Nucl. Technol. 180, 437–442.
Transf. 116, 743–750. Agbodemegbe, V.Y., Cheng, X., Akaho, E.H.K., Allotey, F.K.A., 2016. An investigation
Yuan, F., Chen, Q., 2011. A global optimization method for evaporative cooling of the effect of split-type mixing vane on extent of crossflow between
systems based on the entransy theory. Energy 42, 181–191. subchannels through the fuel rod gaps. Ann. Nucl. Energy 88, 174–185.
Jia, H., Liu, W., Liu, Z., 2012. Enhancing convective heat transfer based on minimum Krauss, T., Meyer, L., 1996. Characteristics of turbulent velocity and temperature in
power consumption principle. Chem. Eng. Sci. 69, 225–230. a wall channel of a heated rod bundle. Exp. Therm Fluid Sci. 12, 75–86.
C. Q., W. M., P. N., 2009. Optimization principle for variable viscosity fluid flow and Xiong, J., Qu, W., Zhang, T., Chai, X., Liu, X., Yang, Y., 2020. Experimental
its application to heavy oil flow drag reduction. Energy Fuels 23, 4470–4478. investigation on split-mixing-vane forced mixing in pressurized water reactor
Li, Q., Yuan, X., P. N., G. F., L. Luo, 2014. A novel optimization approach to convective fuel assembly. Ann. Nucl. Energy 143, 1074–1080.
heat transfer enhancement for solar receiver. Chem. Eng. Sci. 116, 806–816. Holloway, M.H.L., 2002. Development of swirling flow in a rod bundle subchannel. J.
M. J. A., L. X. G., L. Z. X., 2005. Field synergy optimization andenhanced heat transfer Fluids Eng. 124, 123–134.
by multi-longitudinal vortexes flow in tube. Int. J. Heat Mass Transf. 48, 3331– Shen, Y., Cao, Z., Lu, Q., 1991. An investigation of crossflow mixing effect caused by
3337. grid spacer with mixing blades in a rod bundle. Nucl. Eng. Des. 125, 111–119.
K. Guo, Transfer process enhancement method based on flow field structure I. K, M. H, 2006. Development of rod-embedded fiber LDV to measure velocity in
construction and its application, Tianjin University, 2015. fuel rod bundles. J. Nucl. Sci. Technol. 43, 150–158.
W. Zan, Research and prediction model of micro-scale enhanced heat transfer C. T. M., 2007. Thermal-hydraulic analysis of cross-shaped spiral fuel in high power
mechanism of boiling two-phase flow, Zhejiang University, 2013. density BWRs. Massachusetts Institute of Technology, Cambridge, MA.
Kandlikar, S.G., 2002. Fundamental issues related to flow boiling in minichannels Lee, J., Kim, J., Song, C., 2014. Synthesis of the turbulent mixing in a rod bundle with
and microchannels. Exp. Therm Fluid Sci. 26, 389–407. vaned spacer grids based on the OECD-KAERI CFD benchmark exercise. Nucl.
Mehendale, S.S., Jacobi, A.M., Shah, R.K., 2000. Fluid flow and heat transfe at micro- Eng. Des. 279, 3–18.
and meso-scales with application to heat exchanger design. Appl. Mech. Rev. 53, McKrell, C.T.M., Kazimi, T.J., M. S, 2013. Experimental investigation of hydraulics
175–193. and lateral mixing for helical-cruciform fuel rod assemblies. Nucl. Technol. 182,
Chen, L., Tian, Y.S., Karayiannis, T.G., 2006. The effect of tube diameter on vetical 259–273.
two-phase flow regimes in small tubes. Int. J. Heat Mass Transf. 49, 4220–4230. Koroush, S., 2016. Numerical investigation of the boiling crisis for helical cruciform-
Ribatski, G., Wojtan, L., Thome, J.R., 2006. An analysis of experimental data and shaped rods at high pressures. Int. J. Multiph. Flow 83, 51–61.
prediction methods for two-phase frictional pressure drop and flow boiling heat K. S, 2016. Numerical investigation of the boiling crisis for helical cruciform-shaped
transfer in micro-scale channels. Exp. Therm Fluid Sci. 31, 1–19. rods at high pressures. Int. J. Multiph. Flow 83, 51–61.
Thome, J.R., 2004. Review on two-phase flow instabilities in narrow spaces. Int. J. Yan, B.H., 2019. The thermal hydraulic phenomenon in tight lattice bundles: a
Heat Fluid Flow 25, 128–139. review. Ann. Nucl. Energy 126, 330–349.
Kew, P.A., Cornwell, K., 1997. Correlations for the prediction of boiling heat transfer B. V. V., B. S. M., K. L. L., 2007. Experimental study of burnout in channels with
in small-diameter channels. Appl. Therm. Eng. 17, 705–715. twisted fuel rod. J Therm. Eng. 54, 386–389.
Cheng, P., Wu, H.Y., 2006. Mesoscale and microscale phase-change heat transfer. Ma, J., Huang, Y., Huang, J., 2011. Single-phase water longitudinal vortex
Advances in Heat Transfer 39, 461–563. enhancement heat transfer experiment in narrow-gap rectangular channels.
Lu, Q., Chen, D., Wang, Q., 2015. Visual investigation on the interface morphology of Nucl. Power Eng. 32, 89–94.
Taylor bubble and the characteristics of two-phase flow in mini-channel. Chem. Huang, J., Huang, Y., Ma, J., 2011. Heat transfer characteristics of vapor-liquid two-
Eng. Sci. 134, 96–107. phase flow in narrow-gap rectangular channels with longitudinal vortex
Charnay, R., Revellin, R., Bonjour, J., 2015. Flow boiling heat transfer in minichannels generators. Science and Technology Annual Report of China Nuclear Power.
at high saturation temperatures: Part I - Experimental investigation and Research and Design Institute.
analysis of the heat transfer mechanisms. Int. J. Heat Mass Transf. 87, 636–652. Huang, J., Huang, Y., Ma, J., 2013. Longitudinal vortex enhanced vapor-liquid two-
Tran, T.N., Wambsganss, M.W., France, D.M., 1996. Small circular- and rctangular- phase heat transfer characteristics, China Nuclear Power Research and Design
channel boiling with two refrigerants. Int. J. Multiph. Flow 22, 485–498. Institute Science and Technology. Annual Report.
Vlasie, C., Macchi, H., Guilpart, J., Agostini, B., 2004. Flow boiling in small diameter Wang, L., Chen, Q., Zeng, M., Zhou, Y., Wang, Q., Huang, Y., Xiao, Z., 2007.
channels. Int. J. Refrig 27, 191–201. Experimental study of heat transfer enhancement in narrow rectangular
Collier, J.G., Thome, J.R., 1996. Convective Boiling and Condensation. Clarendon channel with longitudinal vortex generators. Nucl. Eng. Des. 237.
Press. R. L. W., 2009. Single-phase heat transfer, friction, and fouling characteristic of
Bao, Z.Y., Fletcher, D.F., Haynes, B.S., 2000. Flow boiling heat transfer of Freon R11 three-dimensional cone roughness in tube flow. Int. J. Heat Mass Transf. 52,
and HCFC123 in narrow passages. Int. J. Heat Mass Transf. 43, 3347–3358. 2624–2631.
Ali, R., Palm, B., Maqbool, M.H., 2011. Flow Boiling Heat Transfer Characteristics of a Xu, W., Wang, S., Zhang, Q., Wang, Q., Lu, H., Tan, H., 2016. Experimental and
Minichannel up to Dryout Condition. J. Heat Transfer 133, 1051–1053. numerical studies of heat transfer and friction factor of therminol liquid phase
Lin, S., Lin, P.A., Kew, P.A., Cornwell, K., 2001. Flow boiling of refrigerant R141B in heat transfer fluid in a ribbed tube. Appl. Therm. Eng. 95, 165–177.
small tubes. Chem. Eng. Res. Des. 79, 417–424. W. W., Z. Y., L. B., 2018. Numerical investigation of tube-side fully developed
Huo, X., Chen, L., Tian, Y.S., Karayiannis, T.G., 2004. Flow boiling and flow regimes in turbulent flow and heat transfer in outward corrugated tubes. International
small diameter tubes. Appl. Therm. Eng. 24, 1225–1239. Journal ofHeat and Mass Transfer 116, 115–126.
Karayiannis, T.G., Mahmoud, M.M., Kenning, D.B.R., 2012. A study of discrepancies S. M, Z. M, 2018. Investigation on turbulent flow and heat transfer characteristics
in flow boiling results in small to microdiameter metallic tubes. Exp. Therm and technical economy of corrugated tube. Appl. Therm. Eng. 129, 1–11.
Fluid Sci. 36, 126–142. P. S, M. S. F, M. S. G, 2016. Computational fluid dynamics simulation of heat
Sumith, B., Kaminaga, F., Matsumura, K., 2003. Saturated flow boiling of water in a enhancement in internally helical grooved tubes. Int. Commun. Heat Mass
vertical small diameter tube. Exp. Therm Fluid Sci. 27, 789–801. Transfer 73, 25–32.
K. K. E, D. T, W. M. W, 1997. Microscale flow visualization of nucleate boiling in B. K, C. M, G. H, 2009. The investigation of groove geometry effect on heat transfer
small channels: mechanisms influencing heat transfer, Proceeding of for internally grooved tubes. Appl. Therm. Eng. 29, 753–761.
international conference on compact heat exchanges for the process A. K, J. C, 2013. Heat transfer and single-phase flow in internally grooved tubes. Int.
industries. Begell House Inc, New York, pp. 343–352. Commun. Heat Mass Transfer 42, 62–68.
B. P, K.S. G, 2005. Experimental study of flow patterns, pressure drop, and flow M. J. A., L. X. G., C. Z. J., 2005. Experimental study on convective heat transfer in
instabilities in parallel rectangular minichannels. Heat Transfer Eng. 26, 20–27. alternating elliptical axis tubes. Exp. Therm Fluid Sci. 29, 457–465.
Owhaib, W., C. Martı́n-Callizo, B. Palm, 2004. Evaporative heat transfer in vertical S. A. R, Y. D. S, A. D, 2014. Experimental and numerical study on heat transfer and
circular microchannels. Appl. Therm. Eng. 24, 1241–1253. flow resistance of oil flow in alternating elliptical axis tubes. Int. J. Heat Mass
Yan, Y., Lin, T., 1998. Evaporation heat transfer and pressure drop of refrigerant R- Transf. 77, 124–130.
134a in a small pipe. Int. J. Heat Mass Transf. 41, 4183–4196. T. X. H, Z. D. S, Z. G. Y, 2013. Heat transfer and pressure drop performance of twisted
Shen, X., Miyazaki, K., Xu, J., 2001. Study on boiling heat transfer in a vertical oval tube heat exchanger. Appl. Therm. Eng. 50, 374–383.
annular narrow channel. Nucl. Sci. Eng. 21, 244–251. J. Z. J, C. F. Q, G. Z. X, 2017. Effects of pitch and corrugation depth on heat transfer
Sobierska, E., Kulenovic, R., Mertz, R., 2007. Heat transfer mechanism and flow characteristics in six-start spirally corrugated tube. Int. J. Heat Mass Transf. 108,
pattern during flow boiling of water in a vertical narrow channel—experimental 1011–1025.
results. Int. J. Therm. Sci. 46, 1172–1181. E. E, A. H, N. A, 2014. Study on heat transfer and friction factor characteritics of c-
L., M., C. He, 2002. Subcooled flow boiling heat transfer in vertical rectangular Al2O3-water through circular tube with twisted tape inserts with different
narrow gap. Therm. Sci. Technol. 1, 185–188. thicknesses. Int. J. Therm. Sci. 82, 72–83.
J. C. P., F. R. H., L. F. B. P, M. B., M. M, 2004. Confined boiling of FC72 and FC87 on a E.-A. S, W. K, E.-A. P, 2010. Heat transfer enhancement in a tube using delta-winglet
downward facing heating copper disk. Int. J. Heat Fluid Flow 25, 313–319. twisted tape inserts. Appl. Therm. Eng. 30, 310–318.
M. M, P. A. G. G, 2011. Nucleate boiling and critical heat flux of HFE7100 in E.-A. S, T. C, E.-A. P, 2010. Turbulent heat transfer enhancement by ounterrco-
horizontal narrow spaces. Exp. Therm Fluid Sci. 35, 772–779. swirling flow in a tube fitted with twin twisted tapes. Exp. Therm Fluid Sci. 34,
Xia, C., Hu, W., Guo, Z., 1996. Natural Convective Boiling in Vertical Rectangular 53–62.
Narrow Channels. Exp. Therm Fluid Sci. 12, 313–324. E.-A. S, Y. K, N. K, 2012. Heat transfer augmentation by helically twisted tapes as
McKrell, C.T.M., Kazimi, T.J., M. S, 2010. Assessment of helical-cruciform fuel rods swirl and turbulence promoters. Chem. Eng. Process. Process Intensif. 60, 42–
for high power density LWRs, MIT-NFC-TR-118. CANES Stewart. 48.
31
Q. Lu, Y. Liu, J. Deng et al. Annals of Nuclear Energy 159 (2021) 108302
S. P, C. W, T. C, 2018. Influence of regularly spaced quadruple twisted tape elements Kousalya, A.S., Singh, K.P., Fisher, T.S., 2015. Heterogeneous wetting surfaces with
on thermal enhancement characteristics. Chem. Eng. Process. Process Intensif. graphitic petal-decorated carbon nanotubes for enhanced flow boiling. Int. J.
128, 114–123. Heat Mass Transf. 87, 380–389.
B.M.M. K, C.M.S. U, S. M, 2013. Heat transfer and friction factor characteristics in E. E, C. R. G, F. S, R. S. A, M. N, 2016. Enhanced Heat Transfer Properties of Magnetite
turbulent flow through a tube fitted with perforated twisted tape inserts. Int. Nano-fluids due to Neel and Brownian Relaxation Mechanisms. Chem. Eng.
Commun. Heat Mass Transfer 46, 49–57. Commun. 203, 1157–1164.
E.-A. S, P. P, 2005. Enhancement of heat transfer in a tube with regularly-spaced K. E, R. M, A. N, 2018. Convective heat transfer enhancement in a pitted
helical tape swirl generators. Sol. Energy 78, 483–494. microchannel by stimulation of magnetic nanoparticles. Chem. Eng.
G. S, O. V, B. O, 2010. Heat transfer enhancement in a tube with equilateral triangle Processing-Process Intens. 126, 156–167.
cross sectioned coiled wire inserts. Exp. Therm Fluid Sci. 34, 684–691. Bahiraei, M., Hangi, M., 2015. Flow and heat transfer characteristics of magnetic
Z. C, W. D, Z. Y, 2015. Numerical study on heat transfer and flow characteristics of a nano-fluids: a review. J Magn Magn Mater 374, 125–138.
tube fitted with double spiral spring. Int. J. Therm. Sci. 94, 18–27. S. L. S, S. M. K, S. A. C. M, 2013. Thermal conductivity of ethylene glycol and water
Ahn, H.S., 2010. Pool boiling CHF enhancement by micro/nanoscale modification of mixture based Fe3O4 nano-fluid. Int Commun Heat Mass, International
zircaloy-4 surface. Nucl. Eng. Des. 240, 3350–3360. Communication. Heat Mass Transf. 49, 17–24.
Haas, C., Kaiser, F., Schulenberg, T., 2018. Critical heat flux for flow boiling of water A. R, D. E, M. T, B. J, H. LW, M. B, 2014. Effect of magnetic field on laminar convective
on micro-structured Zircaloy tube surfaces. International Journalof Heat and heat transfer of magnetite nano-fluids. Int. J. Heat Mass Transf. 68, 94–109.
Mass Transfer 120, 793–806. Z. X, H. H, W. H, 2014. Flow and heat transfer past row of magnetic obstacles for
Fong, R.W.L., Nitheanandan, T., Bullock, C.D., 2001. Correlation Between the Critical various separation ratios. Int. Commun. Heat Mass Transf. 59, 178–187.
Heat Flux and the Fractal Surface Roughness of Zirconium Alloy Tubes. A. M, M. B, K. A, K. S, P. PK, 2016. Effect of magnetic field on laminar convective heat
Enhanced Heat Transfer 8, 137–146. transfer characteristics of ferrofluid flowing through a circular stainless steel
L. C.Y., Z. B.J., K. K.J., 2012. Morphological change of plain and nano-porous surfaces tube. Int. J. Heat Fluid Flow 59, 74–86.
during boiling and its effect on nucleate pool boiling heat transfer, H. M, N.-T. N, 2017. Magnetofluidics for manipulation of convective heat transfer,
ExperimentalThermal and Fluid. Science 40, 150–158. International Communication. Heat Mass Transf. 81, 149–154.
Yu, C., Lu, D.C., Cheng, T., 2006. Pool boiling heat transfer onartificial micro-cavity K. M, T. A, O. H, S. K, 2015. Heat transfer enhancement by external magnetic field for
surfaces in dielectric fluid FC-72. J. Micromech. Microeng. 16, 2092–2099. paramagnetic laminar pipe flow. Int. J. Heat Mass Transf. 90, 388–395.
Zheng, X., 2017. Research on Boiling Heat Transfer Characteristics of Multi-scale and Z. S. A, K. R, S. H, A. H, T. Y, M. M, 2018. Experimental investigation of the flow and
Modified Structures. North China Electric Power University. heat transfer of magnetic nano-fluid in a vertical tube in the presence of
Liu, B., Liu, J., Zhou, J., 2019aa. Experimental study of subcooled boiling pool heat magnetic quadrupole field, Experimental Thermal Fluid. Science 91, 155–165.
transfer and its‘‘hookback”phenomenon on micro/nanostructured surfaces. Int. Wang, J., Li, G., Zhu, H., Luo, J., Sunden, B., 2019. Experimental investigation on
Commun. Heat Mass Transfer 100, 73–82. convective heat transfer of ferrofluids inside a pipe under various magnet
Liu, B., Liu, J., Zhang, Y., 2019bb. Experimental and theoretical study of pool boiling orientations. Int J Heat Mass Tran 132, 407–419.
heat transfer and itsCHF mechanism on femtosecond laser processed surfaces. K. M, N. K, F. H, W. K, S. K, 2019. Effect of magnetic field on natural convection inside
Int. J. Heat Mass Transf. 132, 259–270. a partially-heated vertical duct: Experimental study. Int. J. Heat Mass Transf.
Zhang, Y., Zhou, J., Zhou, W., 2018. CHF correlation of boiling in FC-72 with micro- 132, 1231–1238.
pin-fins for electronics cooling. Appl. Therm. Eng. 138, 494–500. N. P., 2015. Effect of Magnetic Fields on the Boiling Heat Transfer Characteristics of
Zhou, J., Liu, B., Qi, B., 2019. Experimental investigations of bubble behaviors and Nano-fluids. International Journal of Thermophys 36, 281–290.
heat transfer performance on micro nanostructure surfaces. Int. J. Therm. Sci. L. S-Y., J. W-T., Z. C-Y., Z. H., T. W-Q., 2019. Effects of magnetic field on the pool
135, 133–147. boiling heat transfer of water-based alpha-Fe2O3 and gamma-Fe2O3 nano-fluid.
Wu, D., 2017. Research on Flow Boiling Heat Transfer Characteristics in Rectangular Int. J. Heat Mass Transf. 128, 762–772.
Flow Channels with Extremely Wettable Micro-nano Structures. South China A. H, M. M, M. R, 2014. Experimental study on the effect of magnetic field on critical
University of Technology. heat flux of ferrofluid flow boiling in a vertical annulus, Experimental Thermal
L. Dong, XiaojunQuan, P. Cheng, An experimental investigation of enhanced pool Fluid. Science 58, 156–169.
boiling heat transferfrom surfaces with micro/nano-structures, International M. J. G, K. S, F. J. D, W. K. B, B. M. J, T. S. M, 2019. On the energy harvesting and heat
Journal of Heat and Mass Transfer, 71 (2014) 189-196. transfer ability of a ferro-nano-fluid oscillating heat pipe. Int. J. Heat Mass
Kim, S., Lee, G., Kang, J., Moriyama, K., Park, H., Kim, M., 2017. The role of surface Transf. 132, 162–171.
energy inheterogeneous bubble growth on idealsurface. Int. J. Heat Mass Transf. Chen, Y., 2016. Experimental and Numerical Simulation Research on the Heat
108, 1901–1909. Transfer Characteristics of Nano-fluids Enhanced by Electric Field. Shanghai
Takata, Y., Hidaka, S., Uraguchi, T., 2006. Boiling Feature on a SuperWater-Repellent Jiaotong University.
Surface. Heat Transfer Eng. 27, 25–30. Lin, Z., Wang, J., Li, R., 2007. Enhanced Heat Transfer Technology, Chemical Industry
Phan, H.T., Caney, N., Marty, P., Colasson, S., Gavillet, J., 2009. Surface wettability Press.
control by nanocoating: The effects on pool boiling heat transfer B. E., V. J., W. L., 1962. Electroconvection effects on heat transfer. Chem. Eng. Prog. 58.
ndnucleationmechanism. Int. J. Heat Mass Transf. 52, 5459–5471. Yabe, Y.K., 1997. Electrohydrodynamic enhancement of falling film evaporation
S. J. Kim, I. C. Bang, J. Buongiorno, L. W. Hu, Surface wettabilitychange during pool heat exchangers. Trans ASME, J Heat Transfer 119, 339–343.
boiling of nano-fluids and its effect on critical heat flux, International Journal of SunadaK., Y. A., T. T., Experimental study of EHD pseudo-drop-wise condensation,
Heat and Mass Transfer. Proc 3rd ASME/JSME Thermal Engineering Joint Conf, 1991, pp. 47-54.
H. OHanley, C. Coyle, J. Buongiorno, T. McKrell, L. Hu, M. Rubner, R. Cohen, Separate Didkovsky, BologaM, Int J heat MassTransfer, Vapor film condensation heat transfer
effects of surface roughness, wettability, and porosity on the boiling critical and hydrodynamics under the influence of an electric field, 24 (1981) 811~819.
heat flux, Applied Physics Letters, 103 (2013). P. J. S, P. D. M, 2000. Influence of fluid properties on electrohydrod-ynamic heat
Wang, C.H., Dhir, V.K., 1993. Effect of surface wettability on active nucleation site transfer enhancement in liquids under viscous and electrically do-minated flow
density during pool boiling of water on a vertical surface. J. Heat Transfer 115, conditions. Exp. Therm Fluid Sci. 21, 187–197.
659–669. Wang, F., Li, R., Yu, H., 2005. Experimental study on enhanced natural convection
Betz, A.R., Xu, J., Qiu, H., Attinger, D., 2010. Surfaces with mixed hydrophilic and and boiling heat transfer of benzene by external electric field. J. Solar Energy 26,
hydrophobic areas enhance pool boiling. Appl. Phys. Lett. 97, 1409–1419. 277–280.
A. R., Betz, J. Jenkins, C.-J., C. Kim, D. Attinger, Boiling heat transfer on L. Z. H, Y. J, 2004. Augmentation of laminar forced convective heat transfer of an oil
superhydrophilic, superhydrophobic, and superbiphilic surfaces, International flow in an enhanced tube by EHD effect. J. Heat Transfer 126, 121–138.
Journal of Heat and Mass Transfer, 13 (2013) 733-741. M. J. Y., J. S. P., C. S., U. S., 2005. Motion of nanoparticles in nano-fluids under an
Suroto, B.J., Tashiro, M., Hirabayashi, S., Hidaka, S., Kohno, M., Takata, Y., 2013. electric field. ASME 2005 International Mechanical Engineering Congress and
Effects of hydrophobic-spot periphery and subcooling on nucleate pool boiling Exposition, 497–551.
from a mixed-wettability surface. J. Therm. Sci. Technol. 8, 294–308. Liu, Y., Yu, H., Dong, W., 2008. Experimental study on heat transfer enhancement by
Jo, H., Ahn, H.S., Kang, S., Kim, M.H., 2011. A study of nucleate boiling heat transfer electric field using nano-fluid. J. Shanghai Univ. Technol. 30, 439–442.
on hydrophilic, hydrophobic and heterogeneous wetting surfaces. Int. J. Heat Jia, J., 2012. Research on enhanced heat transfer performance of nano-fluid
Mass Transf. 54, 5643–5652. electrophoresis. Yanshan University.
Gong, S., 2015. Lattice Boltzmann method for the hydrophobicity influence on the Huang, H., 2012. Research on Surface Modification and Electrophoresis Behavior of
pool boiling heat transfer. Shanghai Jiaotong University. Nanodiamond and Graphene. Yanshan University.
32