Abstract:With the aim of efficiently simulating three-dimensional multiphase turbulent flows with a phase-field method, we propose a new discretization scheme for the biharmonic term (the 4th-order derivative term) of the Cahn-Hilliard equation. This novel scheme can significantly reduce the computational cost while retaining the same accuracy as the original procedure. Our phase-field method is built on top of a direct numerical simulation solver, named AFiD (www.afid.eu) and open-sourced by our research group. It relies on a pencil distributed parallel strategy and a FFT-based Poisson solver. To deal with large density ratios between the two phases, a pressure split method [1] has been applied to the Poisson solver. To further reduce computational costs, we implement a multiple-resolution algorithm which decouples the discretizations for the Navier-Stokes equations and the scalar equation: while a stretched wall-resolving grid is used for the Navier-Stokes equations, for the Cahn-Hilliard equation we use a fine uniform mesh. The present method shows excellent computational performance for large-scale computation: on meshes up to 8 billion nodes and 3072 CPU cores, a multiphase flow needs only slightly less than 1.5 times the CPU time of the single-phase flow solver on the same grid. The present method is validated by comparing the results to previous studies for the cases of drop deformation in shear flow, including the convergence test with mesh refinement, and breakup of a rising buoyant bubble with density ratio up to 1000. Finally, we simulate the breakup of a big drop and the coalescence of O(103) drops in turbulent Rayleigh-Bénard convection at a Rayleigh number of 108, observing good agreement with theoretical results.

An Efficient Dimension Splitting P-Adaptive Method for the Binary Fluid Surfactant Phase Field Model

Efficient numerical simulation of the conserved Allen-Cahn type flow-coupled binary fluid-surfactant model by a dimension splitting method

An Adaptive Time-Stepping Method for the Binary Fluid-Surfactant Phase Field Model on Evolving Surfaces

Numerical Approximation of a Binary Fluid-Surfactant Phase Field Model of Two-Phase Incompressible Flow

Decoupled and Energy Stable Time-Marching Scheme for the Interfacial Flow with Soluble Surfactants.

Numerical Approximation of a Phase-Field Surfactant Model with Fluid Flow

An efficiently linear and totally decoupled variant of SAV approach for the binary phase-field surfactant fluid model

Fast and stable dimension splitting simulations for the hydrodynamically coupled three-component conserved Allen-Cahn phase field model

Numerical simulation for the conserved Allen–Cahn phase field model of two-phase incompressible flows by an efficient dimension splitting method

Decoupled and energy stable schemes for phase-field surfactant model based on mobility operator splitting technique

SPH simulations of three-dimensional non-Newtonian free surface flows

An adaptive phase field method for the mixture of two incompressible fluids

Convergence and Stability Analysis of Energy Stable and Bound‐preserving Numerical Schemes for Binary Fluid‐surfactant Phase‐field Equations

Modeling and numerical simulation of surfactant systems with incompressible fluid flows on surfaces

Numerical Approximations for the Cahn–Hilliard Phase Field Model of the Binary Fluid-Surfactant System

Interfacial dynamics with soluble surfactants: A phase-field two-phase flow model with variable densities

An efficient phase-field method for turbulent multiphase flows

Efficient Numerical Simulation of Cahn-Hilliard Type Models by a Dimension Splitting Method

Efficiently linear and unconditionally energy-stable time-marching schemes with energy relaxation for the phase-field surfactant model

Decoupled, Energy Stable Schemes for a Phase-Field Surfactant Model.

Long Time Numerical Simulations for Phase-Field Problems Using P-Adaptive Spectral Deferred Correction Methods.