簡述數(shù)值模擬中離散化的方法

董磊磊，崔之健，孫明龍

（1.西安石油大學(xué)，陜西 西安 710065；2.中國石油集團西部鉆探工程有限公司試油公司，新疆 克拉瑪依 834000）

伴隨著全球經(jīng)濟和現(xiàn)代化管理工業(yè)的飛速發(fā)展，相應(yīng)而來的是計算機技術(shù)在各行各業(yè)的普及。這大大減少了人力的消耗，且使有些原本比較復(fù)雜的問題求解起來變得容易。計算機技術(shù)的相對成熟，再結(jié)合一些專業(yè)領(lǐng)域的知識，一些專業(yè)人士開發(fā)出一些軟件，模擬工程領(lǐng)域和自然領(lǐng)域的各種場景，用來解決一些實際問題，其處理問題的過程叫做數(shù)值模擬。經(jīng)過幾十年的發(fā)展，數(shù)值模擬已經(jīng)可以較為真實的模擬現(xiàn)場情況。通過模擬結(jié)果進而調(diào)整各個設(shè)計參數(shù)，以使設(shè)備在最優(yōu)設(shè)計結(jié)構(gòu)下達到最佳的工況。數(shù)值模擬把數(shù)學(xué)和計算機聯(lián)系了起來，在許多工程領(lǐng)域扮演著不可或缺的角色。在數(shù)值模擬過程中，遇到的都是一些連續(xù)域內(nèi)的非線性偏微分的問題，直接求解十分困難，而離散化的出現(xiàn)使這個問題簡單化，變得易于求解。

1 數(shù)值模擬

在早期的工業(yè)領(lǐng)域，許多問題因為其實際情況的復(fù)雜性，得到其解析解非常困難。1953年，Bruce G.H和PeacemanD.W對一維氣相不穩(wěn)定徑向和線形流進行了模擬，這是數(shù)值模擬技術(shù)初次面世。起初數(shù)值模擬的解法相對較少，相關(guān)的計算機技術(shù)不夠成熟，所以只能應(yīng)用于模擬一維問題。1954年，兩相流動模型應(yīng)運而生，West W.J和Garvin W.W針對油藏中存在的問題，模擬了不穩(wěn)定兩相流的流動情況［1］。

數(shù)值模擬，顧名思義就是對現(xiàn)實中某些工況數(shù)值的模擬，那么就必須借助電子計算機來進行［2］。既然是模擬，就會要求提供一些參數(shù)以作為原始數(shù)據(jù)，并且要根據(jù)特定問題選用適當(dāng)?shù)哪M方法。通過對有限元和有限容積有關(guān)知識的學(xué)習(xí)，熟悉了數(shù)值模擬相關(guān)理論后，接著進行大量數(shù)值模擬計算，最后的結(jié)果再以圖像顯示的方法展示出來，可以直觀的解釋某些特定問題［3］。這對學(xué)者了解發(fā)生某種現(xiàn)象的內(nèi)部原因有很大的幫助。數(shù)值模擬在研究生涯里使用非常廣泛，同時因為它需要用到計算機，因此可以叫做計算機模擬。

數(shù)值模擬其實就是用計算機做試驗，并且可以較為真實的展現(xiàn)特定介質(zhì)的工況。例如，某一特定尺寸形狀的機翼在空氣中的繞流情況，建立模型后，通過在計算機上模擬計算，最后將結(jié)果以圖片的形式展示出來，直觀明了，可以看到流場中的許多細節(jié)。計算機模擬結(jié)果與實驗結(jié)果沒有什么較大區(qū)別，而且簡化了許多步驟，耗時較短，有利于模擬較多的不同條件下設(shè)備的運轉(zhuǎn)。最后根據(jù)不同條件下模擬的結(jié)果，對比選出該設(shè)備在某些特定條件下的最佳工況的相關(guān)參數(shù)。數(shù)值模擬在處理一些幾何形狀復(fù)雜，或包含某些非線性特征的問題時，數(shù)值模擬可以簡單地分為以下幾個步驟：

1）針對特定工程問題，建立數(shù)學(xué)模型；2）模型建立完成后，就是確定計算方法；3）編制程序并進行計算；4）通過圖像展示的結(jié)果進行分析，得出結(jié)論［4］。

2 離散化

離散化，就是把無限空間中有限的個體通過映射放到有限的空間中去，以此來提高算法的時空效率［5］，是在有限個體原始數(shù)據(jù)相對大小不變的情況下，按比例進行縮小后來建立相應(yīng)的模型。離散化的應(yīng)用在程序設(shè)計中非常普遍，在諸多可能的情況下，離散化只需考慮要使用到的值。這就大大簡化了實驗，同時也有效地降低了計算機處理過程中的時間復(fù)雜度。離散化可以提高一個算法的質(zhì)量，而且有可能實現(xiàn)原本不會實現(xiàn)的某些算法。

區(qū)域離散化是將要計算的空間區(qū)域，劃分成大量互不重疊的子區(qū)域，每個子區(qū)域有一個節(jié)點，以及此節(jié)點代表的一個控制體積，最后需要確定上述節(jié)點在空間中所處的位置和其代表的控制體積。根據(jù)選用的研究方法不同，控制體積可能是流場中一個固定的區(qū)域，也可能會隨流體流動位置發(fā)生變化。網(wǎng)格節(jié)點一般會被看成是控制體積的代表，因為每一個節(jié)點必然對應(yīng)一個控制體積，不會出現(xiàn)一對多或多對一的情況，且每個節(jié)點代表的位置不會重合［6］。控制體積作為一個特定的區(qū)域，并不能全部與最開始劃分出來的子區(qū)域重合，這種情況也不會影響離散質(zhì)量。

