-convergencia del método de elementos finitos cuando el lado derecho está solo en(Poisson eqn)

Sé que la aproximación de elementos finitos lineales por partes de satisface siempre que U sea ​​lo suficientemente suave y f \ en L ^ 2 (U) .$u_h$

$$\Delta u(x)=f(x)\quad\text{in }U\\ u(x)=0\quad\text{on }\partial U$$ $$\|u-u_h\|_{H^1_0(U)}\leq Ch\|f\|_{L^2(U)}$$ $U$$f\in L^2(U)$

Pregunta: Si $f\in H^{-1}(U)\setminus L^2(U)$ , ¿tenemos la siguiente estimación análoga, en la que se quita una derivada en ambos lados:

$$\|u-u_{h}\|_{L^2(U)}\leq Ch\|f\|_{H^{-1}(U)}\quad?$$

puedes proporcionar referencias?

Pensamientos: Como todavía tenemos $u\in H^1_0(U)$ , debería ser posible obtener convergencia en $L^2(U)$ . Intuitivamente, esto debería ser posible incluso con funciones constantes por partes.

Sí , este es el truco estándar de Aubin-Nitsche (o dualidad ). La idea es utilizar el hecho de que es su propio espacio dual para escribir la forma como norma de operador Por lo tanto, tenemos que estimar para arbitrario . Para hacer eso, " " a considerando primero para arbitraria la solución del problema dual $L^2$$L^2$

$$\|u\|_{L^2} = \sup_{\phi\in L^2\setminus\{0\}} \frac{(u,\phi)}{\|\phi\|_{L^2}}.$$ $(u-u_h,\phi)$$\phi\in L^2$$u-u_h$$H^1_0$$\phi\in L^2$$w_\phi\in H^1_0$ $$\label{eq1}\tag{1} (\nabla w_\phi,\nabla v) =(\phi,v) \qquad\text{for all }v\in H^1_0.$$ Usando la regularidad estándar de la ecuación de Poisson, sabemos que $$\|w_\phi\|_{H^2}\leq C \|\phi\|_{L^2}.$$

Insertar en y usar la ortogonalidad de Galerkin para cualquier elemento finito (en su caso, función lineal por partes) produce la estimación Dado que esto es válido para todo , la desigualdad sigue siendo cierta si tomamos el mínimo sobre todos los lineales por . Por lo tanto obtenemos $v=u-u_h\in H^1_0$$\eqref{eq1}$$w_h$

$$\begin{aligned} (\phi,u-u_h) &= (\nabla w_\phi,\nabla (u-u_h))\\ &= (\nabla w_\phi-\nabla w_h,\nabla (u-u_h))\\ &\leq C \|u-u_h\|_{H^1}\|w_\phi-w_h\|_{H^1}. \end{aligned}$$ $w_h$$w_h$ $$\label{eq2}\tag{2} \|u-u_h\|_{L^2} = \sup_{\phi\in L^2\setminus\{0\}} \frac{(u-u_h,\phi)}{\|\phi\|_{L^2}} \leq C \|u-u_h\|_{H^1} \sup_{\phi\in L^2\setminus\{0\}} \frac{\inf_{w_h}\|w_\phi-w_h\|_{H^1}}{\|\phi\|_{L^2}}.$$ Este es el Aubin-Nitsche-Lemma .

El siguiente paso es utilizar estimaciones de error estándar para la mejor aproximación de elementos finitos de soluciones a la ecuación de Poisson. Como solo está en , no obtenemos una estimación mejor que Pero afortunadamente, podemos usar el hecho de que tiene una mayor regularidad desde el lado derecho lugar de . En este caso, tenemos Insertar y en$u$$H^1$

$$\label{eq3}\tag{3} \|u-u_h\|_{H^1} \leq \inf_{v_h}\|u-v_h\|_{H^1} \leq c \|u\|_{H^1} \leq C \|f\|_{H^{-1}}.$$ $w_\phi$$\phi\in L^2$$H^{-1}$ $$\label{eq4}\tag{4} \inf_{w_h}\|w_\phi-w_h\|_{H^1} \leq c h \|w_\phi\|_{H^2} \leq C h \|\phi\|_{L^2}$$ $\eqref{eq3}$$\eqref{eq4}$$\eqref{eq2}$ ahora produce la estimación deseada.

(Tenga en cuenta que las estimaciones estándar requieren que el grado polinómico de la aproximación de elementos finitos y el exponente de Sobolev de la solución verdadera satisfagan , por lo que este argumento no funciona para la aproximación constante por partes ( ). También hemos usado que , es decir, que tenemos una aproximación conforme , lo que no es cierto para las constantes por partes).$k$$m$$m<k+1$$k=0$$u-u_h\in H^1_0$

Como solicitó una referencia: puede encontrar una declaración (incluso para espacios negativos de Sobolev lugar de ) en el Teorema 5.8.3 (junto con el Teorema 5.4.8) en$H^{-s}$$L^2$

Susanne C. Brenner y L. Ridgway Scott , MR 2373954 La teoría matemática de los métodos de elementos finitos , Textos en Matemática Aplicada ISBN: 978-0-387-75933-3.