Upper bound of high-order derivatives for Wachspress coordinates on polytopes

Résumé technique : Majoration des dérivées d'ordre élevé pour les coordonnées de Wachspress sur les polytopes

Énoncé du problème
Les coordonnées barycentriques généralisées (GBC) sont fondamentales pour la construction de méthodes d'éléments finis (MEF) sur des maillages polytopaux. Bien que les bornes de gradient (dérivées d'ordre un) des GBC soient bien établies et suffisantes pour les estimations d'erreur en norme $H^1$ dans le cas d'équations elliptiques d'ordre deux, les applications MEF d'ordre supérieur (par exemple, les équations elliptiques d'ordre quatre comme la flexion de plaques) nécessitent des majorations pour les dérivées d'ordre élevé (normes $H^k$, $k > 1$). Avant ce travail, de telles bornes n'étaient pas disponibles dans la littérature, créant un vide théorique dans l'analyse d'erreur des MEF non conformes basées sur les GBC, telles que l'élément de type Morley proposé dans [34].

Méthodologie
Les auteurs se concentrent sur les coordonnées de Wachspress définies sur des polytopes convexes, simples et non dégénérés de dimension $d$. La dérivation des majorations pour les dérivées d'ordre arbitraire ($D^\alpha \phi_v$) implique les étapes techniques suivantes :

1. Décomposition algébrique : La coordonnée de Wachspress $\phi_v$ est exprimée comme un rapport $w_v / W$, où $w_v$ est une fonction de poids polynomiale et $W$ est la somme de tous les poids. Les auteurs utilisent la formule de Leibniz généralisée et la formule multivariée de Faà di Bruno pour développer les dérivées de ce quotient. Cela réduit le problème à l'estimation des dérivées du numérateur ($w_v$) et du dénominateur ($W$).
2. Paramètres géométriques : Les bornes sont dérivées en fonction de deux grandeurs géométriques clés :
   • $h_K$ : Le diamètre du polytope $K$.
   • $h^*$ : La distance minimale entre les sommets et les facettes non incidents de $K$.
3. Estimation des dérivées :
   • Numérateur ($w_v$) : Puisque $w_v$ est un produit de fonctions de distance linéaires, ses dérivées sont majorées à l'aide d'arguments combinatoires impliquant le nombre de facettes et le diamètre $h_K$.
   • Dénominateur ($W$) : Une étape critique consiste à établir une borne inférieure strictement positive pour $W$. Les auteurs démontrent que $W(x) \geq C \frac{(h^*)^{|F|}}{h_K^d}$ en analysant la distribution des facettes par rapport à un point $x$, en utilisant la propriété de précision linéaire des GBC et la géométrie des polytopes simples.
4. Synthèse : En combinant les majorations des dérivées de $w_v$ et $W$ avec la minoration de $W$, les auteurs déduisent une majoration explicite pour $|D^\alpha \phi_v|$.

Contributions clés

1. Dérivation de bornes d'ordre élevé : L'article fournit les premières majorations explicites pour les dérivées d'ordre arbitraire des coordonnées de Wachspress sur des polytopes convexes simples. Le résultat principal (Théorème 3.7) énonce que lorsque $h^* = O(h_K)$ (une condition impliquant que le polytope n'est pas « mince » ou dégénéré d'une manière spécifique), les dérivées satisfont :
   $$|D^\alpha \phi_v| \leq C h_K^{-|\alpha|}$$
   où $C$ dépend de la dimension $d$, du nombre de sommets $|V|$, du nombre de facettes $|F|$ et de l'indice multiple $\alpha$.
2. Analyse de la régularité géométrique : L'article compare systématiquement diverses conditions de régularité de forme pour les polytopes convexes (par exemple, longueur minimale des arêtes, nombre maximal de sommets, paramètre d'épaisseur). Il établit des implications logiques entre ces conditions, montrant que l'hypothèse $h^* = O(h_K)$ (Hypothèse H1) est pratiquement raisonnable et implique d'autres conditions de régularité standard utilisées en MEF.
3. Précision de la netteté : Les auteurs clarifient que si les bornes dérivées sont optimales en ordre asymptotique lorsque $h^* = O(h_K)$, elles peuvent ne pas l'être lorsque $h^* \ll h_K$. Des expériences numériques suggèrent que dans ce dernier cas, le comportement des dérivées dépend de la direction des arêtes « courtes » par rapport à la direction de la dérivée.

Résultats numériques
Les auteurs valident leurs résultats théoriques par des expériences numériques en 2D :

• Qualité du maillage : Ils démontrent que les tessellations de Voronoï centrées (CVT), après élimination des arêtes « courtes », satisfont naturellement la condition $h^* = O(h_K)$.
• Vérification des bornes : Pour des polygones aléatoires satisfaisant $h^* = O(h_K)$, les valeurs calculées de $\Lambda_\alpha = \max \sum |D^\alpha \phi_i|$ s'alignent avec la prédiction théorique $O(h_K^{-|\alpha|})$.
• Cas dégénérés : Des expériences sur des polygones où $h^* \ll h_K$ révèlent que la majoration théorique est conservative. La croissance réelle des dérivées dépend de l'orientation des arêtes courtes, certaines dérivées directionnelles présentant des taux de croissance plus lents que ce que la borne générale suggère.
• Application : L'article inclut des résultats issus d'une MEF basée sur les GBC de type Morley appliquée à un problème de flexion de plaque encastrée (équation d'ordre quatre). Les taux de convergence numériques ($O(h)$ en semi-norme $H^2$ et $O(h^2)$ en semi-norme $H^1$) confirment que les bornes de dérivées dérivées sont suffisantes pour prouver la convergence optimale pour de tels problèmes.

Importance et affirmations
L'article prétend combler un vide théorique critique dans l'analyse des méthodes d'éléments finis sur polytopes. En fournissant les bornes de dérivées d'ordre élevé nécessaires, il permet la preuve rigoureuse de la convergence optimale pour les MEF basées sur Wachspress appliquées aux équations elliptiques d'ordre quatre. Les auteurs soulignent que leur approche évite les techniques hautement techniques et spécifiques au gradient des travaux précédents (comme [14]) en utilisant une stratégie algébrique différente adaptée aux ordres arbitraires. De plus, ce travail contribue à la compréhension de la géométrie des polytopes en clarifiant les relations entre diverses hypothèses de régularité de forme, affirmant que la condition $h^* = O(h_K)$ est une hypothèse robuste et pratique pour générer des maillages de haute qualité dans les applications.