3 離散化方法及特點

3.1 有限差分法

有限差分法，是數(shù)值解法中發(fā)展時間較長，相對來說較為完善，也最為經(jīng)典的一種離散方法［7］。差分過程是將需要求解的空間區(qū)域劃分為大量的網(wǎng)格，根據(jù)特定的幾何形狀選擇適當(dāng)?shù)膭澐址绞剑@樣有助于簡化轉(zhuǎn)化方程的步驟。劃分網(wǎng)格后，用有限個網(wǎng)格上的節(jié)點代替原來連續(xù)的求解域［8］，選用適當(dāng)?shù)牟钌绦问剑言酒⒎址匠?（控制方程）中的偏導(dǎo)數(shù)項全部用差商表示出來，再經(jīng)過一些處理就能得到對應(yīng)的差分方程。

有限差分法將原來復(fù)雜的偏微分方程問題直接轉(zhuǎn)換為代數(shù)問題的求解，不需要構(gòu)造特定的函數(shù)，建模簡單，編程也相對容易。適用曲線型和拋物線型相關(guān)問題的求解，橢圓形型問題應(yīng)用較少，對于邊界條件復(fù)雜的問題也表現(xiàn)出相對劣勢［9］。

3.2 有限元法

有限元法，是將空間里一個連續(xù)的求解域，劃分許多合適大小、形狀的單元，分別在每個單元上構(gòu)造插值函數(shù)，再對這些插值函數(shù)使用相應(yīng)的極值原理，然后在這些所有的微小單元上都能找到一個特定的有限元方程，最后用這些有限元方程把初始問題的控制方程所代替。經(jīng)過如此轉(zhuǎn)換后，則所有微小單元極值的和就代表了該系統(tǒng)總體的極值，即將系統(tǒng)劃分為局部單元，分別處理各局部單元后，再進行總體合成，于是使原本問題的初始控制方程就變成了含有特定邊界條件的代數(shù)方程組。

有限元法在求解橢圓型問題時，有較好的適應(yīng)性。相比其它方法，處理問題的過程較慢，對復(fù)雜問題分析時，占用內(nèi)存較大，耗費的計算資源也很驚人，在使用時對操作者的經(jīng)驗要求較高。所以它在商用CFD軟件中應(yīng)用的不多，F(xiàn)IDAP軟件是基于有限元方法建立的，在這個軟件里較為常用［10］。

3.3 有限體積法

有限體積法，是將計算區(qū)域劃分為大量的小網(wǎng)格，每個網(wǎng)格點所處的位置不同，所以其代表的一個在它周圍的控制體積肯定也不同。這些控制體積不會相互重合，但是在求解的時候必須應(yīng)用到他們之間的關(guān)系。用初始微分方程的形式，在這些控制體積上應(yīng)用并進行積分，其中的偏導(dǎo)數(shù)項就會消失，對應(yīng)每個網(wǎng)格點會得到一個相應(yīng)的離散方程，聯(lián)立這些方程，就會得到一個與所劃分網(wǎng)格形式對應(yīng)的離散方程組。離散方程組中的未知數(shù)會因節(jié)點所處位置不同而發(fā)生變化。

有限體積法在處理過程中需要尋求節(jié)點值，這點十分類似于前面所說的有限差分法；在積分過程中時，有限體積法需要根據(jù)網(wǎng)格特點，對物理量值的分布規(guī)律做出假設(shè)，這一點又與有限元法雷同。

有限體積法在使用過程中始終遵循積分守恒原理，因此不管針對微小計算單元還是結(jié)構(gòu)整體，都有良好的守恒性；其插值函數(shù)可以根據(jù)具體問題靈活假設(shè)，而且對于泰勒展開造成的離散能有效的避免；對于有些結(jié)構(gòu)復(fù)雜的工程問題，需要分區(qū)域劃分網(wǎng)格時，能夠很好地適應(yīng)；與有限元法融合，綜合彼此的優(yōu)點，更好的進行流固耦合分析［11］。在對控制體積積分進行計算時，有限積分法需要使用到插值函數(shù)。離散方程在這個插值函數(shù)的作用下，比較容易得出，而且有限積分法不要求微分方程里的所有項必須采取同一種插值函數(shù)。這也是其優(yōu)于其他離散方法的一個方面。

4 結(jié)語

在數(shù)值模擬過程中，常用的離散化方法有三個，都有各自的特點：

1）對于有限差分法，現(xiàn)在發(fā)展的已經(jīng)非常成熟，利用差分直接將偏微分問題簡化為代數(shù)問題，處理結(jié)構(gòu)性網(wǎng)格有較大的優(yōu)勢，對于非結(jié)構(gòu)網(wǎng)格差分效果較差［12］。

2）對于有限元法，是一種具有高效的處理能力、且使用較多的計算方法。適用于一些幾何和物理條件比較復(fù)雜的情況，尤其是在以拉普拉斯和泊松方程所描述的物理場中，而且和程序的契合相對較好，這樣就有利于程序的標(biāo)準(zhǔn)化設(shè)計。

3）對于有限體積法，其原理非常簡單，積分守恒在這個方法中得到了充分的體現(xiàn)，不管是在無限小控制體積還是整個計算區(qū)域，所有因變量都遵循積分守恒原理。