Abstract

High-fidelity flow boiling simulations are conducted in a vertical mini channel with offset strip fins (OSF) using R113 as a working fluid. Finite-element code PHASTA coupled with level set method for interface capturing is employed to model multiple sequential bubble nucleation using transient three-dimensional approach. The code performance is validated against experiments for a single nucleation site in a vertical rectangular channel. To assess code performance, a study on the bubble departure from the wall in a mini channel with OSF is carried out first. Contributions from the microlayer are not considered due to low heat flux values applied to the channel (1 kW/m2). The influence of surface characteristics, such as contact angle and liquid superheat on bubble dynamics, is also analyzed as well as the local two-phase heat transfer coefficient. For higher void fractions, two conical nucleation cavities are introduced in the same channel with OSF. Observed bubble characteristics (departure diameter, bubble departure frequency) are evaluated and bubble trajectories are presented and analyzed. The local heat transfer coefficient is then evaluated for each simulation. The results show approximately a 2.5 time increase in the local heat transfer coefficient when the individual bubbles approach the wall. With smaller bubble nucleation diameters, the heat transfer coefficient can increase by up to a factor of two. Thus, the current work demonstrates the flow modeling capability of the boiling phenomena in complex geometry with OSF.

1 Introduction

Forced convection combined with boiling permits effective heat removal from surfaces. High heat transfer coefficients are vital for flow boiling phenomenon especially in nuclear reactors but have relevant applications in steam generators, microchannels, compact heat exchangers, etc. In the latter examples, different fin shapes (rectangular, triangular, wavy, offset strip fin (OSF), etc.) can be added to increase heat removal from the surface. In addition, refrigerant coolant (i.e., R113) as a working fluid has also been used in conjunction with fins to increase heat transfer rate, especially in compact heat exchanges for low temperature/pressure applications [1,2]. However, accurately predicting the two-phase heat transfer coefficient of novel configurations can be challenging.

From an experimental perspective, several studies have been conducted to predict two-phase heat transfer coefficient in mini channels with OSF. In Ref. [3], R113 behavior is examined in a vertical channel with OSF over a wide range of Reynolds numbers and applied heat fluxes. The authors proposed a correlation to find the two-phase forced convective part of the local two-phase heat transfer coefficient. The experimental results showed that the heat transfer coefficient values may lie within ±25% of their predictions. More recently, in Ref. [4], an experimental study was performed with R134a in a brazed heat exchanger with vapor quality up to 0.7. A correlation for Reynonds number factor F was proposed (used for two-phase forced convective heat transfer coefficient evaluation) as well as two-phase frictional multiplier ϕf (applied in two-phase frictional pressure drop correlations). In Ref. [5], the boiling heat transfer characteristerics of HFE-7100 were investigated in a vertical channel with OSF and without them. The results showed that for low heat loads (less than 80 kW/m2), convective boiling was the dominant process in two-phase heat transfer for which OSF may have had significant impact. For higher heat loads (more than 120 kW/m2), the influence of inner channel geometry is negligible as nucleate boiling governs the two-phase heat transfer. In Ref. [6], the authors conducted an experimental study with serrated fins in a horizontal rectangular channel using R134a as the working fluid. They varied heat flux and captured flow patterns from dispersed bubbly to annular flow followed by intermittent dry-out. The results showed that for low vapor quality (less than 0.6) and under high heat fluxes, the mass flux influence on heat transfer coefficient is almost negligible. Conversely, for high vapor quality flow, the heat transfer coefficient increases with mass flux growth.

Due to the high computational costs, numerical investigations are mostly conducted to determine the heat transfer coefficient for single-phase flows. In Ref. [7], the SST kω turbulence model is used to predict Colburn factor j for a 3D geometry with OSF. The obtained results lie within ±20% range comparing to a well-known correlation proposed by Manglik and Bergles [8]. Later, in Ref. [9], the researchers proposed a new correlation for j that includes the Prandtl number. This was required as the existing correlations were limited to a working fluid with a low Pr (e.g., for air Pr= 0.7). In Ref. [10], it was discovered that the laminar flow heat transfer coefficient is mainly influenced by the growth of the boundary layer. In turbulent flow, vortex shedding also becomes meaningful. 3D simulations using Fluent were performed in Ref. [11] in which the standard kϵ model was used for “far from wall” calculations, and the SST kω turbulence model was applied for “near wall” calculations. The authors discovered a noticeable difference in results when compared to published work in Ref. [8]. Specifically, differences (∼15%) were discovered for the laminar regime (Re < 120) and a significant difference (∼35%) in the transitional (500 < Re < 1250) and turbulent (Re > 1250) regimes. In both Refs. [12] and [13], 3D simulations in Fluent were conducted and compared with the existing correlation [8]. More geometry parameters were analyzed in Ref. [12], a wider range of Re was investigated in Ref. [13].

To model boiling flows with interface capturing, researchers employ various interface-tracking methods (ITM). For example, a volume of fluid (VOF) based flow solver was proposed in Ref. [14] that reduces computational cost by keeping interfacial temperature at saturation and explicitly coupling it with source terms. This solver was validated, and its capability was demonstrated on 3D flow boiling simulation in a rectangular microchannel. In work done by Ref. [15], a new ITM is used that couples both VOF and level set methods (VOSET). Researchers implemented a microlayer model to take into consideration the evaporation in the microlayer underneath a bubble and conjugate heat transfer between the wall and fluid region. They investigated subcooled boiling of water in a vertical rectangular channel in 2D. To ensure multiple bubble nucleation, 15 nucleation cavities are introduced in the domain. Varying heat fluxes from 100 to 500 kW/m2 allowed them to observe flow regime transition from bubbly to elongated bubbly flow. In Ref. [16], the influence of surface characteristics on bubble dynamics is investigated in a vertical rectangular channel. Bubble departure from a solid surface is modeled in 2D. The authors concluded that a hydrophilic surface increased the heat transfer from the wall due to the evaporation of the microlayer. On amphiphilic surfaces, the waiting period (between the departure of one bubble and nucleation of a subsequent bubble) is essentially eliminated which results in an effective contact line heat transfer. Numerical simulations conducted alongside flow boiling experiments are presented in Ref. [17]. In that study, FC-72 is used as the working fluid in a vertical rectangular channel in 2D. The presented computational results showed good agreement with experimental data (including flow visualizations, cross section averaged void fraction, velocity profiles, wall temperature, and heat transfer coefficient). This prediction diverges from known solutions when the mass velocity increases, in part due to 2D modeling limitations. In Ref. [18], the authors investigated flow boiling of FC-72 in a vertical rectangular microchannel in 3D. The results revealed good agreement of numerical predictions for wall temperature variations and overall flow patterns with experimental results. This study also obtained fluid temperature distributions that cannot be measured in small channel experiments.

