汽車空氣動力學(xué) 教學(xué)課件 第八章 汽車空氣動力學(xué)數(shù)值計算

KQ4主編KQ4主編1第八章汽車空氣動力學(xué)數(shù)值計算第一節(jié)空氣動力學(xué)數(shù)值計算概述

第二節(jié)CFD在汽車空氣動力學(xué)中的應(yīng)用

第三節(jié)非粘性流方法

第四節(jié)N—S方程的簡化

第五節(jié)N—S方程的解法

第六節(jié)結(jié)束語第八章汽車空氣動力學(xué)數(shù)值計算第一節(jié)空氣動力學(xué)數(shù)值計算概述2第一節(jié)空氣動力學(xué)數(shù)值計算概述一、空氣動力學(xué)的研究方法

二、CFD的一般方法

三、計算流體力學(xué)的發(fā)展過程

四、CFD的理論基礎(chǔ)及現(xiàn)狀第一節(jié)空氣動力學(xué)數(shù)值計算概述一、空氣動力學(xué)的研究方法

二、3一、空氣動力學(xué)的研究方法CFD是建立在全Navier—Strokes(以下稱N—S)方程近似解基礎(chǔ)上的計算技術(shù)。根據(jù)近似解的精度等級，把N—S方程的解法分成以下四類：

1)線性非粘性流方法；

2)非線性非粘性流方法；

3)平均雷諾數(shù)基礎(chǔ)上的N—S方程；

4)全N—S方程。一、空氣動力學(xué)的研究方法CFD是建立在全Navier—Str4CFD有多種多樣的方法，最常用的方法有以下三種：1)有限差分法(finite-difference

2)有限要素法(finite

3)邊界積分法(boundaryintegrlmethod，BIM)。二、CFD的一般方法CFD有多種多樣的方法，最常用的方法有以下三種：二、CFD的5圖8-1CFD系統(tǒng)圖二、CFD的一般方法圖8-1CFD系統(tǒng)圖二、CFD的一般方法6二、CFD的一般方法二、CFD的一般方法7二、CFD的一般方法二、CFD的一般方法8三、計算流體力學(xué)的發(fā)展過程圖8-2CFD的發(fā)展過程三、計算流體力學(xué)的發(fā)展過程圖8-2CFD的發(fā)展過程9四、CFD的理論基礎(chǔ)及現(xiàn)狀(一)無粘流理論

(二)粘性流理論

(三)CFD存在的問題及當(dāng)前研究的重點(diǎn)四、CFD的理論基礎(chǔ)及現(xiàn)狀(一)無粘流理論

(二)粘性流理論10(一)無粘流理論1.計算結(jié)果與試驗(yàn)結(jié)果的差異

2.不同計算方法的計算結(jié)果的差異(一)無粘流理論1.計算結(jié)果與試驗(yàn)結(jié)果的差異

2.不同計算方111.計算結(jié)果與試驗(yàn)結(jié)果的差異1977年，曾有人用有限要素法計算了兩輛小客車相鄰行駛時的氣動力。計算結(jié)果表明，往前超車的一輛車由于另一輛車的干擾，產(chǎn)生一個使其偏離鄰車的橫擺力矩，這與實(shí)車試驗(yàn)結(jié)果定性一致。1982年又有人用該方法計算了小轎車表面的壓力分布(β＝0°)狀態(tài)。計算時，沒有考慮車身底部外形以及車輪的影響，作了平滑車底的近似。1.計算結(jié)果與試驗(yàn)結(jié)果的差異1977年，曾有人用有限要素法計122.不同計算方法的計算結(jié)果的差異放棄無旋條件則無粘運(yùn)動方程成為非線性的歐拉方程。歐拉法與有限要素法的不同之處在于：有限要素法是把物面分成許多單元，而歐拉法是將物體周圍的流場進(jìn)行單元劃分。2.不同計算方法的計算結(jié)果的差異放棄無旋條件則無粘運(yùn)動方程成13(二)粘性流理論1.計算結(jié)果與試驗(yàn)結(jié)果在分離區(qū)的差異

2.計算結(jié)果與試驗(yàn)結(jié)果雖然定性相似，但定量差異較大(二)粘性流理論1.計算結(jié)果與試驗(yàn)結(jié)果在分離區(qū)的差異

