計算流體力學(xué)（Computational Fluid 動力學(xué)）20世紀(jì)50年代以來，隨著計算機的發(fā)展而產(chǎn)生的一個介于數(shù)學(xué)、流體力學(xué)和計算機之間的交叉學(xué)科，主要研究內(nèi)容是通過計算機和數(shù)值方法來求解流體力學(xué)的控制方程，對流體力學(xué)問題進(jìn)行模擬和分析。
正文
用電子計算機和離散化的數(shù)值方法對流體力學(xué)問題進(jìn)行數(shù)值模擬和分析的一個新分支。
簡史? 流體力學(xué)和其他學(xué)科一樣，是通過理論分析和實驗研究兩種手段發(fā)展起來的。很早就已有理論流體力學(xué)和實驗流體力學(xué)兩大分支。理論分析是用數(shù)學(xué)方法求出問題的定量結(jié)果。但能用這種方法求出結(jié)果的問題畢竟是少數(shù)，計算流體力學(xué)正是為彌補分析方法的不足而發(fā)展起來的。
早在20世紀(jì)初，理查德就已提出用數(shù)值方法來解流體力學(xué)問題的思想。但是由于這種問題本身的復(fù)雜性和當(dāng)時計算工具的落后，這一思想并未引起人們重視。自從40年代中期電子計算機問世以來，用電子計算機進(jìn)行數(shù)值模擬和計算才成為現(xiàn)實。1963年美國的F.H.哈洛和J.E.艾瑞克·弗洛姆用當(dāng)時的IBM7090計算機，成功地解決了二維長方形柱體的繞流問題并給出尾流渦街的形成和演變過程，受到普遍重視。1965年，哈洛和弗羅姆發(fā)表“流體動力學(xué)的計算機實驗”一文，對計算機在流體力學(xué)中的巨大作用作了引人注目的介紹。從此，人們把60年代中期看成是計算流體力學(xué)興起的標(biāo)志。
計算流體力學(xué)的歷史雖然不長，但已廣泛深入到流體力學(xué)的各個領(lǐng)域，相應(yīng)地也形成了各種不同的數(shù)值解法。就目前情況看，主要是有限差分方法和有限元法。有限差分方法在流體力學(xué)中已得到廣泛應(yīng)用。而有限元法是從求解固體力學(xué)問題發(fā)展起來的。近年來在處理低速流體問題中，已有相當(dāng)多的應(yīng)用，而且還在迅速發(fā)展中。
流體力學(xué)基本方程? 為了說明計算流體力學(xué)主要方法，需先了解流體力學(xué)運動的基本方程的性質(zhì)和分類。流體力學(xué)的基本方程是在19 世紀(jì)上半葉由C.-L.-M.-H.納維和G.G.喬治·斯托克斯等人建立的，稱為納維－斯托克斯方程，簡稱N-S方程。二維非定常不可壓縮流體的N-S方程為：
式中u、v為沿著x、y方向上的速度分量；t為時間；p為壓力；ρ為密度；ν為運動粘性系數(shù)。在不同條件下,N-S方程的數(shù)學(xué)性質(zhì)也不一樣。
①N-S方程描述粘性流體隨時間而變的非定常運動。時間項和方程右邊的高階導(dǎo)數(shù)項決定方程的性質(zhì)。它同二維熱傳導(dǎo)方程類似，屬于拋物型方程。
②粘性流體的定常運動是將原方程中的時間項省去。此時N-S方程的性質(zhì)，取決于它的高階導(dǎo)數(shù)項，和拉普拉斯方程一樣，為橢圓型方程。
③無粘流的歐拉方程是將 N-S方程的右邊粘性項略去而得。它也適用于可壓縮流體。從形式上不容易判斷歐拉方程的性質(zhì)。因多數(shù)無粘流動皆為無旋流動，故如將歐拉方程改用速度勢嗞表示，則二維定?？蓧嚎s氣流的方程為：
式中c為聲速。此式是二階偏微分方程
的一般形式，其性質(zhì)要看而定。在超聲速區(qū)，上式類似于波動方程，為雙曲型；在亞聲速區(qū)，上式便與拉普拉斯方程相同，為橢圓型?？傊?，流體力學(xué)的運動方程是極其復(fù)雜的非線性偏微分方程，具有各種不同的類型，而且往往還是混合型的。要全面描述流體的運動，還必須同時考慮其他方程，如連續(xù)性方程、能量方程和狀態(tài)方程等。所以計算流體力學(xué)在很大程度上就是針對不同性質(zhì)的偏微分方程采用和發(fā)展相應(yīng)的數(shù)值解方法。
低速無粘流動數(shù)值解? 在無旋條件下，低速流動的速度勢滿足拉普拉斯方程或泊松方程。很多平面問題利用復(fù)變函數(shù)和保角映射可以求得解析解，這是經(jīng)典流體力學(xué)的重要內(nèi)容。但對幾何形狀比較復(fù)雜的物體，必須用下述的數(shù)值解法。
①迭代解法 這是用逐步近似求解聯(lián)立方程的方法，也是橢圓型微分方程的主要數(shù)值解法。此法程序簡單，存儲量與運算量均比較小，一般先假定一組初值，然后求每個網(wǎng)點上的新值。以五點格式為例，網(wǎng)點上的新值是鄰近四點初值的平均。新值求出后，舊值還要保留，以便計算其他各點的新值。這種簡單迭代收斂很慢，現(xiàn)已很少使用。但若稍加改進(jìn)，用算出的新值沖掉舊值，并引進(jìn)一個松弛因子，以加速收斂，將算出來的新值與原來的舊值加權(quán)平均，就成為50年代發(fā)展起來的逐次超松弛法。
②時間相關(guān)法 這是用非定常方程求解定常問題的方法，常用于求解N-S方程和歐拉方程等。雖然用的是非定常方程，但所解的并不是非定常問題。根據(jù)給定的初始條件以及隨時間改變的約束條件，非定常問題是研究流動隨時間的演變過程。這種非定常行為和給出的初值很有關(guān)系。然而時間相關(guān)法的初值，原則上是隨意選取的，只是須滿足定常問題所規(guī)定的邊界條件。在求解過程中，流動隨時間的變化并不代表真實的物理過程。當(dāng)時間足夠長后，未知函數(shù)值逐步與時間無關(guān)，便漸近趨于定常解。所以時間相關(guān)法實際上也是一種迭代法，時間變量只不過是用來記錄迭代的次數(shù)而已。
③交替方向隱式法 流體力學(xué)的應(yīng)用問題，往往是二維和三維的空間問題。由于穩(wěn)定性的要求，時間步長集團(tuán)受維數(shù)的限制，維數(shù)愈高，要求時間步長愈小，計算工作量也愈大。50年代中期D.W.裴斯曼和J.道格拉斯等人提出所謂交替方向隱式法，以加快計算速度。如在二維非定常方程中，第一步先對x的導(dǎo)數(shù)用隱式差分，而y方向的導(dǎo)數(shù)則用前一個的數(shù)值。第二步對 y的導(dǎo)數(shù)用隱式差分,x方向的導(dǎo)數(shù)則用第一步算出來的數(shù)值。這一方法的優(yōu)點是穩(wěn)定性好，有足夠的二階精度，所產(chǎn)生的差分方程是三對角矩陣方程，便于求解。
④有限基本解法 解位勢流動的一種數(shù)值方法。中國航空工業(yè)集團(tuán)有限公司中的低速飛機設(shè)計采用位勢理論計算各種氣動力參數(shù)，就是求解二維或三維拉普拉斯方程。在經(jīng)典流體力學(xué)中，用基本解的疊加來解拉普拉斯方程的做法是很成功的。這種方法的要點是，用源、匯、偶極子的分布代替機翼和機身對流場的影響。它們的強度由邊界條件確定，結(jié)果需要求解積分方程。對一些簡單情況可以求解，對一般情況則比較困難。高速電子計算機的出現(xiàn)使這種積分方程的數(shù)值解法也有了突破。其主要思想是把積分方程離散化，積分方程代表源、匯等奇點在空間連續(xù)分布的總和。例如，若把機翼和機身表面，分割成若干個小單元，每個單元上的奇點強度取平均值。把這些奇點的總和疊加起來，就得出流場總的效應(yīng)。因此，它用有限項的求和來代替積分，而最后要解的是一組代數(shù)方程。由于基本解都是具有奇點的函數(shù)，所以這種方法又稱為有限奇點法或鱗片法。（見有限基本解法）跨聲速流動數(shù)值解? 跨聲速流動的流場是既有亞聲速區(qū)又有超聲速區(qū)的一種混合流場。在不考慮粘性影響和小擾動的情況下，定常二維速度勢方程是混合型的，即V嗞xx+嗞yy＝0，式中V是來流馬赫數(shù)Ma∞與嗞x的復(fù)雜函數(shù)。V＞0是亞聲速區(qū)(橢圓型)，而V＜0為超聲速區(qū)（雙曲型）。美國的E.M.穆曼和J.D.科爾在1971年首先采用混合差分格式，并運用松弛法成功地解出定常小擾動速度勢方程。混合差分格式就是在亞聲速區(qū)用中心差分格式，所有鄰近網(wǎng)點上的條件都會影響計算點，而在超聲速區(qū)，則用迎風(fēng)格式，因為上游迎風(fēng)網(wǎng)點正好是雙曲型波動方程的依賴區(qū)。（見跨聲速流數(shù)值計算）
超聲速流動數(shù)值解? 在超聲速流動中，主要問題是如何處理激波。用數(shù)值方法處理超聲速流場中的激波現(xiàn)有兩種方法。一是激波捕捉法，另一是激波裝配法。激波捕捉法對激波本身并不需作任何特殊處理，只是在計算公式中，直接或間接地引進(jìn)“粘性”項，自動算出激波的位置和強度，以“捕捉”激波。其中又有所謂人工粘性和格式粘性兩種方法。人工粘性方法是J.von諾伊曼和R.D.里希特邁爾于1950年首先提出的，它是以真實粘性流體的物理理論為基礎(chǔ)的一種自動處理激波近似方法。該法是在激波層內(nèi)，人為地加入粘性項，使激波間斷變成光滑的過渡區(qū)。近年來，在超聲速流動中得到廣泛的應(yīng)用。格式粘性是通過某種差分格式間接地引入粘性項拉克斯格式。拉克斯-文德霍夫格式和麥克馬克格式都具有類似的效果。激波裝配法是把激波仍當(dāng)作間斷面來處理，激波前后要滿足激波跳躍條件。但是在普通坐標(biāo)中，它的實現(xiàn)很困難。一般采用坐標(biāo)變換，使激波位置（此時是未知的）和一個坐標(biāo)軸重合，然后把激波看作內(nèi)邊界。這種處理是比較精確的，但也是很麻煩和不方便的。最好的辦法是把激波捕捉法和激波裝配法結(jié)合起來。例如在流場外圍的離體激波用激波裝配法，在流場內(nèi)的激波用激波捕捉法。（見超聲速無粘繞流數(shù)值解）
粘性流動數(shù)值解法? 可參見納維－斯托克斯方程數(shù)值解、邊界層方程數(shù)值解法和湍流數(shù)值計算等。
參考書目
朱幼蘭等著：《初邊值問題差分方法及繞流》，科學(xué)出版社，北京市，1980。
P.J.羅奇著，鐘錫昌、劉學(xué)宗譯：《計算流體力學(xué)》，科學(xué)出版社，北京，1983年。P.J.Roache,Computation Fluid 動力學(xué),Hermosa Pub., 阿爾伯克基,1972.
參考資料 >