“…Numerical methods of Navier-Stokes/Darcy have attracted a lot of attention. So far, a great deal of numerical methods are proposed to solve this model by virtue of different ways, such as finite element methods [12], discontinuous Galerkin finite element methods [5], two-grid methods [1,15,16], modified two-grid methods [6], partitioned time stepping method [7], characteristic stabilized finite element methods [8], mortar finite element methods [2], grad-div stabilized projection finite element method [14], modular grad-div method [11] and so on. The grad-div stabilized method is first introduced in [4], which can penalize mass conservation and improve the solution quality efficiently.…”