Alternativas al análisis de estabilidad von Neumann para métodos de diferencia finita

13

Estoy trabajando en resolver las ecuaciones de poroelasticidad unidimensionales acopladas (modelo de biot), dadas como:

(λ+2μ)2ux2+px=0
t[γp+ux]κη[2px2]=q(x,t)
en el dominio y con las condiciones de contorno: Ω=(0,1)

p=0,(λ+2μ)ux=u0 en y en .u = 0 , px=0x=1u=0,px=0x=1

Discreté estas ecuaciones usando un esquema de diferencia finita centrada:

γp t + 1 i -p t i

(λ+2μ)ui+1t+12uit+1+ui1t+1Δx2+pi+1t+1pi1t+12Δx=0
γpit+1pitΔt+ui+1t+1ui1t+12ΔxΔt[ui+1tui1t2ΔxΔt]κη[pi+1t+12pit+1+pi1t+1Δx2]=qit+1

Actualmente estoy trabajando en los detalles de la convergencia del esquema analizando su consistencia y estabilidad. La parte de consistencia me parece bastante sencilla, pero ya preveo algunas dificultades con el análisis de estabilidad. En primer lugar, hay dos variables y dos ecuaciones. En segundo lugar, también hay un término derivado espaciotemporal mixto en la segunda ecuación. Estoy familiarizado con el análisis de estabilidad de von neumann y puedo ver que será muy difícil establecer la estabilidad con este método. ¿Hay alguna alternativa al análisis de von neumann que pueda usar?

Paul
fuente
1
Si no se siente cómodo haciendo el análisis con un sistema de ecuaciones, simplemente diferencie la primera ecuación con respecto a y la segunda con respecto a x . Luego usa la igualdad de las derivadas parciales mixtas para eliminar u . txu
David Ketcheson
@DavidKetcheson: Interesante. En esencia, ¿está sugiriendo que podría reducir el sistema a una sola variable y realizar el análisis estándar de von Neumann en sin ninguna pérdida de generalidad para usted ? pu
Paul
Es el mismo problema, ya sea que lo escriba como un sistema o un PDE escalar.
David Ketcheson

Respuestas:

7

Si sustituye, al menos para su análisis, porux, puede escribir su sistema como [ 0 0 I I ]duxux donde todas las constantes se establecen en 1 y donde el subíndice

[00II]ddt[ph(t)ux,h(t)]+[hhΔh0][ph(t)ux,h(t)]=[qh(t)0]()
1 refiere a la discretización del espacio tanto de las variables como de los operadores diferenciales. Su esquema se obtiene aproximando dh través de Euler implícito.ddt

Ahora la estructura diferencial-algebraica (DAE) es evidente. Para las variables hay ecuaciones diferenciales (en el tiempo) y algebraicas.

Si puede demostrar que es invertible, cf. esta preimpresión [p. 3] y la edición a continuación, que el DAE es de índice 1 o libre de extrañeza y se sabe que Euler es convergente, ver Teorema 5.12 en este libro . (Descargo de responsabilidad: este libro no está disponible gratuitamente y está escrito por mi supervisor de doctorado)[hhII]

Con este enfoque, tal vez evite el análisis de estabilidad.

L2()Δhh

()uux

APÉNDICE: Se dice que un DAE es el índice 1, si se puede transformar en un ODE sin diferenciar las ecuaciones.

[E10]y˙+[A1A2]y=f.
[E1A2]y~y[E1A2][E~11E~12A~21A~22]A~22A2A~11E~12A~221A~21

()A2:=[h h]y~2(ph,ux,h)ddty~2()

ene
fuente
[hhII]
@Paul No encontré un teorema de referencia, por lo que insertaré los argumentos en mi respuesta ...
enero
4

No estoy familiarizado con las ecuaciones dadas aquí, pero recuerdo haber aprendido otro método para verificar la estabilidad de un esquema numérico en mi curso. Se conoce como análisis de ecuaciones modificadas.

Aquí hay una buena referencia para eso,

http://193.146.160.29/gtb/sod/usu/$UBUG/repositorio/10291890_Warming.pdf

En la referencia anterior, se establece la conexión entre la teoría de la estabilidad basada en el análisis de ecuaciones modificadas y el análisis de estabilidad de Von Neumann.

Después de un poco de búsqueda en línea, me encontré con las siguientes referencias,

Este artículo analiza el modelo de diferencias finitas de las ecuaciones poroelásticas de Biot a frecuencias sísmicas. Tiene una sección sobre la estabilidad del esquema numérico también.

Este artículo presenta una estrategia de solución para desacoplar el sistema acoplado y verificar la estabilidad del esquema numérico.

Subodh
fuente
No he realizado el análisis de ecuaciones modificadas en las ecuaciones anteriores, pero como la pregunta pedía alternativas al análisis de Von Neumann, escribí la respuesta anterior. Es muy posible que no responda la pregunta. Pero alguien podría encontrar útiles las referencias enumeradas en su trabajo.
Subodh
¡Gracias por la referencia! Puedo ver que la forma necesaria en su documento de Análisis de ecuaciones modificado no se ajusta a las ecuaciones que estoy usando, ¡pero es bastante interesante aprender nuevas técnicas de análisis!
Paul