Arbitrary Lagrangian-Eulerian finite element method for curved and deforming surfaces I. General theory and application to fluid interfaces

JOURNAL OF COMPUTATIONAL PHYSICS

卷:407

Article

Sahu, Amaresh¹  Omar, Yannick A. D.¹  Sauer, Roger A.²  Mandadapu, Kranthi K.^1,3 

[1] Univ Calif Berkeley, Dept Chem & Biomol Engn, Berkeley, CA 94720 USA

[2] Rhein Westfal TH Aachen, Aachen Inst Adv Study Computat Engn Sci AICES, Templergraben 55, D-52056 Aachen, Germany

[3] Lawrence Berkeley Natl Lab, Chem Sci Div, Berkeley, CA 94720 USA

关键词: Interfacial flows;    Fluid film;    Arbitrary Lagrangian-Eulerian;    Finite element method;    Differential geometry;    Lipid membrane;   

DOI  :  10.1016/j.jcp.2020.109253

来源: Elsevier

【 摘 要 】

An arbitrary Lagrangian-Eulerian (ALE) finite element method for arbitrarily curved and deforming two-dimensional materials and interfaces is presented here. An ALE theory is developed by endowing the surface with a mesh whose in-plane velocity need not depend on the in-plane material velocity, and can be specified arbitrarily. A finite element implementation of the theory is formulated and applied to curved and deforming surfaces with in-plane incompressible flows. Numerical inf-sup instabilities associated with inplane incompressibility are removed by locally projecting the surface tension onto a discontinuous space of piecewise linear functions. The general isoparametric finite element method, based on an arbitrary surface parametrization with curvilinear coordinates, is tested and validated against several numerical benchmarks. A new physical insight is obtained by applying the ALE developments to cylindrical fluid films, which are computationally and analytically found to be stable to non-axisymmetric perturbations, and unstable with respect to long-wavelength axisymmetric perturbations when their length exceeds their circumference. A Lagrangian scheme is attained as a special case of the ALE formulation. Though unable to model fluid films with sustained shear flows, the Lagrangian scheme is validated by reproducing the cylindrical instability. However, relative to the ALE results, the Lagrangian simulations are found to have spatially unresolved regions with few nodes, and thus larger errors. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free