Abstract
We develop an accurate and robust numerical scheme for solving the incompressible hydrodynamically coupled Cahn–Hilliard system of the two-phase fluid flow system on complex surfaces. Our algorithm leverages a number of efficient techniques, including the subdivision-based isogeometric analysis (IGA) method for spatial discretization, the explicit Invariant Energy Quadratization (EIEQ) method for linearizing nonlinear potentials, the Zero-Energy-Contribution (ZEC) method for decoupling, and the projection method for the Navier–Stokes equation to facilitate fully decoupled type implementations. The integration of these methodologies results in a fully discrete scheme with desired properties such as linearity, second-order temporal accuracy, full decoupling, and unconditional energy stability. The implementation of the scheme is straightforward, requiring the solution of a few elliptic equations with constant coefficients at each time step. The rigorous stability proof of unconditional energy stability and the implementation procedure are given in detail. Numerous numerical simulations on complex curved surfaces are carried out to verify the effectiveness of the proposed numerical scheme.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call