計(jì)算流體力學(xué)常用的五大類(lèi)數(shù)值方法

 流體力學(xué)數(shù)值方法有很多種，其數(shù)學(xué)原理各不相同，但有二點(diǎn)是所有方法都具備的，即離散化和代數(shù)化?？偟膩?lái)說(shuō)其基本思想是：將原來(lái)連續(xù)的求解區(qū)域劃分成網(wǎng)格或單元子區(qū)域，在其中設(shè)置有限個(gè)離散點(diǎn)(稱為節(jié)點(diǎn))，將求解區(qū)域中的連續(xù)函數(shù)離散為這些節(jié)點(diǎn)上的函數(shù)值；通過(guò)某種數(shù)學(xué)原理，將作為控制方程的偏微分方程轉(zhuǎn)化為聯(lián)系節(jié)點(diǎn)上待求函數(shù)值之間關(guān)系的代數(shù)方程(離散方程)，求解所建立起來(lái)的代表方程以獲得求解函數(shù)的節(jié)點(diǎn)值。

 不同的數(shù)值方法，其主要區(qū)別在于求解區(qū)域的離散方式和控制方程的離散方式上。在流體力學(xué)數(shù)值方法中，應(yīng)用比較廣泛的是有限差分法、有限元法、邊界元法、有限體積法和有限分析法，現(xiàn)簡(jiǎn)述如下。

 深圳市有限元科技有限公司是一家有十年有限元分析項(xiàng)目經(jīng)驗(yàn)的高科技企業(yè)，公司代理國(guó)外多款著名有限元軟件，并提供軟件的銷(xiāo)售與培訓(xùn)業(yè)務(wù)，另研發(fā)出多款行業(yè)有限元通用軟件，并秉承以最高質(zhì)量的產(chǎn)品和最高質(zhì)量的服務(wù)滿足客戶的各種需求的服務(wù)理念，致力于為客戶提供一站式有限元整體解決方案，目前已為全國(guó)超過(guò)500家企業(yè)提供有限元分析服務(wù)。如需購(gòu)買(mǎi)有限元軟件或咨詢服務(wù)請(qǐng)聯(lián)系電話：13632683051，咨詢QQ：2039363860/4006046636。

 有限差分法

 這是最早采用的數(shù)值方法，它是將求解區(qū)域劃分為矩形或正交曲線網(wǎng)格，在網(wǎng)格線交點(diǎn)(即節(jié)點(diǎn))上，將控制方程中的每一個(gè)微商用差商來(lái)代替，從而將連續(xù)函數(shù)的微分方程離散為網(wǎng)格節(jié)點(diǎn)上定義的差分方程，每個(gè)方程中包含了本節(jié)點(diǎn)及其附近一些節(jié)點(diǎn)上的待求函數(shù)值，通過(guò)求解這些代數(shù)方程就可獲得所需的數(shù)值解。

 有限差分法的優(yōu)點(diǎn)是它建立在經(jīng)典的數(shù)學(xué)逼近理論的基礎(chǔ)上，容易為人們理解和接受；有限差分法的主要缺點(diǎn)是對(duì)于復(fù)雜流體區(qū)域的邊界形狀處理不方便，處理得不好將影響計(jì)算精度。

 有限元法

 有限元法的基本原理是把適定的微分問(wèn)題的解域進(jìn)行離散化，將其剖分成相連結(jié)又互不重疊的具有一定規(guī)則幾何形狀的有限個(gè)子區(qū)域(如：在二維問(wèn)題中可以劃分為三角形或四邊形；在三維問(wèn)題中可以劃分為四面體或六面體等)，這些子區(qū)域稱之為單元，單元之間以節(jié)點(diǎn)相聯(lián)結(jié)。函數(shù)值被定義在節(jié)點(diǎn)上，在單元中選擇基函數(shù)(又稱插值函數(shù))，以節(jié)點(diǎn)函數(shù)值與基函數(shù)的乘積的線性組合成單元的近似解來(lái)逼近單元中的真解。

 利用古典變分方法(里茲法或伽遼金法)由單元分析建立單元的有限元方程，然后組合成總體有限元方程，考慮邊界條件后進(jìn)而求解。由于單元的幾何形狀是規(guī)則的，因此在單元上構(gòu)造基函數(shù)可以遵循相同的法則，每個(gè)單元的有限元方程都具有相同的形式，可以用標(biāo)準(zhǔn)化的格式表示，其求解步驟也就變得很規(guī)范，即使是求解域剖分各單元的尺寸大小不一樣，其求解步驟也不用改變，這就為利用計(jì)算機(jī)編制通用程序進(jìn)行求解帶來(lái)了方便。

 有限元法的主要優(yōu)點(diǎn)是對(duì)于求解區(qū)域的單元剖分沒(méi)有特別的限制，因此特別適合處理具有復(fù)雜邊界流場(chǎng)的區(qū)域。

 邊界元法

 邊界元法是在經(jīng)典積分方程和有限元法基礎(chǔ)上發(fā)展起來(lái)的求解微分方程的數(shù)值方法，其基本思想是：將微分方程相應(yīng)的基本解作為權(quán)函數(shù)，應(yīng)用加權(quán)余量法并應(yīng)用格林函數(shù)導(dǎo)出聯(lián)系解域中待求函數(shù)值與邊界上的函數(shù)值與法向?qū)?shù)值之間關(guān)系的積分方程；令積分方程在邊界上成立，m.tyrp.net獲得邊界積分方程，該方程表述了函數(shù)值和法向?qū)?shù)值在邊界上的積分關(guān)系，而在這些邊界值中，一部份是在邊界條件中給定的，另一部份是待求的未知量，邊界元法就是以邊界積分方程作為求解的出發(fā)點(diǎn)，求出邊界上的未知量；在所導(dǎo)出的邊界積分方程基礎(chǔ)上利用有限元的離散化思想，把邊界離散化，建立邊界元代數(shù)方程組，求解后可獲得邊界上全部節(jié)點(diǎn)的函數(shù)值和法向?qū)?shù)值；將全部邊界值代入積分方程中，即可獲得內(nèi)點(diǎn)函數(shù)值的計(jì)算表達(dá)式，它可以表示成邊界節(jié)點(diǎn)值的線性組合。

 邊界元法的優(yōu)點(diǎn)是：

 （1）將全解域的計(jì)算化為解域邊界上的計(jì)算，使求解問(wèn)題的維數(shù)降低了一維，減少了計(jì)算工作量；

 （2）能夠方便地處理無(wú)界區(qū)域問(wèn)題。例如對(duì)于勢(shì)流等的無(wú)限區(qū)域問(wèn)題，使用邊界元法求解時(shí)由于基本解滿足無(wú)窮遠(yuǎn)處邊界條件，在無(wú)窮遠(yuǎn)處邊界上的積分恒等于零。因此對(duì)于無(wú)限區(qū)域問(wèn)題，例如具有無(wú)窮遠(yuǎn)邊界的勢(shì)流問(wèn)題，無(wú)需確定外邊界，只需在內(nèi)邊界上進(jìn)行離散即可；

 （3）邊界元法的精度一般高于有限元法。邊界元法的主要缺點(diǎn)是邊界元方程組的系數(shù)矩陣是不對(duì)稱的滿陣，該方法目前只適用于線性問(wèn)題。

 有限體積法

 有限體積法又稱為控制體積法，其導(dǎo)出離散方程的基本思路是：

 （1）將計(jì)算區(qū)域劃分為一系列不重復(fù)的控制體積，每一個(gè)控制體積都有一個(gè)節(jié)點(diǎn)作代表，將待求的守恒型微分方程在任一控制體積及一定時(shí)間間隔內(nèi)對(duì)空間與時(shí)間作積分；

 （2）對(duì)待求函數(shù)及其導(dǎo)數(shù)對(duì)時(shí)間及空間的變化型線或插值方式作出假設(shè)；

 （3）對(duì)步驟1中各項(xiàng)按選定的型線作出積分并整理成一組關(guān)于節(jié)點(diǎn)上未知量的離散方程。有限體積法著重從物理觀點(diǎn)來(lái)構(gòu)造離散方程，每一個(gè)離散方程都是有限大小體積上某種物理量守恒的表示式，推導(dǎo)過(guò)程物理概念清晰，離散方程系數(shù)具有一定的物理意義，并可保證離散方程具有守恒特性，這是有限體積法的主要優(yōu)點(diǎn)。

 就離散方法而言，有限體積法可視作有限元法和有限差分法的中間物，該方法的主要缺點(diǎn)是不便對(duì)離散方程進(jìn)行數(shù)學(xué)特性分析。

 有限分析法

 有限分析法在某種意義上說(shuō)是在有限元法基礎(chǔ)上發(fā)展起來(lái)的一種數(shù)值方法，其基本思想是：將求解區(qū)域劃分成矩形網(wǎng)格，網(wǎng)格線的交點(diǎn)為計(jì)算節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)與相鄰的四個(gè)網(wǎng)格組成一個(gè)計(jì)算單元，即一個(gè)計(jì)算單元由一個(gè)中心節(jié)點(diǎn)與8個(gè)相鄰節(jié)點(diǎn)組成；在每個(gè)單元中函數(shù)的近似解不是象有限元方法那樣采用單元基函數(shù)的線性組合來(lái)表達(dá)，而是以單元中未知函數(shù)的分析解來(lái)表達(dá)；為了獲得單元中的分析解，m.tyrp.net單元邊界條件采用插值函數(shù)來(lái)逼近，在單元中把控制方程中非線性項(xiàng)局部線性化(如N-S方程中的對(duì)流項(xiàng)中認(rèn)為其流速為已知，并對(duì)單元中待求函數(shù)的組合形式作出假設(shè)，找出其系數(shù)用單元邊界節(jié)點(diǎn)上待求函數(shù)值表達(dá)的分析解；利用單元分析解確定單元中心節(jié)點(diǎn)與8個(gè)相鄰節(jié)點(diǎn)間待求函數(shù)值之間關(guān)系的一個(gè)代數(shù)方程，稱為單元有限分析方程；將所有內(nèi)點(diǎn)上的單元有限分析方程聯(lián)立，就構(gòu)成總體有限分析方程，通過(guò)代數(shù)方程組求解，即可獲得求解區(qū)域中全部離散點(diǎn)的函數(shù)值。

 雖然有限分析解獲得的是求解區(qū)域中離散點(diǎn)的函數(shù)值，但是由于每個(gè)單元內(nèi)部都有與其中心節(jié)點(diǎn)對(duì)應(yīng)的分析解表達(dá)式，因此有限分析解在每一個(gè)節(jié)點(diǎn)的局部區(qū)域內(nèi)都是連續(xù)可微的，這對(duì)于需要計(jì)算求解函數(shù)導(dǎo)數(shù)的計(jì)算流體力學(xué)問(wèn)題具有明顯的優(yōu)勢(shì)。

 該計(jì)算方法與有限元、有限差分法比較具有較高的精度。此外，有限分析法具有自動(dòng)迎風(fēng)特性，能準(zhǔn)確地模擬對(duì)流項(xiàng)，同時(shí)不存在數(shù)值振蕩失真問(wèn)題。有限分析法的缺點(diǎn)是對(duì)復(fù)雜形狀的求解區(qū)域適應(yīng)性較差。

 本文出自深圳有限元科技有限公司官網(wǎng)：m.tyrp.net 轉(zhuǎn)載請(qǐng)注明