偏微分方程邊值問題差分方法

偏微分方程邊值問題差分方法

物理學中的平衡態或定常態問題,例如彈性膜的平衡、彈性柱的扭轉、定常態熱傳導、電場、磁場、滲流、亞聲速流及不可壓縮無旋流等等,通常都可歸結為橢圓型偏微分方程邊值問題。

基本介紹

  • 中文名:偏微分方程邊值問題差分方法
  • 特點:直觀、簡易、普遍
  • 典型例子:泊松方程的邊值問題
  • 屬於:學術
正文,

正文

典型例子是泊松方程
(1)
的邊值問題。即要求定出未知函式=(,),使之在某個區域內滿足(1),並在區域邊界嬠上滿足一定的邊界條件,通常有如下三類:
第一類:=,
第二類:  第三類:  式中表示外法嚮導數,、和是定義在嬠上的已知函式。附加第一類和第二類邊界條件的問題分別稱為狄利克雷問題和馮·諾伊曼問題;有些問題在邊界不同區段滿足不同類的邊界條件,稱為混合邊值問題。
橢圓型邊值問題的求解,只在很特殊情況下才能用解析方法,一般情況下實際有效的途徑是數值方法,差分法是其中一類。
差分法的思想和做法是,把定解區域剖分為格線,在格線結點上以差商代替微商或用某種插值方式,把微分方程化為包含有限個未知數的差分方程組。差分法直觀、簡易、能普遍用於各種類型的微分方程和任意形狀的區域。因為它包含巨大的運算量,所以只在電子計算機問世之後,才得到廣泛的套用和發展。
從微分方程出發的差分化 格線剖分的一種最簡單又常用的做法是取平行於坐標軸的直線作為格線線,例如取=,=,、為步長,、取一切整數,這時格線結點為(,)。對方程(1)進行差分化、以表示差分近似解、表示在格線結點(,)上的分量。如果(,)是內結點,即鄰近四個格線結點都在上,則用中心二階差商代替二階微商代入(1),即得相應的差分方程
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
(2)
對第一類邊值問題,如果邊界結點正好落在嬠上,則取相應的邊界值即可。一般的情況如圖 1所示,嬠與格線線交於和,於是靠近邊界的結點可利用偏心差商建立差分方程:
式中()、()取邊界值,,是線段、的長度。
對第二、第三類邊值問題,可取外最靠近嬠的一層格線結點為邊界結點,相應的差分方程可建立如下:設邊界結點及其鄰近的邊界如圖2所示,過結點作嬠的法線,它與嬠和格線線分別交於和,用差商代替的外法嚮導數,其中()利用格線結點和上的值作線性插值,即
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
代入第二類、第三類邊界條件,就得到邊界結點的差分方程。
還可以用其他插值方法作邊界處理。但是,這種對微分方程及其邊界條件分開處理的方法,對自共軛邊值問題,包括現在討論的最簡單的典型例子,所得差分方程組的係數矩陣一般都不具有對稱性。
從積分守恆原理出發的差分化 與平衡態或定常態緊密聯繫的橢圓型邊值問題,在物理上表示某種守恆規律,在數學上表現為某種積分守恆形式。例如與方程(1)等價的積分守恆形式為
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
(3)
式中是區域的任一子區域。假設定解區域為多邊形,對作任意三角剖分,然後過三角形的邊作中垂線,如圖3中的虛線所示。對應每個格線結點,都存在一個以中垂線為邊的多邊形。若以每個這樣的多邊形作為(3)中的積分區域,建立格線結點的差分方程。若是內結點,對應的如圖4所示,則(3)中左端的環路積分可作如下逼近:
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
(4)
若是邊界結點,對應的如圖5所示,其中線段與落在邊界嬠上,這時有
  
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
(5)
對於(5)中的可用第二、第三類邊界條件代入。對於第一類邊界條件,須作某種插值處理。
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法對(3)的右端,作逼近
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
(6)
上述各式中的|·|表示線段的長度或區域的面積。把(4)、(5)、(6)的逼近公式代入(3),就得到差分方程組。對第二、第三類邊值問題,這種差分化途徑在處理上統一,所得差分方程組的係數矩陣具有對稱性。
從變分原理出發的差分化 平衡態或定常態的物理問題,往往可用變分原理表達,即表示為一個極小值問題。例如微分方程(1)的第一邊值問題,就等價於泛函
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
(7)
在滿足第一類邊界條件的容許函式集上的極小解問題。對於第二、第三類邊值問題,對應的泛函為
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
,  (8)
而容許函式不必滿足任何邊界條件,第二、第三類邊界條件是極小解自然滿足的。因此,也稱之為自然邊界或自由邊界問題。
從變分原理出發進行差分化,其步驟是先對區域作格線剖分,然後對積分(7)或(8)進行差分逼近,得到一個有限和式,它是定義在格線結點上的差分解的二次函式,它的極小解可歸結為解線性方程組,此方程組的係數矩陣恆具有對稱性。
差分方程組的求解 隨著差分法的實際套用,產生了在計算機上求解高階稀疏矩陣問題的種種方法,其中最簡單而且常用的是點鬆弛法。對代數方程組
點鬆弛法的一般疊代格式是
是疊代序號,是任意給定的初始值,是疊代參數。當係數矩陣對稱正定時,取0<1時,稱超鬆弛。對於五點格式的差分方程組(2),存在一個最優的疊代參數>1, 與 ω =1相比, 達到同樣精度的運算量從()階降為()階,是格線點的總數。
差分方程組的求解,還有各種直接法和其他疊代法。直接法大多是高斯消去法的變形,其中心問題是如何採取適當的消去順序,使得在不影響解的精度的前提下,儘可能在運算量、存貯量及程式複雜性等方面得到好處或達到某種平衡。在疊代法方面,則還有切比雪夫疊代和共軛斜量法,它們也常作為加速手段與點鬆弛法結合使用。對於特殊形狀區域(如矩形域),則有高效的快速傅立葉變換方法和交替方向法。特別引人注目的是近年發展起來的多重格線法,其運量可達到()階。
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法
偏微分方程邊值問題差分方法偏微分方程邊值問題差分方法

相關詞條

熱門詞條

聯絡我們