Isogeometric analysis of the advective CahnHilliard equation: Spinodal decomposition under shear flow
We present a numerical study of the spinodal decomposition of a binary fluid undergoing shear flow using the advective CahnHilliard equation, a stiff, nonlinear, parabolic equation characterized by the presence of fourth-order spatial derivatives. Our numerical solution procedure is based on isogeometric analysis, an approximation technique for which basis functions of high-order continuity are employed. These basis functions allow us to directly discretize the advective CahnHilliard equation without resorting to a mixed formulation. We present steady state solutions for rectangular domains in two-dimensions and, for the first time, in three-dimensions. We also present steady state solutions for the two-dimensional TaylorCouette cell. To enforce periodic boundary conditions in this curved domain, we derive and utilize a new periodic Bzier extraction operator. We present an extensive numerical study showing the effects of shear rate, surface tension, and the geometry of the domain on the phase evolution of the binary fluid. Theoretical and experimental results are compared with our simulations.
isogeometric analysis cahnhilliard equation spinodal decomposition shear flow steady state bzier extraction