Well-Posedness and Discretization for a Class of Models for Mixed-Dimensional Problems with High-Dimensional Gap

In this work, we illustrate the underlying mathematical structure of mixed-dimensional models arising from the composition of graphs and continuous domains. Such models are becoming popular in applications, in particular, to model the human vasculature. We first discuss the model equations in the strong form which describes the conservation of mass and Darcy's law in the continuum and network as well as the coupling between them. By introducing proper scaling, we propose a weak form that avoids degeneracy. Well-posedness of the weak form is shown through standard Babuska-Brezzi theory. We also develop the mixed formulation finite-element method and prove its well-posedness. A mass-lumping technique is introduced to derive the two-point flux approximation type discretization as well, due to its importance in applications. Based on the Babuska-Brezzi theory, error estimates can be obtained for both the finite-element scheme and the TPFA scheme. We also discuss efficient linear solvers for discrete problems. Finally, we present some numerical examples to verify the theoretical results and demonstrate the robustness of our proposed discretization schemes.

[1]  Jinchao Xu,et al.  Robust preconditioners for incompressible MHD models , 2015, J. Comput. Phys..

[2]  Barbara Wohlmuth,et al.  Mathematical modelling, analysis and numerical approximation of second order elliptic problems with inclusions , 2018 .

[3]  Axel R Pries,et al.  A hybrid discrete-continuum approach for modelling microcirculatory blood flow. , 2019, Mathematical medicine and biology : a journal of the IMA.

[4]  Jan M. Nordbotten,et al.  Robust Discretization of Flow in Fractured Porous Media , 2016, SIAM J. Numer. Anal..

[5]  Anders M. Dale,et al.  Cortical Surface-Based Analysis I. Segmentation and Surface Reconstruction , 1999, NeuroImage.

[6]  Andrew J. Wathen,et al.  Analysis of Preconditioners for Saddle-Point Problems , 2004, SIAM J. Sci. Comput..

[7]  Patrick Jenny,et al.  Modeling tissue perfusion in terms of 1d-3d embedded mixed-dimension coupled problems with distributed sources , 2020, J. Comput. Phys..

[8]  Michel Fortin,et al.  Error analysis of piecewise constant pressure approximations of Darcy’s law , 2006 .

[9]  Jan Mar,et al.  Mixed-dimensional models for real-world applications , 2019 .

[10]  Todd Arbogast,et al.  Mixed Methods for Two-Phase Darcy-Stokes Mixtures of Partially Melted Materials with Regions of Zero Porosity , 2017, SIAM J. Sci. Comput..

[11]  Jean E. Roberts,et al.  Mixed and hybrid methods , 1991 .

[12]  Erlend Hodneland,et al.  A new framework for assessing subject-specific whole brain circulation and perfusion using MRI-based measurements and a multi-scale continuous flow model , 2019, PLoS Comput. Biol..

[13]  Nikos Stergiopulos,et al.  Pulse Wave Propagation in the Arterial Tree , 2011 .

[14]  Barbara I. Wohlmuth,et al.  A Local Error Estimate for the Poisson Equation with a Line Source Term , 2015, ENUMATH.

[15]  Ingeborg G. Gjerde,et al.  Splitting method for elliptic equations with line sources , 2018, ESAIM: Mathematical Modelling and Numerical Analysis.

[16]  Xiaozhe Hu,et al.  Algebraic Multigrid for Least Squares Problems on Graphs with Applications to HodgeRank , 2017, 2017 IEEE International Parallel and Distributed Processing Symposium Workshops (IPDPSW).

[17]  M. Fortin,et al.  Mixed Finite Element Methods and Applications , 2013 .

[18]  Radim Blaheta,et al.  A multilevel method with correction by aggregation for solving discrete elliptic problems , 1986 .

[19]  Ludmil T. Zikatanov,et al.  Algebraic multigrid methods * , 2016, Acta Numerica.

[20]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[21]  L. Zikatanov,et al.  Robust Block Preconditioners for Biot's Model , 2017, CSE 2017.

[22]  Jinchao Xu,et al.  Comparative Convergence Analysis of Nonlinear AMLI-Cycle Multigrid , 2012, SIAM J. Numer. Anal..

[23]  Ludmil T. Zikatanov,et al.  An Adaptive Multigrid Method Based on Path Cover , 2018, SIAM J. Sci. Comput..

[24]  A. Quarteroni,et al.  On the coupling of 1D and 3D diffusion-reaction equations. Applications to tissue perfusion problems , 2008 .

[25]  J. Maître,et al.  Connection between finite volume and mixed finite element methods , 1996 .

[26]  J. Nordbotten,et al.  Functional analysis and exterior calculus on mixed-dimensional geometries , 2017, Annali di Matematica Pura ed Applicata (1923 -).

[27]  Daniele Bertaccini,et al.  Non-local network dynamics via fractional graph Laplacians , 2019, J. Complex Networks.

[28]  Paolo Zunino,et al.  Derivation and analysis of coupled PDEs on manifolds with high dimensionality gap arising from topological model reduction , 2019, ESAIM: Mathematical Modelling and Numerical Analysis.

[29]  Ludmil T. Zikatanov,et al.  A multigrid method based on graph matching for convection–diffusion equations , 2003, Numer. Linear Algebra Appl..

[30]  D. Owen Handbook of Mathematical Functions with Formulas , 1965 .

[31]  Lenore Cowen,et al.  Computing the diffusion state distance on graphs via algebraic multigrid and random projections , 2018, Numer. Linear Algebra Appl..

[32]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[33]  Jean E. Roberts,et al.  Mixed and hybrid finite element methods , 1987 .

[34]  R. Hiptmair,et al.  MULTIGRID METHOD FORH(DIV) IN THREE DIMENSIONS , 1997 .