Hostname: page-component-78c5997874-t5tsf Total loading time: 0 Render date: 2024-11-10T14:12:04.565Z Has data issue: false hasContentIssue false

A Coupled Immersed Interface and Level Set Method for Three-Dimensional Interfacial Flows with Insoluble Surfactant

Published online by Cambridge University Press:  03 June 2015

Jian-Jun Xu*
Affiliation:
School of Mathematical and Computational Sciences, Xiangtan University, Hunan 410005, China Hunan Key Lab for Computation and Simulation in Science and Engineering, Hunan 410005, China
Yunqing Huang*
Affiliation:
School of Mathematical and Computational Sciences, Xiangtan University, Hunan 410005, China Hunan Key Lab for Computation and Simulation in Science and Engineering, Hunan 410005, China
Ming-Chih Lai*
Affiliation:
Department of Applied Mathematics, National Chiao Tung University, Hsinchu 30050, Taiwan
Zhilin Li*
Affiliation:
Department of Mathematics, North Carolina State University, Raleigh, NC, 27695, USA School of Mathematical Sciences, Nanjing Normal University, Nanjing, China
Get access

Abstract

In this paper, a numerical method is presented for simulating the 3D interfacial flows with insoluble surfactant. The numerical scheme consists of a 3D immersed interface method (IIM) for solving Stokes equations with jumps across the interface and a 3D level-set method for solving the surfactant convection-diffusion equation along a moving and deforming interface. The 3D IIM Poisson solver modifies the one in the literature by assuming that the jump conditions of the solution and the flux are implicitly given at the grid points in a small neighborhood of the interface. This assumption is convenient in conjunction with the level-set techniques. It allows standard Lagrangian interpolation for quantities at the projection points on the interface. The interface jump relations are re-derived accordingly. A novel rotational procedure is given to generate smooth local coordinate systems and make effective interpolation. Numerical examples demonstrate that the IIM Poisson solver and the Stokes solver achieve second-order accuracy. A 3D drop with insoluble surfactant under shear flow is investigated numerically by studying the influences of different physical parameters on the drop deformation.

