Algorithms for Solving Darcian Flow in Structured Porous Media
Type of document
articlePeer-reviewed
publishedVersion
Author
Kuráž , Michal
Mayer , Petr
Rights
Creative Commons Attribution 4.0 International Licensehttp://creativecommons.org/licenses/by/4.0/
openAccess
Metadata
Show full item recordAbstract
This paper presents several algorithms that were implemented in DRUtES [1], a new open source project. DRUtES is a finite element solver for coupled nonlinear parabolic problems, namely the Richards equation with the dual porosity approach (modeling the flow of liquids in a porous medium). Mass balance consistency is crucial in any hydrological balance and contaminant transportation evaluations. An incorrect approximation of the mass term greatly depreciates the results that are obtained. An algorithm for automatic time step selection is presented, as the proper time step length is crucial for achieving accuracy of the Euler time integration method. Various problems arise with poor conditioning of the Richards equation: the computational domain is clustered into subregions separated by a wetting front, and the nonlinear constitutive functions cover a high range of values, while a very simple diagonal preconditioning method greatly improves the matrix properties. The results are presented here, together with an analysis.
Collections
The following license files are associated with this item:
Except where otherwise noted, this item's license is described as Creative Commons Attribution 4.0 International License