2.計141.計算結(jié)果與試驗(yàn)結(jié)果在分離區(qū)的差異用邊界法將物體周圍的流場分為粘性流區(qū)和非粘性流區(qū)兩個區(qū)域。在物體附近很薄的一層稱為邊界層，考慮了粘性的影響。邊界層外為非粘性流，用上述無粘流方法進(jìn)行處理。邊界層內(nèi)為粘性流。最近有人用這種方法計算了小轎車的壓力分布，計算結(jié)果表明，由于在分離區(qū)沒有考慮分離的影響，分離區(qū)的計算結(jié)果與試驗(yàn)結(jié)果有較大差異。1.計算結(jié)果與試驗(yàn)結(jié)果在分離區(qū)的差異用邊界法將物體周圍的流場152.計算結(jié)果與試驗(yàn)結(jié)果雖然定性相似，但定量差異較大目前進(jìn)行了關(guān)于分離和渦流的CFD研究并積極地利用CFD對風(fēng)洞試驗(yàn)進(jìn)行補(bǔ)充。現(xiàn)在用CFD可解決下列問題：1)用CFD預(yù)測分離和渦流較少的車身前部的表面壓力分布。

2)對流態(tài)顯示比較困難的車身底部，可用CFD得到的數(shù)據(jù)對風(fēng)洞試驗(yàn)進(jìn)行補(bǔ)充。

3)用CFD方法進(jìn)行渦流和分離少的低阻車的設(shè)計和開發(fā)，從而縮短設(shè)計、開發(fā)的時間。2.計算結(jié)果與試驗(yàn)結(jié)果雖然定性相似，但定量差異較大目前進(jìn)行了16(三)CFD存在的問題及當(dāng)前研究的重點(diǎn)1)目前尚無成熟的非線性偏微分方程數(shù)值方法的數(shù)學(xué)理論，尚未建立嚴(yán)格的誤差估計及嚴(yán)格的穩(wěn)定性分析方法。

2)數(shù)值方法只能給出離散數(shù)據(jù)或單元近似函數(shù)，因而不能完全代替解析解。

3)數(shù)值方法只能分辨相對網(wǎng)格尺寸而言的長波現(xiàn)象，而為提高分辨率減小網(wǎng)格，又受到計算機(jī)容量、速度的限制。

(三)CFD存在的問題及當(dāng)前研究的重點(diǎn)1)目前尚無成熟的非線17目前，CFD主要研究下述問題：1)基本方程和湍流模型的處理；

2)數(shù)值方法和計算格式的發(fā)展；

3)計算坐標(biāo)的選取、物面描述和網(wǎng)格生成；

4)邊界條件等影響計算精度的因素的處理；

5)CFD計算結(jié)果與試驗(yàn)結(jié)果的對比分析以及采用PIV法等流態(tài)顯示技術(shù)的發(fā)展。(三)CFD存在的問題及當(dāng)前研究的重點(diǎn)目前，CFD主要研究下述問題：(三)CFD存在的問題及當(dāng)前18第二節(jié)CFD在汽車空氣動力學(xué)中的應(yīng)用一、汽車流場數(shù)值模擬的特點(diǎn)和難點(diǎn)

二、流場中作用的力

三、理想流體的計算準(zhǔn)則

四、空氣的兩種近似法第二節(jié)CFD在汽車空氣動力學(xué)中的應(yīng)用一、汽車流場數(shù)值模擬的19一、汽車流場數(shù)值模擬的特點(diǎn)和難點(diǎn)汽車周圍的流場是非常復(fù)雜的，它存在著顯著的分離現(xiàn)象，如：具有大曲面的二維或全三維流動產(chǎn)生空泡形分離渦，與車身形狀有關(guān)的尾流結(jié)構(gòu)以及頂部、底部、側(cè)面氣流之間產(chǎn)生強(qiáng)烈的相互影響。一、汽車流場數(shù)值模擬的特點(diǎn)和難點(diǎn)汽車周圍的流場是非常復(fù)雜的，20二、流場中作用的力圖8-3作用于流體上的力二、流場中作用的力圖8-3作用于流體上的力21二、流場中作用的力（8-1）（8-2）（8-3）（8-4）二、流場中作用的力（8-1）22二、流場中作用的力（8-5）（8-6）（8-7）二、流場中作用的力（8-5）23（8-8）（8-9）（8-10）（8-11）三、理想流體的計算準(zhǔn)則（8-8）三、理想流體的計算準(zhǔn)則24（8-12）（8-13）（8-14）（8-15）三、理想流體的計算準(zhǔn)則（8-12）三、理想流體的計算準(zhǔn)則25四、空氣的兩種近似法前面已分析了汽車?yán)@流研究的難點(diǎn)，為了使汽車空氣動力學(xué)問題易于處理，必須對空氣的物理特性進(jìn)行近似和簡化。汽車空氣動力學(xué)把空氣分為兩種類型，即完全流體和粘性不可壓縮流體。四、空氣的兩種近似法前面已分析了汽車?yán)@流研究的難點(diǎn)，為了使汽26第三節(jié)非粘性流方法一、控制方程

