Simulating two-phase flows with thermodynamically consistent energy stable Cahn-Hilliard Navier-Stokes equations on parallel adaptive octree based meshes

JOURNAL OF COMPUTATIONAL PHYSICS

卷:419

Article

Khanwale, Makrand A.¹  Lofquist, Alec D.^1,4  Sundar, Hari³  Rossmanith, James A.²  Ganapathysubramanian, Baskar¹ 

[1] Iowa State Univ, Dept Mech Engn, Ames, IA 50011 USA

[2] Iowa State Univ, Dept Math, Ames, IA 50011 USA

[3] Univ Utah, Sch Comp, Salt Lake City, UT 84112 USA

[4] Affirm Inc, San Francisco, CA USA

关键词: Two-phase flows;    Energy stable;    Adaptive finite elements;    Octrees;    Scalable;   

DOI  :  10.1016/j.jcp.2020.109674

来源: Elsevier

【 摘 要 】

We report on simulations of two-phase flows with deforming interfaces at various density contrasts by solving thermodynamically consistent Cahn-Hilliard Navier-Stokes equations. An (essentially) unconditionally energy-stable Crank-Nicolson-type time integration scheme is used. Detailed proofs of energy stability of the semi-discrete scheme and for the existence of solutions of the advective-diffusive Cahn-Hilliard operator are provided. We discretize spatial terms with a conforming continuous Galerkin finite element method in conjunction with a residual-based variational multi-scale (VMS) approach in order to provide pressure stabilization. We deploy this approach on a massively parallel numerical implementation using fast octree-based adaptive meshes. A detailed scaling analysis of the solver is presented. Numerical experiments showing convergence and validation with experimental results from the literature are presented for a large range of density ratios. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free