Discrete Optimization of Robin Transmission Conditions for Anisotropic Diffusion with Discrete Duality Finite Volume Methods
- Others:
- Section de mathématiques [Genève] ; Université de Genève = University of Geneva (UNIGE)
- Laboratoire Analyse, Géométrie et Applications (LAGA) ; Université Paris 8 Vincennes-Saint-Denis (UP8)-Centre National de la Recherche Scientifique (CNRS)-Université Sorbonne Paris Nord
- Institut de Mathématiques de Marseille (I2M) ; Aix Marseille Université (AMU)-École Centrale de Marseille (ECM)-Centre National de la Recherche Scientifique (CNRS)
- COmplex Flows For Energy and Environment (COFFEE) ; Inria Sophia Antipolis - Méditerranée (CRISAM) ; Institut National de Recherche en Informatique et en Automatique (Inria)-Institut National de Recherche en Informatique et en Automatique (Inria)-Laboratoire Jean Alexandre Dieudonné (JAD) ; Université Nice Sophia Antipolis (1965 - 2019) (UNS) ; COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Centre National de la Recherche Scientifique (CNRS)-Université Côte d'Azur (UCA)-Université Nice Sophia Antipolis (1965 - 2019) (UNS) ; COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Centre National de la Recherche Scientifique (CNRS)-Université Côte d'Azur (UCA)
Description
Discrete Duality Finite Volume (DDFV) methods are very well suited to discretize anisotropic diffusion problems, even on meshes with low mesh quality. Their performance stems from an accurate reconstruction of the gradients between mesh cell boundaries, which comes however at the cost of using both a primal (cell centered) and a dual (vertex centered) mesh, and thus leads to larger system sizes. To solve these systems, we propose to use non-overlapping optimized Schwarz methods with Robin transmission conditions, which can also well take into account anisotropic diffusion across subdomain interfaces. We study these methods here directly at the discrete level, and prove convergence using energy estimates for general decompositions including cross points and fully anisotropic diffusion. Our analysis reveals that primal and dual meshes might be coupled using different optimized Robin parameters in the optimized Schwarz methods. We present both the separate and coupled optimization of Robin transmission conditions and derive parameters which lead to the fastest possibe convergence in each case. We illustrate our analysis with numerical experiments for the model problem covered by our analysis, and also in situations that go beyond our analysis, with an application to anisotropic image reconstruction.
Abstract
International audience
Additional details
- URL
- https://hal.archives-ouvertes.fr/hal-02539124
- URN
- urn:oai:HAL:hal-02539124v2
- Origin repository
- UNICA