二、作為歐拉方程解的無旋流場

三、無旋流

四、渦格法和面元法的基本概念

五、渦格法

六、面元法

七、非粘性無粘流方法第三節(jié)非粘性流方法一、控制方程

二、作為歐拉方程解的無旋流27一、控制方程（8-16）（8-17）一、控制方程（8-16）28二、作為歐拉方程解的無旋流場（8-18）（8-19）（8-20）二、作為歐拉方程解的無旋流場（8-18）29三、無旋流1.速度場與壓力場

2.流線函數(shù)三、無旋流1.速度場與壓力場

2.流線函數(shù)301.速度場與壓力場（8-21）（8-22）（8-23）（8-24）1.速度場與壓力場（8-21）311.速度場與壓力場（8-25）（8-26）（8-27）1.速度場與壓力場（8-25）322.流線函數(shù)（8-29）（8-30）（8-31）（8-32）2.流線函數(shù)（8-29）332.流線函數(shù)（8-33）（8-34）（8-35）2.流線函數(shù)（8-33）342.流線函數(shù)圖8-4流線和等勢線2.流線函數(shù)圖8-4流線和等勢線35四、渦格法和面元法的基本概念用渦格法和面元法模擬繞任意形狀物體的勢流，允許將流動視為兩個流場疊加而成的流場，其一為定常v∞的均勻“原始”流場，另一個是具有速度為vp(x，y，z)的“擾動”流場。均勻原始流定義為沒有物體時的流動，擾動流則使流場的每一點(diǎn)實(shí)際速度為用?、?0和?p分別表示結(jié)果流動、原始流動和干擾流動的勢，上述兩個邊界條件可表示為（8-36）（8-37）四、渦格法和面元法的基本概念用渦格法和面元法模擬繞任意形狀物36五、渦格法圖8-5渦格法

a)渦流網(wǎng)格表面分割b)分離表面的馬蹄渦的排列五、渦格法圖8-5渦格法

a)渦流網(wǎng)格表面分割b)分離表37五、渦格法（8-38）（8-39）（8-40）五、渦格法（8-38）38五、渦格法圖8-6沿轎車上表面中心線的壓力分布五、渦格法圖8-6沿轎車上表面中心線的壓力分布39六、面元法圖8-7用面元法進(jìn)行汽車模型數(shù)值模擬的基本特征六、面元法圖8-7用面元法進(jìn)行汽車模型數(shù)值模擬的基本特征40六、面元法（8-41）（8-42）（8-43）（8-44）（8-45）六、面元法（8-41）41六、面元法面元法是用被稱為面元的平面四邊形或三角形的面積單元代替連續(xù)的車身和尾流表面，假定面元上的點(diǎn)源或偶極子強(qiáng)度為常數(shù)。這種方法的特點(diǎn)是：1)通過表面離散化，把式(8-44)中的面積分簡化為每個面元上的有限積分的和。

2)設(shè)面元上σ和μ為常數(shù)，則將它們從積分中提出，使方程線性化。若車身面元數(shù)為m，尾流面元數(shù)為n，則式(8-45)中的積分可表示為（8-46）（8-47）六、面元法面元法是用被稱為面元的平面四邊形或三角形的面積單元42六、面元法（8-49）（8-50）六、面元法（8-49）43七、非粘性無粘流方法1.不可壓流非線性方程的解

