“…The model that belongs to the class of diffuse interface or phase-field methods, has been used in physics, chemistry, biology, and engineering fields. In recent years, driven by the major developments of numerical algorithms and by increased availability of computational resources, direct numerical simulation of Cahn-Hilliard-Navier-Stokes equations has become increasingly popular [1,14,16,25,27,30].…”