The instability of a moving viscous drop

1990 ◽  
Vol 210 ◽  
pp. 1-21 ◽  
Author(s):  
C. Pozrikidis
Keyword(s):  
Surface Tension ◽  
Thin Filament ◽  
Viscosity Ratio ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Fredholm Integral Equations ◽  
Spherical Drop ◽  
Viscous Drop ◽  
Secondary Mechanism ◽  
Axisymmetric Perturbations

The deformation of a moving spherical viscous drop subject to axisymmetric perturbations is considered. The problem is formulated using two different variations of the boundary integral method for Stokes flow, one due to Rallison & Acrivos, and the other based on an interfacial distribution of Stokeslets. An iterative method for solving the resulting Fredholm integral equations of the second kind is developed, and is implemented for the case of axisymmetric motion. It is shown that in the absence of surface tension, a moving spherical drop is unstable. Prolate perturbations cause the ejection of a tail from the rear of the drop, and the entrainment of a thin filament of ambient fluid into the drop. Oblate perturbations cause the drop to develop into a nearly steady ring. The viscosity ratio plays an important role in determining the timescale and the detailed pattern of deformation. Filamentation of the drop emerges as a persistent but secondary mechanism of evolution for both prolate and oblate perturbations. Surface tension is not capable of suppressing the growth of perturbations of sufficiently large amplitude, but is capable of preventing filamentation.

scholarly journals Nonlinear two-dimensional free surface solutions of flow exiting a pipe and impacting a wedge

2021 ◽  
Vol 126 (1) ◽  
Author(s):  
Alex Doak ◽  
Jean-Marc Vanden-Broeck
Keyword(s):  
Surface Tension ◽  
Integral Equation ◽  
Free Surface ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Two Dimensional ◽  
Numerical Schemes ◽  
Additional Advantage ◽  
Infinite Wedge ◽  
Good Agreement

AbstractThis paper concerns the flow of fluid exiting a two-dimensional pipe and impacting an infinite wedge. Where the flow leaves the pipe there is a free surface between the fluid and a passive gas. The model is a generalisation of both plane bubbles and flow impacting a flat plate. In the absence of gravity and surface tension, an exact free streamline solution is derived. We also construct two numerical schemes to compute solutions with the inclusion of surface tension and gravity. The first method involves mapping the flow to the lower half-plane, where an integral equation concerning only boundary values is derived. This integral equation is solved numerically. The second method involves conformally mapping the flow domain onto a unit disc in the s-plane. The unknowns are then expressed as a power series in s. The series is truncated, and the coefficients are solved numerically. The boundary integral method has the additional advantage that it allows for solutions with waves in the far-field, as discussed later. Good agreement between the two numerical methods and the exact free streamline solution provides a check on the numerical schemes.

Theory and experiment on the low-Reynolds-number expansion and contraction of a bubble pinned at a submerged tube tip

1998 ◽  
Vol 356 ◽  
pp. 93-124 ◽  
Author(s):  
HARRIS WONG ◽  
DAVID RUMSCHITZKI ◽  
CHARLES MALDARELLI
Keyword(s):  
Surface Tension ◽  
Reynolds Number ◽  
Numerical Solutions ◽  
Full Range ◽  
Boundary Integral ◽  
Gas Pressure ◽  
Boundary Integral Method ◽  
Liquid Density ◽  
Gravitational Acceleration ◽  
Dynamic Surface

The expansion and contraction of a bubble pinned at a submerged tube tip and driven by constant gas flow rate Q are studied both theoretically and experimentally for Reynolds number Re[Lt ]1. Bubble shape, gas pressure, surface velocities, and extrapolated detached bubble volume are determined by a boundary integral method for various Bond (Bo=ρga2/σ) and capillary (Ca=μQ/σa2) numbers, where a is the capillary radius, ρ and μ are the liquid density and viscosity, σ is the surface tension, and g is the gravitational acceleration.Bubble expansion from a flat interface to near detachment is simulated for a full range of Ca (0.01–100) and Bo (0.01–0.5). The maximum gas pressure is found to vary almost linearly with Ca for 0.01[les ]Ca[les ]100. This correlation allows the maximum bubble pressure method for measuring dynamic surface tension to be extended to viscous liquids. Simulated detached bubble volumes approach static values for Ca[Lt ]1, and asymptote as Q3/4 for Ca[Gt ]1, in agreement with analytic predictions. In the limit Ca→0, two singular time domains are identified near the beginning and the end of bubble growth during which viscous and capillary forces become comparable.Expansion and contraction experiments were conducted using a viscous silicone oil. Digitized video images of deforming bubbles compare well with numerical solutions. It is observed that a bubble contracting at high Ca snaps off.

