On a fully implicit, monolithic finite element method-multigrid solution approach for dynamic porous media problems