Conforming Discretizations of Mixed-Dimensional Partial Differential Equations
MetadataVis full innførsel
Mixed-dimensional partial differential equations (PDEs) are coupled equations defined on connected manifolds of different dimensionalities. Two main examples of mixed-dimensional PDEs are considered in this dissertation, namely flow in fractured porous media and mechanics of composite materials. We focus on the discretization of these examples using hierarchical finite elements defined on coupled manifolds of codimension one, successively. By uncovering their underlying structure, we use the corresponding tools to define, analyze and discretize mixed-dimensional partial differential equations. Our first example concerning mixed-dimensional PDEs arises in the context of fracture flow. Here, the planar fractures, intersection lines, as well as intersection points are represented as lower-dimensional manifolds. In turn, the entire embedded fracture network forms a mixed-dimensional geometry. We continue by defining the conservation and constitutive laws on the mixed-dimensional geometry, leading to a hierarchically coupled system of partial differential equations. Next, we extend these concepts from flow to derive the governing equations concerning mechanics of materials with thin inclusions in an analogous manner. Together, the embedded features and their surroundings form the mixed-dimensional geometry and the behavior of the system can be captured by prescribing significantly different material parameters. The analysis of these systems introduces several new concepts including mixed-dimensional function spaces and semi-discrete differential operators. With the aim of discretization, we use finite element exterior calculus to construct mixed finite element schemes on the mixed-dimensional geometry. We focus on two families of mixed-dimensional finite elements, hierarchically ordered by dimensionality. We refer to these families as the first and second kind and show that both are of interest in the context of fracture flow, with different behavior in terms of convergence and computational cost. On the other hand, the mixed formulation of the mechanics equations requires the family of elements of the second kind. For fracture flow, stability and optimal convergence of the discretization method are shown with the use of weighted, mixed-dimensional Sobolev spaces. A novel way of incorporating the fracture aperture leads to a scheme capable of handling arbitrarily small and spatially varying apertures. In case of fractures pinching out, the degeneration of the equations eliminates the possibility for flow resulting in a natural termination of fractures. In a benchmark study concerning flow through fractured porous media, the proposed scheme is compared to various other numerical methods. Four two-dimensional test cases of varying complexity are considered, specifically designed to highlight the typical difficulties with modeling flow through fracture networks.With eight participating numerical schemes, a clear view is given of the performance and limitations of the state-of-the-art numerical schemes. Finally, we consider the evolution of the water table and identify the water table itself as a lower-dimensional manifold. Its location is governed by a partial differential equation which is coupled to the underlying saturated region. To solve this problem, a numerical scheme is proposed which maps the problem to a stationary reference domain. We analyze the properties of this scheme and successfully apply it to a real world problem concerning ground flow patterns surrounding meandering streams.
Består avPaper I: W.M. Boon, J.M.Nordbotten, I.Yotov (2018), Robust Discretization of Flow in Fractured Porous Media, SIAM Journal of Numerical Analysis, 56(4), 2203–2233. The article is available at: http://hdl.handle.net/1956/18157
Paper II: W.M. Boon, J.M. Nordbotten, Stable Mixed Finite Elements for Linear Elasticity with Thin Inclusions. The article is not available in BORA.
Paper III: J.M. Nordbotten, W.M. Boon, Modeling, Structure and Discretization of Mixed-dimensional Partial Differential Equations, Domain Decomposition Methods in Science and Engineering XXIV, Lecture Notes in Computational Science and Engineering. The article is available in the main thesis.
Paper IV: B. Flemisch, I. Berre, W.M. Boon, A. Fumagalli, N. Schwenck, A. Scotti, I. Stefansson, A. Tatomir (2018), Benchmarks for Single-phase Flow in Fractured Porous Media, Advances in Water Resources 111, 239–258. The article is available in the main thesis. The article is also available at: https://doi.org/10.1016/j.advwatres.2017.10.036
Paper V: J.M. Nordbotten, W.M. Boon, A. Fumagalli, E. Keilegavlen, Unified Approach to Discretization of Flow in Fractured Porous Media, Computational Geosciences. The article is available in the main thesis.
Paper VI: W.M. Boon, N. Balbarini, P.J. Binning, J.M. Nordbotten (2017), Efficient Water Table Evolution Discretization using Domain Transformation, Computational Geosciences 21(1), 3–11. The article is available at: http://hdl.handle.net/1956/18158
Paper VII: N. Balbarini, W.M. Boon, P.L. Bjerg, J.M. Nordbotten, P.J. Binning (2017), A 3-D Model of the Influence of Meanders on Groundwater Discharge to a Gaining Stream in an Unconfined Sandy Aquifer, Journal of Hydrology 552, 168–181. The article is available in the main thesis. The article is also available at: https://doi.org/10.1016/j.jhydrol.2017.06.042
Paper VIII: W.M. Boon, J.M. Nordbotten, An Adaptive Penalty Method for Constrained Minimization Problems, Optimization Letters. Full text not available in BORA.