2.可壓流歐拉方程的解七、非粘性無粘流方法1.不可壓流非線性方程的解

2.可壓流歐441.不可壓流非線性方程的解基于解不可壓、無粘、非定常流的非線性方程，將模擬尾流分離區(qū)的形成視為時間相關(guān)過程［式(8-16)］，通過微分消除方程中的壓力項，得出僅含速度項u、v、w的方程，引入渦量W，即（8-51）

（8-52）（8-53）（8-54）（8-55）1.不可壓流非線性方程的解基于解不可壓、無粘、非定常流的非線451.不可壓流非線性方程的解圖8-8各種時間步長后，在車尾下底邊處，渦的邊緣和剪切層上卷的發(fā)展

a)5個時間步長b)10個時間步長c)15個時間步長1.不可壓流非線性方程的解圖8-8各種時間步長后，在車尾下46圖8-9CFD結(jié)果與試驗(yàn)結(jié)果的比較2.可壓流歐拉方程的解圖8-9CFD結(jié)果與試驗(yàn)結(jié)果的比較2.可壓流歐拉方程的解47圖8-10二維歐拉方程解的一例2.可壓流歐拉方程的解圖8-10二維歐拉方程解的一例2.可壓流歐拉方程的解48第四節(jié)N—S方程的簡化1.不同計算方法的計算結(jié)果的差異

2.邊界條件

3.面積分割第四節(jié)N—S方程的簡化1.不同計算方法的計算結(jié)果的差異

2491.不同計算方法的計算結(jié)果的差異圖8-11快背式車身表面二維計算的例子

a)CFD各種方法的比較b)VW公司對Audi2000轎車的壓力分布計算結(jié)果與試驗(yàn)結(jié)果的比較1.不同計算方法的計算結(jié)果的差異圖8-11快背式車身表面二50圖8-12用差分法和有限體積法計算車身表面壓力分布的差異1.不同計算方法的計算結(jié)果的差異圖8-12用差分法和有限體積法計算車身表面壓力分布的差異1511.不同計算方法的計算結(jié)果的差異由于網(wǎng)格分割的密度對各種計算方法的不同影響，使不同計算方法的計算結(jié)果不同。圖8-12表明用差分法和有限體積法對車身表面壓力分布計算結(jié)果的影響。計算采用k-ε湍流模型。由圖可見，計算結(jié)果雖定性相似，但定量分析存在如下差別：1)采用有限體積法的計算結(jié)果表明，在車身頂蓋前、后端棱角部位，Cp有較大的峰值。

2)采用差分法計算車身下部的壓力分布，得出較低的Cp值。1.不同計算方法的計算結(jié)果的差異由于網(wǎng)格分割的密度對各種計算522)采用差分法計算車身下部的壓力分布，得出較低的Cp值。圖8-13用差分法和有限體積法計算車身尾流的差異

a)數(shù)值解法對尾流速度分布的影響b)數(shù)值解法對尾流能量分布的影響2)采用差分法計算車身下部的壓力分布，得出較低的Cp值。圖8532.邊界條件2.邊界條件54圖8-14采用不同邊界條件計算車身上、下表面壓力分布的結(jié)果

a)邊界條件對車身表面壓力分布的影響b)固定壁對尾流的影響c)移動壁對尾流的影響2.邊界條件圖8-14采用不同邊界條件計算車身上、下表面壓力分布的結(jié)果553.面積分割圖8-15網(wǎng)格分割密度對車身表面壓力分布計算結(jié)果的影響

1—采用160×80等間距(車身為25×12等間距)網(wǎng)格分割

2—采用149×69不等間距(車身為51×29等間距)網(wǎng)格分割3.面積分割圖8-15網(wǎng)格分割密度對車身表面壓力分布計算結(jié)563.面積分割圖8-16網(wǎng)格分割密度對車身表面三維3.面積分割圖8-16網(wǎng)格分割密度對車身表面三維573.面積分割817.TIF3.面積分割817.TIF583.面積分割3.面積分割593.面積分割圖8-18不同Re數(shù)下的CFD計算結(jié)果3.面積分割圖8-18不同Re數(shù)下的CFD計算結(jié)果60第五節(jié)N—S方程的解法圖8-19用差分法解層流定常N—S方程示例1第五節(jié)N—S方程的解法圖8-19用差分法解層流定常N—S61第五節(jié)N—S方程的解法圖8-20用差分法解層流定常N—S方程示例2第五節(jié)N—S方程的解法圖8-20用差分法解層流定常N—S62第五節(jié)N—S方程的解法圖8-21車身表面的網(wǎng)格劃分第五節(jié)N—S方程的解法圖8-21車身表面的網(wǎng)格劃分63第五節(jié)N—S方程的解法圖8-22車身表面的網(wǎng)格劃分

