Alternativen zur von-Neumann-Stabilitätsanalyse für Finite-Differenzen-Methoden

13

Ich arbeite an der Lösung der gekoppelten eindimensionalen Poroelastizitätsgleichungen (Biot-Modell), gegeben als:

(λ+2μ)2ux2+px=0
t[γp+ux]κη[2px2]=q(x,t)
in der Domäne und mit den Randbedingungen: Ω=(0,1)

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

Ich habe diese Gleichungen mit einem zentrierten Finite-Differenzen-Schema diskretisiert:

&ggr;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

Ich arbeite derzeit an den Details der Konvergenz des Schemas, indem ich dessen Konsistenz und Stabilität analysiere. Der Konsistenzteil scheint mir ziemlich einfach zu sein, aber ich sehe bereits einige Schwierigkeiten mit der Stabilitätsanalyse. Zunächst gibt es zwei Variablen und zwei Gleichungen. Zweitens gibt es in der zweiten Gleichung auch einen gemischten raumzeitlichen abgeleiteten Term. Ich kenne die Stabilitätsanalyse von Neumann und kann sehen, dass es sehr schwierig sein wird, mit dieser Methode Stabilität herzustellen. Gibt es Alternativen zur von neumann-Analyse, die ich verwenden kann?

Paul
quelle
1
Wenn Sie sich bei der Analyse mit einem Gleichungssystem nicht wohl fühlen, differenzieren Sie einfach die erste Gleichung in Bezug auf und die zweite in Bezug auf x . Verwenden Sie dann die Gleichheit gemischter partieller Derivate, um u zu eliminieren . txu
David Ketcheson
@ DavidKetcheson: Interessant. Im Wesentlichen schlagen Sie vor, dass ich das System auf eine einzige Variable reduzieren und die Standard-von-Neumann-Analyse auf ohne die Allgemeingültigkeit für u ? pu
Paul
Es ist das gleiche Problem, ob Sie es als System oder als skalare PDE schreiben.
David Ketcheson

Antworten:

7

Wenn Sie, zumindest für Ihre Analyse, ersetzen, vonuxkönnen Sie Ihr System als[ 0 0 I I ]dschreiben uxux wo alle Konstanten auf 1 gesetzt sind und wo der Index

[00ichich]ddt[ph(t)ux,h(t)]+[-hh-Δh0][ph(t)ux,h(t)]=[qh(t)0]()
1 bezieht sich auf die Raumdiskretisierung sowohl der Variablen als auch der Differentialoperatoren. Ihr Schema wird dann erhalten, indem man d approximierth über implizite Euler.ddt

Nun ist die differentiell-algebraische (DAE) Struktur offensichtlich. Für die Variablen gibt es sowohl Differentialgleichungen (in der Zeit) als auch algebraische Gleichungen.

Wenn Sie zeigen können, dass invertierbar ist, vgl. dieser Vorabdruck [S. 3] und die folgende Bearbeitung, als der DAE Index 1 oder Fremdheit-frei ist und implizite Euler bekannt ist, konvergent zu sein, siehe Theorem 5.12 in diesem Buch . (Haftungsausschluss: Dieses Buch ist nicht frei verfügbar und wird von meinem Doktorvater verfasst.)[-hhichich]

Mit diesem Ansatz können Sie die Stabilitätsanalyse umgehen.

L2()Δhh

()uux

ANHANG: Eine DAE wird als Index 1 bezeichnet, wenn sie ohne Differenzierung der Gleichungen in eine ODE umgewandelt werden kann.

[E10]y˙+[EIN1EIN2]y=f.
[E1EIN2]y~y[E1EIN2][E~11E~12EIN~21EIN~22]EIN~22EIN2EIN~11-E~12EIN~22-1EIN~21

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

Jan
quelle
[-hhichich]
@ Paul Ich habe keinen Satz Referenz finden, so dass ich die Argumente in meine Antwort einfügen wird ...
Jan
4

Die hier angegebenen Gleichungen sind mir nicht bekannt, aber ich erinnere mich, dass ich in meinem Kurs eine andere Methode zur Überprüfung der Stabilität eines numerischen Schemas gelernt habe. Es ist als Modified Equation Analysis bekannt.

Hier ist eine gute Referenz dafür,

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

In der obigen Literaturstelle wird der Zusammenhang zwischen der auf der Analyse der modifizierten Gleichung basierenden Stabilitätstheorie und der Stabilitätsanalyse nach Von Neumann hergestellt.

Nach einer kurzen Online-Suche stieß ich auf folgende Referenzen:

In diesem Artikel wird die Finite-Differenz-Modellierung von Biots poroelastischen Gleichungen bei seismischen Frequenzen erörtert. Es gibt auch einen Abschnitt zur Stabilität des numerischen Schemas.

Dieses Papier stellt eine Lösungsstrategie des gekoppelten Systems zu entkoppeln, und die Stabilität des numerischen Verfahrens zu überprüfen.

Subodh
quelle
Ich habe die modifizierte Gleichungsanalyse nicht für die obigen Gleichungen durchgeführt, aber als die Frage nach Alternativen zur Von-Neumann-Analyse gestellt wurde, schrieb ich die obige Antwort. Es ist durchaus möglich, dass es die Frage nicht beantwortet. Aber jemand könnte die aufgelisteten Referenzen für seine Arbeit nützlich finden.
Subodh
Vielen Dank für den Hinweis! Ich kann sehen, dass das Formular, das in Ihrem Artikel zur Analyse modifizierter Gleichungen benötigt wird, nicht ganz zu den von mir verwendeten Gleichungen passt, aber es ist faszinierend, neue Analysetechniken zu erlernen!
Paul