Mathematical models and numerical methods for water management in soils
Duration: 02/2022 – 02/2024
Funding scheme: Funding for Mobility Cooperation with Germany
Funder: Academy of Finland
Budget: 18,000 €
My role: Finnish principal investigator
Sibling project: Discontinuous Galerkin methods and parameter estimation for microstructure models in porous media
German principal investigator: Nadja Ray
Public description
Nowadays, soils should meet more and more requirements. For example, they should be stable to build houses on them safely. In flood-prone areas, they should be able to absorb a lot of water, and in arid regions, they should hold water as well as possible. Soils also play an essential role in the greenhouse effect, as they (and not just the plants growing on their surface) can sequester enormous amounts of carbon and release it if not treated properly.
However, if we want to treat soils to meet all or at least some of these requirements, we need to understand the basic building blocks (i.e., microaggregates; see the picture) of soils and how their development affects soil properties. To this end, we develop two models and numerical methods to simulate them. First, we describe the formation of soil microaggregates, and second, we consider the flow of fluids through consolidated soil. Our goal is to apply the results to support water management.
The image is taken from Microaggregates in soils by Kai Uwe Totsche et al., licensed under CC BY 4.0.
Scientific abstract
This exchange program develops new mathematical tools for porous media applications. The main goal is to support water management by providing models for the consolidation of soils and the flow and transport through them as well as numerical methods to simulate them on the computer. We use these models and techniques to optimize soil parameters and ultimately serve the goal of treating the soil in such a way that, for example, plants find optimal growing conditions or icroorganisms in the ground can sequester as much CO2 as possible.
The first research focuses on parameter estimation for cellular automaton (CA) based models. Such models have been successfully developed and applied in soil science research by Andreas Rupp (Lappeenranta-Lahti University of Technology, abbreviated as LUT) and the group at the Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU). However, physically meaningful parameters that enter the CA model are not experimentally accessible. Consequently, the parameterization of the underlying models remains a challenge. Heikki Haario’s group (LUT) has published research articles on determining relevant parameters in partial differential equations. We use adapted versions of these methods to find appropriate parameters for the CA models.
In the second part of the exchange program, we address discontinuous Galerkin methods for models used to simulate soil applications with two porosities. Dual porosity models accurately describe flow through a porous medium whose solid matrix is composed of a porous medium. Here, we characterize the flow within the large pores by a Stokes equation while formulating the flow through the (micropores of the) porous matrix itself with Darcy’s law. Coupling strategies between the different flow models are constructed and analyzed using techniques developed in the context of petroleum applications or the interactions of lakes with their bottoms. Ultimately, our research may facilitate the design of more efficient batteries and accumulators. Thus, this research topic can be a starting point for collaborating with Prof. Bernardo Barbiellini (LUT) to optimize electric batteries further.
Keywords
cellular automaton model, coupled Stokes-Darcy flow, coupling condition, discontinuous Galerkin, dual-porosity flow, inverse problem, mathematical modeling, microaggregate, numerical methods, parameter estimation, soil science, water management
Related publications
Lu, Peipei; Wang, Wei; Kanschat, Guido; Rupp, Andreas
Homogeneous multigrid for HDG applied to the Stokes equation Journal Article Forthcoming
In: IMA Journal of Numerical Analysis, Forthcoming, ISBN: 0272-4979.
@article{LuWKR23,
title = {Homogeneous multigrid for HDG applied to the Stokes equation},
author = {Peipei Lu and Wei Wang and Guido Kanschat and Andreas Rupp},
url = {https://academic.oup.com/imajna/advance-article-pdf/doi/10.1093/imanum/drad079/52422702/drad079.pdf},
doi = {10.1093/imanum/drad079},
isbn = {0272-4979},
year = {2023},
date = {2023-10-23},
urldate = {2023-10-23},
journal = {IMA Journal of Numerical Analysis},
abstract = {We propose a multigrid method to solve the linear system of equations arising from a hybrid discontinuous Galerkin (in particular, a single face hybridizable, a hybrid Raviart–Thomas, or a hybrid Brezzi–Douglas–Marini) discretization of a Stokes problem. Our analysis is centered around the augmented Lagrangian approach and we prove uniform convergence in this setting. Beyond this, we establish relations, which resemble those in Cockburn and Gopalakrishnan (2008, Error analysis of variable degree mixed methods for elliptic problems via hybridization. Math. Comput., 74, 1653–1677) for elliptic problems, between the approximates that are obtained by the single-face hybridizable, hybrid Raviart–Thomas and hybrid Brezzi–Douglas–Marini methods. Numerical experiments underline our analytical findings.},
keywords = {},
pubstate = {forthcoming},
tppubtype = {article}
}
We propose a multigrid method to solve the linear system of equations arising from a hybrid discontinuous Galerkin (in particular, a single face hybridizable, a hybrid Raviart–Thomas, or a hybrid Brezzi–Douglas–Marini) discretization of a Stokes problem. Our analysis is centered around the augmented Lagrangian approach and we prove uniform convergence in this setting. Beyond this, we establish relations, which resemble those in Cockburn and Gopalakrishnan (2008, Error analysis of variable degree mixed methods for elliptic problems via hybridization. Math. Comput., 74, 1653–1677) for elliptic problems, between the approximates that are obtained by the single-face hybridizable, hybrid Raviart–Thomas and hybrid Brezzi–Douglas–Marini methods. Numerical experiments underline our analytical findings.
Lu, Peipei; Maier, Roland; Rupp, Andreas
A localized orthogonal decomposition strategy for hybrid discontinuous Galerkin methods Online Forthcoming
Forthcoming, visited: 28.07.2023.
@online{LuMR23,
title = {A localized orthogonal decomposition strategy for hybrid discontinuous Galerkin methods},
author = {Peipei Lu and Roland Maier and Andreas Rupp},
url = {https://arxiv.org/abs/2307.14961},
year = {2023},
date = {2023-07-28},
urldate = {2023-07-28},
abstract = {We formulate and analyze a multiscale method for an elliptic problem with an oscillatory coefficient based on a skeletal (hybrid) formulation. More precisely, we employ hybrid discontinuous Galerkin approaches and combine them with the localized orthogonal decomposition methodology to obtain a coarse-scale skeletal method that effectively includes fine-scale information. This work is a first step to reliably merge hybrid skeletal formulations and localized orthogonal decomposition and unite the advantages of both strategies. Numerical experiments are presented to illustrate the theoretical findings. },
keywords = {},
pubstate = {forthcoming},
tppubtype = {online}
}
We formulate and analyze a multiscale method for an elliptic problem with an oscillatory coefficient based on a skeletal (hybrid) formulation. More precisely, we employ hybrid discontinuous Galerkin approaches and combine them with the localized orthogonal decomposition methodology to obtain a coarse-scale skeletal method that effectively includes fine-scale information. This work is a first step to reliably merge hybrid skeletal formulations and localized orthogonal decomposition and unite the advantages of both strategies. Numerical experiments are presented to illustrate the theoretical findings.
Musch, Markus; Rupp, Andreas; Aizinger, Vadym; Knabner, Peter
Hybridizable discontinuous Galerkin method with mixed-order spaces for non-linear diffusion equations with internal jumps Journal Article
In: GEM - International Journal on Geomathematics, vol. 14, no. 18, pp. 25, 2023, ISSN: 1869-2680.
@article{MuschRAK23,
title = {Hybridizable discontinuous Galerkin method with mixed-order spaces for non-linear diffusion equations with internal jumps},
author = {Markus Musch and Andreas Rupp and Vadym Aizinger and Peter Knabner},
url = {https://link.springer.com/article/10.1007/s13137-023-00228-7},
doi = {10.1007/s13137-023-00228-7},
issn = {1869-2680},
year = {2023},
date = {2023-07-13},
urldate = {2023-07-13},
journal = {GEM - International Journal on Geomathematics},
volume = {14},
number = {18},
pages = {25},
abstract = {We formulate a hybridizable discontinuous Galerkin method for parabolic equations with non-linear tensor-valued coefficients and jump conditions (Henry’s law). The analysis of the proposed scheme indicates the optimal convergence order for mildly non-linear problems. The same order is also obtained in our numerical studies for simplified settings. A series of numerical experiments investigate the effect of choosing different order approximation spaces for various unknowns. },
keywords = {},
pubstate = {published},
tppubtype = {article}
}
We formulate a hybridizable discontinuous Galerkin method for parabolic equations with non-linear tensor-valued coefficients and jump conditions (Henry’s law). The analysis of the proposed scheme indicates the optimal convergence order for mildly non-linear problems. The same order is also obtained in our numerical studies for simplified settings. A series of numerical experiments investigate the effect of choosing different order approximation spaces for various unknowns.
Kazarnikov, Alexey; Ray, Nadja; Haario, Heikki; Lappalainen, Joona; Rupp, Andreas
Parameter estimation for cellular automata Online Forthcoming
Forthcoming, visited: 01.02.2023.
@online{KazarnikovRHLR23,
title = {Parameter estimation for cellular automata},
author = {Alexey Kazarnikov and Nadja Ray and Heikki Haario and Joona Lappalainen and Andreas Rupp},
url = {https://arxiv.org/abs/2301.13320},
year = {2023},
date = {2023-02-01},
urldate = {2023-02-01},
abstract = {Self organizing complex systems can be modeled using cellular automaton models. However, the parametrization of these models is crucial and significantly determines the resulting structural pattern. In this research, we introduce and successfully apply a sound statistical method to estimate these parameters. The method is based on constructing Gaussian likelihoods using characteristics of the structures such as the mean particle size. We show that our approach is robust with respect to the method parameters, domain size of patterns, or CA iterations. },
keywords = {},
pubstate = {forthcoming},
tppubtype = {online}
}
Self organizing complex systems can be modeled using cellular automaton models. However, the parametrization of these models is crucial and significantly determines the resulting structural pattern. In this research, we introduce and successfully apply a sound statistical method to estimate these parameters. The method is based on constructing Gaussian likelihoods using characteristics of the structures such as the mean particle size. We show that our approach is robust with respect to the method parameters, domain size of patterns, or CA iterations.
Kaarnioja, Vesa; Rupp, Andreas
Quasi-Monte Carlo and discontinuous Galerkin Online Forthcoming
Forthcoming, visited: 15.07.2022.
@online{KaarniojaR22,
title = {Quasi-Monte Carlo and discontinuous Galerkin},
author = {Vesa Kaarnioja and Andreas Rupp},
url = {https://arxiv.org/abs/2207.07698},
year = {2022},
date = {2022-07-15},
urldate = {2022-07-15},
abstract = {In this study, we design and develop Quasi-Monte Carlo (QMC) cubatures for non-conforming discontinuous Galerkin (DG) approximations of elliptic partial differential equations (PDEs) with random coefficient. That is, we consider the affine and uniform, and the lognormal models for the input random field, and investigate the QMC cubatures to compute the response statistics (expectation and variance) of the discretized PDE. In particular, we prove that the resulting QMC convergence rate for DG approximations behaves in the same way as if continuous finite elements were chosen. Numerical results underline our analytical findings. Moreover, we present a novel analysis for the parametric regularity in the lognormal setting.},
keywords = {},
pubstate = {forthcoming},
tppubtype = {online}
}
In this study, we design and develop Quasi-Monte Carlo (QMC) cubatures for non-conforming discontinuous Galerkin (DG) approximations of elliptic partial differential equations (PDEs) with random coefficient. That is, we consider the affine and uniform, and the lognormal models for the input random field, and investigate the QMC cubatures to compute the response statistics (expectation and variance) of the discretized PDE. In particular, we prove that the resulting QMC convergence rate for DG approximations behaves in the same way as if continuous finite elements were chosen. Numerical results underline our analytical findings. Moreover, we present a novel analysis for the parametric regularity in the lognormal setting.