a)表面均勻網(wǎng)格劃分b)外形急劇變化處用密網(wǎng)格劃分第五節(jié)N—S方程的解法圖8-22車身表面的網(wǎng)格劃分

a)64第五節(jié)N—S方程的解法圖8-23車身周圍的壓力分布第五節(jié)N—S方程的解法圖8-23車身周圍的壓力分布65第五節(jié)N—S方程的解法圖8-24地板下部的流速矢量第五節(jié)N—S方程的解法圖8-24地板下部的流速矢量66第六節(jié)結(jié)束語本章介紹的許多數(shù)值計算方法，在航空工業(yè)領(lǐng)域已是常規(guī)的設(shè)計工具，而到目前為止，還沒有一種計算方法能對汽車的真實(shí)繞流進(jìn)行既定性又定量的描述。其原因之一是汽車的設(shè)計和使用都沒有像航空界那樣對CFD有迫切的需要；另一方面，在飛機(jī)設(shè)計初期就充分考慮了空氣動力學(xué)，飛機(jī)的外形使其周圍的流場沒有或很少出現(xiàn)分離。因此，大部分為航空工業(yè)開發(fā)的CFD技術(shù)是建立在附著流基礎(chǔ)上的。第六節(jié)結(jié)束語本章介紹的許多數(shù)值計算方法，在航空工業(yè)領(lǐng)域67KQ4主編KQ4主編68第八章汽車空氣動力學(xué)數(shù)值計算第一節(jié)空氣動力學(xué)數(shù)值計算概述

第二節(jié)CFD在汽車空氣動力學(xué)中的應(yīng)用

第三節(jié)非粘性流方法

第四節(jié)N—S方程的簡化

第五節(jié)N—S方程的解法

第六節(jié)結(jié)束語第八章汽車空氣動力學(xué)數(shù)值計算第一節(jié)空氣動力學(xué)數(shù)值計算概述69第一節(jié)空氣動力學(xué)數(shù)值計算概述一、空氣動力學(xué)的研究方法

二、CFD的一般方法

三、計算流體力學(xué)的發(fā)展過程

四、CFD的理論基礎(chǔ)及現(xiàn)狀第一節(jié)空氣動力學(xué)數(shù)值計算概述一、空氣動力學(xué)的研究方法

二、70一、空氣動力學(xué)的研究方法CFD是建立在全Navier—Strokes(以下稱N—S)方程近似解基礎(chǔ)上的計算技術(shù)。根據(jù)近似解的精度等級，把N—S方程的解法分成以下四類：

1)線性非粘性流方法；

2)非線性非粘性流方法；

3)平均雷諾數(shù)基礎(chǔ)上的N—S方程；

4)全N—S方程。一、空氣動力學(xué)的研究方法CFD是建立在全Navier—Str71CFD有多種多樣的方法，最常用的方法有以下三種：1)有限差分法(finite-difference

2)有限要素法(finite

3)邊界積分法(boundaryintegrlmethod，BIM)。二、CFD的一般方法CFD有多種多樣的方法，最常用的方法有以下三種：二、CFD的72圖8-1CFD系統(tǒng)圖二、CFD的一般方法圖8-1CFD系統(tǒng)圖二、CFD的一般方法73二、CFD的一般方法二、CFD的一般方法74二、CFD的一般方法二、CFD的一般方法75三、計算流體力學(xué)的發(fā)展過程圖8-2CFD的發(fā)展過程三、計算流體力學(xué)的發(fā)展過程圖8-2CFD的發(fā)展過程76四、CFD的理論基礎(chǔ)及現(xiàn)狀(一)無粘流理論

(二)粘性流理論

(三)CFD存在的問題及當(dāng)前研究的重點(diǎn)四、CFD的理論基礎(chǔ)及現(xiàn)狀(一)無粘流理論

