Steady flow and interfacial shapes of a highly viscous dispersed phase

A perturbation theory for the steady flow of immiscible liquids is develope d when the dispersed phase is much more viscous than the continuous phase, as is the case in emulsions of highly viscous bitumen in water and in water lubricated pipelines of heavy crude. The perturbation is nonsingular, but nonstandard; the partitioning of the boundary conditions at different order s is not conventional. At zero-th order the dispersed phase moves as a rigi d solid with an as yet unknown, to-be-determined, pressure. The flow of the continuous phase at zero-th order is determined by a Dirichlet problem wit h prescribed velocities on a to-be-iterated interfacial boundary. The first order problem in the dispersed phase is determined from the solution of a Stokes flow problem driven by the previously determined shear strain on the as yet undetermined interfacial boundary. This Stokes problem determines t he unknown, to-be-determined, lowest order pressure distribution. At this p oint we have enough information to test the balance of normal stresses at l owest order; by iterating the interface shapes we may now complete the desc ription of the lowest order problems. The perturbation sequence in powers o f the viscosity ratio has a similar structure at every order and all the pr oblems may be solved sequentially with the caveat that interface shape must be determined iteratively in each perturbation loop. A perturbation solution for the wavy interfacial shapes on core-annular flo ws of very viscous oil is presented and the results are compared with exper iments and a simpler approximation in which the core moves as rigid, but de formable body with no secondary motions. The perturbation theory gives rise to an accurate description of the bamboo waves observed in experiments whe n the holdup ratio measured in the experiments is assumed in the theoretica l calculation. The perturbation solution and the rigid body approximation a re in a relatively good agreement with errors of the order 10% in the flow curves and wave shapes; the error is associated with the neglect of the sec ondary motion in the rigid-deformable core approximation. (C) 2000 Elsevier Science Ltd. All rights reserved.