CFD analysis in effects of spacer‑vane on thermal–hydraulic performance of two‑phase flow in annular channel

Numerous systems that rely on the phenomenon of two-phase flow are prevalent across various industrial sectors, including chemical, petrochemical, and

and bubble size).Therefore, the examination of thermal-hydraulic performance in twophase flow plays an essential part in the effective design of rod bundles equipped with grid spacer-vane and in the safe functioning of the reactor core.Many researchers have used numerical and experimental methods to predict twophase thermal-hydraulic performance in nuclear rod applications.Serizawa et al. [1] investigated experimentally the turbulent phenomena of water-air bubbly flow in a vertical pipe.A tracer technique with a concentration of helium gas was used for the analysis of heat and transport of bubbles in two-phase flow.Numerical and experimental analyses of the flow behavior of subcooled boiling in vertical annulus were conducted by Lee et al. [2].Water and vapor were considered as fluid of two-phase flow, and analyses were performed for different mass fluxes range 476-1061 kg/m 2 s.Ustinenko et al. [3] developed the methodology using CFD (computational fluid dynamics) to validate the various experimental work performed in two-phase flow by past researchers.Thermal-hydraulic performance was analyzed with constant heat flux in annular channels and pipes.
Ala et al. [4] studied the impacts of partial and complete blockage of the flow passage of the subchannels in nuclear fuel rod bundle for two-phase fluid flow.The authors investigated the blockage effects on void fraction, turbulent intensity, pressure drop, and velocity distributions inside the subchannels at different locations downstream.Chalgeri and Jeong [5] studied the flow regimes transition of two-phase fluid flow in rectangular channels with narrow gaps.Air and water were employed as two-phase fluids inside the domain.They proposed a model for transition criteria and showed good agreement against measurement data.Sadatomi et al. [6,7] performed experiments on rod bundles having triangular and square arrays for two-phase (air-water) flow.The fluctuation in differences in static pressure between the subchannels of the rod bundle was measured.Turbulent mixing rates were compared for the different arrays used in the experiments and concluded that array configurations of subchannels play an important part in examining the turbulent mixing rate.
Kawahara et al. [8][9][10] conducted experiments on the rectangular-shaped subchannels for two-phase fluid flow.Water and air flow characteristics visualization was accomplished with the help of acrylic glass.Authors reported that large bubbles have a significant effect on turbulent mixing rate as compared to that of small bubbles, and also, they developed a model for slug flow to predict the effect on turbulent mixing rate.
Dhurandhar et al. [11] presented a thorough review of mixing vane impacts on thermal-hydraulic characteristics for different channel geometries.A detailed review was provided by considering supercritical fluids (as single phase) and different two-phase fluids.Dhurandhar et al. [12][13][14][15] performed computer simulations to examine the impacts of spacers and vanes on thermal-hydraulic characteristics of various fluids in a vertical flow annular channel and concluded that spacers and the vanes are significant channel structures for improving heat transfer.
Arai et al. [16] executed an experimental analysis of void fraction in a boiling water reactor for non-heated provisions of the rods.Pressure drop and density distributions along the bundle length were disused.Williams et al. [17] examined two-phase flow patterns experimentally.Bubbly, annular, and slug flow were analyzed at constant pressure in a vertical flow rod bundle.Ren et al. [18] studied the impact of grid spacer and mixing vane on the phase distribution within rod bundle assembly.Air and water were employed as two-phase fluids.Also, experimental results explained the void fraction distributions in the subchannels.
Drucker et al. [19] explained the heat transfer mechanism in two-phase fluid flow in a tube and fuel rod assembly.Analyses of heat transfer were performed for different ranges of Reynolds numbers.Air and water were employed as two-phase fluids inside the domain, and void fraction distributions were discussed in detail.Li et al. [20] conducted experiments to analyze the heat transfer performance of two-phase fluid flow in a vertical rod bundle assembly.They found that for their test conditions, increasing the system pressure or mass flux will reduce the wall surface superheating degree.Mixing vane effects were studied numerically on the fuel rod assembly for two-phase fluid flow by Xiao et al. [21].CFD results validation against experimental data was performed.They analyzed the flow and thermal characteristics within the sub-channels of the rod bundle.The authors found that the cross-flow plays an important role in the accumulation of vapor on the fuel rod in the domain.
Yang et al. [22] investigated the impacts of spacer grids, including plastic protrusions, on air-water (two-phase) flow within assembly of an 8 × 8 array.The downstream turbulence behavior is stronger, and the low-pressure zone is visible behind the grid spacer.Tian et al. [23] conducted tests to measure the pressure drop within rod bundle assembly for a two-phase (air-water) flow.They discussed the turbulent intensity and velocity distributions in the domain.Also, they developed a correlation for a two-phase friction multiplier and tested it successfully.
Hosokawa et al. [24] examined flow behaviors of two phases like liquid velocity, void fraction, and bubble velocity in rod bundle subchannels.The work was done on a lattice with a square shape assembly of 4 × 4 array having 10-mm rod diameter, 9.1-mm hydraulic diameter, and 1.25 pitch-to-rod diameter ratio, and the authors reported that the acting lift force across bubbles in rod bundle subchannels affects void fraction distributions significantly.Inoue et al. [25] experimentally analyzed the distributions of void fraction within rod bundle.Measurements were performed for normal operating conditions of boiling water reactor rod bundles.An X-ray CT scanner was used to measure the void fractions in the domain.
Based on the aforementioned study of literature, it is inferred that the inclusion of a spacer or grid spacer with vanes has a positive impact on the rate of turbulence mixing and heat transfer improvement.This effect is particularly notable at the spacer location and in the downstream zone of the spacer.Despite the existence of several studies on spacer analysis, the precise forecast of thermal-hydraulic characteristics using computational analysis remains an ongoing area of research for two-phase.This is mainly significant for the design of spacers with varying vane shapes.Moreover, the detailed comparative analysis for the impacts of spacer-vane and spacer alone on thermalhydraulic performance of two-phase flow in annular channel has received little attention in the literature.Most investigations have been conducted in adiabatic flow conditions.This study presents a novel approach to predicting the thermal-hydraulic performance of two-phase annulus flow with a spacer-vane.The present study examines the impacts of spacer-vane on thermal-hydraulic performance of two-phase flow in annular channels using computational CFD analysis.The Eulerian model was employed to solve the two-phase fluid flow within a three-dimensional annular channel model.The SST k-ω mixture turbulence model was utilized for two-phase turbulent flow simulations.Simulations were achieved in ANSYS Fluent for the flow considerations of mass flux 714.4 kg/ m 2 s, liquid temperature 330 K, volume fraction 0.35, heat flux 197.2 kW/m 2 , and atmospheric pressure condition.A detailed investigation was conducted to examine the impacts of spacer-vane on thermal-hydraulic performance of two-phase annulus flow.

