ポアソン方程式における境界条件(例:周期的)の役割


8

与えられた3次元ポアソン方程式 と右辺とドメイン、Iは、機能上の任意の境界条件(BC)を課すことが解放午前φ、または実行彼らはどういうわけか右側と一致する必要がありますか?特に、定期的なBCを課す場合、右側のソリューションは1つだけでしょうか?

∇2ϕ(x,y,z)=f(x,y,z)
ϕ

たとえば、聞かせて とIボックスに解決する(0 、1 )× (0 、1 )× (0 、1 )。今、すべてのソリューションは、の合計でなければなりませんφ 0 + φ 1: φ

f(x,y,z)=−3π2sin⁡(πx)sin⁡(πy)sin⁡(πz)
(0,1)×(0,1)×(0,1)ϕ0+ϕ1
ϕ0(x,y,z)=sin⁡(πx)sin⁡(πy)sin⁡(πz),
なぜなら∇2ϕ0=f、およびϕ1任意であり、調和関数(すなわち、∇2ϕ1=0)。正しい?

ゼロのディリクレBCを課す場合、BCを満たし、方程式を満たし、解は一意でなければならないため、ϕ0が唯一の解です。正しい?

ϕ1ϕ0+ϕ1ϕ1


技術的には、この質問はmath.SEに適しています。
— David Ketcheson 2013

質問では明示的には触れませんでしたが、有限要素ベースの数値ソルバーのコンテキストで暗黙的に想定しています(そのため、私は答えに興味を持っています)が、質問自体は確かに一般的です。
— オンドレジ・セティク

回答:


10

数値的な観点からは、離散化について直接議論するのがおそらく最も簡単です。

Ax=bAbfAbf

fff

Kx≡CC

ARange⁡A=K⊥Ax=b1⋅b=01

ϕ∈C2(Ω)Ω

∫ΩΔϕdx=∫∂Ω∂ϕ∂ndS.
ϕ
∫ΩΔϕdx=0.
ϕΔϕ=f
∫Ωfdx=0.
1⋅b=0fb

ジョン、この素晴らしい答えを本当にありがとう!それがどのように機能するかについての明確な洞察を提供します。周期的な状態には落とし穴があるように感じましたが、なかなか指を向けられませんでした。再度、感謝します。
— オンドレジ・セティク

6

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

ステートメント:周期解の場合、∫ρdv= 0条件(上記のBenによって説明)はシステムを正味中立に設定します。これは、周期的なボックスでのfの平均値を考えるのに良い(物理的な)方法です。

次に、このステートメントをテストしてみましょう。周期的な境界条件の場合、電界の積分はボックスの表面上でゼロでなければなりません(ボックスの表面上のポイントのペアは常に積分で互いに相殺されます)。

次に、発散定理(ガウスの法則)により、ボックスは中立である必要があり、これは私たちが証明しようと試みたことです。

∮SEndA=1ε0Qinside=0

ここでは物理学について少し説明しますが、ここでの議論に素晴らしい補強を提供すると思います。


0

ポアソン方程式の製造ソリューションの生成方法には特に注意が必要です。ソース項の定義は、PDEの弱い形式だけでなく、PDEの強い形式も満たさなければならないためです。上記の派生は、基本的にPDEの弱い形式を使用しています。言い換えると、境界での解の勾配とドメインでのソース項積分との間に関係があります。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.