A computational method for approximating a Darcy–Stokes system governing a vuggy porous medium Academic Article