(二)粘性流理論77(一)無粘流理論1.計算結(jié)果與試驗(yàn)結(jié)果的差異

2.不同計算方法的計算結(jié)果的差異(一)無粘流理論1.計算結(jié)果與試驗(yàn)結(jié)果的差異

2.不同計算方781.計算結(jié)果與試驗(yàn)結(jié)果的差異1977年，曾有人用有限要素法計算了兩輛小客車相鄰行駛時的氣動力。計算結(jié)果表明，往前超車的一輛車由于另一輛車的干擾，產(chǎn)生一個使其偏離鄰車的橫擺力矩，這與實(shí)車試驗(yàn)結(jié)果定性一致。1982年又有人用該方法計算了小轎車表面的壓力分布(β＝0°)狀態(tài)。計算時，沒有考慮車身底部外形以及車輪的影響，作了平滑車底的近似。1.計算結(jié)果與試驗(yàn)結(jié)果的差異1977年，曾有人用有限要素法計792.不同計算方法的計算結(jié)果的差異放棄無旋條件則無粘運(yùn)動方程成為非線性的歐拉方程。歐拉法與有限要素法的不同之處在于：有限要素法是把物面分成許多單元，而歐拉法是將物體周圍的流場進(jìn)行單元劃分。2.不同計算方法的計算結(jié)果的差異放棄無旋條件則無粘運(yùn)動方程成80(二)粘性流理論1.計算結(jié)果與試驗(yàn)結(jié)果在分離區(qū)的差異

2.計算結(jié)果與試驗(yàn)結(jié)果雖然定性相似，但定量差異較大(二)粘性流理論1.計算結(jié)果與試驗(yàn)結(jié)果在分離區(qū)的差異

2.計811.計算結(jié)果與試驗(yàn)結(jié)果在分離區(qū)的差異用邊界法將物體周圍的流場分為粘性流區(qū)和非粘性流區(qū)兩個區(qū)域。在物體附近很薄的一層稱為邊界層，考慮了粘性的影響。邊界層外為非粘性流，用上述無粘流方法進(jìn)行處理。邊界層內(nèi)為粘性流。最近有人用這種方法計算了小轎車的壓力分布，計算結(jié)果表明，由于在分離區(qū)沒有考慮分離的影響，分離區(qū)的計算結(jié)果與試驗(yàn)結(jié)果有較大差異。1.計算結(jié)果與試驗(yàn)結(jié)果在分離區(qū)的差異用邊界法將物體周圍的流場822.計算結(jié)果與試驗(yàn)結(jié)果雖然定性相似，但定量差異較大目前進(jìn)行了關(guān)于分離和渦流的CFD研究并積極地利用CFD對風(fēng)洞試驗(yàn)進(jìn)行補(bǔ)充。現(xiàn)在用CFD可解決下列問題：1)用CFD預(yù)測分離和渦流較少的車身前部的表面壓力分布。

2)對流態(tài)顯示比較困難的車身底部，可用CFD得到的數(shù)據(jù)對風(fēng)洞試驗(yàn)進(jìn)行補(bǔ)充。

3)用CFD方法進(jìn)行渦流和分離少的低阻車的設(shè)計和開發(fā)，從而縮短設(shè)計、開發(fā)的時間。2.計算結(jié)果與試驗(yàn)結(jié)果雖然定性相似，但定量差異較大目前進(jìn)行了83(三)CFD存在的問題及當(dāng)前研究的重點(diǎn)1)目前尚無成熟的非線性偏微分方程數(shù)值方法的數(shù)學(xué)理論，尚未建立嚴(yán)格的誤差估計及嚴(yán)格的穩(wěn)定性分析方法。

2)數(shù)值方法只能給出離散數(shù)據(jù)或單元近似函數(shù)，因而不能完全代替解析解。

3)數(shù)值方法只能分辨相對網(wǎng)格尺寸而言的長波現(xiàn)象，而為提高分辨率減小網(wǎng)格，又受到計算機(jī)容量、速度的限制。

(三)CFD存在的問題及當(dāng)前研究的重點(diǎn)1)目前尚無成熟的非線84目前，CFD主要研究下述問題：1)基本方程和湍流模型的處理；

