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

Wenn die 3D-Poisson-Gleichung

$$\nabla^2 \phi(x, y, z) = f(x, y, z)$$ und die rechte Seite und die Domäne gegeben sind, kann ich der Funktion Randbedingungen (BC) auferlegen $\phi$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 \pi^{2} \sin{\left (\pi x \right )} \sin{\left (\pi y \right )} \sin{\left (\pi z \right )}$$ und ich löse auf einer Box $(0, 1)\times (0, 1) \times (0, 1)$ . Jetzt muss jede Lösung eine Summe von $\phi_0+\phi_1$ wobei: $$\phi_0(x, y, z) = \sin{\left (\pi x \right )} \sin{\left (\pi y \right )} \sin{\left (\pi z \right )}\,,$$ weil$\nabla^2\phi_0 = f$ ist und$\phi_1$ eine beliebigeharmonische Funktion ist(dh$\nabla^2\phi_1 = 0$ ). Richtig?

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

$\phi_1$$\phi_0+\phi_1$$\phi_1$

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

$Ax = b$$A$$b$$f$$A$$b$$f$

$f$$f$$f$

$K$$x \equiv C$$C$

$A$$\operatorname{Range} A = K^\perp$$Ax = b$$\mathbf{1} \cdot b = 0$$\mathbf{1}$

$\phi \in C^2(\Omega)$$\Omega$

$$\int_\Omega \Delta \phi \, dx = \int_{\partial \Omega} \frac{\partial \phi}{\partial \mathbf{n}}\, dS.$$ $\phi$ $$\int_\Omega \Delta\phi\, dx = 0.$$ $\phi$$\Delta \phi = f$ $$\int_\Omega f\, dx = 0.$$ $\mathbf{1} \cdot b = 0$$f$$b$

$$\nabla^2\phi =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.

$$\oint_S {E_n dA = \frac{1}{{\varepsilon _0 }}} Q_{inside} =0$$

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

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.