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

Resumen Técnico: Cota Superior de Derivadas de Alto Orden para Coordenadas de Wachspress en Polítopos

Enunciado del Problema
Las coordenadas baricéntricas generalizadas (GBC) son fundamentales para construir métodos de elementos finitos (MEF) en mallas polipolares. Si bien los límites de gradiente (derivadas de primer orden) de las GBC están bien establecidos y son suficientes para estimaciones de error en norma $H^1$ en ecuaciones elípticas de segundo orden, las aplicaciones de MEF de orden superior (por ejemplo, ecuaciones elípticas de cuarto orden como la flexión de placas) requieren cotas superiores para derivadas de alto orden (normas $H^k$, $k > 1$). Antes de este trabajo, tales cotas no estaban disponibles en la literatura, creando una brecha teórica en el análisis de error de MEF no conformes basados en GBC, como el elemento tipo Morley propuesto en [34].

Metodología
Los autores se centran en las coordenadas de Wachspress definidas en polítopos convexos, simples y no degenerados de dimensión $d$. La derivación de cotas superiores para derivadas de orden arbitrario ($D^\alpha \phi_v$) implica los siguientes pasos técnicos:

1. Descomposición Algebraica: La coordenada de Wachspress $\phi_v$ se expresa como un cociente $w_v / W$, donde $w_v$ es una función de peso polinómica y $W$ es la suma de todos los pesos. Los autores utilizan la fórmula general de Leibniz y la fórmula multivariante de Faà di Bruno para expandir las derivadas de este cociente. Esto reduce el problema a estimar las derivadas del numerador ($w_v$) y del denominador ($W$).
2. Parámetros Geométricos: Las cotas se derivan en términos de dos cantidades geométricas clave:
   • $h_K$: El diámetro del polítopo $K$.
   • $h^*$: La distancia mínima entre vértices no incidentes y facetas de $K$.
3. Estimación de Derivadas:
   • Numerador ($w_v$): Dado que $w_v$ es un producto de funciones de distancia lineales, sus derivadas se acotan utilizando argumentos combinatorios que involucran el número de facetas y el diámetro $h_K$.
   • Denominador ($W$): Un paso crítico consiste en establecer una cota inferior estrictamente positiva para $W$. Los autores demuestran que $W(x) \geq C \frac{(h^*)^{|F|}}{h_K^d}$ analizando la distribución de las facetas con respecto a un punto $x$, utilizando la propiedad de precisión lineal de las GBC y la geometría de los polítopos simples.
4. Síntesis: Combinando las cotas superiores de las derivadas de $w_v$ y $W$ con la cota inferior de $W$, los autores derivan una cota superior explícita para $|D^\alpha \phi_v|$.

Contribuciones Clave

1. Derivación de Cotas de Alto Orden: El artículo proporciona las primeras cotas superiores explícitas para derivadas de orden arbitrario de las coordenadas de Wachspress en polítopos convexos simples. El resultado principal (Teorema 3.7) establece que cuando $h^* = O(h_K)$ (una condición que implica que el polítopo no es "delgado" o degenerado de una manera específica), las derivadas satisfacen:
   $$|D^\alpha \phi_v| \leq C h_K^{-|\alpha|}$$
   donde $C$ depende de la dimensión $d$, el número de vértices $|V|$, el número de facetas $|F|$ y el multiíndice $\alpha$.
2. Análisis de Regularidad Geométrica: El artículo compara sistemáticamente varias condiciones de regularidad de forma para polítopos convexos (por ejemplo, longitud mínima de arista, conteo máximo de vértices, parámetro de robustez). Establece implicaciones lógicas entre estas condiciones, mostrando que la suposición $h^* = O(h_K)$ (Suposición H1) es prácticamente razonable e implica otras condiciones de regularidad estándar utilizadas en MEF.
3. Clarificación de la Agudeza: Los autores aclaran que, si bien las cotas derivadas son agudas en orden asintótico cuando $h^* = O(h_K)$, pueden no ser agudas cuando $h^* \ll h_K$. Los experimentos numéricos sugieren que en este último caso, el comportamiento de las derivadas depende de la dirección de las aristas "cortas" con respecto a la dirección de la derivada.

Resultados Numéricos
Los autores validan sus hallazgos teóricos mediante experimentos numéricos en 2D:

• Calidad de la Malla: Demuestran que las Teselaciones de Voronoi Centroidales (CVT), tras eliminar las aristas "cortas", satisfacen naturalmente la condición $h^* = O(h_K)$.
• Verificación de Cotas: Para polígonos aleatorios que satisfacen $h^* = O(h_K)$, los valores calculados de $\Lambda_\alpha = \max \sum |D^\alpha \phi_i|$ se alinean con la predicción teórica $O(h_K^{-|\alpha|})$.
• Casos Degenerados: Los experimentos en polígonos donde $h^* \ll h_K$ revelan que la cota superior teórica es conservadora. El crecimiento real de las derivadas depende de la orientación de las aristas cortas, con algunas derivadas direccionales que exhiben tasas de crecimiento más lentas de lo que sugiere la cota general.
• Aplicación: El artículo incluye resultados de un MEF basado en GBC tipo Morley aplicado a un problema de flexión de placa empotrada (ecuación de cuarto orden). Las tasas de convergencia numérica ($O(h)$ en la seminorma $H^2$ y $O(h^2)$ en la seminorma $H^1$) confirman que las cotas de derivadas derivadas son suficientes para probar la convergencia óptima para tales problemas.

Significado y Afirmaciones
El artículo afirma llenar una brecha teórica crítica en el análisis de métodos de elementos finitos polipolares. Al proporcionar las cotas de derivadas de alto orden necesarias, permite la demostración rigurosa de la convergencia óptima para MEF basados en Wachspress aplicados a ecuaciones elípticas de cuarto orden. Los autores enfatizan que su enfoque evita las técnicas altamente técnicas y específicas del gradiente de trabajos anteriores (como [14]) al utilizar una estrategia algebraica diferente adecuada para órdenes arbitrarios. Además, el trabajo contribuye a la comprensión de la geometría de polítopos al clarificar las relaciones entre varias suposiciones de regularidad de forma, afirmando que la condición $h^* = O(h_K)$ es una suposición robusta y práctica para generar mallas de alta calidad en aplicaciones.