  A discretization is proposed for models coupling free flow with anisotropic
porous medium flow. Our approach employs a staggered grid finite volume method
for the Navier-Stokes equations in the free flow subdomain and a MPFA finite
volume method to solve Darcy flow in the porous medium. After appropriate
spatial refinement in the free flow domain, the degrees of freedom are
conveniently located to allow for a natural coupling of the two discretization
schemes. In turn, we automatically obtain a more accurate description of the
flow field surrounding the porous medium. Numerical experiments highlight the
stability and applicability of the scheme in the presence of anisotropy and
show good agreement with existing methods, verifying our approach.
