“…This method, often called penalization approach (see, e.g., [9,34,35,39]), consists in considering in the whole computational domain a modified formulation of the Navier-Stokes equations which reduce to their classical form in the fluid region while they include additional resistance terms in the porous region. This approach is similar to the so-called fictitious domain method [1,39].…”