In this paper, we propose a phase-field-based spectral element method by solving the Navier–Stokes/Cahn–Hilliard equations for incompressible two-phase flows. With use of Newton–Raphson Cahn–Hilliard equation and time-stepping scheme Navier–Stokes equation, construct three constant (time-independent) coefficient matrixes solutions velocity, pressure, phase variable. Moreover, invoke modified bu...