有限元素法與有限體積法

本質與自然邊界條件(essential vs natural boundary conditions)

解一條 PDE 時,你還必須說明區域邊緣會發生什麼事——邊界條件。在有限元素法中,邊界依照它們「如何」進入弱形式而有兩種非常不同的型態。有些你必須在函數空間上手動施加(本質);有些則自行打理,自動地烘進積分裡(自然)。這些名稱描述的是方法與它們各自的關係,而非物理本身。

「本質(Dirichlet)」條件指定解在部分邊界上的「值」,例如把一片膜沿著夾緊的邊固定在高度零:那裡 u = g。弱形式無法透過積分施加這一點,所以它必須被直接強加——你把試驗空間限制在那段邊界上已等於 g 的函數,並要求測試函數在那裡「為零」(值已知之處沒有可測試的東西)。「自然(Neumann)」條件指定「通量」或導數,例如已知的熱流或零牽引力:du/dn = h。這一個免費地進入弱形式——它恰好就是分部積分產生的那個邊界項。你只要把那個邊界項「留著」(而不是用為零的測試函數消掉它),條件就在弱的平均意義下自動被滿足。

這種劃分是有限元素法最優雅的特徵之一,也是初學者常見的困惑來源。要記住的誠實:自然條件只在「弱」意義下被滿足——是在積分意義下近似地、而非逐點地滿足——並且只在網格細化時才趨於精確;而本質條件則是精確地施加在離散解上。零通量(齊次 Neumann)邊界字面上完全不需要任何程式碼:什麼都不做就施加了它。而把兩者搞混——試圖手動「強加」一個 Neumann 值,或忘了約束一個 Dirichlet 節點——會導致奇異的方程組或錯誤的答案。

桿中的熱:把左端固定在 100 度(u(0) = 100)是本質的——你直接設定那個節點值,並把它的列從未知數中移除。把右端絕熱(du/dn = 0,無熱通量)是自然的——你字面上什麼都不做,弱形式中為零的邊界項就替你施加了它。

值的條件要由人手施加;通量的條件搭著邊界項免費進入。

自然(Neumann)條件只在弱意義下成立,並隨 h 縮小趨於精確;本質(Dirichlet)條件則是精確施加。齊次 Neumann 條件完全不需程式碼——但純 Neumann 問題會讓解在差一個常數的範圍內無法確定,產生奇異的剛度矩陣。

又稱
Dirichlet vs Neumann conditionsimposed vs natural BCs狄利克雷與諾伊曼邊界條件強制與自然邊界條件