Mathematical equations
The ANSYS Fluent multiphase Eulerian model was used to solve two-phase fluid flow in three-dimensional models of annular channels.Water (liquid) and water vapor were considered as primary and secondary phases, respectively.The governing equations [26] were solved for the two-phase fluid flow during the simulations, which are given as follows:

Continuity equation of phase q
where − → V q = velocity of q th phase ṁpq = mass transfer from p th to q th phase ṁqp = mass transfer from q th to p th phase α q = volume fraction of q th phase S q = source term

Momentum equation of phase q
where τ q = stress-strain tensor of q th phase where µ q = shear viscosity of phase q q = bulk viscosity of phase q (1)

Energy equation of phase q
where h q = specific enthalpy of q th phase − → q q = heat flux Q pq = heat exchange intensity between p th and q th phases h pq = interphase enthalpy

Turbulence model
The equations of SST k-ω mixture turbulence model [26] are given as follows: where mixture density, The production of turbulence kinetic energy; G k,m

Drag force
The interphase momentum transfer between two phases due to the drag force can be obtained via the drag model of Schiller and Naumann model [27].The drag function f is given as follows: where coefficient of drag (C D ) is given as follows: The relative Reynolds number for primary phase q and secondary phase p is obtained as follows: (4) where d p is the bubble diameter of phase p.

Lift force
The lift force acting on a secondary phase p in a primary phase q can be calculated as follows: where ρ q = primary phase density α p = secondary phase volume fraction − → V q and − → V p = velocity of primary and secondary phase respectively C l = lift coefficient Tomiyama lift force model [28] has been used in the simulations to solve the lift coefficient (C l ) , which is given as follows: where Eo′ is a modified Eotvos number based on the long axis of the deformable bubble, d h where σ = surface tension g = gravity d p = bubble diameter

