“…Over the last few decades, a great deal of effort has been devoted to developing an appropriate approximate solution to study the Stokes-Darcy model. Many different techniques and numerical methods were applied to investigate the Stokes-Darcy fluid flow model, such as coupled finite elements methods [15][16][17][18][19][20][21], two-grid/multi-grid methods [22][23][24][25][26][27][28], discontinuous Galerkin finite element methods [29][30][31][32], partitioned time stepping methods [33][34][35][36], least squares methods [37][38][39], domain decomposition methods [1,2,4,[40][41][42][43][44][45][46][47][48], local-parallel finite element methods [49], interface relaxation method [50,51], motar finite element methods [52,53], Lagrange multiplier methods [11,[54][55]…”