Numerical solution of porous media flow with a dual-permeability model
Numerická simulace proudění v porézním prostředí s duální permeabilitou
diploma thesis (DEFENDED)
View/ Open
Permanent link
http://hdl.handle.net/20.500.11956/185518Identifiers
Study Information System: 249823
Collections
- Kvalifikační práce [10932]
Author
Advisor
Referee
Congreve, Scott
Faculty / Institute
Faculty of Mathematics and Physics
Discipline
Computational Mathematics
Department
Department of Numerical Mathematics
Date of defense
14. 9. 2023
Publisher
Univerzita Karlova, Matematicko-fyzikální fakultaLanguage
English
Grade
Excellent
Keywords (Czech)
proudění v porézním prostředí|duální permebilita|numerická řešeníKeywords (English)
porous media flow|dual permeability|numerical solutionProudě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
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