2)數(shù)值方法和計算格式的發(fā)展；

3)計算坐標(biāo)的選取、物面描述和網(wǎng)格生成；

4)邊界條件等影響計算精度的因素的處理；

5)CFD計算結(jié)果與試驗(yàn)結(jié)果的對比分析以及采用PIV法等流態(tài)顯示技術(shù)的發(fā)展。(三)CFD存在的問題及當(dāng)前研究的重點(diǎn)目前，CFD主要研究下述問題：(三)CFD存在的問題及當(dāng)前85第二節(jié)CFD在汽車空氣動力學(xué)中的應(yīng)用一、汽車流場數(shù)值模擬的特點(diǎn)和難點(diǎn)

二、流場中作用的力

三、理想流體的計算準(zhǔn)則

四、空氣的兩種近似法第二節(jié)CFD在汽車空氣動力學(xué)中的應(yīng)用一、汽車流場數(shù)值模擬的86一、汽車流場數(shù)值模擬的特點(diǎn)和難點(diǎn)汽車周圍的流場是非常復(fù)雜的，它存在著顯著的分離現(xiàn)象，如：具有大曲面的二維或全三維流動產(chǎn)生空泡形分離渦，與車身形狀有關(guān)的尾流結(jié)構(gòu)以及頂部、底部、側(cè)面氣流之間產(chǎn)生強(qiáng)烈的相互影響。一、汽車流場數(shù)值模擬的特點(diǎn)和難點(diǎn)汽車周圍的流場是非常復(fù)雜的，87二、流場中作用的力圖8-3作用于流體上的力二、流場中作用的力圖8-3作用于流體上的力88二、流場中作用的力（8-1）（8-2）（8-3）（8-4）二、流場中作用的力（8-1）89二、流場中作用的力（8-5）（8-6）（8-7）二、流場中作用的力（8-5）90（8-8）（8-9）（8-10）（8-11）三、理想流體的計算準(zhǔn)則（8-8）三、理想流體的計算準(zhǔn)則91（8-12）（8-13）（8-14）（8-15）三、理想流體的計算準(zhǔn)則（8-12）三、理想流體的計算準(zhǔn)則92四、空氣的兩種近似法前面已分析了汽車?yán)@流研究的難點(diǎn)，為了使汽車空氣動力學(xué)問題易于處理，必須對空氣的物理特性進(jìn)行近似和簡化。汽車空氣動力學(xué)把空氣分為兩種類型，即完全流體和粘性不可壓縮流體。四、空氣的兩種近似法前面已分析了汽車?yán)@流研究的難點(diǎn)，為了使汽93第三節(jié)非粘性流方法一、控制方程

二、作為歐拉方程解的無旋流場

三、無旋流

四、渦格法和面元法的基本概念

五、渦格法

六、面元法

七、非粘性無粘流方法第三節(jié)非粘性流方法一、控制方程

二、作為歐拉方程解的無旋流94一、控制方程（8-16）（8-17）一、控制方程（8-16）95二、作為歐拉方程解的無旋流場（8-18）（8-19）（8-20）二、作為歐拉方程解的無旋流場（8-18）96三、無旋流1.速度場與壓力場

2.流線函數(shù)三、無旋流1.速度場與壓力場

2.流線函數(shù)971.速度場與壓力場（8-21）（8-22）（8-23）（8-24）1.速度場與壓力場（8-21）981.速度場與壓力場（8-25）（8-26）（8-27）1.速度場與壓力場（8-25）992.流線函數(shù)（8-29）（8-30）（8-31）（8-32）2.流線函數(shù)（8-29）1002.流線函數(shù)（8-33）（8-34）（8-35）2.流線函數(shù)（8-33）1012.流線函數(shù)圖8-4流線和等勢線2.流線函數(shù)圖8-4流線和等勢線102四、渦格法和面元法的基本概念用渦格法和面元法模擬繞任意形狀物體的勢流，允許將流動視為兩個流場疊加而成的流場，其一為定常v∞的均勻“原始”流場，另一個是具有速度為vp(x，y，z)的“擾動”流場。均勻原始流定義為沒有物體時的流動，擾動流則使流場的每一點(diǎn)實(shí)際速度為用?、?0和?p分別表示結(jié)果流動、原始流動和干擾流動的勢，上述兩個邊界條件可表示為（8-36）（8-37）四、渦格法和面元法的基本概念用渦格法和面元法模擬繞任意形狀物103五、渦格法圖8-5渦格法