From the literature review it can be seen that due to high complexity of the flow boiling phenomenon, the majority of high-resolution simulations are conducted in 2D. Only a few cases are conducted in 3D, which impedes comprehensive understanding. In addition, of the studies published, most consider simple rectangular channel geometries which are not necessarily realistic. Thus, we present high fidelity 3D simulations to provide a deeper understanding of the boiling in a vertical mini channels with OSF. We use finite-element code PHASTA to solve Navier–Stokes equations (Sec. 2). The code performance is validated against experiments for a single nucleation site in a vertical rectangular channel [19,20] (Sec. 3). Afterwards, studies for a bubble departing from the wall in a mini channel with OSF are also carried out (Sec. 4). Influence of surface characteristics such as contact angle and liquid superheat on bubble dynamics are analyzed. The results show that there is a small difference in bubble departure diameters (∼4%) for contact angles 30 deg and 45 deg. Local two-phase heat transfer coefficient is also investigated. To achieve higher void fractions, two nucleation cavities are introduced in the same channel with OSF (Sec. 5). For these studies, liquid superheat influence on bubble dynamics is also evaluated. Heat transfer enhancement due to bubble presence is quantified on planes along the bubbles' path.

2 Numerical Tools

2.1 PHASTA Description.

Finite element flow solver PHASTA (parallel, hierarchic, higher-order accurate, adaptive, stabilized, transient analysis) is used for the current simulations. The code solves incompressible Navier–Stokes equations in the strong form shown as follows:
uixi=0
(1)
ρuit+ρujuixj=pxi+τijxj+fi
(2)
ρcp(Tt+(u)T)=(kT)+q
(3)

where ui is the fluid velocity, ρ is the density, p is the static pressure, τij=μ(uixj+ujxi) is the viscous stress tensor, fi is the body force, cp is the specific heat at the constant pressure, T is the temperature, k is the thermal conductivity of a fluid, and q is the dissipation function. For surface tension, the continuum surface force model is used introduced by Ref. [21]. Discretization schemes that are implemented in PHASTA are described in more details in Refs. [22,23].

To capture the bubble interface, the level set method is employed [24]. The method introduces a distance function, φ, that defines the shortest distance to the interface from any point in the domain. At the interface, φ=0, whereas in the vapor phase, φ<0 and in the liquid phase φ>0. The distance function is advected by the flow velocity as
φt+uiφxi=0
(4)
To decrease a jump in properties across the interface, a smoothed Heaviside function Hε is applied in which interface thickness ε is used as
Hε(φ)=0,φ<ε12[1+φε+1πsin(πφε)]1,φ>ε,|φ|<ε
(5)

The Navier–Stokes equations are solved using the “one-fluid” approach for which fluid properties (density, viscosity, specific heat, and thermal conductivity) are found using Hε (e.g., ρ(φ)=ρlHε(φ)+ρg(1Hε(φ))ρl and ρg are the liquid and vapor phase densities respectively).

2.2 Contact Angle Control Algorithm.

The contact angle control algorithm [25] is employed in the current simulations to represent surface characteristics. For this algorithm, a target contact angle needs to be specified based on experimental data. Spatially varying subgrid force FCA is applied at the wall-adjacent cells to ensure that the current contact angle is equal to (or in a range of) the desired value.
FCA=F1cos(φTCAπ2)ρ(HCAdwall)2
(6)

In the above expression, F1 is the parameter that changes FCA with respect to the desired contact angle value [26], TCA and HCA are the thickness and height of FCA application region, respectively, anddwall is the distance from the wall. Current contact angle is calculated using interface normal vector (which is found using distance function as φ) and wall normal vector in wall-adjacent elements. More details on contact angle algorithm implementation in PHASTA and its verification could be found in Ref. [25].

2.3 Evaporation and Condensation Algorithm.

In the presented flow boiling simulations, the evaporation and condensation algorithm is used (previously introduced in PHASTA [27]). In this algorithm, it is assumed that the vapor phase temperature inside a bubble is constant (and is equal to the saturation temperature for the current conditions). Average temperature gradient is found for each bubble using the bubble tracking algorithm [28]. Local heat flux across the interface is calculated as
q=(R0R1)2klTn
(7)
In the above equation, R0 is the bubble radius, R1 is the radius of temperature gradient collection shell [27], kl is the liquid phase thermal conductivity, and Tn is the average temperature gradient across the interface. Volume change that is later applied to the RHS of the continuity Eq. (1) is found as
dV=4πR0q¯hfg(ρlρg)
(8)

Here, q¯ is the volume-averaged heat flux, and hfg is the latent heat of evaporation. This algorithm was verified against analytical bubble growth model [29] and experimental data for pool and flow boiling [30].

3 Model Validation

In order to validate the performance of the boiling and contact angle algorithms for flow boiling conditions, a separate study is conducted. Experiments in a vertical rectangular channel with HFE-301 were chosen as the references [19,20]. These experiments were conducted for the single nucleation site and such parameters as bubble departure diameters and departure frequencies are available for comparison.

The experiments [19,20] were carried out in a wide range of parameters. However, to match the desired flow conditions, the following flow characteristics are considered: mass flux 140 kg/m2·s (Re = 3125), wall heat flux 9.7 kW/m2, inlet subcooling 13.5 K, and atmospheric pressure.

The simulation domain constructed for this study is shown in Fig. 1. The channel dimensions are 324 × 10.1 × 10 mm. It has a total of 264 mm before the test section (including 50 mm unheated and 214 mm heated parts of the channel). The test section is 10 mm long and has a nucleation site of a trimmed conic shape placed in the middle of the section. An additional 50 mm are added after the test section for outlet treatment (to avoid backflow issue and dissipate bubbles since mesh is comparatively coarse there). The heated area in the simulations matches the experimental one.

Fig. 1
Computational domain for validation studies
Fig. 1
Computational domain for validation studies
Close modal

The real cavity size is of the order of nm which is computationally expensive to model. Instead, three bigger cavity sizes are considered: 0.288, 0.18, and 0.1125 mm of the cavity diameter (the depth is the same for all three cases 0.128 mm) (see Fig. 2). The mesh used near cavities will be reused for OSF simulations (discussed in Secs. 4 and 5).

Fig. 2
Different cavity sizes considered: baseline (top), 38% smaller (middle), and 61% smaller (bottom)
Fig. 2
Different cavity sizes considered: baseline (top), 38% smaller (middle), and 61% smaller (bottom)
Close modal

Exact surface characteristics are unknown [19,20]. However, using information from measured contact angles for refrigerant with similar properties [31] as well as magnified views of high-resolution camera [20], a contact angle 15 deg is set for these validation studies.

Bubble departure diameter and departure frequency are compared for three validation cases in Fig. 3, top. The waiting time is part of the fully resolved simulation and not specified by any modeling in our results. However, the limited number of nucleation sites results in large superheated area around the nucleation sites, and thus, the resolved waiting time is fairly small. Departure frequency is calculated as
fi=1tditni
(9)

where fi bubble departure frequency for ithbubble, tdi is the time of departure of ithbubble, and tni is the time of nucleation of ithbubble.

Fig. 3
Bubble statistics are compared with experiments [20] (dash lines on the left plot define uncertainty in departure diameter measurements)
Fig. 3
Bubble statistics are compared with experiments [20] (dash lines on the left plot define uncertainty in departure diameter measurements)
Close modal

