Abstract:In this paper, we propose a spectral element-based phase field method by solving the Navier–Stokes/Cahn–Hilliard equations for incompressible two-phase flows. The high-order differential nonlinear term of the Cahn–Hilliard equation poses a great challenge for obtaining numerical solutions, so the Newton–Raphson method is adopted to tackle this conundrum. Furthermore, we employ the time-stepping scheme to decouple the Navier–Stokes equations to favor the computations with large density and viscosity contrast, in terms of both stability and convergence efficiency. In addition, the continuum surface tension model is used to account for the tangential component of the interfacial force, and thus our numerical method has the ability to simulate thermocapillary flows. We present four examples to demonstrate the interface capture accuracy of the proposed method. The shape of the rotating Zalesak's disk is accurately preserved by the present method even for two periods, which implies less dissipation and higher accuracy at long time numerical simulation. It is also noted that the present method with fourth-order element can achieve similar accuracy with the consistent scheme by evaluating the convective fluxes with the fifth-order weighted essentially non-oscillatory scheme. Moreover, the proposed method appears to comply well with mass conservation. And the results of thermocapillary flow test show good agreement with theoretical prediction. Finally, the rising bubble and Rayleigh–Taylor instability are considered to verify the proposed method for complex changes in interfacial topology, as well as its performance under large density and viscosity contrasts and high Reynolds number conditions.

New high accuracy three-point finite spectral scheme

High order method for solving unsteady incompressible Navier-Stokes equations

Accuracy and Stability Enhancements in the Incompressible Finite-Volume-particle Method for Multiphase Flow Simulations.

IMPROVEMENT OF FINITE-ANALYTIC DISCRETIZATION OF INCOMPRESSIBLE THREE-DIMENSIONAL NAVIER-STOKES EQUATIONS

High Order Compact Generalized Finite Difference Methods for Solving Inviscid Compressible Flows.

High-order Compact Finite Volume Methods on Unstructured Grids with Adaptive Mesh Refinement for Solving Inviscid and Viscous Flows

A new family of high-order compact upwind difference schemes with good spectral resolution

Convergence analysis of three semi-discrete numerical schemes for nonlocal geometric flows including perimeter terms

An Efficient Fourth-Order Three-Point Scheme for Solving Some Nonlinear Dispersive Wave Equations

Finite Spectral Essentially Non-Oscillatory Scheme

3D Implicit-Explicit Hybrid Finite Difference/Spectral Element/Finite Element Time Domain Method Without A Buffer Zone

A New Higher-Order Super Compact Finite Difference Scheme to Study Three-Dimensional Non-Newtonian Flows

A family of efficient high-order hybrid finite difference schemes based on WENO schemes

L2 APPROXIMATION HIGH-ACCURATE SCHEME: BASIC PRINCIPLE AND APPLICATION

A spectral element-based phase field method for incompressible two-phase flows

High-accuracy three-dimensional surface detection in smoothed particle hydrodynamics for free-surface flows

Improved SIP scheme for numerical solutions of 3-D Euler and N-S equations and its application

High-order finite volume methods for solving compressible multicomponent flows

A novel central compact finite-difference scheme for third derivatives with high spectral resolution

Finite Spectral Semi-Lagrangian Method for Incompressible Flows

Efficient energy-stable parametric finite element methods for surface diffusion flow and applications in solid-state dewetting