The subdivision-based IGA-EIEQ numerical scheme for the Cahn-Hilliard-Darcy system of two-phase Hele-Shaw flow on complex curved surfaces

This paper investigates the numerical approximation of a two-phase Hele-Shaw fluid flow system, which consists of a phase-field Cahn-Hilliard model coupled with the Darcy flow, confined on a complex curved surface. To avoid geometric errors caused by the surface approximation and additional approximation errors introduced by numerical methods, for spatial discretization, we consider the recently developed strategy of subdivision-based isogeometric analysis (IGA), where we can accurately represent complex surfaces of arbitrary topology due to its high smoothness and hierarchical refinement properties. By incorporating multiple methods of temporal discretizations, including the Explicit-Invariant-Energy-Quadratization (EIEQ) approach for linearizing nonlinear potentials, the Zero-Energy-Contribution (ZEC) method for decoupling of nonlinear couplings, and the projection method for decoupling of linear couplings between the fluid velocity and pressure, and the finite-element method based on IGA for spatial discretizations, we arrive at a fully discrete scheme that possesses desirable properties such as decoupling, linearity, second-order accuracy in time, and unconditional energy stability. At each time step, only several elliptic equations with constant coefficients are needed to be solved. We also present the rigorous proof of unconditional energy stability, accompanied by numerous numerical examples to verify the stability and accuracy of the numerical scheme, such as the Saffman-Taylor fingering instability occurring on various complex curved surfaces induced by rotational forces.