偏微分方程數值方法是數值分析的一個分支,研究如何得到偏微分方程(PDE) 的數值解。 [1] [2]

一般來說,對於雙曲型方程[3]拋物型方程[4]橢圓方程[5]都有專門的數值方法。 [6] [7]

方法概述

有限差分法

  在這種方法中,函數由它們在某些網格點處的值表示,並通過這些值的差分來近似導數。

有限元法

有限元法 (FEM)是一種數值技術,用於尋找微分方程邊值問題的近似解。它使用變分法,以最小化誤差函數和得到穩定的解。類似於連接許多短線可以逼近一個更大的圓的想法,FEM 是指通過在許多小的子區域(稱為有限元)上構建許多簡單的有限元方程以在更大的區域上逼近更複雜方程的一類方法。

將小區域上的方程聯立後,有限元法一般會得到一個大的代數方程組。

有限體積法

有限體積法是一種以代數方程的形式表示和計算偏微分方程的方法 [LeVeque, 2002;托羅,1999]。類似於有限差分法有限元法,取值是在網格上的離散位置計算的。 「有限體積」是指網格上每個節點周圍的小體積,函數在這個節點的離散值被視為函數在這個小區域內的平均值。有限體積方法中,使用散度定理包含散度項的偏微分方程中的體積積分轉換為曲面積分。然後將這些項作為每個有限體積表面處的通量進行求值。因為進入給定體積的通量與離開相鄰體積的通量相同,所以這些方法是守恆的。有限體積法的另一個優點是它易於定製以使用於非結構化網格。該方法用於許多計算流體力學軟件包。

譜方法

譜方法英語Spectral method是在應用數學科學計算中用於對某些微分方程進行數值求解的技術,通常會涉及到使用快速傅立葉變換。這個想法是將微分方程的解寫為某些「基函數」的和(例如,在譜方法中常用的傅立葉級數,它是三角函數的和),然後選擇出最符合微分方程的係數。

譜方法和有限元方法密切相關並且建立在相同的思想之上;它們之間的主要區別在於譜方法使用在整個域上非零的基函數,而有限元方法使用僅在小子域上非零的基函數。換句話說,譜方法採用全局逼近,而有限元方法使用局部逼近。部分出於這個原因,譜方法具有出色的誤差特性,當解是光滑的時,所謂的「指數收斂」是最快的。然而,沒有已知的三維單區域的譜方法激波捕獲結果。 [8]在有限元社區,有限元的階非常高,或隨當網格參數 h 減小趨於零而增加的方法有時稱為譜元法。

無網格方法

無網格方法英語Meshfree methods不需要連接模擬域數據點的網格。無網格方法可以模擬一些其他困難類型的問題,但需要額外的計算時間和編程工作。

區域分解方法

域分解方法英語Domain_decomposition_methods通過將邊界值問題拆分為子域上的較小邊界值問題并迭代以協調相鄰子域之間的解決方案來解決邊界值問題。每個子域一個或幾個未知數的粗略問題用於進一步協調全局子域之間的解決方案。子域上的問題是獨立的,這使得域分解方法適用於並行計算。域分解方法通常用作Krylov 空間迭代方法的預處理子,例如共軛梯度方法GMRES

多重網格方法

數值分析中的多重網格方法英語Multigrid method(MG)是一組使用不同層次離散化求解微分方程的算法它們是稱為多分辨率方法的一類技術的示例,在(但不限於)表現出多種行為尺度的問題中非常有用。例如,許多基本鬆弛方法對短波長和長波長分量表現出不同的收斂速度,這表明對這些不同尺度的處理方式不同,就像在多重網格的傅立葉分析方法中一樣。 [9] MG 方法可以用作求解器和預條件子。

比較

有限差分法通常被認為是最容易學習和使用的方法。有限元法和有限體積法廣泛應用於工程計算流體動力學中,非常適用於複雜幾何體中的問題。如果解足夠光滑,譜方法通常是最精確的。

參見

參考文獻

外部連結

Wikiwand in your browser!

Seamless Wikipedia browsing. On steroids.

Every time you click a link to Wikipedia, Wiktionary or Wikiquote in your browser's search results, it will show the modern Wikiwand interface.

Wikiwand extension is a five stars, simple, with minimum permission required to keep your browsing private, safe and transparent.