## Interciencia

##
*versión impresa* ISSN 0378-1844

### INCI v.34 n.8 Caracas ago. 2009

A simple model to describe dimple dynamics

Alejandro J. Acevedo-Malavé, Eloy Sira and Máximo García-Sucre

**Alejandro José Acevedo-Malavé**. Ph.D. in Physics, Venezuelan Institute for Scientific Research (IVIC), Venezuela. Researcher, IVIC, Venezuela. Address: Centro de Física Aplicada, Carretera Panamericana vía Jají, Finca El Tucuche, Mérida 5101, Venezuela. e-mail: alejanacev@gmail.com

**Eloy Sira**. Ph.D. in Physics, IVIC, Venezuela. Researcher, IVIC, Venezuela. e-mail: easira@cida.ve

**Máximo García-Sucre**. Máximo García-Sucre. Ph.D. in Physics, Paris University, France. Researcher, IVIC, Physics Center, Venezuela. e-mail: mgs@ivic.ve

**SUMMARY**

A model based on the hydrodynamics equations that allows to describe the dynamics of a dimple, once it has formed, is proposed. The Navier-Stokes equations are considered, and two fundamental approaches are used to simplify the mathematical treatment of the hydrodynamics equations. Certain conditions are considered that must be fulfilled at the interface, which serve to close the system of differential equations and lead to an evolution equation that describes the interfacial film dynamics. With the intention of solving this equation, the method of finite differences has been used.

**Un modelo sencillo para la descripción de la dinámica de un dimple **

**RESUMEN**

Se propone un modelo que permite la descripción de la dinámica de una depresión superficial (dimple) una vez que se ha formado. Las ecuaciones de Navier-Stokes son consideradas, y dos enfoques fundamentales son utilizados para simplificar el tratamiento matemático de las ecuaciones hidrodinámicas. Se consideraron ciertas condiciones que deben cumplirse en la interfase, las cuales sirven para completar el sistema de ecuaciones diferenciales y llevan a una ecuación de evolución que describe la dinámica de la película interfacial. A fin de resolver la ecuación se utilizó el método de diferencias finitas.

**Um modelo simples para a descrição da dinâmica de um dimple **

**RESUMO**

Propõe-se um modelo que permite a descrição da dinâmica de uma depressão superficial (dimple) após ter se formado. As equações de Navier-Stokes são consideradas, e dois enfoques fundamentais são utilizados para simplificar o tratamento matemático das equações hidrodinâmicas. Consideraram-se certas condições que devem cumprir-se na interfase, as quais servem para completar o sistema de equações diferenciais e conduzem a uma equação de evolução que descreve a dinâmica da película interfacial. A fim de resolver a equação se utilizou o método de diferenças finitas.

**Keywords / ** Coalescence / Dimple / Drops / Fluid Mechanics / Interfacial Film /

Received: 10/13/2008. Accepted: 07/10/2009.

Very different mechanisms are known to take place in the problem of emulsion stability (Kashchiev and Exerowa, 1980; Exerowa *et al*., 1983; Ivanov, 1988; Bibette, 1992; Bibette *et al*., 1992; Sonin *et al*., 1994; Kabalnov and Wennerström, 1996). Depending on the thickness of the interfacial film, diverse factors must be taken into account that have a marked influence on the physics of the problem. When thickness of the interfacial film is <300nm, the electrostatic interactions must be taken into account. Mathematically is necessary to include a disjoining pressure P¹0.

It has been reported that one of the most important factors within the stability problem of emulsions is the coalescence process, which is related to the stability of the interfacial film (Vrij, 1966; Vrij and Overbeeck, 1968; Ivanov *et al*., 1969; Sharma and Ruckenstein, 1987).

When drops approach each other in an emulsified system, the coalescence process begins. In this process, the liquid between the drops drains off until the two drops hit and a new one of a greater volume forms (Denkov *et al*., 1991; Kralchevsky *et al*., 1991; Tsekov and Radoev, 1992; Danov *et al*., 1993; Jaeger *et al*., 1994; Ivanov and Kralchevsky, 1997). This process can be divided into two stages: i) the drops approach each other without deformation until, at a certain distance between them a flat circular film appears, and ii) the thickness of this film begins to diminish from a certain separation, until it arrives to the critical separation and the two drops form a larger one.

Before coalescence occurs, a protuberance in the shape of a hole forms at the interfacial film, which depending on its dynamics can oscillate until disappearance or can coalesce, if the two drop surfaces are superposed (Ivanov, 1988; Velev *et al*., 1993; Hartland and Jeelani, 1994; Chesters and Bazhlekov, 2000; Yeo *et al*., 2001; Yeo *et al*., 2003). This protuberance is known as a dimple (Figure 1).

In this paper, a hydrodynamic model based on the quasi-static and lubrication approaches is proposed, allowing the simulation the dimple evolution.

Mathematical formulation of the model

