A numerical study of the applicability of the Boussinesq approximation for a fluid‐saturated porous medium