Analysis and approximation of mixed-dimensional PDEs on 3D-1D domains coupled with Lagrange multipliers
Preprint
Abstract
Coupled partial differential equations defined on domains with different dimensionality are usually called mixed dimensional PDEs. We address mixed dimensional PDEs on three-dimensional (3D) and one-dimensional domains, giving rise to a 3D-1D coupled problem. Such problem poses several challenges from the standpoint of existence of solutions and numerical approximation. For the coupling conditions across dimensions, we consider the combination of essential and natural conditions, basically the combination of Dirichlet and Neumann conditions. To ensure a meaningful formulation of such conditions, we use the Lagrange multiplier method, suitably adapted to the mixed dimensional case. The well posedness of the resulting saddle point problem is analyzed. Then, we address the numerical approximation of the problem in the framework of the finite element method. The discretization of the Lagrange multiplier space is the main challenge. Several options are proposed, analyzed and compared, with the purpose to determine a good balance between the mathematical properties of the discrete problem and flexibility of implementation of the numerical scheme. The results are supported by evidence based on numerical experiments.Keywords
All Related Versions
This publication has 22 references indexed in Scilit:
- A singularity removal method for coupled 1D–3D flow modelsComputational Geosciences, 2019
- Modelling solute transport in the brain microcirculation: is it really well mixed inside the blood vessels?Journal of Fluid Mechanics, 2019
- Mathematical analysis, finite element approximation and numerical solvers for the interaction of 3D reservoirs with 1D wellsGEM - International Journal on Geomathematics, 2019
- Well‐Posedness and Global Behavior of the Peskin Problem of an Immersed Elastic Filament in Stokes FlowCommunications on Pure and Applied Mathematics, 2018
- Preconditioning trace coupled 3d‐1d systems using fractional LaplacianNumerical Methods for Partial Differential Equations, 2018
- Local error estimates of the finite element method for an elliptic problem with a Dirac source termNumerical Methods for Partial Differential Equations, 2017
- Projection stabilization of Lagrange multipliers for the imposition of constraints on interfaces and boundariesNumerical Methods for Partial Differential Equations, 2013
- Preconditioning discretizations of systems of partial differential equationsNumerical Linear Algebra with Applications, 2010
- Self-intersection avoidance and integral properties of generalized cylindersComputer Aided Geometric Design, 2002
- The finite element method with Lagrangian multipliersNumerische Mathematik, 1973