The fundamental object of this study is to construct a model that allows to simulate the dimple dynamics. With this purpose, the Navier-Stokes equations for the interfacial film and the disperse phase were considered in the spirit of the quasi-static and the lubrication approximations.

The lubrication approximation is very useful to simplify the hydrodynamics equations, being applicable under the following conditions:

- The space between the two surfaces is small in comparison with the radius of the interface film (h(r,t) <<R, where h(r,t): thickness of the interface film and R: its radius).

- The inertial forces that act on the interface film are smaller than viscous forces (small Reynold number).

- The "z" component of the velocity is smaller than its radial component.

- The dependencies on angular velocity are very small.

- The variation of v

_{r}with r is much smaller than its variation with z (dv_{r}/dr<<dv_{r}/dz).

In agreement with the previous suppositions, the Navier-Stokes equations in cylindrical coordinates take the form

(1)

(2)

where u: interface velocity, and t: radial component of the stress tensor.

The equations for the flow inside the drops are given by the equation of continuity for an incompressible fluid and the Navier-Stokes equation within the quasi-static approximation, so that

Ñ·v= 0 (3)

-ÑP_{d}+µ_{d}Ñ^{2}v= 0 (4)

where P_{d}: pressure at the dispersed phase, and μ_{d} its viscosity.

On the interfacial film the following conditions must be satisfied:

u= v_{r }(5)

and the sum of the shear stress of the dispersed phase and the interfacial film must be zero, so that:

t + t= 0 (6)

In order to obtain an adimensional system of equations, the variables of the system are scaled according to the following relationships:

(7)

where R_{eq}: drop radius, and V: approach velocity between the droplets.

Now, the equation system has only one parameter, which is the capillary number (Ca= μ_{d}V/s, where s: interfacial tension). This parameter can be eliminated carrying out the following scaling on the system variables:

(8)

The adimensional system of equations is:

(9)

(10)

(11)

(12)

_{ }(13)

_{ }(14)

(15)

where D*: adimensional Laplacian operator acting over h*(r*,t*) in cylindrical coordinates.

The central idea of this section is to obtain from the equation system (Eqs. 9-15) one equation for the surface h*(r*,t*) without the unknown variable u*. For this purpose, Eq. 10 is inserted in Eq. 14, so that

(16)

where

(17)

If, additionally, the fact that at the interface u*= v_{r}^{*} is taken into account, then

(18)

and the pressure gradient can be calculated from Eq. 15. The resulting equation is inserted in Eq. 18, thus giving

(19)

Now, manipulating Eq. 9,

(20)

On Eq. 20, the term in parenthesis is zero (continuity equation), so that

(21)

The derivative of Eq. 21 with respect to r* is

(22)

Next, u* from Eq. 21 is introduced in Eq. 22, to obtain

(23)

On the other hand, if the surfactant is absent there are no Marangoni stresses on the drops surfaces, so that Eq. 12 takes the form

(24)

but the radial component of the vector field in Eq. 24 is

(25)

Taking into account that at the interface u*= v^{*}_{r} and inserting Eqs. 19, 21 and 23 in Eq. 25, it is obtained

(26)

and in turn, Eq. 26 can be manipulated so as to lead to

(27)

The initial and boundary conditions of the problem are

(28)

where R: interfacial film radius, and m: a real number with value of 0.5.

Eq. 27 is the evolution equation for the interfacial film of thickness h*, which was solved numerically using the finite differences method.

Results

The resolution of the evolution equation allows to know the front-wave dynamics at the interfacial film. As initial condition, a Gaussian form disturbance was used in order to simulate the dimple dynamics.

The evolution of the thickness h* of the interfacial film with respect to position and time is shown in Figure 2. It can be seen that the initial disturbance in Gaussian form decreases with time, until it reaches a maximum amplitude and finally returns without the interfacial film broken (the rupture of the interfacial film is the last stage of the coalescence process). Also, it can be observed that the zones near the edges of the interfacial film remain approximately stable. This is a consequence of the approach employed in the development of the model, where an approximately laminar flow draining towards the outside of the film has been assumed.

In the drainage process there are changes in the pressure field on the interfacial film. In Figure 3 it can be observed that the pressure field remains constant for each t up to r* =5, and thereafter the pressure field decays down to zero. Also, it can be seen that the pressure distribution tends to grow with time, which means that the pressure grows when the liquid interface drains towards the outside.

At the center of the interfacial film pressure remains constant with a value of 2, which means that in this region the term

of D*h* is annulled with its second derivative .

The behavior of the radial component of the viscous stress tensor can be appreciated in Figure 4. It is observed that for values of r* <4.5 the stress is zero. For values of r* >4.5 stress grows up to a maximum value and finally it decreases until being annulled.

When the liquid interface drains towards the outside it is observed that stress decreases, until the dimple drains completely out of the drop surface.

Conclusions

A model based on Navier-Stokes equations that describes dimple dynamics was elaborated. To this end, the hydrodynamic equations were considered in the spirit of two fundamental approaches, the quasi-static and the lubrication approximations.