Wall lubrication force
Wall lubrication force proposed by Antal et al. [29] is given as follows: (10) where − → V q − − → V p = relative velocity between phases − → n w = unit normal from wall C wl = wall lubrication coefficient Wall lubrication coefficient (C wl ) proposed by Antal et al. [29] is given as follows: where C w1 = −0.01 and C w2 = 0.05 are dimensionless coefficients.
y w = nearest wall distance

Turbulent dispersion force
The turbulence-induced dispersion force proposed by Lopez de Bertodano model [30] is given as follows: where k q = primary phase turbulent kinetic energy ∇α p = secondary phase volume fraction gradient

Virtual mass force:
The equation for virtual or added mass force [26] is given as follows: where d q dt = time derivative of phase material

Geometry development and boundary details
CFD model of an upward flow annular channel with spacer-vane has been developed.Additionally, an annular channel model with a spacer but without vane has been generated and modeled to study the impacts of spacer-vane in a two-phase flow.Figure 1 a, b, c, and d shows the geometry details of a three-dimensional annular channel with spacer-vane.The blockage ratio of spacer located at the center of the annular channel is assumed to be 0.4.Vanes have been positioned on the spacer tip in the flow direction, and the vane dimensions are shown in Fig. 2. The vane angle is taken as 7° with the flow direction.Table 1 represents the dimensions of simulation geometry.
In the present work, an Eulerian model was used for two-phase fluid (water and water vapor mixture) formulation.Transient simulations were performed with time step size of 0.005 s and total number of time steps of 600.Temperature, hydraulic diameter, volume fraction, mass flux, and turbulent intensity were applied at the domain inlet.The saturation temperature of water has been set as 373.15K at atmospheric pressure.Constant pressure was employed at the domain outlet.Adiabatic and constant heat flux criteria were employed on the inner and outer pipes, respectively.( 18) Smooth wall with no-slip criteria was employed on the walls for momentum equation.Table 2 provides the details of boundary condition.

Mesh model and mesh independence
Mesh independence analysis was executed employing three distinct mesh models.Meshes 1, 2, and 3 possess element counts of 1.9, 1.6, and 1.3 million, respectively.The unstructured tetrahedron mesh models were generated employing ICEM CFD.Mesh was generated by specifying the size of global element with a fixed scaling factor.Mesh refinement was accomplished by retaining distinct sizes of elements on the vane, spacer, and inner-outer pipe surfaces.In order to effectively represent the heat transfer phenomenon, prismatic layers of elements were built close to pipe walls.The initial layer's height on the pipe surface was kept low enough; thus, y + is smaller than 1.An expansion ratio of 1.2 was employed among layers.Figure 3    surface and volume mesh over the inner pipe with spacer-vane and at the annulus gap, respectively.Mesh independence analysis was executed in an annular channel with spacer-vane.The dimensions of the annular channel and flow parameters for mesh independence analysis are given in Tables 1 and 2, respectively.Eulerian model was used to solve the two-phase fluid (water and water vapor) flow.The simulation was accomplished using SST k-ω mixture turbulence model.The tangential velocity distributions of liquid in radial coordinates for three mesh models are displayed in Fig. 4. The tangential velocities were taken at 460 mm from channel inlet, which is over 40 times of the hydraulic diameter to avoid entrance effect [31].Hydraulic diameter of the annular channel is 6 mm.Based on the comparison made for the tangential velocity distributions (Fig. 4), mesh 2 was adopted for further numerical analysis as the variation of tangential velocities is less than 0.3% compared to the other two mesh models, i.e., meshes 1 and 3 (it is per the authors consideration).Table 3 shows the numerical percentage deviation of meshes 1 and 3 with respect to mesh 2.

Methods
The multiphase Eulerian model was applied to solve two-phase (water -primary phase and vapor -secondary phase) fluid flow in a three-dimensional model of annular channel.The time step size was taken as 0.005 s for the present transient simulations.The phase-coupled SIMPLE algorithm was employed to couple the velocity and pressure.The finite volume method is the basis of this algorithm.The value of the residuals was specified as 10 −6 to get converged results for the turbulence model and governing equa- tions.The residuals were achieved using a first-order upwind discretization scheme.The equation for volume fraction was solved by the QUICK scheme.The precision of the results is confirmed by these algorithms and the criteria for residuals to determine the  convergence of the solution.During simulation, the values for under-relaxation were reduced for better convergence of solutions.The convergence was verified with the flux report of net mass flow rate from inlet and outlet boundaries of the domain.Also, the gravity influence was specified for the accuracy of numerical results.

