A two-scale finite element formulation of Stokes flow in porous media

Seepage through saturated porous material with an open pore system is modeled as a non-linear Stokes flow through a rigid matrix. Based on variationally consistent homogenization, the resulting macroscale problem becomes a Darcy-type flow. The prolongation of the Darcy flow fulfills a macrohomogeneity condition, which in a Galerkin context implies a symmetric macroscale problem. The homogenization is of 1st order and periodic boundary conditions are adopted on a Representative Volume Element. A nonlinear nested multiscale technique, in which the subscale problem is used as a constitutive model, is devised. In the presented numerical investigation, the effects of varying physical parameters as well as of the discretization are considered. In particular, it is shown that the two-scale results agree well with those of the fully resolved fine-scale problem.