 Solving the nonlinear and nonstationary Richards equation with two-level adaptive domain decomposition (dd-adaptivity)Michal Kuraz, Petr Mayer and Pavel Pech Applied Mathematics and Computation, 2015, vol. 267, issue C, 207-222 Abstract: Modeling the transport processes in the vadose zone, e.g. modeling contaminant transport, the effect of the soil water regime on changes in soil structure and composition, plays an important role in predicting the reactions of soil biotopes to anthropogenic activities. Water flow is governed by the quasilinear Richards equation, while the constitutive laws are typically supplied by the van Genuchten model, which can be understood as a pore size distribution function. Certain materials with dominantly uniform pore sizes (e.g. coarse-grained materials) can exhibit ranges of constitutive function values within several orders of magnitude. Numerical approximation of the Richards equation requires sequential solutions of systems of linear equations arising from discretization and linearization of the problem. Typically, in the case of two- and three-dimensional problems, it is necessary to solve huge systems of linear equations to obtain only a few local updates of the solution. Since the Richards equation is nonstationary, the local updates typically represent local disturbances originating e.g. from a moving wetting front. An efficient method for splitting the solutions of nonlinear systems should therefore be time adaptive. Aspects of this problem have already been discussed in Kuraz et al. (2013, 2014), where the subdomain split resulted in a high number of subdomains. This paper extends our work in the following directions. The global data (matrix, b- and x-vectors) are now handled as a scattered adaptively changing progression of local data, so we can effectively deactivate parts of the computational domain. This can be understood as a sort of multi-time-step approach. Secondly, we can now handle overlapping domain decomposition together with multi-level treatment of our time-adaptive Schwarz–Picard iterative method. Keywords: Computational hydropedology (hydrology); Spatial adaptivity; Nonlinear operator treatment; Schwarz domain decomposition with coarse level; dd-adaptivity; Multi-time-step methods (search for similar items in EconPapers) Downloads: (external link) Related works: Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text Persistent link: https://EconPapers.repec.org/RePEc:eee:apmaco:v:267:y:2015:i:c:p:207-222 DOI: 10.1016/j.amc.2015.03.130 Access Statistics for this article Applied Mathematics and Computation is currently edited by Theodore Simos More articles in Applied Mathematics and Computation from Elsevier |
Is your work missing from RePEc? Here is how to contribute.
Questions or problems? Check the EconPapers FAQ or send mail to .
EconPapers is hosted by the
School of Business at Örebro University.