a)渦流網(wǎng)格表面分割b)分離表面的馬蹄渦的排列五、渦格法圖8-5渦格法

a)渦流網(wǎng)格表面分割b)分離表104五、渦格法（8-38）（8-39）（8-40）五、渦格法（8-38）105五、渦格法圖8-6沿轎車上表面中心線的壓力分布五、渦格法圖8-6沿轎車上表面中心線的壓力分布106六、面元法圖8-7用面元法進(jìn)行汽車模型數(shù)值模擬的基本特征六、面元法圖8-7用面元法進(jìn)行汽車模型數(shù)值模擬的基本特征107六、面元法（8-41）（8-42）（8-43）（8-44）（8-45）六、面元法（8-41）108六、面元法面元法是用被稱為面元的平面四邊形或三角形的面積單元代替連續(xù)的車身和尾流表面，假定面元上的點(diǎn)源或偶極子強(qiáng)度為常數(shù)。這種方法的特點(diǎn)是：1)通過表面離散化，把式(8-44)中的面積分簡化為每個面元上的有限積分的和。

2)設(shè)面元上σ和μ為常數(shù)，則將它們從積分中提出，使方程線性化。若車身面元數(shù)為m，尾流面元數(shù)為n，則式(8-45)中的積分可表示為（8-46）（8-47）六、面元法面元法是用被稱為面元的平面四邊形或三角形的面積單元109六、面元法（8-49）（8-50）六、面元法（8-49）110七、非粘性無粘流方法1.不可壓流非線性方程的解

2.可壓流歐拉方程的解七、非粘性無粘流方法1.不可壓流非線性方程的解

2.可壓流歐1111.不可壓流非線性方程的解基于解不可壓、無粘、非定常流的非線性方程，將模擬尾流分離區(qū)的形成視為時間相關(guān)過程［式(8-16)］，通過微分消除方程中的壓力項，得出僅含速度項u、v、w的方程，引入渦量W，即（8-51）

（8-52）（8-53）（8-54）（8-55）1.不可壓流非線性方程的解基于解不可壓、無粘、非定常流的非線1121.不可壓流非線性方程的解圖8-8各種時間步長后，在車尾下底邊處，渦的邊緣和剪切層上卷的發(fā)展

a)5個時間步長b)10個時間步長c)15個時間步長1.不可壓流非線性方程的解圖8-8各種時間步長后，在車尾下113圖8-9CFD結(jié)果與試驗(yàn)結(jié)果的比較2.可壓流歐拉方程的解圖8-9CFD結(jié)果與試驗(yàn)結(jié)果的比較2.可壓流歐拉方程的解114圖8-10二維歐拉方程解的一例2.可壓流歐拉方程的解圖8-10二維歐拉方程解的一例2.可壓流歐拉方程的解115第四節(jié)N—S方程的簡化1.不同計算方法的計算結(jié)果的差異

2.邊界條件

3.面積分割第四節(jié)N—S方程的簡化1.不同計算方法的計算結(jié)果的差異

21161.不同計算方法的計算結(jié)果的差異圖8-11快背式車身表面二維計算的例子

a)CFD各種方法的比較b)VW公司對Audi2000轎車的壓力分布計算結(jié)果與試驗(yàn)結(jié)果的比較1.不同計算方法的計算結(jié)果的差異圖8-11快背式車身表面二117圖8-12用差分法和有限體積法計算車身表面壓力分布的差異1.不同計算方法的計算結(jié)果的差異圖8-12用差分法和有限體積法計算車身表面壓力分布的差異11181.不同計算方法的計算結(jié)果的差異由于網(wǎng)格分割的密度對各種計算方法的不同影響，使不同計算方法的計算結(jié)果不同。圖8-12表明用差分法和有限體積法對車身表面壓力分布計算結(jié)果的影響。計算采用k-ε湍流模型。由圖可見，計算結(jié)果雖定性相似，但定量分析存在如下差別：1)采用有限體積法的計算結(jié)果表明，在車身頂蓋前、后端棱角部位，