Numerical solution of porous media flow with a dual-permeability model

Abstract

The flow in porous media can be described by the Richards equation. However, porous media often exhibit a variety of heterogeneities, thus treat- ing a porous medium as homogeneous does not often fit the reality well. Therefore, we describe the flow in the porous medium using the Richards equation with the dual-permeability model, which assumes that the porous medium can be separated into two different media. This thesis is con- cerned with the numerical solution of the Richards equation with the dual- permeability model. We present the derivation of the dual-permeability model, and for the numerical solution, we use the space-time discontinu- ous Galerkin method. This produces a system of nonlinear algebraic equa- tions that need to be linearized. We perform a 1D experiment to verify the method and, finally, we present a 2D single-ring experiment to demonstrate the method. 1

Proudění v porézním prostředí lze popsat Richardsovou rovnicí. Porézní prostředí je často různorodé a homogenní popis nemusí odpovídat realitě. Proto budeme uvažovat model Richardsovy rovnice s duální permeabilitou. Tento model předpokládá, že můžeme rozdělit porézní prostředí na dvě různé na sobě ležící propojené části. V této práci se budeme zabývat numerickým řešením Richardsovy rovnice s duální permeabilitou. Popíšeme odvození to- hoto modelu a následně problém diskretizujeme nespojitou Galerkinovou metodou vzhledem k prostoru a času. Tato diskretizace vede na systém nelineárních algebraických rovnic. Na závěr provedeme dva numerické ex- perimenty. 1