Vesicle electrohydrodynamic simulations by coupling immersed boundary and immersed interface method

JOURNAL OF COMPUTATIONAL PHYSICS

卷:317

Article

Hu, Wei-Fan¹  Lai, Ming-Chih^2,3  Seol, Yunchang⁴  Young, Yuan-Nan⁵ 

[1] Natl Chung Hsing Univ, Dept Appl Math, 145 Xingda Rd, Taichung 402, Taiwan

[2] Natl Chiao Tung Univ, Ctr Math Modeling & Sci Comp, 1001 Ta Hsueh Rd, Hsinchu 300, Taiwan

[3] Natl Chiao Tung Univ, Dept Appl Math, 1001 Ta Hsueh Rd, Hsinchu 300, Taiwan

[4] Natl Taiwan Univ, Natl Ctr Theoret Sci, 1,Sec 4,Rd Roosevelt, Taipei 10617, Taiwan

[5] New Jersey Inst Technol, Dept Math Sci, Newark, NJ 07102 USA

关键词: Immersed boundary method;    Immersed interface method;    Vesicle electrohydrodynamics;    Leaky dielectric model;    Navier-Stokes equations;   

DOI  :  10.1016/j.jcp.2016.04.035

来源: Elsevier

【 摘 要 】

In this paper, we develop a coupled immersed boundary (IB) and immersed interface method (IIM) to simulate the electrodeformation and electrohydrodynamics of a vesicle in Navier-Stokes leaky dielectric fluids under a DC electric field. The vesicle membrane is modeled as an inextensible elastic interface with an electric capacitance and an electric conductance. Within the leaky dielectric framework and the piecewise constant electric properties in each fluid, the electric stress can be treated as an interfacial force so that both the membrane electric and mechanical forces can be formulated in a unified immersed boundary method. The electric potential and transmembrane potential are solved simultaneously via an efficient immersed interface method. The fluid variables in Navier-Stokes equations are solved using a projection method on a staggered MAC grid while the electric potential is solved at the cell center. A series of numerical tests have been carefully conducted to illustrate the accuracy and applicability of the present method to simulate vesicle electrohydrodynamics. In particular, we investigate the prolateoblate-prolate (POP) transition and the effect of electric field and shear flow on vesicle electrohydrodynamics. Our numerical results are in good agreement with those obtained in previous work using different numerical algorithms. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free