CFD results validation and adoption of turbulence model
The simulation was accomplished in ANSYS Fluent for CFD results validation.Data from an experiment executed by Serizawa et al. [1] was used for CFD results validation.Model dimensions and boundary parameters for CFD validation are identical to an experiment executed by Serizawa et al. [1].The CFD model used for the validation is shown in Fig. 5.It is a vertical pipe, and the model dimensions are as follows: the diameter and length are 0.06 m and 2.15 m, respectively.The simulation was performed for two-phase fluid (water-air) flow with flow parameters of liquid superficial velocity as 1.03 m/s, gas velocity as 0.151 m/s, bubble diameter as 4 mm, and atmospheric pressure condition.The turbulence models, such as standard k-ε and SST k-ω, were examined to analyze flow performance of this CFD model.Figure 6a and b shows the comparative results for radial distributions of void fraction and liquid velocity obtained from CFD and the experiment, respectively.The radial distributions of void fraction and liquid velocity were taken on line L1, which is at 1.8 m from pipe inlet (Fig. 5). Figure 6a shows that the peak of void fraction is close to the pipe wall because of the lift force acting on gaseous phase and the higher bubble diameter.Figure 6b shows that the liquid velocity is maximum at the center of the pipe and zero at the pipe wall.Tables 4 and 5 represent the uncertainty in CFD results for void fraction and liquid velocity distributions compared with experiment, respectively.It is obvious from Fig. 6a and b that the CFD results of void fraction and liquid velocity determined by turbulence model of SST k-ω were a better match to experimental data than the standard k-ε.Thus, the impacts of spacer-vane were analyzed by adopting the SST k-ω for subsequent simulations of CFD in the annular channel.Moreover; several researchers have found that the SST k-ω is more appropriate for thermal-hydraulic analysis within nuclear rod bundle assemblies compared to other turbulence models [32][33][34][35][36][37][38].Based on the previous research works conducted by several researchers [32][33][34][35][36][37][38] for flow thermal characteristics analysis in nuclear fuel channels, the authors examined the SST k-ω and the standard k-ε in the present work.
The SST k-ω turbulence model predicted more accurately than the standard k-ε model.

Results and discussion
In this section, the thermal-hydraulic performance of two-phase flow in annular channels with spacer-vane was discussed.Simulation results were produced for the flow considerations of mass flux 714.4 kg/m 2 .s,liquid temperature 330 K, volume fraction 0.35, heat flux 197.2 kW/m 2 , and at atmospheric pressure condition.The parametric plot for convergence of the present simulation in an annular channel with spacer-vane is shown in Fig. 7.

