tft每日頭條

 > 生活

 > 微分方程數值解知識點總結

微分方程數值解知識點總結

生活 更新时间:2025-01-31 23:01:35

微分方程數值解知識點總結(微分方程數值計算中的有限元法)1

一.有限差分方法(Finite Differential Method)

有限差分方法是計算機數值模拟最早采用的方法,至今仍被廣泛運用。該方法将求解域劃分為差分網格,用有限個網格節點代替連續的求解域。有限差分法以泰勒級數展開等方法,把控制方程中的導數用網格節點上的函數值的差商代替進行離散,從而建立以網格節點上的值為未知數的代數方程組。該方法是一種直接将微分問題變為代數問題的近似數值解法,數學概念直觀,表達簡單,是發展較早且比較成熟的數值方法。

對于有限差分格式,從格式的精度來劃分,有一階格式、二階格式和高階格式。從差分的空間形式來考慮,可分為中心格式和逆風格式。考慮時間因子的影響,差分格式還可以分為顯格式、隐格式、顯隐交替格式等。目前常見的差分格式,主要是上述幾種形式的組合,不同的組合構成不同的差分格式。差分方法主要适用于有結構網格,網格的步長一般根據實際地形的情況和柯朗穩定條件來決定。

構造差分的方法有多種形式,目前主要采用的是泰勒級數展開方法。其基本的差分表達 式主要有三種形式:一階向前差分、一階向後差分、一階中心差分和二階中心差分等, 其中前兩種格式為一階計算精度,後兩種格式為二階計算精度。通過對時間和空間這幾種不同差分格式的組合,可以組合成不同的差分計算格式。

二.有限元法(Finite Element Method)

有限元法的基礎是變分原理和加權餘量法,其基本求解思想是把計算域劃分為有限個互不重疊的單元,在每個單元内,選擇一些合适的節點作為求解函數的插值點,将微分 方程中的變量改寫成由各變量或其導數的節點值與所選用的插值函數組成的線性表達式 ,借助于變分原理或加權餘量法,将微分方程離散求解。采用不同的權函數和插值函數形式,便構成不同的有限元方法。有限元方法最早應用于結構力學,後來随着計算機的發展慢慢用于流體力學的數值模拟。在有限元方法中,把計算域離散剖分為有限個互不重疊且相互連接的單元,在每個單元内選擇基函數,用單元基函數的線形組合來逼近單元中的真解,整個計算域上總體的基函數可以看為由每個單元基函數組成的,則整個計算域内的解可以看作是由所有單元上的近似解構成。

根據所采用的權函數和插值函數的不同,有限元方法也分為多種計算格式。從權函數的選擇來說,有配置法、矩量法、最小二乘法和伽遼金法。從計算單元網格的形狀來劃分,有三角形網格、四邊形網格和多邊形網格,從插值函數的精度來劃分,又分為線性插值函數和高次插值函數等。不同的組合同樣構成不同的有限元計算格式。對于權函數,伽遼金(Galerkin)法是将權函數取為逼近函數中的基函數 ;最小二乘法是令權函數等于餘量本身,而内積的極小值則為對代求系數的平方誤差最小;在配置法中,先在計算域内選取N個配置點 。令近似解在選定的N個配置點上嚴格滿足微分方程,即在配置點上令方程餘量為0。插值函數一般由不同次幂的多項式組成,但也有采用三角函數或指數函數組成的乘積表示,但最常用的多項式插值函數。有限元插值函數分為兩大類,一類隻要求插值多項式本身在插值點取已知值,稱為拉格朗日(Lagrange)多項式插值;另一種不僅要求插值多項式本身,還要求它的導數值在插值點取已知值,稱為哈密特(Hermite)多項式插值。單元坐标有笛卡爾直角坐标系和無因次自然坐标,有對稱和不對稱等。常采用的無因次坐标是一種局部坐标系,它的定義取決于單元的幾何形狀,一維看作長度比,二維看作面積比,三維看作體積比。在二維有限元中,三角形單元應用的最早,近來四邊形等參元的應用也越來越廣。對于二維三角形和四邊形電源單元,常采用的插值函數為有Lagrange插值直角坐标系中的線性插值函數及二階或更高階插值函數、面積坐标系中的線性插值函數、二階或更高階插值函數等。

三.有限體積法(Finite Volume Method)

有限體積法又稱為控制體積法。其基本思路是:将計算區域劃分為一系列不重複的控制體積,并使每個網格點周圍有一個控制體積;将待解的微分方程對每一個控制體積積分,便得出一組離散方程。其中的未知數是網格點上的因變量的數值。為了求出控制體積的積分,必須假定值在網格點之間的變化規律,即假設值的分段的分布的分布剖面。從積分區域的選取方法看來,有限體積法屬于加權剩餘法中的子區域法;從未知解的近似方法看來,有限體積法屬于采用局部近似的離散方法。簡言之,子區域法屬于有限體積發的基本方法。

有限體積法的基本思路易于理解,并能得出直接的物理解釋。離散方程的物理意義,就是因變量在有限大小的控制體積中的守恒原理,如同微分方程表示因變量在無限小的控制體積中的守恒原理一樣。 有限體積法得出的離散方程,要求因變量的積分守恒對任意一組控制體積都得到滿足,對整個計算區域,自然也得到滿足。這是有限體積法吸引人的優點。有一些離散方法,例如有限差分法,僅當網格極其細密時,離散方程才滿足積分守恒;而有限體積法即使在粗網格情況下,也顯示出準确的積分守恒。

就離散方法而言,有限體積法可視作有限單元法和有限差分法的中間物。有限單元法必須假定值在網格點之間的變化規律(既插值函數),并将其作為近似解。有限差分法隻考慮網格點上的數值而不考慮值在網格點之間如何變化。有限體積法隻尋求的結點值,這與有限差分法相類似;但有限體積法在尋求控制體積的積分時,必須假定值在網格點之間的分布,這又與有限單元法相類似。在有限體積法中,插值函數隻用于計算控制體積的積分,得出離散方程之後,便可忘掉插值函數;如果需要的話,可以對微分方程中不同的項采取不同的插值函數。

,

更多精彩资讯请关注tft每日頭條,我们将持续为您更新最新资讯!

查看全部

相关生活资讯推荐

热门生活资讯推荐

网友关注

Copyright 2023-2025 - www.tftnews.com All Rights Reserved