Imponer las condiciones de compatibilidad para el método de elementos finitos mixtos en la ecuación de Stokes

13

Supongamos que tenemos la siguiente ecuación del modelo de flujo de Stokes:

{-reyov(νtu)+pag=Freyovtu=0 0
donde la viscosidad ν(X) es una función, para el elemento finito mixto estándar, digamos que usamos el par estable: espacio Crouzeix-Raviart Vh para la velocidad tu y el espacio constante en función del elemento Sh para la presión pag , tenemos la siguiente forma variacional:

L([tu,pag],[v,q])=Ωνtu:v-Ωqreyovtu-Ωpagreyovv=ΩFvv×qVh×Sh

y sabemos que dado que el multiplicador de Lagrange pag puede determinarse hasta una constante, la matriz finalmente ensamblada debe tener un espacio nulo 1 , para evitar esto podríamos hacer que la presión pag sobre algún elemento determinado sea cero, de modo que no tengamos que resolver un sistema singular

Así que aquí está mi pregunta 1:

  • (P1) ¿Hay otra manera que hacer cumplir pag=0 0 en algún elemento para eliminar el núcleo para el elemento finito mixto estándar? o digamos, ¿algún solucionador que pueda resolver el sistema singular para obtener una solución compatible? (o algunas referencias son bienvenidas)

Y sobre la compatibilidad, para (1) debería ser y el pequeño truco es calcular ser la que obtuvimos de la solución de el sistema lineal se resta por su promedio ponderado: ˜ p p ˜ p = p - ν

Ων-1pag=0 0
pag~pag
(2)pag~=pag-νEl |ΩEl |Ων-1pag

Sin embargo, recientemente acabo de implementar un elemento finito mixto estabilizado para la ecuación de Stokes por Bochev, Dohrmann y GunzbergerPAG1-PAG0 0˜ L ([u,p],[v,q])= L ([u,p],[v,q])-Ω (p- Π 1 p)(q- Π 1 q)=Ω fv , en el que agregaron un término estabilizado a la formulación variacional (1):

L~([tu,pag],[v,q])=L([tu,pag],[v,q])-Ω(pag-Π1pag)(q-Π1q)=ΩFvv×qVh×Sh
donde es la proyección desde el espacio constante por partes hasta el continuo por partes , y el núcleo constante del elemento finito mixto original desapareció, sin embargo, sucedieron cosas extrañas, (2) no ya no funciona, acuñé el problema de prueba deΠ1PAG0 0PAG1un problema de interfaz para la ecuación de difusión , esto es lo que obtuve para la presión , la derecha es la solución verdadera y la izquierda es la aproximación numérica:pag

Prueba de Stokes 1

sin embargo, si es una constante, el problema de la prueba funciona bien: νPrueba de Stokes 2

Supongo que se debe a la forma en que estoy imponiendo la condición de compatibilidad, ya que está vinculada con la estabilidad de inf-sup de todo el sistema, aquí está mi segunda pregunta:

  • (P2): ¿hay alguna otra forma que no sea (2) para imponer la compatibilidad para la presión ? o al acuñar el problema de la prueba, ¿qué tipo de debo usar?pagpag
Shuhao Cao
fuente
MathML no funciona?
Shuhao Cao
Usamos MathJaX en StackExchange, todo lo que publicaste se muestra muy bien, gracias por la pregunta detallada.
Aron Ahmadia

Respuestas:

8

La condición de compatibilidad se refiere a la velocidad, no a la presión. Establece que si solo tiene condiciones de límite de Dirichlet para la velocidad, entonces estas deberían ser compatibles con la restricción libre de divergencia, es decir, conΩtunorte=0 0Ω el límite de la dominio computacional (no la celda).

En este caso, no se puede distinguir de con una constante arbitraria porque no tiene una condición límite en para fijar la constante. Por lo tanto, hay infinitas soluciones para la presión y, para comparar soluciones, se necesita una convención. Los matemáticos prefieren elegir modo que (porque pueden integrarse) mientras que el físico prefiere (porque pueden medir en un punto). Si es su equivalente discreto de tiene un espacio nulo que consiste en el vector de identidad.pag(pag+C)CpagCpag¯=pagrmiFpag(XrmiF)=pagrmiFsipagpag , implica quesi

Los métodos del subespacio de Krylov pueden resolver un sistema singular al eliminar el espacio nulo del subespacio de Krylov en el que buscan la solución. Sin embargo, eso no significa que obtendrá la solución que coincida con una convención dada, siempre tendrá que determinar la constante usted mismo en un paso posterior al procesamiento, ningún solucionador puede hacerlo por usted.pag

Aquí hay algunas sugerencias para abordar su problema:

  • La ecuación (2) parece extraña. Si es una función de ¿cómo puede estar fuera de la integral?νX
  • ¿Su campo de velocidad satisface la restricción de compatibilidad?
  • Trate de no hacer nada por la presión, solo deje que el solucionador tenga una , luego mire . ¿Es una constante?pagpag-pagmiXunCt
  • Si no, ¿estás seguro de que el espacio nulo de es de hecho el vector de identidad y nada más? ¿Tanto en papel como en el código? El problema parece lo suficientemente pequeño como para calcular realmente el espacio nulo.si
Chris
fuente
2

En cuanto a (Q1), puede elegir un solucionador de problemas de punto de silla que calcule una solución de mínimos cuadrados para su sistema. Entonces, se puede imponer una condición adicional al multiplicador, como establecer un grado específico de libertad, imponer un promedio específico.

En general, y creo que esto responde (Q1), puede usar una restricción lineal que pueda distinguir diferentes constantes.

Esta restricción puede imponerse en un paso posterior al procesamiento, o por una elección adecuada del espacio de prueba (por ejemplo, si deja de lado un grado de libertad).

shuhalo
fuente