A two-dimensional elliptical droplet (density ratio 1/1000) is released in a large domain. Under the effect of surface-tension the shape of the droplet oscillates around its (circular) equilibrium shape. The fluids inside and outside the droplet are inviscid so ideally no damping of the oscillations should occur. As illustrated on Figure 89 some damping occurs in the simulation due to numerical dissipation.
This simulation is also a stringent test case of the accuracy of the surface tension representation as no explicit viscosity can damp eventual parasitic currents.
The initial shape of the droplet is given by:
|r(θ) = r0 + αcos(nθ)|
The oscillation frequency is then :
A comparison between the theoretical and numerical values of the frequency is given in Figure 90.
The amount of numerical damping can be estimated by computing an equivalent viscosity. With viscosity, kinetic energy is expected to decrease as:
where C is a constant, ν the viscosity and D the droplet diameter. Using curve fitting the damping coefficient b=Cν/D2 can be estimated (black curves on Figure 89). An equivalent Laplace number can then be computed as:
The equivalent Laplace number depends on spatial resolution as illustrated in Figure 91.