Abstract
In this work, we present a novel isogeometric analysis discretization for the Navier-Stokes-Cahn-Hilliard equation, which uses divergence-conforming spaces. Basis functions generated with this method can have higher-order continuity, and allow to directly discretize the higher-order operators present in the equation. The discretization is implemented in PetIGA-MF, a high-performance framework for discrete differential forms. We present solutions in a two-dimensional annulus, and model spinodal decomposition under shear flow.