The resolution of the evolution equation allows to know how does the dimple dynamics develop.

By solving the evolution equation, the dependence of the interfacial film thickness on position and time was obtained. An initial disturbance in Gaussian form decreased with time until reaching maximum amplitude and no coalescence was observed.

The pressure field diminished with time, from the center of the interfacial film to the barrier ring. This behavior was inverted when the dimple changed its concavity.

When the liquid interface drained towards the outside, the stress on the interfacial film diminished from the barrier ring to the center, and this tendency was inverted when the dimple returned to the equilibrium position.

References

1. Bibette J (1992) Stability of thin films in concentrated emulsions. *Langmuir 8*: 3178-3182. [ Links ]

2. Bibette J, Morse DC, Witten TA, Weitz DA (1992) Stability criteria for emulsions. *Phys. Rev. Lett. 69*: 2439-2443. [ Links ]

3. Chesters AK, Bazhlekov IB (2000) Effect of insoluble surfactants on drainage and rupture of a film between drops interacting under a constant force. *J. Coll. Interf. Sci. 230*: 229-243. [ Links ]

4. Danov K, Petsev D, Denkov N (1993) Pair interaction energy between deformable drops and bubbles. *J. Chem. Phys. 99*: 9-15. [ Links ]

5. Denkov N, Kralchevsky P, Ivanov I, Vassilieff C (1991) Effect of droplet deformation on the interactions in microemulsions. *J. Coll. Interf. Sci. 143*: 157-173. [ Links ]

6. Exerowa D, Balinov B, Kashchiev D (1983) Nucleation mechanism of rupture of newtonian black films : II. Experiment. *J. Coll. Interf. Sci. 94*: 45-53. [ Links ]

7. Hartland S, Jeelani SAK (1994) Effect of interfacial tension gradients on emulsion stability. *Coll. Surf. A 88*: 289-302. [ Links ]

8. Ivanov IB (1988) *Thin Liquid Films. Fundamentals and Applications*. Dekker. New York, USA. 1160 pp. [ Links ]

9. Ivanov I, Kralchevsky PA (1997) Stability of emulsions under equilibrium and dynamic conditions. *Coll. Surf. A 128*: 155-175. [ Links ]

10. Ivanov IB, Radoev B, Manev E, Scheludko A (1969) Theory of the critical thickness of rupture of thin liquid films. *Disc. Farad. Soc. 1262*: 13-21. [ Links ]

11. Jaeger PhT, Janssen JJM, Groeneweg F, Agterof WGM (1994) Coalescence in emulsions containing inviscid drops with high interfacial mobility. *Coll. Surf. A 85*: 255-264. [ Links ]

12. Kabalnov A, Wennerström H (1996) Macroemulsion stability: The oriented wedge theory revisited. *Langmuir 12*: 276-292. [ Links ]

13. Kashchiev D, Exerowa D (1980) Nucleation mechanism of rupture of newtonian black films. I. Theory. *J. Coll. Interf. Sci. 77*: 501-511. [ Links ]

14. Kralchevsky P, Ivanov I, Dimitrov S (1991) Interfacial fluctuation-dissipation processes and contact angles of thin liquid films. *Chem. Phys. Lett. 187*: 129-136. [ Links ]

15. Sharma A, Ruckenstein E (1987) Stability, critical thickness, and the time of rupture of thinning foam and emulsion films. *Langmuir 3*: 760-768. [ Links ]

16. Sonin AA, Bonfillon A, Langevin D (1994) Thinning of soap films: The role of surface viscoelasticity. *J. Coll. Interf. Sci. 162*: 323-330. [ Links ]

17. Tsekov R, Radoev B (1992) Rupture of thinning liquid films. *J. Chem. Soc. Farad. Trans. 88*: 251-253. [ Links ]

18. Velev OD, Gurkov TD, Borwankar RP (1993) Spontaneous cyclic dimpling in emulsion films due to surfactant mass transfer between the phases. *J. Coll. Interf. Sci. 159*: 497-501. [ Links ]

19. Vrij A (1966) Possible mechanism for the spontaneous rupture of thin, free liquid films. *Discuss. Farad. Soc. 42*: 23-33. [ Links ]

20. Vrij A, Overbeek JThG (1968) Rupture of thin liquid films due to spontaneous fluctuations in thickness. *J. Am. Chem. Soc. 90*: 3074-3078. [ Links ]

21. Yeo LY, Matar OK, Pérez de Ortiz ES, Hewitt GF (2001) The dynamics of Marangoni-driven local film drainage between two drops. *J. Coll. Interf Sci. 241*: 233-247. [ Links ]

22. Yeo LY, Matar OK, Pérez de Ortiz ES, Hewitt GF (2003) Film drainage between two surfactant-coated drops colliding at constant approach velocity. *J. Coll. Interf. Sci. 257*: 93-107. [ Links ]