Science and Technology Archive

Lift Forces in Bubbly Flows

The dynamics of two-phase (gas/liquid) bubbly flows are complex: bubbles deform and disperse; large latent heats and heat capacity differentials influence local boiling; and relatively small changes in heated surface temperatures yield order of magnitude changes in boiling complexity.  Because the local void volume has a direct feedback effect on reactor neutron flux and fuel rod power production, prediction of local boiling rates and bulk boiling effects in nuclear reactors is key in achieving a higher fidelity prediction of the fuel system temperature distribution, local and  global power production, local fuel rod burnup, and CRUD deposition.

When the fuel rod surface temperature is high enough, steam bubbles nucleate and grow on the surface, detaching when the lift forces acting on the bubble are large     enough to overcome the suction force of the bubble on the rod.  Shear lift forces are reaction forces that result as the flowing water passes over the surface of the bubble. They evolve as the shape and size of the bubble evolve and are also influenced by asymmetries in the upstream flow field.

Figure 1 provides a general illustration of a spherical bubble in flow. The general equation describing the forces acting on the bubble can be written as [1]:

Where FI is the inertial force, FA is the added mass term, FL is the lift force, FD is the drag force, FB is the buoyancy force.  When these terms are expanded, the resolved forces that are the target of the work are more easily identified:

Where Ub is the bubble velocity, v is the fluid velocity, CM is the added mass coefficient, CL is the lift coefficient, CD is the drag coefficient, γ is the ratio of the densities (ρb/ρ), g denotes gravity and ω is vorticity.

In general, the coefficients are not constant, although their functional dependence has not been explicitly stated.  A few effects have been left out, including those due to time history, density gradients, and temperature gradients.  One of the lesser understood momentum closures is that due to the transverse forces.  These transverse forces are very critical for LWR applications as transverse motion of bubbles in the narrow and long flow channels between the fuel rods can have a huge impact on the heat transfer and neutron moderation.

Many widely used two-fluid Comp utational Multiphase Fluid Dynamics (CMFD) codes, such as Ansys' CFX & Fluent, CD-adapco’s Star-CCM+,  and RPI’s NPHAS E-CMFD, implement Auton's shear lift model [5] and Antal's wall force model [7] as the two primary lift forces.  All allow the user to tweak the model lift coefficients (CL), although by default the coefficients are constants; some codes include alternative correlations for the lift coefficients. While the base form of the Auton lift force with a constant lift coefficient is an acceptable approximation at high Reynolds number (Re) and low Strouhal number (Sr) flows about roughly spherical bubbles, for different bubble shapes and orientations, the Auton closure should be modified. Additionally, Auton's model derivation assumes a Sr number of 1, and this assumption isn’t always acceptable near the wall; when Sr is not small Auton's model overpredicts the shear lift force significantly. Substitute correlations presented in the past do include some bubble deformation and alternate orientations, but they don’t appear to be sufficient for the range of PWR flow conditions.

Thus, for CASL, the immediate need from a CMFD perspective is to include the effect of the wall on shear lift forces. To address this need, ORNL researchers Thomas Daly and Sreekanth Pannala and UT-K professor Art Ruggles collaborated in developing a new closure for lift force that includes the wall effects.  The team tested the closure using a comparison with fully developed turbulent pipe flow, and the new closure displays the correct asymptotic behavior. 

Proposed CASL Closure Form

In order to accumulate a large body of data for use in constructing the new closure, the team first did a survey of existing work and extracted the information necessary to calculate the primary dimensionless parameters (see Table 1 for a listing; note that the characteristic length is the bubble diameter, d).  The reference data and closures were plotted against the dimensionless parameters in Figure 6.  As can be seen upon inspection of the plots, the existing closures do not capture the effect of all the different parameters and are limited in range. For example , Rastello's model [3, 4] is good at moderate to large Re and away from the wall, but is not as accurate at very low Re.

Using the referenced closures as the starting point, the research team worked to develop a new closure that could    capture the empirically representative behavior of the lift coefficient versus Sr and Re. In order to better capture the influence of the wall, the proposed new closure includes an additional parameter, E, that renders the distance from the wall non-dimensional. The form was selected such that E=1 as the bubble touches the wall and approaches 0 as the bubble moves away from the wall.  Eight different variations of the closure were investigated, with the best match (based on a least-squares approach), called “DRP,” proposed for CASL adoption and future augmentation.

The new DRP lift closure was plotted with various independent parameters in Figures 6 and 7, and is shown with the reference data set.

The terms in the DRP lift closure that do not include the distance to the wall are based on Legendre & Magnaudet's model [2].  The new closure appears to capture both low and high Re behavior well. This was achieved by modifying Legendre's high Re model to allow for a local maximum near Re = 50. Like the other models it approaches the theoretical result of CL = 1/2 as Re approaches infinity. Unlike the other models discussed in this report, the proposed model also includes a dependence on distance from the wall. This allows it to remain physically reasonable as a bubble approaches the wall.

Next steps

The correlation developed represents the initial formulation of the closure.  Additional work is needed to provide a stronger theoretical footing and to include additional effects encountered in a PWR. The critical effects in increasing order of importance and difficulty that need to be accounted for in two-phase PWR conditions are:

¨ Blockage (reduced area on the side of the wall for flow)
¨ No-slip condition and the resulting velocity/vorticity profile
¨ Turbulence (especially for small bubbles when fluctuations can influence bubble transport) 
¨  Deformation (for large bubbles)
¨ Swarming (high bubble density, polydispersity, break-up and coalescence
¨ Boiling effects.