It was observed that as the cavity diameter decreases, the departure diameters of bubbles also become smaller. This happens as the cavity is always filled with vapor and thus, the bubble interface is attached to the cavity's edges (constant boundary conditions are set for level set function on cavity walls). These boundary conditions are needed to ensure multiple bubble nucleation. Otherwise, there would be no vapor phase left in the cavity. Due to this, a smaller cavity diameter leads to faster bubble departure.

The results for the smallest cavity agree well with experimental data [20]. From bubble departure frequency plot (Fig. 3, top right) it can be seen that the data is consistent (i.e., smaller bubbles depart faster from the cavity). It is also noted that the modeled values are somewhat lower than observed experimentally (e.g., MSE for the smallest cavity is 146.2 Hz).

To better understand this, the amount of steam released from a cavity can be found with
rs=Vtb
(10)

where V is bubble volume (πD36) in cubic meters and tb is the nucleation period in seconds. The results are shown in the bottom of Fig. 3. From this plot it could be seen that the amount of steam is consistently lower for the smallest cavity (as bubble departure frequency results discussed before). This could be attributed to microlayer evaporation that is not considered in this study. In future work, additional source terms will be considered to account for the amount of steam generated in this tiny liquid film close to the wall.

4 Bubble Departure From the Wall

4.1 Model Setup.

To model the flow through OSF, a part of a vertical rectangular channel used in experiments [3] is considered (Fig. 4, left). In the domain, only the flow region is modeled. Therefore, there is a “step” of 0.2 mm corresponding to the fin thickness (Fig. 4, right) that is present in the y-direction. The flow is resolved for the first 12 mm of the domain, other parts are created to increase the code robustness (coarse meshes are used to prevent backflow and dissipate bubbles, the outlet treatment box includes last 6 mm (not shown in Fig. 4, right)).

Fig. 4
Geometry of OSF (left), corresponding computational domain (right, dimensions are shown in m)
Fig. 4
Geometry of OSF (left), corresponding computational domain (right, dimensions are shown in m)
Close modal

A no-slip boundary condition is imposed at the top and bottom of the domain where the flow is periodic in the z-direction. Inflow/outflow boundary conditions are set for velocity and temperature profiles in x-direction. Gravity is applied in the opposite direction of the flow. There is an uncertainty in the velocity and temperature profiles in the channel. Therefore, it is assumed that the flow is fully developed in the modeled part of the channel. Also, based on experimental data from Ref. [3], Re = 200 is chosen for these simulations. Parabolic velocity and temperature profiles are considered to be representative for this flow. Constant temperature Tw = 323.3 K (2 deg superheat) is set at the right and left walls of the domain (which represent inner channel walls assuming highly conductive fin material). According to FLUENT simulations [32], this value corresponds to heat flux of 1000 W/m2, which was used in the experiment in Ref. [3]. Since there is an uncertainty in exact temperature values at the inner walls, superheat studies were performed. Fluid properties for R113 at 102.4 kPa are presented in Table 1. Initial time-step is set 5 × 10−6 s and was automatically adjusted during the simulation to keep Courant–Friedrichs–Lewy (CFL) number = 1.0.

Table 1

Fluid properties

NameLiquid phase valuesVapor phase values
Dynamic viscosity (Pa·s)5×1041.1×105
Density (kg/m3)15077.5
Specific heat J/(kg·K)941692
Thermal conductivity W/(m·K)67.4×1039.51×103
Surface tension coefficient (N/m)0.01456
NameLiquid phase valuesVapor phase values
Dynamic viscosity (Pa·s)5×1041.1×105
Density (kg/m3)15077.5
Specific heat J/(kg·K)941692
Thermal conductivity W/(m·K)67.4×1039.51×103
Surface tension coefficient (N/m)0.01456

4.2 Mesh Design.

In such a complex geometry, considerable friction losses are present when the flow passes through the OSF. To properly resolve these losses, boundary layer meshing structure (BLs) is created on wall surfaces (Fig. 5). To accurately represent a bubble (around 25–30 elements across the bubble diameter are needed [33]) as well as to minimize the computational cost, several mesh refinement regions are introduced (Table 2). The bubble with an initial radius of 0.3 mm is placed at the location (0.00825; 0.000135; and 0.004775) m. This size was chosen to save computational resources (as smaller bubble would require higher mesh resolutions). The bubble is located at the bottom step of the fin. Due to stagnant liquid at the step corners, more superheated liquid is accumulated here. Therefore, bubble departure diameter is expected to be bigger rather than when the bubble is initialized at the top of the fin. In addition, the target contact angle of 30 deg is set. This value was also a source of uncertainty (since exact surface characteristics are unknown), therefore, a contact angle study was performed. For the contact angle algorithm to perform as designed, BLs are needed on the wall (Fig. 5, bottom).

Fig. 5
Mesh refinement regions (top), magnified view of the mesh near the initialized bubble (bottom)
Fig. 5
Mesh refinement regions (top), magnified view of the mesh near the initialized bubble (bottom)
Close modal
Table 2

Mesh sizes and number of elements for each mesh

#Mesh resolution, mPresence of BLsNumber of elements
11.6 × 10−4Yes707,092
21.28 × 10−3No
31.0 × 10−5No410,502
42.0 × 10−5Yes793,797
54.0 × 10−5no1,984,200
68.0 × 10−5No472,760
#Mesh resolution, mPresence of BLsNumber of elements
11.6 × 10−4Yes707,092
21.28 × 10−3No
31.0 × 10−5No410,502
42.0 × 10−5Yes793,797
54.0 × 10−5no1,984,200
68.0 × 10−5No472,760

4.3 Initial Conditions for Boiling Simulations.

It was mentioned that parabolic velocity and temperature profiles are imposed at the channel inlet to achieve laminar (Re = 200) flow. However, due to the staggered geometry of the channel (see Fig. 5, top), initially these profiles are somewhat different near the right and left walls. To overcome this discrepancy, single-phase simulations are conducted first to make sure that the temperature and velocity gradients are consistent near the walls. The time step size was adjusted automatically keeping maximum CFL number 20.0. The resulted velocity and temperature profiles are shown in Fig. 6. Single-phase simulations are carried out for all three bulk superheats considered (0, 2, and 4 deg). Then, the solutions at the last time step are set as initial conditions (ICs) for flow boiling simulations.

Fig. 6
Velocity and temperature profiles applied for three superheat cases (Color version online)
Fig. 6
Velocity and temperature profiles applied for three superheat cases (Color version online)
Close modal

4.4 Contact Angle Parametric Studies.

Contact angle is an important parameter for boiling simulations. To investigate the influence on bubble growth and departure, two contact angles are considered: 30 deg and 45 deg. Bubbles at the initial time step are presented in Fig. 7, top. The resulting bubble characteristics are shown in Table 3. Using the bubble tracking algorithm [28], bubble growth statistics are obtained (see Supplemental Figure S2 available in Supplemental Materials on the ASME Digital Collection).

