Rolle der Randbedingungen (z. B. periodisch) in der Poisson-Gleichung

8

Wenn die 3D-Poisson-Gleichung

2ϕ(x,y,z)=f(x,y,z)
und die rechte Seite und die Domäne gegeben sind, kann ich der Funktion Randbedingungen (BC) auferlegen ϕoder tun müssen sie irgendwie mit der rechten Seite übereinstimmen? Insbesondere wenn ich periodische BC auferlege, gibt es genau eine Lösung für eine rechte Seite?

Zum Beispiel sei:

f(x,y,z)=3π2sin(πx)sin(πy)sin(πz)
und ich löse auf einer Box (0,1)×(0,1)×(0,1) . Jetzt muss jede Lösung eine Summe von ϕ0+ϕ1 wobei:
ϕ0(x,y,z)=sin(πx)sin(πy)sin(πz),
weil2ϕ0=f ist undϕ1 eine beliebigeharmonische Funktion ist(dh2ϕ1=0 ). Richtig?

Wenn ich Dirichlet BC Null auferlege, ist ϕ0 die einzige Lösung, da es das BC erfüllt, die Gleichung erfüllt und die Lösung eindeutig sein muss. Richtig?

ϕ1ϕ0+ϕ1ϕ1

Ondřej Čertík
quelle
Technisch ist diese Frage eher für math.SE geeignet.
David Ketcheson
Ich habe es in der Frage nicht explizit erwähnt, aber ich gehe implizit im Kontext eines auf finiten Elementen basierenden numerischen Lösers davon aus (deshalb interessiert mich die Antwort), aber die Frage selbst ist in der Tat allgemein.
Ondřej Čertík

Antworten:

10

Aus numerischer Sicht ist es vielleicht am einfachsten, die Diskretisierungen direkt zu diskutieren.

Ax=bAbfAbf

fff

KxCC

ARangeA=KAx=b1b=01

ϕC2(Ω)Ω

ΩΔϕdx=ΩϕndS.
ϕ
ΩΔϕdx=0.
ϕΔϕ=f
Ωfdx=0.
1b=0fb
Ben
quelle
John, vielen Dank für diese erstaunliche Antwort! Es bietet einen klaren Einblick in die Funktionsweise. Ich hatte das Gefühl, dass es einen Haken für den periodischen Zustand gibt, konnte aber keinen Finger darauf zeigen. Danke noch einmal.
Ondřej Čertík
6

2ϕ=f
f=(ρ/εo)

Aussage: Für die periodische Lösung setzt die Bedingung ∫ρdv = 0 (oben von Ben beschrieben) das System auf netzneutral. Dies ist eine gute (physikalische) Möglichkeit, über den Durchschnittswert von f über die periodische Box nachzudenken.

Testen wir nun diese Aussage. Für periodische Randbedingungen muss das Integral des elektrischen Feldes über der Oberfläche der Box Null sein (Sie können immer Punktepaare auf der Boxoberfläche finden, die sich im Integral gegenseitig aufheben).

Dann sollte nach dem Divergenzsatz (Gauß-Gesetz) die Box neutral sein, und das wollen wir beweisen.

SEndA=1ε0Qinside=0

Ein bisschen Physik hier, aber ich denke, es ist eine schöne Verstärkung für die Diskussion hier.

Ajay Muralidharan
quelle
0

Seien Sie vorsichtig, wenn Sie ein Verfahren zur Herstellung von Lösungen für die Poisson-Gleichung erstellen. Da die Definition des Quellterms sowohl die starke Form der PDE als auch die schwache Form der PDE erfüllen muss. Die oben angegebene Ableitung verwendet grundsätzlich die schwache Form der PDE. Mit anderen Worten, es gibt eine Beziehung zwischen dem Gradienten der Lösung an der Grenze und dem Quellterm Integral in der Domäne.

Hasbestein
quelle