Mortaring Method for Non-Conforming Meshes
Hi Christophe and the GetDP community,
We are currently looking into mortaring methods in order to handle non-conforming meshes on interface curves/surfaces for scalar problems. To achieve the interface matching, we would need to compute surface/line integrals of nodal basis/test functions of the form
I_{12} = \int_{\Gamma_1} v_1 \, \lambda_2' \, \text{d}\Gamma,
where v_1
is supported on \Gamma_1
while \lambda_2'
is supported on \Gamma_2
and the meshes discretizing \Gamma_1
and \Gamma_2
(which may be coincident) are non-matching. Hence, it is not possible (?) to use a Link
constraint. To compute I_{12}
, we would thus need an appropriate quadrature rule. The interested reader can find more information in Section 2.5.1 (in particular Remark 2.5.1) in [1].
Do you know if someone has done mortaring in GetDP in the past without using iterative methods? Would you have some insights on how to tackle integrals such as I_{12}
?
Please let me know if you need more details, many thanks for all your help!
[1]: A. Quarteroni and A. Valli, “Domain Decomposition Methods for Partial Differential Equations,” Oxford Science Publications, Oxford, 1999.