Numerical studies of cusp formation at fluid interfaces in Stokes flow

1998 ◽  
Vol 357 ◽  
pp. 29-57 ◽  
Author(s):  
C. POZRIKIDIS
Keyword(s):  
Surface Tension ◽  
Stokes Flow ◽  
Boundary Integral ◽  
Numerical Results ◽  
Boundary Integral Method ◽  
Fluid Interfaces ◽  
Two Dimensional ◽  
Numerical Studies ◽  
Insoluble Surfactant ◽  
Cusp Formation

Numerical studies are performed addressing the development of regions of high curvature and the spontaneous occurrence of cusped interfacial shapes in two-dimensional and axisymmetric Stokes flow. In the numerical simulations, the velocity field is computed using a boundary-integral method, and the evolution of the concentration of an insoluble surfactant over an evolving interface is computed using an implicit finite-volume method. Three configurations are considered in detail, and the results are used to elucidate three different aspects of cusp formation. In the first series, the deformation of a two-dimensional bubble immersed in a family of straining flows devised by Antanovskii, and of an axisymmetric bubble immersed in an analogous family of flows devised by Sherwood, are examined. The numerical results indicate that highly elongated and cusped two-dimensional shapes, and pointed or cusped axisymmetric shapes, are unstable and should not be expected to occur in practice. In the second series of studies, the role of an insoluble surfactant on the transient deformation of bubbles subject to the Antanovskii or Sherwood flow is investigated. Under certain conditions, the reduced surface tension at the tips raises the local curvature to high values and causes the ejection of a sheet or column of gas by means of tip streaming. In the third series of studies, the coalescence of a polygonal formation of five viscous columns of a fluid placed in an arrangement that differs only slightly from one proposed recently by Richardson is examined. The numerical results confirm Richardson's predictions that transient cusps may occur at a finite time in the presence of surface tension. The underlying physical mechanism is discussed on the basis of reversibility of surface-driven Stokes flow and with reference to the regularity of the motion driven by negative surface tension. Replacing the inviscid ambient gas with a slightly viscous fluid whose viscosity is as low as one hundredth the viscosity of the cylinders suppresses the cusp formation.

Highly nonlinear standing water waves with small capillary effect

1998 ◽  
Vol 369 ◽  
pp. 253-272 ◽  
Author(s):  
WILLIAM W. SCHULTZ ◽  
JEAN-MARC VANDEN-BROECK ◽  
LEI JIANG ◽  
MARC PERLIN
Keyword(s):  
Surface Tension ◽  
Water Waves ◽  
Standing Waves ◽  
High Amplitude ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Ripple Formation ◽  
Node Distribution ◽  
Highly Nonlinear ◽  
Periodic Standing Waves

We calculate spatially and temporally periodic standing waves using a spectral boundary integral method combined with Newton iteration. When surface tension is neglected, the non-monotonic behaviour of global wave properties agrees with previous computations by Mercer & Roberts (1992). New accurate results near the limiting form of gravity waves are obtained by using a non-uniform node distribution. It is shown that the crest angle is smaller than 90° at the largest calculated crest curvature. When a small amount of surface tension is included, the crest form is changed significantly. It is necessary to include surface tension to numerically reproduce the steep standing waves in Taylor's (1953) experiments. Faraday-wave experiments in a large-aspect-ratio rectangular container agree with our computations. This is the first time such high-amplitude, periodic waves appear to have been observed in laboratory conditions. Ripple formation and temporal symmetry breaking in the experiments are discussed.

The flow of a liquid film along a periodic wall

