A projection-based, semi-implicit time-stepping approach for the Cahn-Hilliard Navier-Stokes equations on adaptive octree meshes
arXiv:2107.05123 · doi:10.1016/j.jcp.2022.111874
Abstract
The Cahn-Hilliard Navier-Stokes (CHNS) system provides a computationally tractable model that can be used to effectively capture interfacial dynamics in two-phase fluid flows. In this work, we present a semi-implicit, projection-based finite element framework for solving the CHNS system. We use a projection-based semi-implicit time discretization for the Navier-Stokes equation and a fully-implicit time discretization for the Cahn-Hilliard equation. We use a conforming continuous Galerkin (cG) finite element method in space equipped with a residual-based variational multiscale (RBVMS) formulation. Pressure is decoupled using a projection step, which results in two linear positive semi-definite systems for velocity and pressure, instead of the saddle point system of a pressure-stabilized method. All the linear systems are solved using an efficient and scalable algebraic multigrid (AMG) method. We deploy this approach on a massively parallel numerical implementation using parallel octree-based adaptive meshes. The overall approach allows the use of relatively large time steps with much faster time-to-solve than similar fully-implicit methods. We present comprehensive numerical experiments showing detailed comparisons with results from the literature for canonical cases, including the single bubble rise and Rayleigh-Taylor instability.
arXiv comment51 pages, 18 Figures, version accepted for publication Journal of Computational Physics. arXiv admin note: text overlap with arXiv:2009.06628
References in corpus (13)
- Machine Learning for Fluid Mechanics
- Perspectives on Machine Learning-augmented Reynolds-averaged and Large Eddy Simulation Models of Turbulence
- A second order in time, uniquely solvable, unconditionally stable numerical scheme for Cahn-Hilliard-Navier-Stokes equation
- A thermodynamically consistent phase-field model for two-phase flows with thermocapillary effects
- Mass Conservative and Energy Stable Finite Difference Methods for the Quasi-incompressible Navier-Stokes-Cahn-Hilliard system: Primitive Variable and Projection-Type Schemes
- A conservative diffuse-interface method for compressible two-phase flows
- Wall-bounded multiphase flows of N immiscible incompressible fluids: consistency and contact-angle boundary condition
- Massively Parallel Simulations of Binary Black Hole Intermediate-Mass-Ratio Inspirals
- Industrial scale large eddy simulations (LES) with adaptive octree meshes using immersogeometric analysis
- Simulating two-phase flows with thermodynamically consistent energy stable Cahn-Hilliard Navier-Stokes equations on parallel adaptive octree based meshes
- A fully-coupled framework for solving Cahn-Hilliard Navier-Stokes equations: Second-order, energy-stable numerical methods on adaptive octree based meshes
- Scalable adaptive PDE solvers in arbitrary domains
- Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows
Cited by in corpus (3)
- Efficient finite element schemes for a phase field model of two-phase incompressible flows with different densities
- On the Energy Analysis of Two-phase Flows Simulated with the Diffuse Interface Method
- A simple, fully-discrete, unconditionally energy-stable method for the two-phase Navier-Stokes Cahn-Hilliard model with arbitrary density ratios