“…where Ω ⊂ R 2 is a bounded and convex domain with the boundary Ω, u = (u 1 , u 2 ) and p are the fluid velocity and pressure, respectively, > 0 is the viscosity coefficient, 2 ≤ r < ∞ and > 0 are two damping parameters, f is a given external force, and |u| ∶= placement equations, 34 the Cahn-Hilliard equations, 33 and the nonlinear parabolic equation. 35 In Huang et al, 36 a novel two-step method was proposed for the Navier-Stokes equations, which uses a lower-order element to solve one nonlinear system as an initial approximation and then uses a higher-order element to solve one linear system. In fact, the two-step method can be cast in the framework of the two-level method, but compared with the two-level method, the two-step method only needs one mesh size, which can avoid the discussion on the relation of the coarse and fine meshes.…”