1988 ◽  
Vol 188 ◽  
pp. 275-300 ◽  
Author(s):  
C. Pozrikidis
Keyword(s):  
Surface Tension ◽  
Liquid Film ◽  
Numerical Procedure ◽  
Surface Profile ◽  
Rotating Disk ◽  
Boundary Integral ◽  
Creeping Flow ◽  
Boundary Integral Method ◽  
Free Surface Profile ◽  
Gravity Driven Flow

The creeping flow of a liquid film along an inclined periodic wall of arbitrary geometry is considered. The problem is formulated using the boundary-integral method for Stokes flow. This method is extended to two-dimensional flows involving free surfaces, and is implemented in an iterative numerical procedure. Detailed calculations for flow along a sinusoidal wall are perfomed. The free-surface profile is studied as a function of flow rate, inclination angle, wave amplitude, and surface tension, and is compared with previous asymptotic solutions. The results include streamline patterns, velocity profiles and wall-shear-stress distributions, and establish criteria for flow reversal. For specified wall geometry, the asymptotic behaviour for very small flow rates is shown to be a strong function of surface tension. It is demonstrated that these results are valid in a qualitative sense for general wall geometries. The analogy between gravity-driven flow and the flow of a liquid layer on a rotating disk (spin coating) is also discussed.

The effect of surfactant on bursting gas bubbles

1995 ◽  
Vol 302 ◽  
pp. 231-257 ◽  
Author(s):  
Jeremy M. Boulton-Stone
Keyword(s):  
Surface Tension ◽  
Free Surface ◽  
Stress Condition ◽  
Numerical Technique ◽  
Gas Bubbles ◽  
Free Surface Flows ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Surface Shear ◽  
Surface Compression

A numerical technique, based on the boundary integral method, is developed to allow the modelling of unsteady free-surface flows at large Reynolds numbers in cases where the surface is contaminated by some surface-active compound. This requires the method to take account of the tangential stress condition at the interface and is achieved through a boundary-layer analysis. The constitutive relation that forms the surface stress condition is assumed to be of the Boussinesq type and allows the incorporation of surface shear and dilatational viscous forces as well as Marangoni effects due to gradients in surface tension. Sorption kinetics can be included in the model, allowing calculations for both soluble and insolube surfactants. Application of the numerical model to the problem of bursting gas bubbles at a free surface shows the greatest effect to be due to surface dilatational viscosity which drastically reduces the amount of surface compression and can slow and even prevent the information of a liquid jet. Surface tension gradients give dilatational elasticity to the surface and thus also significantly prevent surface compression. Surface shear viscosity has a smaller effect on the interface motion but results in initially increased surface concentrations due to the sweeping up of surface particles ahead of the inward-moving surface wave.

Influence of internal viscosity on the large deformation and buckling of a spherical capsule in a simple shear flow

2011 ◽  
Vol 672 ◽  
pp. 477-486 ◽  
Author(s):  
É. FOESSEL ◽  
J. WALTER ◽  
A.-V. SALSAC ◽  
D. BARTHÈS-BIESEL
Keyword(s):  
Shear Flow ◽  
Simple Shear ◽  
Viscosity Ratio ◽  
Fluid Motion ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Simple Shear Flow ◽  
Flow Strength ◽  
Capsule Deformation ◽  
Internal Viscosity

The motion and deformation of a spherical elastic capsule freely suspended in a simple shear flow is studied numerically, focusing on the effect of the internal-to-external viscosity ratio. The three-dimensional fluid–structure interactions are modelled coupling a boundary integral method (for the internal and external fluid motion) with a finite element method (for the membrane deformation). For low viscosity ratios, the internal viscosity affect the capsule deformation. Conversely, for large viscosity ratios, the slowing effect of the internal motion lowers the overall capsule deformation; the deformation is asymptotically independent of the flow strength and membrane behaviour. An important result is that increasing the internal viscosity leads to membrane compression and possibly buckling. Above a critical value of the viscosity ratio, compression zones are found on the capsule membrane for all flow strengths. This shows that very viscous capsules tend to buckle easily.

Sign in / Sign up

Export Citation Format

Share Document