Fig. 7
Bubbles at initial time step (top) and at the departure moments (bottom)
Fig. 7
Bubbles at initial time step (top) and at the departure moments (bottom)
Close modal
Table 3

Contact angle studies parameters

Case #Contact angle valueSuperheat at the wall (in the channel center), KBubble departure diameter, mDeparture time (ms)
130 deg2 (0)4.113 × 10−45.16
245 deg2 (0)3.964 × 10−44.37
Case #Contact angle valueSuperheat at the wall (in the channel center), KBubble departure diameter, mDeparture time (ms)
130 deg2 (0)4.113 × 10−45.16
245 deg2 (0)3.964 × 10−44.37

From Fig. 7, one may see that the bubble with a higher contact angle reaches smaller departure diameter and departs faster. Overall, the difference in bubble departure diameters for two contact angles considered is ∼4%. The results are obtained for the single bubble departure event only. To verify this observation, multiple bubble nucleation events are needed (to have more reliable bubble statistics). For that purpose, Sec. 5 discusses flow boiling simulations conducted in the same geometry with OSF, but with two nucleation sites.

4.5 Superheat Studies.

Parabolic temperature profile as well as wall temperature are the sources of uncertainty in the current study. To investigate the influence of wall temperature values on bubble growth, three superheat studies are conducted (Table 4). The simulations were conducted until a bubble reached the coarser mesh region (mesh #6, Fig. 5). The average distance from the nucleation site that the bubble covered is around 0.00375 m. Bulk superheat is chosen to keep the same temperature difference between the wall and the channel center. Mesh design and other model conditions are the same as in the contact angle study. For the current simulation, the target contact angle is 30 deg, advancing contact angle is 25 deg and receding contact angle is 35 deg. The resulted bubble growth statistics are shown in Supplemental Figure S3 available in Supplemental Materials, and bubble positions at the departure moments are presented in Supplemental Figure S4 available in Supplemental Materials on the ASME Digital Collection.

Table 4

Superheat studies parameters

Case #Contact angle valueSuperheat at the wall (in the channel center), KBubble departure diameter, mDeparture time (ms)
330 deg2 (0)5.237 × 10−410.62
44 (2)8.126 × 10−417.28
56 (4)1.410 × 10−327.09
Case #Contact angle valueSuperheat at the wall (in the channel center), KBubble departure diameter, mDeparture time (ms)
330 deg2 (0)5.237 × 10−410.62
44 (2)8.126 × 10−417.28
56 (4)1.410 × 10−327.09

Higher superheat values at the surface lead to a more intensive evaporation process. Because of that, bubbles reach bigger sizes for cases #4–5 compared to case #3. However, since vapor prevents liquid from coming closer to the surface, it exacerbates the heat transfer for case #4–5 relative to case #3. On the other hand, the small-size bubble (as in case #3) departs much sooner from the wall and, thereby, heat transfer coefficient increases. Another observation is that a larger bubble (cases #4–5) slides longer along the wall (see bubble departure locations in Supplemental Figure S4 available in Supplemental Materials on the ASME Digital Collection: in case #4 the bubble slides to the fin corner and then departs; in case #5 the bubble departs only in the middle of the next fin since it was sliding along the surface before that). Such bubble behavior tends to disturb thermal boundary layer, hence, increasing heat transfer coefficient [34,35]. More details on heat transfer coefficient evaluation will be presented in Sec. 5.3.3.

4.6 Local Void Fraction Calculation.

Since only one bubble is present in each simulation, a smaller part of the channel is considered (Fig. 8, left) for local void fraction evaluation. The resulting void fraction statistics for each case are shown in Fig. 8, right. It is noted that since for case #5 the bubble size is larger, its corresponding void fraction has higher values compared to cases #4 and 3. The bubble growth is somewhat similar for each case: as the bubble attaches to the surface, its volume grows nonlinearly. However, when the bubble departs and starts to move toward the channel center, the growth rate becomes more linear. Since for case #3, there is no superheat at the channel center, bubble growth rate in this region is considerably smaller than, for example, in case #4. For case #5, there is 4 deg superheat at the channel center which, consequently, leads to a much higher bubble growth rate compared to cases #3 and 4. The final local void fraction achieved for cases #3–5 are 0.57%, 1.42%, and 4.78%, respectively.

Fig. 8
Local void fraction evaluation region (left) and void fraction statistics for cases #3–5 (right) (Color version online)
Fig. 8
Local void fraction evaluation region (left) and void fraction statistics for cases #3–5 (right) (Color version online)
Close modal

Bulk superheat studies show that the temperature distribution in the channel plays a crucial role in heat transfer evaluation. The difference in 2 deg (cases #3 and 4) results in ∼50% increase of the bubble departure diameter and the change in 4 deg (cases #3 and 5) — in nearly 170% growth.

4.7 Local Heat Transfer Evaluation.

Local two-phase heat transfer coefficient is found for the cases #3–5 (see Fig. 9) with
h=qΔT=kΔT1SSTdS
(11)

In the equation above, values are found from Fourier's Law using area-averaged normal to the wall temperature gradient values (1SSTdS in Eq. (11)), thermal conductivity of R113 liquid phase k (see Table 1 for reference), and known temperature differences at the walls and the centerline (ΔT = 2 deg). More details on that could be found in Ref. [36]. These temperature gradients were obtained at 0.04 mm from the wall. This value is chosen to make sure that temperature changes linearly over this distance. The calculations are performed for each plane on the bubble path (see Supplemental Figure S6 available in Supplemental Materials) at four time intervals: 0, 9, 18, and 27 ms. Points at t= 0 s correspond to single-phase heat transfer coefficients.

Bubbles with low density vapor move faster than liquid. Correspondingly, a bubble of a bigger size travels further downstream (see Supplemental Figure S1 available in Supplemental Materials). The results show that local heat transfer coefficient for case #5 is generally higher than for cases #3–4. This could imply that the nucleate boiling part of the two-phase heat transfer coefficient that depends on the liquid superheat plays an important role (for case #5 the superheat rate at the wall is the highest). At the same time, convective boiling heat transfer part is also significant for case #5: from snapshots of bubbles' location, one could notice a larger disturbance of the thermal boundary layer caused by the biggest bubble sliding along the wall (compared to cases #3–4).

By analyzing the change in the heat transfer coefficient on plane #1, it was concluded that as the bubble grows, it acquires more energy due to phase change. Then as the bubble starts to move downstream following the flow velocity, it leaves the plane, and the heat transfer coefficient correspondingly decreases for all cases. When the bubble comes closer to plane #2, the heat transfer coefficient is increased significantly. It should be noted that since for cases #3–4 bubbles traveled smaller distance over this plane, the growth of the heat transfer coefficient is not as significant as for case #5. For plane #3, the heat transfer coefficient starts to increase as the bubble would closer approach this surface. For planes #4 and #5 (side walls), the heat transfer coefficient trend can be explained as follows: when the bubble is relatively far away from the side surface (as, for example, in cases #3 and 4), due to phase change, it consumes the amount of liquid in the vicinity of the wall, thereby, degrading convective heat transfer. However, when the bubble comes closer to the wall (or when its size increases), the boiling heat transfer prevails. Observed local heat transfer coefficient values are obtained for the section of the channel where the bubble is nucleated. Overall, the heat transfer enhancement effect is more pronounced for plane #1 where the bubble is nucleated. Once it starts to grow (at t= 9 ms), the heat transfer coefficient is increased by 44.2% (case #3), 86.1% (case #4), and 315.7% (case #5). To validate the results, the authors intend to compare it with the experiments [3] as the part of the future work. However, for this a higher void fraction should be achieved which imposes additional numerical challenges.

5 Bubble Departure From Nucleation Cavities

5.1 Model Setup.

To increase void fraction, two nucleation cavities are introduced in the same channel on the right and left walls (see Supplemental Figure S7 available in Supplemental Materials on the ASME Digital Collection). Cavities have conic shape (diameter 0.288 mm, depth 0.128 mm) and are trimmed at the bottom. The two cavities are located at (8.25; 0; 4.675) and (9.75; 2.8; 4.675) mm. Exact locations of nucleation sites in experiments [3] are unknown as they depend on surface imperfections and may change for different heat flux conditions. Information on mesh used for these studies is available in Supplemental Materials on the ASME Digital Collection.

As for the previous runs, since the exact temperature values at the inner channel walls likely depend on the streamwise location in the full heat exchanger design (in experiments the size of the channel is 750 × 3 × 100 mm), several simulations were conducted with three different bulk superheat values: 0, 2, and 4 deg. These simulations are conducted with Re = 200. For these simulations the target contact angle is set to 15 deg, advancing contact angle is 10 deg, and receding contact angle is 20 deg. The same working fluid is considered R113 (Table 1). Simulations are carried out with the initial time step 2.1 ms, and it was adjusted during the runs to keep flow CFL number = 0.5.

It has been observed that when a bubble grows at the wall, there is a thin liquid film exists between the bubble and the wall called microlayer [37,38]. The contribution of the microlayer to bubble evaporation is substantial (heat fluxes reach 1 MW/m2 [39]) especially at atmospheric pressure. However, as was shown in Ref. [40], heat flux from liquid evaporation in the microlayer region is negligible for low heat fluxes. Since in the current simulations, applied wall temperature corresponds to heat flux 1000 W/m2, the contribution from the microlayer is not considered.

Bulk superheat plays a major role in the bubble growth, and it is important to accurately represent it in the simulations. There are several examples in which conjugate heat transfer between the solid and fluid domains was accounted for in boiling simulations, including [4144]. However, in the current studies, wall thickness was small (0.2 mm), and the material is highly conductive (copper) [3], therefore, conjugate heat transfer is neglected.

5.2 Superheat Studies

5.3.1 Difference in Temperature Distribution.

To quantify the difference in the temperature distribution in the channel, sensitivity studies with three different bulk superheats are conducted. Departure diameters of the first-generation bubbles (e.g., the first bubble produced at each site) are summarized in Table 5. The case with 0 deg superheat is considered as a baseline for this study. From Table 5 it can be seen that a 2 deg variation in the bulk superheat causes a significant difference between departure diameters. This effect is more pronounced for 4 deg superheat difference. In Supplemental Figure S5 available in Supplemental Materials on the ASME Digital Collection, bubble shapes for the first-generation bubbles at the right wall cavity are shown at simulation time 6 ms. It could be seen that for 2 and 4 deg bulk superheat, bubbles readily depart, whereas for the 0 deg bulk superheat case, the bubble is more closely attached to the wall. Compared to the bubble departure from the wall simulations (the results are shown in Table 4), current bubble departure diameters are smaller since a smaller contact angle (i.e., 15 deg) is imposed.

Fig. 9
Local heat transfer coefficient changes for different planes on the bubble path
Fig. 9
Local heat transfer coefficient changes for different planes on the bubble path
Close modal
Table 5

Departure diameter of the first-generation bubbles

Departure diameter, mm
Bulk superheat, KNucleation site #1Nucleation site #2Average difference, %
00.45350.453
20.51260.512213.05
40.63790.632940.19
Departure diameter, mm
Bulk superheat, KNucleation site #1Nucleation site #2Average difference, %
00.45350.453
20.51260.512213.05
40.63790.632940.19

5.3.2 Bubble Statistics.

Some of the bubble statistics are shown in Fig. 10, bubble trajectories are available in Supplemental Material. From these plots it can be seen that departure diameter is generally higher for the higher bulk superheat. For 2 and 0 deg bulk superheat cases, the obtained departure diameter values are more consistent whereas for 4 deg a higher discrepancy was observed (which was attributed to coalescence between just departed bubbles and bubbles left in cavities). Bubble departure frequency is more consistent when there is no coalescence between already departed bubbles and newly nucleated ones.

Fig. 10
Bubble statistics for nucleation site #1 (on the right wall, open symbols) and #2 (on the left wall, solid symbols) for different bulk superheats (Color version online)
Fig. 10
Bubble statistics for nucleation site #1 (on the right wall, open symbols) and #2 (on the left wall, solid symbols) for different bulk superheats (Color version online)
Close modal

Snapshots of simulations for 4 deg bulk superheat are shown in Fig. 11; for the 0 and 2 deg bulk superheats in Fig. 12. One could notice that for higher temperature gradients values near the walls, bubbles grow larger (for example, compare all three cases at 8 ms). For 4 deg bulk superheat, bubbles reach around 1 mm in diameter (more detailed data is presented below) and they extensively coalesce with each other (see snapshot for 32 ms, Fig. 11). At the end of the simulation, the vapor phase occupies more than a half of the channel cross section which could later result in a flow transition (from bubbly to slug flow).

Fig. 11
Bubble dynamics for 4 deg bulk superheat
Fig. 11
Bubble dynamics for 4 deg bulk superheat
Close modal
Fig. 12
Bubble dynamics for 0 and 2 deg bulk superheat
Fig. 12
Bubble dynamics for 0 and 2 deg bulk superheat
Close modal

Results shown in Fig. 12 indicate that for lower bulk superheat values (e.g., 0 deg) bubbles have smaller departure diameters. Due to the size, bubbles coalesce less often and stop growing once they are near the centerline (there is no superheat here). For 2 deg bulk superheat, more bubbles are observed compared to the previous case due to higher temperature values in the channel.

5.3.3 Local Heat Transfer Evaluation.

Heat transfer coefficient is calculated locally on planes along bubbles' path (see Fig. 13) for three considered cases. The same procedure was used to calculate local heat transfer coefficients as was discussed previously in Sec. 4.7. The results are shown in Fig. 14. For all three plots presented, points at t = 0 s correspond to single-phase heat transfer (these are solutions at the last time-step of single-phase simulations which are recycled as initial conditions for the two-phase flow runs).

Fig. 13
Planes used for local heat transfer evaluation
Fig. 13
Planes used for local heat transfer evaluation
Close modal
Fig. 14
Local heat transfer coefficient
Fig. 14
Local heat transfer coefficient
Close modal

If one considers heat transfer evolution for bulk superheat of 0 deg, it may be seen that for plane #1, the value does not change significantly. This implies that conditions for bubble growth and departure stay consistent over time. The plots from Fig. 10 for nucleation site #1 also support this point. The same conclusion can be made for heat transfer coefficient behavior on plane #1 for bulk superheat 2 deg. For 4 deg bulk superheat, there is a decline in heat transfer after ∼10 ms. This happens due to coalescence between the departed bubble and the newly nucleated one. The newly merged bubble stays longer near the wall consuming more energy from the surrounding liquid. For the subsequent bubble (point at ∼16 ms), no coalescence event was observed. At ∼32 ms, an increase in heat transfer for plane #1 attributes to the fact that the cavity is full of liquid (no vapor phase is left at this time). The situation is different for plane #6: increases in heat transfer coefficient (for bulk superheats 0 and 2 deg after ∼60 and 30 ms correspondingly) suggest that small-size bubbles are nucleated (could be seen in Fig. 10 as well). Thus, more liquid occupies the cavity causing the rise of heat transfer. For 4 deg bulk superheat, there is a coalescence event again (at ∼10 ms) happening between the departed and newly nucleated bubbles which contributes to a decrease in heat transfer coefficient.

For plane #2, a decrease in heat transfer coefficient is observed for all three cases (starting from ∼8 ms). This is the result of the numerical assumption from the boiling algorithm — vapor phase is considered to be at saturation temperature. During single-phase simulations, bubbles initially nucleated in the domain work as heat sinks. As the code runs, this saturated temperature affects thermal boundary layer (areas of lower temperature near walls could be seen at 8 ms in Figs. 11 and 12). Although this effect distorts heat transfer evaluation at the beginning, it disappears as boiling simulations are run longer.

For plane #3 one could clearly see the effect of bulk superheat on heat transfer: before bubbles reach this plane, the average heat transfer coefficient values are 114.7 (0 deg), 189.5 (2 deg), and 253.5 (4 deg) W/(m2·K).

For planes #4 and 5, the situation is similar: values are nearly constant until bubbles come closer to these planes causing disturbances in thermal and hydrodynamic boundary layers (e.g., see the growth on plane #4 after ∼60 ms for bulk superheat 0 deg, Fig. 14).

6 Conclusions

Interface capturing flow boiling simulations were conducted for a vertical mini channel with offset strip fins (OSF) using code PHASTA. To establish the credibility of the results, model validation is conducted for a single nucleation site in a vertical rectangular channel. The observed bubble departure diameter agrees well with the experiments for the smallest cavity size considered. The steam generation and bubble departure frequency are underpredicted which could be attributed to microlayer evaporation (not considered in the current study).

Second, a complex geometry with OSF under single-bubble flow boiling was analyzed. Bubble dynamics characteristics are investigated with varying contact angle and bulk superheat. It was observed that temperature distribution in the channel plays a major role in determining heat transfer process for this laminar flow (Re = 200). The increase in bulk superheat by 2 deg leads to a nearly 50% growth of the bubble departure diameter and the increase by 4 deg allows to achieve almost a 170% growth.

Third, to observe multiple consequent bubble nucleation, two nucleation cavities were introduced in the same channel with OSF. Single-phase solutions are used as initial conditions for boiling simulations. Bulk superheat sensitivity tests are conducted to analyze bubble dynamics. Analysis of departure diameter values shows that the 4 deg temperature difference in the bulk fluid may result in a 40% difference in bubble departure diameter value. Local heat transfer coefficient is evaluated on each plane along the bubbles' path. It is shown that for some cases (e.g., for bulk superheat 0 deg, plane #3) bubbles approaching the wall may increase heat transfer coefficient by up to 2.5 times. With a smaller nucleation diameter, more liquid is cooling the cavity which, in turn, results in nearly twice as high heat transfer coefficient (e.g., for bulk superheat 2 deg, plane #6).

Future work will be focused on improving model predictions of bubble departure frequency and amount of steam released from a cavity as well as testing PHASTA capabilities on modeling flows with a higher void fraction.

Acknowledgment

ACUSIM linear algebra solution library was employed which was provided by Altair Engineering Inc. along with mesh and geometry modeling libraries provided by Simmetrix Inc. for simulations. We also acknowledge the computing resources provided on Henry2, a high-performance computing cluster operated by North Carolina State University. Authors greatly appreciate the constructive and useful comments from the reviewer as well as valuable feedback from Edward Chen.

Funding Data

  • Current work is performed under Two-Phase Flow DNS Phase 2 Project supported by Mitsubishi Heavy Industries, Ltd.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

Nomenclature

cp =

specific heat at the constant pressure, J/(kg·K)

dwall =

distance from the wall, m

fi =

body force (including gravity and surface tension), Pa/m

FCA =

contact angle force, H

F1 =

contact angle algorithm parameter

hfg =

latent heat of evaporation, J/kg

HCA =

height of contact angle force application region, m

Hε =

smoothed Heaviside function

k =

thermal conductivity of a fluid, W/(m·K)

p =

static pressure, Pa

q =

dissipation function, J/(m3·s)

q¯ =

volume-averaged heat flux, W/m4

R0 =

bubble radius, m

R1 =

radius of temperature gradient collection outer shell, m

t =

time, s

tb =

nucleation period, s

T =

temperature, K

TCA =

thickness of contact angle force application region, m

ui =

fluid velocity, m/s

V =

bubble volume, m3

ρ =

density, kg/m3

τij =

viscous stress tensor, Pa

ε =

interface thickness, m

ρl =

liquid density, kg/m3

ρg =

vapor density, kg/m3

References

1.
Lee
,
J.
, and
Mudawar
,
I.
,
2005
, “
Two-Phase Flow in High-Heat-Flux Micro-Channel Heat Sink for Refrigeration Cooling Applications: Part II - Heat Transfer Characteristics
,”
Int. J. Heat Mass Transfer.
,
48
(
5
), pp.
941
955
.10.1016/j.ijheatmasstransfer.2004.09.019
2.
Bertsch
,
S.
,
Groll
,
E.
, and
Garimella
,
S.
,
2009
, “
Effects of Heat Flux, Mass Flux, Vapor Quality, and Saturation Temperature on Flow Boiling Heat Transfer in Microchannels
,”
Int. J. Multiphase Flow
,
35
(
2
), pp.
142
154
.10.1016/j.ijmultiphaseflow.2008.10.004
3.
Kim
,
B.
, and
Sohn
,
B.
,
2006
, “
An experimental Study of Flow Boiling in a Rectangular Channel With Offset Strip Fins
,”
Int. J. Heat Fluid Flow
,
27
, pp.
514
521
.10.1016/j.ijheatfluidflow.2005.11.008
4.
Raju
,
M. A.
,
Ashok Babu
,
T.
, and
Ranganayakulu
,
C.
,
2017
, “
Flow boiling Heat Transfer and Pressure Drop Analysis of R134a in a Brazed Heat Exchanger With Offset Strip Fins
,”
Heat Mass Transfer
,
53
, pp.
3167
3180
.10.1007/s00231-017-2060-1
5.
Pulvirenti
,
B.
,
Matalone
,
A.
, and
Barucca
,
U.
,
2010
, “
Boiling heat Transfer in Narrow Channels With Offset Strip Fins: Application to Electronic Chipsets Cooling
,”
Appl. Therm. Eng.
,
30
, pp.
2138
2145
.10.1016/j.applthermaleng.2010.05.026
6.
Wen
,
T.
,
Zhan
,
H.
, and
Zhang
,
D.
,
2019
, “
Flow boiling Heat Transfer in Mini Channel With Serrated Fins: Experimental Investigation and Development of New Correlation
,”
Int. J. Heat Mass Transfer
,
128
, pp.
1081
1094
.10.1016/j.ijheatmasstransfer.2018.09.071
7.
Kim
,
M.
, and
Lee
,
K.
,
2010
, “
The thermoflow Characteristics of an Oscillatory Flow in Offset-Strip Fins
,”
Numer. Heat Transfer, Part A
,
58
, pp.
835
851
.10.1080/10407782.2010.529033
8.
Manglik
,
R. M.
, and
Bergles
,
A. E.
,
1995
, “
Heat transfer and Pressure Drop Correlations for the Rectangular Offset Strip Fin Compact Heat Exchanger
,”
Exp. Therm. Fluid Sci.
,
10
, pp.
171
180
.10.1016/0894-1777(94)00096-Q
9.
Kim
,
M.
,
Lee
,
J.
,
Yook
,
S.
, and
Lee
,
K.
,
2011
, “
Correlations and Optimization of a Heat Exchanger With Offset-Strip Fins
,”
Int. J. Heat Mass Transfer
,
54
, pp.
2073
2079
.10.1016/j.ijheatmasstransfer.2010.11.056
10.
DeJong
,
N.
,
Zhang
,
L.
,
Jacobi
,
A.
,
Balachandar
,
S.
, and
Tafti
,
D.
,
1998
, “
A Complementary Experimental and Numerical Study of the Flow and Heat Transfer in Offset Strip-Fin Heat Exchangers
,”
ASME J. Heat Mass Transfer-Trans. ASME
,
120
, pp.
690
698
.10.1115/1.2824338
11.
Saad
,
S.
,
Clement
,
P.
,
Gentric
,
C.
,
Fourmigue
,
J.
, and
Leclerc
,
J.
,
2011
, “
Experimental distribution of Phases and Pressure Drop in a Two-Phase Offset Strip Fin Type Compact Heat Exchanger
,”
Int. J. Multiphase Flow
,
37
, pp.
576
584
.10.1016/j.ijmultiphaseflow.2011.03.009
12.
Song
,
R.
,
Cui
,
M.
, and
Liu
,
J.
,
2017
, “
A correlation for Heat Transfer and Flow Friction Charactiristics of the Offset Strip Fin Heat Exchanger
,”
Int. J. Heat Mass Transfer
,
115
, pp.
695
705
.10.1016/j.ijheatmasstransfer.2017.08.054
13.
Raju
,
M. A.
,
Ashok Babu
,
T.
, and
Ranganayakulu
,
C.
,
2015
, “
Development of Single Phase Heat Transfer Correlations for Water & R134a in Rectangular Channel With Smooth Wavy Fin
,”
J. Phys. Sci.
,
5
, pp.
199
208.
10.17265/2159-5348/2015.03.006
14.
Pan
,
Z.
,
Weibel
,
J.
, and
Garimella
,
S.
,
2016
, “
A saturated-Interface-Volume Phase Change Model for Simulating Flow Boiling
,”
Int. J. Heat Mass Transfer
,
93
, pp.
945
956
.10.1016/j.ijheatmasstransfer.2015.10.044
15.
Ling
,
K.
,
Zhang
,
S.
,
Liu
,
W.
,
Sui
,
X.
, and
Tao
,
W.
,
2021
, “
Interface Tracking Simulation for Subcooled Flow Boiling Using VOSET Method
,”
Front. Energy Res.
,
8
, p. 526035.10.3389/fenrg.2020.526035
16.
Allred
,
T. P.
,
Weibel
,
J. A.
, and
Garimella
,
S. V.
,
2022
, “
The effect of Dynamic Wetting Behavior on Boiling Heat Transfer Mechnisms During Bubble Growth and Departure
,”
Int. J. Heat Mass Transfer
,
184
, p.
122276
.10.1016/j.ijheatmasstransfer.2021.122276
17.
Lee
,
J.
,
O'Neill
,
L. E.
,
Lee
,
S.
, and
Mudawar
,
I.
,
2019
, “
Experimental and Computational Investigation on Two-Phase Flow and Heat Transfer of Highly Subcooled Flow Boiling in Vertical Upflow
,”
Int. J. Heat Mass Transfer
,
136
, pp.
1199
1216
.10.1016/j.ijheatmasstransfer.2019.03.046
18.
Lei
,
Y.
,
Mudawar
,
I.
, and
Chen
,
Z.
,
2020
, “
Computational and Experimental Investigation of Condensation Flow Patterns and Heat Transfer in Parallel Rectangular Micro-Channels
,”
Int. J. Heat Mass Transfer
,
149
, pp.
119
158
.10.1016/j.ijheatmasstransfer.2019.119158
19.
Yoo
,
J.
,
Estrada-Perez
,
C. E.
, and
Hassan
,
Y. A.
,
2016
, “
Experimental study on Bubble Dynamics and Wall Heat Transfer Arising From a Single Nucleation Site at Subcooled Flow Boiling conditions - Part 1: Experimental Methods and Data Quality Verification
,”
Int. J. Multiphase Flow
,
84
, pp.
315
324
.10.1016/j.ijmultiphaseflow.2016.04.018
20.
Yoo
,
J.
,
Estrada-Perez
,
C. E.
, and
Hassan
,
Y. A.
,
2016
, “
Experimental study on Bubble Dynamics and Wall Heat Transfer Arising From a Single Nucleation Site at Subcooled Flow Boiling conditions - Part 2: Data Analysis on Sliding Bubble Characteristics and Associated Wall Heat Transfer
,”
Int. J. Multiphase Flow
,
84
, pp.
292
314
.10.1016/j.ijmultiphaseflow.2016.04.019
21.
Brackbill
,
J.
,
Kothe
,
D.
, and
Zemach
,
C.
,
1992
, “
A continuum Method for Modeling Surface Tension
,”
J. Comput. Phys.
,
100
(
2
), pp.
335
354
.10.1016/0021-9991(92)90240-Y
22.
Jansen
,
K.
,
1999
, “
A stabilized Finite Element Method for Computing Turbulence
,”
Comput. Methods Appl. Mech. Eng.
,
174
(
3
), pp.
299
317
.10.1016/S0045-7825(98)00301-6
23.
Whiting
,
C.
, and
Jansen
,
K.
,
2001
, “
A stabilized Finite Element Method for the Incompressible Navier-Stokes Equations Using a Hierarchical Basis
,”
Int. J. Numer. Methods Fluids
,
35
(
1
), pp.
93
116
.10.1002/1097-0363(20010115)35:1<93::AID-FLD85>3.0.CO;2-G
24.
Sussman
,
M.
,
Almgren
,
A.
,
Bell
,
J.
,
Colella
,
P.
,
Howell
,
L.
, and
Welcome
,
M.
,
1999
, “
An adaptive Level Set Approach for Incompressible Two-Phase Flows
,”
J Comput. Phys.
,
148
(
1
), pp.
81
124
.10.1006/jcph.1998.6106
25.
Li
,
M.
,
Zeng
,
K.
,
Wonnell
,
L.
, and
Bolotnov
,
I.
,
2019
, “
Development of a New Contact Angle Control Algorithm for Level-Set Method
,”
ASME
J. Fluids Eng.,
141
, p.
061301
.10.1115/1.4041987
26.
Mishra
,
A.
, and
Bolotnov
,
I.
,
2013
, “
Contact Angle Control Algorithm Development for Level-Set Interface Tracking Method
,”
Transactions of 2013 ANS Annual Meeting
, Atlanta, GA, June 16–20, pp.
1017
1019
.
27.
Li
,
M.
, and
Bolotnov
,
I.
,
2016
, “
Interface Tracking Simulation of Phase-Change Phenomena: Boiling and Condensation Verification
,”
ASME
Paper No. FEDSM2016-7701.10.1115/FEDSM2016-7701
28.
Fang
,
J.
, and
Bolotnov
,
I.
,
2017
, “
Bubble tracking Analysis of PWR Two-Phase Flow Simulations Based on the Level Set Method
,”
Nucl. Eng. Des.
,
323
, pp.
68
77
.10.1016/j.nucengdes.2017.07.034
29.
Scriven
,
L.
,
1959
, “
On the Dynamics of Phase Growth
,”
Chem. Eng. Sci.
,
10
(
1–2
), pp.
1
13
.10.1016/0009-2509(59)80019-1
30.
Li
,
M.
,
2019
, “
High Resolution Boiling Simulation Using Interafce Tracking Method
,” Ph.D. dissertation,
North Carolina State University
,
Raleigh, NC
.
31.
Yue Sze
,
V.
,
2020
, “
Flow Boiling of HFE-7200 in Multi-Microchannel Heat Sinks for High-Heat Flux Applications
,” Ph.D. dissertation,
Brunel University
,
UK
.
32.
Kondo
,
Y.
, and
Tanimoto
,
K.
,
2020
,
Porous two-Phase CFD on Offset-Strip Fin Apparatus to Provide Boundary Data for Two-Phase Flow DNS
,
Mitsubishi Heavy Industries, Ltd
.,
Japan
.
33.
Li
,
M.
, and
Bolotnov
,
I.
,
2020
, “
The Evaporation and Condensation Model With Interface Tracking
,”
Int. J. Heat Mass Transfer
,
150
, pp.
119
256
.10.1016/j.ijheatmasstransfer.2019.119256
34.
Zhou
,
K.
,
Coyle
,
C.
,
Li
,
J.
,
Buongiorno
,
J.
, and
Li
,
W.
,
2017
, “
Flow boiling in Vertical Narrow Microchannels of Different Surface Wettability Characteristics
,”
Int. J. Heat Mass Transfer
,
109
, pp.
103
114
.10.1016/j.ijheatmasstransfer.2017.01.111
35.
Ahmadi
,
R.
, and
Okawa
,
T.
,
2015
, “
Influence of Surface Wettability on Bubble Behavior and Void Evolution in Subcooled Flow Boiling
,”
Int. J. Therm. Sci.
,
97
, pp.
114
125
.10.1016/j.ijthermalsci.2015.06.012
36.
Bayazit
,
B.
,
Keith Hollingsworth
,
D.
, and
Witte
,
L.
,
2003
, “
Heat Transfer Enhancement Caused by Sliding Bubbles
,”
ASME J. Heat Mass Transfer-Trans.
,
125
, pp.
503
509
.10.1115/1.1565090
37.
Yuan
,
D.
,
Chen
,
D.
,
Yan
,
X.
,
Xu
,
J.
,
Lu
,
Q.
, and
Huang
,
Y.
,
2018
, “
Bubble behavior and Its Contribution to Heat Transfer of Subcooled Flow Boiling in a Vertical Rectangular Channel
,”
Ann. Nucl. Energy
,
119
, pp.
191
202
.10.1016/j.anucene.2018.05.010
38.
Iskhakova
,
A.
,
Dinh
,
N. T.
,
Bolotnov
,
I. A.
,
Kondo
,
Y.
, and
Tanimoto
,
K.
,
2022
, “
Flow Boiling Heat Transfer Study in a Mini Channel Using DNS Approach
,”
The 19th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-19)
, Brussels, Belgium, Ma. 6–11, (virtual).
39.
Sato
,
Y.
, and
Niceno
,
B.
,
2015
, “
A depletable Micro-Layer Model for Nucleate Pool Boiling
,”
J. Comput. Phys.
,
300
, pp.
20
52
.10.1016/j.jcp.2015.07.046
40.
Dhir
,
V.
,
Abarajith
,
H.
, and
Li
,
D.
,
2007
, “
Bubble Dynamics and Heat Transfer During Pool and Flow Boiling
,”
Heat Transfer Eng.
,
28
(
7
), pp.
608
624
.10.1080/01457630701266421
41.
Yabuki
,
T.
, and
Nakabeppu
,
O.
,
2014
, “
Heat transfer Mechanisms in Isolated Bubble Boiling of Water Observed With MEMS Sensor
,”
Int. J. Heat Mass Transfer
,
76
, pp.
286
297
.10.1016/j.ijheatmasstransfer.2014.04.012
42.
Sato
,
Y.
, and
Niceno
,
B.
,
2017
, “
Nucleate Pool Boiling Simulations Using the Interface Tracking Method: Boiling Regime From Discrete Bubble to Vapor Mushroom Region
,”
Int. J. Heat Mass Transfer
,
105
, pp.
505
524
.10.1016/j.ijheatmasstransfer.2016.10.018
43.
Sato
,
Y.
, and
Niceno
,
B.
,
2018
, “
Pool boiling Simulation Using an Interface Tracking Method: From Nucleate Boiling to Film Boiling Regime Through Critical Heat Flux
,”
Int. J. Heat Mass Transfer
,
125
, pp.
876
890
.10.1016/j.ijheatmasstransfer.2018.04.131
44.
Fuchs
,
T.
,
Kern
,
J.
, and
Stephan
,
P.
,
2006
, “
A Transient Nucleate Boiling Model Including Microscale Effects and Wall Heat Transfer
,”
ASME J. Heat Mass Transfer-Trans.
,
128
, pp.
1257
1265
.10.1115/1.2349502

Supplementary data