AEM Home page> People > Faculty > Daniel D. Joseph> Archive on Irrotational Motions of Viscous and Viscoelastic Fluids

Rise velocity of a spherical cap bubble


University of Minnesota, Aerospace Engineering and Mechanics, Minneapolis, 110 Union St. SE, Minneapolis,
MN 55455, USA

(Received 23 October 2002 and in revised form 26 February 2003)


Potential flows of incompressible fluids admit a pressure (Bernoulli) equation when the divergence of the stress is a gradient as in inviscid fluids, viscous fluids, linear viscoelastic fluids and second-order fluids. We show that in potential flow without boundary layers the equation balancing drag and acceleration is the same for all these fluids, independent of the viscosity or any viscoelastic parameter, and that the drag is zero when the flow is steady. But, if the potential flow is viewed as an approximation to the actual flow field, the unsteady drag on bubbles in a viscous (and possibly in a viscoelastic) fluid may be approximated by evaluating the dissipation integral of the approximating potential flow because the neglected dissipation in the vorticity layer at the traction-free boundary of the bubble gets smaller as the Reynolds number is increased. Using the potential flow approximation, the actual drag on a spherical gas bubble of radius rising with velocity in a linear viscoelastic liquid of density and shear modules is estimated to be


and, in a second-order fluid,


where is the coefficient of the first normal stress and is the viscosity of the fluid. Because is negative, we see from this formula that the unsteady normal stresses oppose inertia; that is, oppose the acceleration reaction. When is slowly varying, the two formulae coincide. For steady flow, we obtain the approximate drag for both viscous and viscoelastic fluids. In the case where the dynamic contribution of the interior flow of the bubble cannot be ignored as in the case of liquid bubbles, the dissipation method gives an estimation of the rate of total kinetic energy of the flows instead of the drag. When the dynamic effect of the interior flow is negligible but the density is important, this formula for the rate of total kinetic energy leads to where is the density of the fluid (or air) inside the bubble and is the volume of the bubble.


Classical theorems of vorticity for potential flow of ideal fluids hold equally for second-order fluid. The drag and lift on two-dimensional bodies of arbitrary cross-section in a potential flow of second-order and linear viscoelastic fluids are the same as in potential flow of an inviscid fluid but the moment in a linear viscoelastic fluid is given by


where is the inviscid moment and is the circulation, and


in a second-order fluid. When is slowly varying, the two formulae for coincide. For steady flow, they reduce to


which is also the expression for in both steady and unsteady potential flow of a viscous fluid. Moreover, when there is no stream, this moment reduces to the actual moment on a rotating rod.


Potential flows of models of a viscoelastic fluid like Maxwell's are studied. These models do not admit potential flows unless the curl of the divergence of the extra stress vanishes. This leads to an over-determined system of equations for the components of the stress. Special potential flow solutions like uniform flow and simple extension satisfy these extra conditions automatically but other special solutions like the potential vortex can satisfy the equations for some models and not for others.