A conservative SPH method for surfactant dynamics

In this paper, a Lagrangian particle method is proposed for the simulation of multiphase flows with surfactant. The model is based on the multiphase smoothed particle hydrodynamics (SPH) framework of Hu and Adams (2006) [1]. Surface-active agents (surfactants) are incorporated into our method by a scalar quantity describing the local concentration of molecules in the bulk phase and on the interface. The surfactant dynamics are written in conservative form, thus global mass of surfactant is conserved exactly. The transport model of the surfactant accounts for advection and diffusion. Within our method, we can simulate insoluble surfactant on an arbitrary interface geometry as well as interfacial transport such as adsorption or desorption. The flow-field dynamics and the surfactant dynamics are coupled through a constitutive equation, which relates the local surfactant concentration to the local surface-tension coefficient. Hence, the surface-tension model includes capillary and Marangoni-forces. The present numerical method is validated by comparison with analytic solutions for diffusion and for surfactant dynamics. More complex simulations of an oscillating bubble, the bubble deformation in a shear flow, and of a Marangoni-force driven bubble show the capabilities of our method to simulate interfacial flows with surfactants.

[1]  Patrick J. Fox,et al.  Smoothed Particle Hydrodynamics Model for Diffusion through Porous Media , 2001 .

[2]  O. Basaran,et al.  Repeated formation of fluid threads in breakup of a surfactant-covered jet. , 2005, Physical review letters.

[3]  K. Stebe,et al.  Marangoni effects on drop deformation in an extensional flow: The role of surfactant physical chemistry. I. Insoluble surfactants , 1996 .

[4]  C. Pozrikidis,et al.  Effect of surfactants on the deformation of drops and bubbles in Navier–Stokes flow , 2006 .

[5]  Hongkai Zhao,et al.  An Eulerian Formulation for Solving Partial Differential Equations Along a Moving Interface , 2003, J. Sci. Comput..

[6]  D. An,et al.  The effects of surfactants on drop deformation and breakup By , 2005 .

[7]  Ying-Chih Liao,et al.  Effects of dynamic surface tension and fluid flow on the oscillations of a supported bubble , 2006 .

[8]  J A CLEMENTS,et al.  Surface Tension of Lung Extracts , 1957, Proceedings of the Society for Experimental Biology and Medicine. Society for Experimental Biology and Medicine.

[9]  J. Lowengrub,et al.  A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant , 2004 .

[10]  J. Li,et al.  Numerical simulation of breakup of a viscous drop in simple shear flow through a volume-of-fluid method , 2000 .

[11]  Nikolaus A. Adams,et al.  A multi-phase SPH method for macroscopic and mesoscopic flows , 2006, J. Comput. Phys..

[12]  R D Kamm,et al.  Dynamic behavior of lung surfactant. , 2001, Journal of biomechanical engineering.

[13]  Darsh T. Wasan,et al.  The kinetics of adsorption of surface active agents at gas-liquid surfaces , 1983 .

[14]  R D Kamm,et al.  Dynamic surface tension of surfactant TA: experiments and theory. , 1994, Journal of applied physiology.

[15]  J. K. Chen,et al.  A corrective smoothed particle method for boundary value problems in heat conduction , 1999 .

[16]  Ming-Chih Lai,et al.  An immersed boundary method for interfacial flows with insoluble surfactant , 2008, J. Comput. Phys..

[17]  S. Schürch,et al.  Surface activity of lipid extract surfactant in relation to film area compression and collapse. , 1994, Journal of applied physiology.

[18]  L. Scriven,et al.  The Marangoni Effects , 1960, Nature.

[19]  Metin Muradoglu,et al.  A front-tracking method for computation of interfacial flows with soluble surfactants , 2008, J. Comput. Phys..

[20]  Zhilin Li,et al.  A level-set method for interfacial flows with surfactant , 2006, J. Comput. Phys..

[21]  Jie Zhang,et al.  A front tracking method for a deformable intravascular bubble in a tube with soluble surfactant transport , 2006, J. Comput. Phys..

[22]  Geoffrey Ingram Taylor,et al.  The formation of emulsions in definable fields of flow , 1934 .

[23]  M. Timmermans,et al.  The effect of surfactant on the stability of a liquid thread , 2002, Journal of Fluid Mechanics.

[24]  K. v. Neergaard,et al.  Neue Auffassungen über einen Grundbegriff der Atemmechanik , 1929 .

[25]  Petros Koumoutsakos,et al.  Simulations of (an)isotropic diffusion on curved biological surfaces. , 2006, Biophysical journal.

[26]  J. Morris,et al.  Modeling Low Reynolds Number Incompressible Flows Using SPH , 1997 .

[27]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[28]  R. E. Pattle,et al.  Properties, Function and Origin of the Alveolar Lining Layer , 1955, Nature.

[29]  Andrew W. Fitzgibbon,et al.  Direct Least Square Fitting of Ellipses , 1999, IEEE Trans. Pattern Anal. Mach. Intell..

[30]  Elias I. Franses,et al.  Adsorption dynamics of surfactants at the air/water interface: a critical review of mathematical models, data, and mechanisms , 1995 .

[31]  C. Pozrikidis,et al.  A Finite-volume/Boundary-element Method for Flow Past Interfaces in the Presence of Surfactants, with Application to Shear Flow Past a Viscous Drop , 1998 .

[32]  W. Bartsch,et al.  The surfactant system of the adult lung: physiology and clinical perspectives , 1992, The clinical investigator.

[33]  O. Basaran,et al.  Effects of insoluble surfactants on the nonlinear deformation and breakup of stretching liquid bridges , 1999 .

[34]  Guillermo Sapiro,et al.  Variational Problems and Partial Differential Equations on Implicit Surfaces: Bye Bye Triangulated Surfaces? , 2003 .

[35]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[36]  L. Gary Leal,et al.  The Influence of Surfactant on the Deformation and Breakup of a Viscous Drop: The Effect of Surfactant Solubility , 1994 .