Effect of spacer-vane on velocity and pressure distributions
The tangential and radial velocity distributions of liquid in radial coordinates have been presented in Figs. 8 and 9, respectively.The velocity distributions and distance in radial coordinates are shown on the ordinate and abscissa, respectively.The velocity distributions have been taken at different sites (X/D = 1, 2, 5, 20, 30, 40, and 50) in downstream from the spacer tip.X and D represent the downstream distance from spacer tip and a hydraulic diameter, respectively.Figure 8 represents that the tangential velocity magnitude is significantly improved near the wall of outer pipe in the downstream area near the spacer tip (at X/D = 1).This improved tangential velocity is because of the spacer influence as it decreases the area of flow.The tangential velocity gradually reduces adjacent to the outer pipe wall further downstream from the spacer tip (i.e., X/D = 2 and 5) due to the recovery of flow area [36].It is obvious from Fig. 8 that from X/D = 20 onwards downstream, the flow has become fully developed since the variations in the flow velocity profiles are negligible along the flow direction and the spacer effects are negligible.Also, the turbulence induced by the spacer-vane in the flow greatly dies down from X/D = 20 onwards due to the axial decay of vortexes, which appeared due to the presence of the spacer-vane in the flow domain.Similarly, Fig. 9 shows that the radial velocity magnitude is maximum near the wall of outer pipe which is around 0.8 m/s for location of X/D = 1.Further downstream up to X/D = 5, the radial velocity reduces gradually near the outer pipe wall due to the recovery of flow area, and X/D = 20 onwards in downstream of the spacer effects disappear [36].
Figure 10 a and b represents the axial velocity contour of liquid for spacer only (without vane) and spacer-vane in annular channel, respectively.The velocity contour has been taken on the midplane of the channel annulus gap.The spacer serves the purpose of decreasing domain flow area by acting as a flow blockage.Decreased domain flow area increases flow velocity.The spacer's vane further increases flow velocity considerably in the downstream area near the spacer.Figure 10 clearly shows that the maximum velocity magnitude exceeds 1.3 m/s with spacer-vane, which is higher than that of spacer alone in the annular channel.Enhanced flow velocity contributes significantly to reducing the channel wall temperature, which enhances heat transfer performance.
Figure 11 shows the axial velocity distribution of liquid in radial coordinates of the annular channel with spacer-vane.Figure 11 shows that the axial velocity magnitude is maximum near the wall of outer pipe which is around 1.1 m/s for location of X/D = 1.This enhanced velocity (for X/D = 1) is because of the spacer influence as it decreases the area of flow.Due to the recovery of flow area, the velocity gradually reduces near the wall of outer pipe in downstream from spacer tip (i.e., X/D = 2 and 5).Further downstream, the flow has become fully developed [36].
Figure 12 shows the dynamic pressure distribution of liquid.The dynamic pressure magnitude is maximum near the outer pipe wall which is about 650 Pa for location    The maximum turbulent intensity is observed in the core zone of the annulus gap, which is over 16%.This enhanced turbulent intensity in the annulus gap is about 11% higher than the turbulent intensity specified at the inlet of channel.The magnitude of turbulent intensity is observed to be greater using spacer-vane compared to spacer alone in downstream of the channel.Additionally, the effect of vane on turbulent intensity is noticed up to the farthest downstream from the spacer tip compared to spacer alone.This prolonged effect of turbulent intensity in the annulus gap enhances the heat transfer rate in the channel.Figure 14 shows the turbulent intensity distribution of the mixture in radial coordinates.The distribution has been taken at different sites (X/D = 1, 2, 5, 20, 30, 40, and 50) in downstream from the spacer tip. Figure 14 clearly shows that near the spacer tip (X/D = 1, 2, and 5) downstream, the turbulent intensity magnitudes are higher at the core zone [39] (at radial distance = 11 mm) of the annulus gap, and the turbulent intensities are again improved significantly near the wall of outer pipe due to spacer influence.The spacer provides a flow blockage and yields vortexes behind the spacer downstream.These vortexes increase the turbulent mixing in the core zone of the annulus gap, and thus, the turbulent intensity magnitudes are higher at the core zone of the annulus.Maximum turbulent intensity is observed about 14% at the core zone for location of X/D = 2. Also, Fig. 13b clearly shows that the turbulent intensity magnitude exceeds 16% in the core zone of the annulus gap for location of X/D = 3.Further, from X/D = 20 onwards downstream, the turbulent intensity variations are noticed as negligible.and spacer-vane in annular channel, respectively.The contour has been taken on the midplane which is created in annulus gap along the channel length.Figure 15 a and b clearly shows that near the spacer tip downstream, the temperature magnitude is less near the wall of outer pipe using spacer-vane compared to spacer alone in the channel.The spacer's vane further enhances the turbulent mixing in the downstream area near spacer; thus, the drop in temperature is observed to be greater using spacer-vane compared to spacer alone [12].In the downstream area near the spacer tip, the temperature magnitude is observed in the range of 365 to 373 K near the wall of outer pipe using spacer-vane.Hence, it is obvious that the spacer's vane further reduces the wall temperature significantly near the spacer downstream.
Figure 16 shows the temperature distribution of liquid in annular channel with spacervane.The temperature magnitude is minimum near the wall of outer pipe, which is about 365 K for location of X/D = 1.Due to the formation of vortexes, the turbulent mixing is predominant in the downstream area near the spacer.Turbulent mixing of fluid enhances the turbulence kinetic energy, which causes a reduction in temperature near the spacer downstream [12].Further downstream, a gradual increase in temperature is observed due to the axial decay of vortexes [40].In downstream, the temperature magnitude increases gradually near the wall of outer pipe and is observed as 375 K for X/D = 5 and about 395 K for X/D = 50.

