A second-order backward difference time-stepping scheme for penalized Navier-Stokes equations modeling filtration through porous media

We investigate the stability and convergence of a fully implicit, linearly extrapolated second‐order backward difference time‐stepping scheme for the penalized Navier–Stokes equations modeling filtration through porous media. In the penalization approach, an extended Navier–Stokes equation is used in the entire computational domain with suitable resistance terms to mimic the presence of porous medium. It is widely used as an alternative to the heterogeneous approach in which different types of partial differential equations (PDEs) are used in fluid and porous subregions along with suitable continuity conditions at the interface. However, the introduction of extra resistance terms makes the penalized Navier–Stokes equations more nonlinear. We prove that the linearly extrapolated scheme is unconditionally stable and derive optimal order error estimates without any stability condition. To show feasibility and applicability of the approach, it is used to numerically solve a passive control problem in which flow around a solid body is controlled by adding porous layers on the surface. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 681–705, 2016