La aproximación de Taylor-Hood del flujo de Stokes es un método mixto de elementos finitos, para el cual las estimaciones de error generalmente tienen la forma
donde es exactamente solución y es la aproximación. Para el flujo de Stokes, y (con valor medio cero). Para el elemento - -Taylor-Hood, consta de polinomios cuadráticos por partes continuas y
∥u−uh∥V+∥p−ph∥M≤C(infwh∈Vh∥u−wh∥V+infqh∈Mh∥p−qh∥M),(1)
(u,p)∈V×M(uh,ph)∈Vh×MhV=H10(Ω)dM=L2(Ω)P2P1VhMhde polinomios lineales por partes continuos, para los cuales ambos términos en el lado derecho pueden estar delimitados por errores de aproximación cuadrática utilizando argumentos estándar (por ejemplo, lema de Bramble-Hilbert y reglas de transformación):
(regla general "número de derivadas a la izquierda potencias de número de derivadas a la derecha"). Insertar esto en produce
(suponiendo que la solución exacta es de hecho lo suficientemente regular).
infwh∈Vh∥u−wh∥H1infqh∈Mh∥p−qh∥L2≤Ch2∥u∥H3≤Ch2∥p∥H2
+h =(1)∥u−uh∥H1+∥p−ph∥L2≤Ch2(∥u∥H3+∥p∥H2).
Dado que tanto el problema de Stokes continuo como el discreto son un sistema lineal acoplado triangular superior de la forma
la estimación del error (que explota que la diferencia de las soluciones satisface un sistema similar para y y usa la invertibilidad de en el núcleo de ) para depende en general del error en .
Au+B∗pBu=f=0
AhBhABu−uhp−ph
Sin embargo, si observa la prueba, hay un agujero de bucle: si el espacio nulo de está contenido en el espacio nulo de , entonces el término de acoplamiento de hecho desaparece, y obtiene una estimación de error que involucra solo a :
Desde allí, puede aplicar el truco estándar de Aubin-Nitsche (si la ecuación adjunta está bien planteada, que es el caso si el dominio es lo suficientemente regular (un polígono convexo en 2D o tiene un límite que puede ser parametrizado por una función diferenciable de Lipschitz) para obtener una tasa de convergencia para el error de un orden superior:
BhBu
∥u−uh∥H1≤Ch2∥u∥H3
ΩL2∥u−uh∥L2≤Ch3∥u∥H3
Puede encontrar estos resultados en Ern, Guermond: Theory and Practice of Finite Elements , Springer, 2004 . (Las estimaciones de error se recopilan en el Teorema 4.26, mientras que la regularidad necesaria para se define en el Lema 4.17; desafortunadamente, las pruebas están dispersas en el libro, y creo que no se verifica explícitamente en ninguna parte).ΩkerBh⊂kerB