Effect of spacer-vane on volume fraction
Figure 17 represents the volume fraction distribution of water vapor in radial coordinates of the annular channel with spacer-vane.Figure 17 shows that the maximum volume fraction of vapor is about 0.44 (near the spacer tip, i.e., at X/D = 1 to 5) at a radial distance 11.7 mm from inner pipe wall and then gradually decreases to a minimum at the outer pipe wall.When the flow has fully developed, i.e., X/D = 20 onwards, the maximum volume fraction (about 0.45 at X/D = 50) is noticed at the center zone of the annulus gap (i.e., at a radial distance of 11 mm).It is observed that in the farthest downstream, most of the water vapor is present at the core of the annulus gap.The peak of the volume fraction is observed at the core of the annulus gap due to the higher water vapor input (0.35) at the inlet [1] and the small size (diameter 0.01 mm) of the water vapor

Conclusions
A detailed investigation was conducted to examine the impacts of spacer-vane on thermal-hydraulic performance of two-phase annulus flow.Additionally, an annular channel model with a spacer but without vane was generated and modeled in order to examine the impact of the spacer-vane in a two-phase flow.The multiphase Eulerian model was applied to solve two-phase (water -primary phase and vapor -secondary phase) fluid flow in a three-dimensional model of annular channel.Simulation results were produced for the flow considerations of mass flux 714.4 kg/m 2 .s,liquid temperature 330 K, volume fraction 0.35, heat flux 197.2 kW/m 2 , and at atmospheric pressure condition.The simulation was accomplished using SST k-ω mixture turbulence model.The impacts of spacer-vane on thermal-hydraulic performance of two-phase annulus flow were studied, and the conclusions are as follows: • Maximum turbulent intensity is observed in the core zone of annulus gap, which is over 16%.The augmented turbulent intensity in the annulus gap is approximately 11% higher than the turbulent intensity indicated at the channel's inlet.• The velocity magnitude near the wall of outer pipe is significantly improved in the downstream area near the spacer tip.This improved velocity is because of the spacer influence as it decreases the area of flow.• Enhanced flow velocity plays a significant role in reducing the wall temperature of channel, and hence, it enhances the heat transfer performance.• Due to the small size of water vapor used in the flow field, most of the water vapor is present at the core of the annulus gap in farthest downstream of the channel.The volume fraction of water vapor is observed as minimum at the outer pipe wall, which is favorable to increase the heat transfer rate from the channel wall.

Page 2 of 26
Dhurandhar et al.Journal of Engineering and Applied Science (2024) 71:127 represents a mesh model for the geometry used in the present analysis.Figures 3 a and b represents the

Fig. 1 Fig. 2
Fig. 1 Geometry details. a A three-dimensional annular channel, b an annular channel with its boundaries, c channel with spacer-vane, and d vane on the spacer tip

Fig. 3
Fig. 3 Mesh model.a Surface mesh on inner pipe with spacer-vane and b volume mesh at annulus gap

Fig. 4
Fig. 4 Mesh independence analysis for annular channel with spacer-vane

Fig. 6 a
Fig. 6 a Void fraction distribution-CFD b Liquid velocity distribution-CFD validation

of X/D = 1 .
The dynamic pressure increased near the wall of outer pipe due to the augmented flow velocity near spacer downstream.The difference in dynamic pressure between the locations X/D = 1 and 2 is observed as 250 Pa at radial distance of 12 mm (near the wall of outer pipe).Further downstream, up to X/D = 5, the dynamic pressure reduces gradually near the outer pipe wall, and then the spacer effects are observed to be negligible (X/D = 20 onwards).

Fig. 7 Fig. 8
Fig. 7 Convergence plot for the present simulation

Fig. 9
Fig. 9 Radial velocity distribution of liquid for spacer-vane annular channel

Fig. 10 Fig. 11
Fig. 10 Axial velocity contour of liquid a spacer only and b spacer-vane

Fig. 12
Fig. 12 Dynamic pressure distribution of liquid for spacer-vane annular channel

Fig. 14
Fig. 14 Turbulent intensity distribution of the mixture for spacer-vane annular channel

Fig. 16 Fig. 17
Fig. 16 Temperature distribution of liquid for spacer-vane annular channel

Table 1
Dimensions of simulation geometry

Table 2
Details of boundary condition

Table 3
Numerical percentage deviation of mesh

Table 4
Uncertainty in CFD results of volume fraction with experiment

Table 5
Uncertainty in CFD results of liquid velocity with experiment