¿Qué pasa con esta simple estimación de error para PDE lineal?

10

Sea un dominio Lipschitz convexo poligonalmente limitado en , deje .R 2 f L 2 ( Ω )ΩR2FL2(Ω)

Entonces, la solución del problema de Dirichlet en , on tiene una solución única en y está bien planteada, es decir, para alguna constante tenemos .Ω traza u = 0 Ω H 2 C u H 2C f L 2Δtu=FΩrastrotu=0 0ΩH2CtuH2CFL2

Para alguna aproximación de elementos finitos tuh , digamos, con elementos nodales en una cuadrícula uniforme, tenemos la estimación del error

tu-tuhH1ChtuH2

Parece (tal vez estoy equivocado con eso) que las personas generalmente no usan la estimación de error obvia

tu-tuhH1ChFL2

que podemos obtener mediante la combinación de las dos desigualdades anteriores. En cambio, los estimadores de error a posteriori se desarrollan en varias formas. La única objeción que puedo imaginar contra la ecuación anterior es que la constante C podría en la práctica ser demasiado pesimista o no confiablemente estimable.

shuhalo
fuente

Respuestas:

8

La razón por la cual las personas prefieren usar la primera estimación, en mi opinión, es que la primera surge naturalmente de la ortogonalidad de Galerkin del FEM, la propiedad de aproximación de interpolación y, lo más importante, la coercitividad de la forma bilineal (para el problema del valor límite de la ecuación de Poisson , es equivalente a la desigualdad de Poincaré / Friedrichs para las funciones ): u - u h 2 H 1 ( Ω )H0 01

tu-tuhH1(Ω)2C1(tu-tuh)L2(Ω)2(tu-tuh)L2(Ω)2=Ω(tu-tuh)(tu-tuh)=Ω(tu-tuh)(tu-yotu)(tu-tuh)L2(Ω)(tu-yotu)L2(Ω)(tu-tuh)L2(Ω)(tu-yotu)L2(Ω)C2htuH2(Ω)
donde depende de la constante en la desigualdad de Poincaré / Friedrichs para las funciones , es la interpolación de en el finito elemento espacio, yC1H0 01yotutuC2 depende de los ángulos mínimos de la malla.

Mientras que la estimación de regularidad elíptica está únicamente en el nivel PDE, no tiene nada que ver con la aproximación, más el argumento anterior, se mantiene incluso cuando es una distribución.tuH2(Ω)CFL2(Ω)FH-1

Ahora pasemos a la razón por la cual las estimaciones de error a posteriori se usan ampliamente, principalmente porque:

  • Es computable, no hay una constante genérica en la expresión de las estimaciones.

  • El estimador tiene su forma local, que podría ser el indicador de error local que se utiliza en el procedimiento de refinación de malla adaptativa. Por lo tanto, el problema con singularidades o geometrías realmente "malas" podría ser tratado.

Las dos estimaciones de tipo a priori que usted enumeró son válidas, nos proporcionan la información de los órdenes de convergencia, sin embargo, ninguna de ellas podría ser un indicador de error local solo para un triángulo / tetraedro, porque ninguna de ellas es computable debido a la constante , ni se definen localmente.

EDITAR: Para obtener una visión más general de la FEM para PDE elípticas, recomiendo leer el Capítulo 0 en el libro de Brenner y Scott: La teoría matemática de los métodos de elementos finitos , que consta de solo 20 páginas y cubre brevemente casi todos los aspectos de los métodos de elementos finitos. , desde la formulación de Galerkin del PDE, hasta la motivación por la que nos gustaría utilizar FEM adaptativo para abordar algún problema. Espero que esto te ayude más.

Shuhao Cao
fuente
1

Su estimación es demasiado pesimista en dos frentes. Ya ha identificado el primero ( ahora no solo incluye la constante de interpolación sino también la constante de estabilidad). El segundo es que la estimación del error realmente lee Tenga en cuenta que el lado derecho tiene el , no la norma. Por supuesto, puede atar los rhs por la norma completa, pero pierde de nuevo de esta manera.C

miL2ChEl |tuEl |H2.
H2
Wolfgang Bangerth
fuente