Type
Research Article
Copyright
Copyright © Global Science Press Limited 2014

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Adams, J., Swarztrauber, P., and Sweet, R., Fishpack: Efficient Fortran subprograms for the solution of separable elliptic partial differential equations, http://www.netlib.org/fishpack/.Google Scholar
[2]Beale, J.T., and Layton, A.T., On the accuracy of finite differencec methods for elliptic problems with interfaces, Commun. Appl. Math. Comput. Sci., 1(2006), 91119.CrossRefGoogle Scholar
[3]Chen, K.-Y., Feng, K.-A., Kim, Y. and Lai, M.-C., A note on pressure accuracy in immersed boundary method for Stokes flow, J. Comput. Phys., 230(2011), 43774383.Google Scholar
[4]Deng, S., Ito, K., and Li, Z., Three dimensional elliptic solvers for interface problems and applications, J. Comput. Phys., 184(2003), 215243.Google Scholar
[5]Ganesan, S. and Tobiska, L., A coupled arbitrary Lagrangian-Eulerian and Lagrangian method for computation of free surface flows with insoluble surfactants, J. Comput. Phys., 228(2009), 28592873.CrossRefGoogle Scholar
[6]Huang, H., and Li, Z., Convergence analysis of the immersed interface method, IMA J. Numer. Anal., 19(1999), 583608.Google Scholar
[7]Ito, K., Li, Z., and Wan, X., Pressure jump conditions for Stokes equations with discontinuous viscosity in 2D and 3D, Methods and Appl. Anal., 19(2006), 229234.Google Scholar
[8]James, A.J. and Lowengrub, J., A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant, J. Comput. Phys., 201(2004), 685722.Google Scholar
[9]Jiang, G.-S. and Peng, D., Weighted ENO schemes for Hamilton-jacobi equations, SIAM J. Sci. Comput., 21(2000), 2126.CrossRefGoogle Scholar
[10]Khatri, S.M. and Tornberg, A.K., A numerical method for two-phase flows with insoluble surfactants, Comput. Fluids, 49(2011), 150165.Google Scholar
[11]Lai, M.-C. and Li, Z., A remark on jump conditions for the three-dimensional Navier-Stokes equations involving an immersed moving membrane, Applied Math. Lett., 14(2001), 149154.Google Scholar
[12]Lai, M.-C., Tseng, Y.-H. and Huang, H., An immersed boundary method for interfacial flows with insoluble surfactant, J. Comput. Phys., 227(2008), 7279.Google Scholar
[13]Lai, M.-C., Tseng, Y.-H. and Huang, H., Numerical simulation of moving contact lines with surfactant by immersed boundary method, Communications in Computational Physics, 8(2010), 735757.CrossRefGoogle Scholar
[14]Lai, M.-C., Huang, C.-Y. and Huang, Y.-M., Simulating the axisymmetric interfacial flows with insoluble surfactant by immersed boundary method, International Journal of Numerical Analysis and Modeling, 8(2011), 105117.Google Scholar
[15]R. J., LeVeque, and Z., LiThe immersed interface method for elliptic equation with discontinuous coefficients and singular sources, SIMA J. Numer. Anal., 31(1994), 10191044Google Scholar
[16]Li, X. and Pozrikidis, C., Effect of surfactants on drop deformation and on the rheology of dilute emulsion in Stokes flow, J. Fluid Mech., 385(1999), 7999.Google Scholar
[17]Li, Z., and Ito, K., The immersed interface method: Numerical for PDEs involving interfaces and irregular domains, SIAM frontiers in applied mathematics, SIAM, Philadelphia, USA, 2006.Google Scholar
[18]Lowengrub, J., Xu, J., and Voigt, A., Surface phase separation and flow in a simple model of multicomponent drops and vesicles, Fluid Dyn. Mater. Proc., 3(2007), 119.Google Scholar
[19]Muradoglu, M. and Tryggvason, G., A front-tracking method for the computation of interfacial flows with soluble surfactants, J. Comput. Phys., 227(2008), 2238.CrossRefGoogle Scholar
[20]Murry, G., Rotation about an arbitrary axis in 3 dimensions, http://inside.mines.edu/∼gmurray/ArbitraryAxisRotation/Google Scholar
[21]Osher, S. and Fedkiw, R.P., Level set methods: An overview and some recent results, J. Comput. Phys., 169(2001), 463.Google Scholar
[22]Peskin, C., Numerical analysis of blood flow in the heart, em J. Comput. Phys., 25(1977), 220252.Google Scholar
[23]Peskin, C., The immersed boundary methods, Acta Numerica,11(2002), 479517.Google Scholar
[24]Sethian, J.A. and Smereka, P., Level set methods for fluid interfaces, Ann. Rev. Fluid Mech., 35(2003), 341.Google Scholar
[25]Shu, C.-W., Total-variation-diminishing time discretization, SIAM J. Sci. Stat. Comput., 9(1988), 1073.Google Scholar
[26]Stone, H.A. and Leal, L.G., The effects of surfactants on drop deformation and breakup, J. Fluid Mech., 220(1990), 161186.Google Scholar
[27]Teigen, K.E., Song, P., Lowengrub, J. and Voigt, A., A diffusive-interface method for two-phase flows with soluble surfactants, J. Comput. Phys., 230(2011), 375393.Google Scholar
[28]Xu, J., Li, Z., Lowengrub, J., and Zhao, H., A level set method for solving interfactial flows with surfactant, J. Comput. Phys., 212(2006), 590616.CrossRefGoogle Scholar
[29]Xu, J., Li, Z., Lowengrub, J., Zhao, H., Numerical study of surfactant-laden drop-drop interactions, Commun. Comput. Phys., 10(2011), 453473.Google Scholar
[30]Xu, J., Yuan, H.-Z. and Huang, Y.-Q., A three dimensional level-set method for solving convection-diffusion equations along moving interfaces (in Chinese). Sci Sin Math, 42(2012), 445454Google Scholar
[31]Xu, J., Yang, Y., and Lowengrub, J.A level-set continuum method for two-phase flows with insoluble surfactants, J. Comput. Phys., 231(2012), 58975909.CrossRefGoogle Scholar
[32]Yang, X. and James, A.J., An arbitrary Lagrangian-Eulerian (ALE) method for interfacial flows with insoluble surfactants, Fluid Dyn. Mater. Proc., 3(2007), 6596.Google Scholar