Tasas de convergencia de elementos finitos para problemas mixtos


8

Codifiqué un problema de Stokes Flow usando elementos finitos y estoy en el proceso de verificar que funcione. Simplemente no estoy seguro de qué tasa de convergencia debería esperar ya que refino globalmente la malla.

Sé que para problemas escalares usando funciones de base lineal esperaría convergencia de orden ( es el tamaño del elemento), y usando funciones de base cuadrática esperaría convergencia de orden en la norma y una potencia menos en el . El problema que tengo ahora es que al codificar el flujo de Stokes utilicé el elemento Taylor-Hood, que utiliza lineales para la presión y cuadrático para los componentes de velocidad. ¿Es tan simple como las velocidades que convergen en y la presión en el orden ?h2hh3L2H1h3h2

Publiqué esto en mathoverflow primero y me dijeron que podría ser más adecuado para este foro.

Respuestas:


7

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

(1)uuhV+pphMC(infwhVhuwhV+infqhMhpqhM),
(u,p)V×M(uh,ph)Vh×MhV=H01(Ω)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).
infwhVhuwhH1Ch2uH3infqhMhpqhL2Ch2pH2
+h =(1)
uuhH1+pphL2Ch2(uH3+pH2).

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+Bp=fBu=0
AhBhABuuhpph

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

uuhH1Ch2uH3
ΩL2
uuhL2Ch3uH3

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).ΩkerBhkerB


Parece que pones mucho trabajo en compilar esta respuesta. Realmente lo aprecio. ¿Puede aclarar de dónde obtuvo la estimación: , ¿o eso se acaba de dar? uuhH1+pphL2Ch2(uH3+pH2)
Lukas Bystricky

1
No hay problema, me alegro si es de alguna ayuda. La estimación proviene de conectar los errores de interpolación estándar en la primera estimación; He añadido ese paso.
Christian Clason

Experimentalmente (usando un dominio de caja que no puede ser parametrizado por una función diferenciable) convergencia en y convergencia en para todo el sistema. h2L2hH1
Lukas Bystricky

1
Hmm, entonces no tengo mucha idea en este momento. Esta es realmente una pregunta nueva, que las personas más familiarizadas con los problemas de Stokes están mejor equipadas para responder. ¿Por qué no hace una nueva pregunta (use el enlace a continuación), donde describe exactamente qué problema está resolviendo (lado derecho, condiciones de contorno, geometría, solución exacta), qué espera ver y qué está obteniendo en su lugar ? (Si el código deal.II es lo suficientemente pequeño, también puede publicarlo; algunos de los principales desarrolladores son colaboradores frecuentes aquí.)
Christian Clason

1
Resulta que estaba mirando los componentes incorrectos en mi solución. Cuando miro los componentes correctos obtengo en y en para las velocidades y 1 potencia menos para la presión. h3L2h2H1
Lukas Bystricky
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.