第06章 fluent邊界條件_第1頁
第06章 fluent邊界條件_第2頁
第06章 fluent邊界條件_第3頁
第06章 fluent邊界條件_第4頁
第06章 fluent邊界條件_第5頁
已閱讀5頁,還剩46頁未讀 繼續(xù)免費閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認(rèn)領(lǐng)

文檔簡介

1、定義壁面的組分邊界條件FLUENT默認(rèn)所有的組分在壁面處具有零梯度條件(除了參加表面化學(xué)反應(yīng)的組分),但是可以指定壁面處的組分質(zhì)量分?jǐn)?shù)。也就是如同在入口處指定的Dirichlet邊界條件,也可以用于壁面。如果你希望保留默認(rèn)的零梯度條件,你就不必輸入任何東西了。如果你希望指定壁面處的組分質(zhì)量分?jǐn)?shù),步驟如下:1.在壁面面板的組分邊界條件中,選擇組分名字右邊的下拉列表指定的質(zhì)量分?jǐn)?shù)(而不是零梯度),此時面板會擴展為包含組分質(zhì)量分?jǐn)?shù)的對話框。Figure 1: 組分邊界條件輸入的壁面面板2.輸入相應(yīng)的組分質(zhì)量分?jǐn)?shù)。每一組分的邊界條件類型是分別指定的,所以對于不同的組分你可以采用不同的方法。注意:如果在

2、湍流流動中你使用組分的Dirichle條件,F(xiàn)LUENT就不會是用壁面函數(shù)來計算壁面處的組分?jǐn)U散流量。定義壁面的反應(yīng)邊界條件如果你在組分模型面板中激活了表面反應(yīng)的模擬,你就可以表明在壁面處表面反應(yīng)是否被激活。激活或關(guān)閉表面反應(yīng),壁面面板就會相應(yīng)地打開或關(guān)閉表面反應(yīng)選項。注意:組分在壁面處是假定為零梯度條件的,它不參加任何表面反應(yīng)。定義壁面的輻射邊界條件如果你打算使用P-1輻射模型、DTRM或者DO模型,你就需要設(shè)定壁面的(內(nèi)部)發(fā)散率以及(可選)黑體溫度。詳情請參閱設(shè)定邊界條件一節(jié)(Rosseland不需要任何邊界條件的輸入,因為FLUENT假定發(fā)射率為1,如果你使用DO模型你也要定義壁面為漫

3、反射、鏡面反射或者半透明,詳情請參閱設(shè)定輻射邊界條件)定義壁面的離散相邊界條件如果你是在模擬粒子的離散相,你就可以在壁面處設(shè)定粒子軌道詳情請參閱離散向模型的邊界設(shè)定。壁面邊界的默認(rèn)設(shè)定默認(rèn)熱邊界條件為固定的熱流為零,壁面默認(rèn)為不移動。壁面處的剪應(yīng)力計算程序?qū)τ诜腔票诿鏃l件,F(xiàn)LUENT使用鄰近壁面或者流體邊界的流動性質(zhì)來預(yù)測壁面處流體的剪應(yīng)力。在層流流動制,這一計算簡單地依賴于壁面處的速度梯度,在湍流流動中則使用壁面限制湍流流動的近壁面處理方法。對于指定剪切的壁面, FLUENT會在邊界處計算切向速度。如果是無粘流動,所有的壁面都使用滑移條件,所以它們是無摩擦的而且對鄰近流體單元不施加剪應(yīng)力

4、。層流中的剪應(yīng)力計算在層流流動中壁面剪應(yīng)力和法向速度梯度的關(guān)系為:當(dāng)壁面處的速度梯度很大時,你必須保證網(wǎng)格足夠精細(xì),這樣才能解出邊界層的精確結(jié)果。層流流動中近壁面節(jié)點放置的指導(dǎo)方針在節(jié)點密度和節(jié)點束中介紹。湍流中的剪應(yīng)力計算湍流流動的壁面處理,在壁面限制的湍流流動的近壁面處理一節(jié)中敘述。壁面邊界的熱傳導(dǎo)計算溫度邊界條件當(dāng)在壁面處應(yīng)用固定溫度條件,從流體單元到壁面的熱傳導(dǎo),由下式計算:其中:h_f=流體邊界當(dāng)?shù)責(zé)醾鲗?dǎo)系數(shù)T_w=壁面表面溫度T_f=當(dāng)?shù)亓黧w溫度q=壁面處傳來的對流熱流量q_rad=輻射熱流量注意:流體邊界熱傳導(dǎo)系數(shù)是基于當(dāng)?shù)亓鲌鰲l件計算得來的(比如說湍流層次、溫度以及速度輪廓),

5、請參閱流體邊界熱傳導(dǎo)計算一節(jié)的方程1,以及標(biāo)準(zhǔn)壁面函數(shù)9。從固體單元到壁面邊界的熱傳導(dǎo)公式為: 其中:k_s=固體的熱傳導(dǎo)率T_s=當(dāng)?shù)毓腆w溫度D n=壁面表面和固體單元中心的距離。熱流邊界條件當(dāng)你在壁面處定義熱流邊界條件時,你需要在壁面表面指定熱流量。FLUENT使用溫度邊界條件中的方程1,然后你就可以輸入熱流量來確定鄰近流體單元的壁面表面溫度:其中,流體邊界熱傳導(dǎo)系數(shù)已經(jīng)在溫度邊界條件中敘述了,它是基于當(dāng)?shù)亓鲌鰲l件計算得到的。當(dāng)壁面和固體區(qū)域交界時,壁面表面的溫度為:上述兩式的變量請參閱溫度邊界條件一節(jié)。對流熱傳導(dǎo)邊界條件當(dāng)你在壁面處指定對流熱傳導(dǎo)系數(shù)作為邊界條件時,F(xiàn)LUENT使用你所輸

6、入的外部熱傳導(dǎo)系數(shù)以及外部熱沉(heat sink)溫度來計算到壁面的熱流量:其中:h_ext=你所定義的外部熱傳導(dǎo)系數(shù)T_ext=你所定義的外部熱沉溫度q_rad=輻射熱流量上述方程假定壁面零厚度。外部輻射邊界條件當(dāng)使用外部輻射條件時,流入壁面的熱流量為:其中:e_ext=你所定義的外部壁面表面的發(fā)射率s=Stefan-Boltzmann常數(shù)T_w=壁面的表面溫度T_?=區(qū)域外部的溫度的輻射源或者消失(sink)處q_rad=從內(nèi)部去向壁面輻射的熱流量Equation 1假定壁面厚度為零。外部對流和輻射結(jié)合的邊界條件當(dāng)你選擇組合的外部熱傳導(dǎo)方程條件時,到壁面的熱流量為:其中的變量已經(jīng)在對流熱

7、傳導(dǎo)邊界條件和外部輻射邊界條件中定義了。Equation 1假定壁面厚度為零。流動邊界熱傳導(dǎo)系數(shù)的計算在層流流動中,壁面處流體邊界熱傳導(dǎo)是用應(yīng)用于壁面的Fourier定律計算得到的,F(xiàn)LUENT使用它的離散格式為: 其中n是垂直于壁面的當(dāng)?shù)刈鴺?biāo)。對于湍流流動,F(xiàn)LUENT對于從熱和動量遷移中類比得到的溫度使用壁面定律93。詳細(xì)內(nèi)容請參閱標(biāo)準(zhǔn)壁面函數(shù)。對稱邊界條件對稱邊界條件用于所計算的物理外形以及所期望的流動/熱解具有鏡像對稱的特征的情況中。也可以用它們來模擬粘性流動的滑移壁面。本節(jié)描述了對稱平面內(nèi)流動的處理,并提供了一些使用對稱邊界的例子。在對稱邊界條件中你不需要定義任何邊界條件,但是你必須

8、謹(jǐn)慎地定義對稱邊界的位置。在對稱外形的中線處,你應(yīng)該使用軸邊界類型而不是對稱邊界類型,如軸邊界條件一節(jié)中的的圖1,詳細(xì)內(nèi)容請參閱軸邊界條件。對稱邊界的計算程序FLUENT假定所有量通過對稱邊界的流量為零。經(jīng)過對稱平面的對流流量為零,因此對稱邊界的法向速度為零。通過對稱平面沒有擴散流量:因此所有流動變量的法向梯度在對稱平面內(nèi)為零。因此對稱邊界條件可以總結(jié)如下:l 對稱平面內(nèi)法向速度為零l 對稱平面內(nèi)所有變量的法向梯度為零如上所述,對稱的定義要求這些條件決定流過對稱平面的流量為零。因為對稱邊界的剪應(yīng)力為零,所以在粘性流動計算中它也可以用滑移壁面來解釋。對稱邊界的例子對稱邊界用于減少計算模擬的范圍,

9、它只需要模擬所有物理系統(tǒng)的一個對稱子集。下面兩個圖是通過該種方法使用對稱邊界的例子。Figure 1:使用對稱邊界模擬三維管道的四分之一Figure 2: 使用對稱邊界模擬圓形截面的四分之一下面的圖則是誤用對稱平面的兩個例子,在這兩個例子中,雖然幾何外形是對稱的,但是流動本身卻不符合對稱邊界條件的要求。在第一個例子中浮力產(chǎn)生了非對稱流動。在第二個例子中,流動中的渦流產(chǎn)生了一個垂直于應(yīng)該是對稱平面的流動。需要注意的是,這兩個粒子都要使用旋轉(zhuǎn)周期性邊界(請參閱周期性邊界一節(jié)的圖一)Figure 3: 對稱的誤用周期性邊界條件周期性邊界條件用來解決,物理模型和所期待的流動的流動/熱解具有周期性重復(fù)的

10、特點。FLUENT提供了兩種類型的周期性邊界條件。第一種類型不允許通過周期性平面具有壓降(對于FLUENT4用戶來說:這一類型的周期性邊界是指FLUENT4中的圓柱形邊界)。第二種類型允許通過平移周期性邊界具有壓降,它是你能夠模擬完全發(fā)展的周期性流動(在FLUENT4中是周期性邊界)。本節(jié)討論了無壓降的周期性邊界條件。在周期性流動和熱傳導(dǎo)一節(jié)中,完全發(fā)展的周期性模擬能力得到了詳盡的描述。周期性邊界的例子周期性邊界條件用于模擬通過計算模型內(nèi)的兩個相反平面的流動是相同的情況。下圖是周期性邊界條件的典型應(yīng)用。在這些例子中,通過周期性平面進入計算模型的流動和通過相反的周期性平面流出流場的流動是相同的。

11、正如這些例子所示,周期性平面通常是成對使用的。Figure 1: 在圓柱容器中使用周期性邊界定義渦流周期性邊界的輸入對于沒有任何壓降的周期性邊界,你只需要輸入一個東西,那就是你的所模擬的幾何外形是旋轉(zhuǎn)性周期還是平移性周期。(對于有周期性壓降的周期流還要輸入其它的東西,請參閱周期性流動和熱傳導(dǎo)一節(jié)。)旋轉(zhuǎn)性周期邊界是指關(guān)于旋轉(zhuǎn)對稱幾何外形中線形成了一個包括的角度。本節(jié)中的圖一就是旋轉(zhuǎn)性周期。平移性周期邊界是指在直線幾何外形內(nèi)形成周期性邊界。下面兩圖是平移性周期邊界:Figure 1: 物理區(qū)域Figure 2: 所模擬的區(qū)域?qū)τ谥芷谛赃吔纾阈枰谥芷谛悦姘澹ㄏ聢D)中指定平移性邊界還是旋轉(zhuǎn)性邊界

12、,該面板是從設(shè)定邊界條件菜單中打開的。Figure 3: 周期性面板(對于耦合解算器,周期性面板中將會有附加的選項,這一選項允許你指定壓力跳躍,詳細(xì)內(nèi)容請參閱周期性流動和熱傳導(dǎo)一節(jié)。)如果區(qū)域是旋轉(zhuǎn)性區(qū)域,請選擇旋轉(zhuǎn)性區(qū)域類型。如果是平移性就選擇平移性區(qū)域類型。對于旋轉(zhuǎn)性區(qū)域,解算器會自動計算通過周期性區(qū)域的旋轉(zhuǎn)角度。旋轉(zhuǎn)軸是為鄰近單元指定的旋轉(zhuǎn)軸。注意:對于使用旋轉(zhuǎn)周期性邊界來說,你不必指定鄰近單元區(qū)域為移動的。例如,你能夠使用具有管的平切片的非旋轉(zhuǎn)坐標(biāo)系來模擬三維管流,管的切片需要具有旋轉(zhuǎn)性周期。你可以使用Grid/Check菜單選項(參閱檢查網(wǎng)格一節(jié))來計算和顯示周期性邊界所有表面的旋轉(zhuǎn)

13、角度的最大值、最小值和平均值。如果最大值、最小值和平均值之間的差別可以忽略 ,那么網(wǎng)格有一個問題:對于指定軸來說網(wǎng)格幾何外形不是周期性的。周期性邊界的默認(rèn)設(shè)定默認(rèn)為平移周期性邊界條件周期性邊界的計算程序FLUENT在周期性邊界處理流動就像反向周期性平面是和前面的周期性邊界直接相鄰一樣,因此,當(dāng)計算流過鄰近流體單元的周期性邊界時,就會使用與反向周期性平面相鄰的流體單元的流動條件。軸邊界的計算程序軸邊界條件軸邊界類型必須使用在對稱幾何外形的中線處(見下圖)。它也可以用在圓柱兩極的四邊形和六面體網(wǎng)格的中線上(比如:像FLUENT4之類的結(jié)構(gòu)網(wǎng)格生成代碼所產(chǎn)生的網(wǎng)格)。在軸邊界處,你不必定義任何邊界條

14、件。Figure 1: 在軸對稱幾何外形的中線處軸邊界條件的使用軸邊界的計算程序要確定軸上特定點的適當(dāng)物理值,F(xiàn)LUENT使用鄰近單元中的單元值。流體條件流體區(qū)域是一組所有現(xiàn)行的方程都被解出的單元。對于流體區(qū)域只需要輸入流體材料類型。你必須指明流體區(qū)域內(nèi)包含哪種材料,以便于使用適當(dāng)?shù)牟牧蠈傩?。如果你模擬組分輸運或者燃燒,你就不必在這里選擇材料屬性,當(dāng)你激活模型時,組分模型面板中會指定混合材料。相似地,對于多相流動你也不必指定材料屬性,當(dāng)你在多相流模型面板中激活模型時,你會選擇它們。可選擇的輸入允許你設(shè)定熱、質(zhì)量、動量、湍流、組分以及其它標(biāo)量屬性的源項。你也可以為流體區(qū)域定義運動。如果鄰近流體區(qū)

15、域內(nèi)具有旋轉(zhuǎn)周期性邊界,你就需要指定旋轉(zhuǎn)軸。如果你使用k-e模型或者Spalart-Allmaras模型來模擬湍流,你可以選擇定義流體區(qū)域為層流區(qū)域。如果你用DO模型模擬輻射,你可以指定流體是否參加輻射。對于多孔區(qū)域的信息,請參閱多孔介質(zhì)條件一節(jié)。流體區(qū)域的輸入在流體面板中(下圖),你需要設(shè)定所有的流體條件,該面板是從設(shè)定邊界條件菜單中打開的。Figure 1: 流體面板定義流體材料要定義流體區(qū)域內(nèi)包含的材料,請在材料名字下拉列表中選擇適當(dāng)?shù)倪x項。這一列表中會包含所有已經(jīng)在使用材料面板中定義的流體材料(或者從材料數(shù)據(jù)庫中加載)。如果你模擬組分輸運或者多相流,在流體面板的下拉列表中不會出現(xiàn)材料名

16、。對于組分計算,所有流體區(qū)域的混合材料將會是你在組分模型面板中所指定的材料。對于多相流,所有流體區(qū)域的材料將會是你在多相流模型面板中所指定的材料。定義源項如果你希望在流體區(qū)域內(nèi)定義熱、質(zhì)量、動量、湍流、組分以及其它標(biāo)量屬性的源項,你可以激活源項選項來實現(xiàn)。詳情請參閱定義質(zhì)量、動量、能量和其它源項一節(jié)。指定層流區(qū)域如果你使用k-e模型或者Spalart-Allmaras模型來模擬湍流,在指定的流體區(qū)域關(guān)掉湍流模擬是可能的(即:使湍流生成和湍流粘性無效,但是湍流性質(zhì)的輸運仍然保持)。如果你知道在某一區(qū)域流動是層流這一功能是很有用的。比方說:如果你知道機翼上的轉(zhuǎn)唳點的位置,你可以在層流單元區(qū)域邊界和

17、湍流區(qū)域邊界創(chuàng)建一個層流/湍流過渡邊界。這一功能允許你模擬機翼上的湍流過渡。要在流體區(qū)域內(nèi)取消湍流模擬,請在流體面板中打開層流區(qū)域選項。 指定旋轉(zhuǎn)軸如果鄰近流體區(qū)域存在旋轉(zhuǎn)性周期邊界,或者區(qū)域是旋轉(zhuǎn)的,你必須指定旋轉(zhuǎn)軸。要定義旋轉(zhuǎn)軸,請設(shè)定旋轉(zhuǎn)軸方向和起點。這個軸和任何鄰近壁面區(qū)域或任何其它單元區(qū)域所使用的旋轉(zhuǎn)軸是獨立的。對于三維問題,旋轉(zhuǎn)軸起點是從旋轉(zhuǎn)軸起點中輸入的起點,方向為旋轉(zhuǎn)軸方向選項中輸入的方向。對于二維非軸對稱問題,你只需要指定旋轉(zhuǎn)軸起點,方向就是通過指定點的z方向。(z向是垂直于幾何外形平面的,這樣才能保證旋轉(zhuǎn)出現(xiàn)在該平面內(nèi))。對于二維軸對稱問題,你不必定義軸,旋轉(zhuǎn)通常就是關(guān)于x

18、軸的,起點為(0,0)。定義區(qū)域運動對于旋轉(zhuǎn)和平移坐標(biāo)系要定義移動區(qū)域,請在運動類型下菜單(如果你用滾動條向右滾動到旋轉(zhuǎn)軸起點和方向,就是可見的了)中選擇運動參考坐標(biāo)系。然后在面板的擴展部分設(shè)定適當(dāng)?shù)膮?shù)。要對移動或者滑移網(wǎng)格定義移動區(qū)域,在移動類型下拉列表中選擇移動網(wǎng)格,然后在擴展面板中設(shè)定適當(dāng)?shù)膮?shù)。詳情請參閱滑動網(wǎng)格。對于包括線性、平移運動的流體區(qū)域問題,通過設(shè)定X, Y,和Z分量來指定平移速度。對于包括旋轉(zhuǎn)運動的問題,在旋轉(zhuǎn)速度中指定旋轉(zhuǎn)速度。旋轉(zhuǎn)軸的定義請參閱指定旋轉(zhuǎn)軸一節(jié)。關(guān)于在移動參考系中模擬流動的詳細(xì)內(nèi)容請參閱移動區(qū)域的流動一節(jié)。定義輻射參數(shù)如果你使用DO輻射模型,你可以用參加

19、輻射選項指定流體區(qū)域是否參加輻射的計算。詳情請參閱輻射邊界條件一節(jié)。固體條件固體區(qū)域是僅用來解決熱傳導(dǎo)問題的一組區(qū)域。作為固體處理的材料可能事實上是流體,但是假定其中沒有對流發(fā)生。固體區(qū)域僅需要輸入材料類型。你必須表明固體區(qū)域包含哪種材料,以便于計算是使用適當(dāng)?shù)牟牧?。可選擇的輸入允許你設(shè)定體積熱生成速度(熱源)。你也可以定義固體區(qū)域的運動。如果在鄰近的固體單元內(nèi)有旋轉(zhuǎn)性周期邊界,你就需要指定旋轉(zhuǎn)軸。如果你模擬DO輻射模型,你可以指定固體材料是否參加輻射的計算。固體區(qū)域的輸入流體區(qū)域的輸入在固體面板中(下圖),你需要設(shè)定所有的固體條件,該面板是從設(shè)定邊界條件菜單中打開的。Figure 1: 固體

20、面板定義流體材料要定義固體區(qū)域內(nèi)包含的材料,請在材料名字下拉列表中選擇適當(dāng)?shù)倪x項。這一列表中會包含所有已經(jīng)在使用材料面板中定義的固體材料(或者從材料數(shù)據(jù)庫中加載)。定義熱源如果你希望在固體區(qū)域內(nèi)定義熱源項,你可以激活源項選項來實現(xiàn)。詳情請參閱定義質(zhì)量、動量、能量和其它源項一節(jié)。指定旋轉(zhuǎn)軸如果鄰近固體區(qū)域存在旋轉(zhuǎn)性周期邊界,或者區(qū)域是旋轉(zhuǎn)的,你必須指定旋轉(zhuǎn)軸。要定義旋轉(zhuǎn)軸,請設(shè)定旋轉(zhuǎn)軸方向和起點。這個軸和任何鄰近壁面區(qū)域或任何其它單元區(qū)域所使用的旋轉(zhuǎn)軸是獨立的。對于三維問題,旋轉(zhuǎn)軸起點是從旋轉(zhuǎn)軸起點中輸入的起點,方向為旋轉(zhuǎn)軸方向選項中輸入的方向。對于二維非軸對稱問題,你只需要指定旋轉(zhuǎn)軸起點,方向

21、就是通過指定點的z方向。(z向是垂直于幾何外形平面的,這樣才能保證旋轉(zhuǎn)出現(xiàn)在該平面內(nèi))。對于二維軸對稱問題,你不必定義軸,旋轉(zhuǎn)通常就是關(guān)于x軸的,起點為(0,0)。定義區(qū)域運動對于旋轉(zhuǎn)和平移坐標(biāo)系要定義移動區(qū)域,請在運動類型下菜單(如果你用滾動條向右滾動到旋轉(zhuǎn)軸起點和方向,就是可見的了)中選擇運動參考坐標(biāo)系。然后在面板的擴展部分設(shè)定適當(dāng)?shù)膮?shù)。要對移動或者滑移網(wǎng)格定義移動區(qū)域,在移動類型下拉列表中選擇移動網(wǎng)格,然后在擴展面板中設(shè)定適當(dāng)?shù)膮?shù)。詳情請參閱滑動網(wǎng)格。對于包括線性、平移運動的流體(?原文是流體,按理說應(yīng)該是固體)區(qū)域問題,通過設(shè)定X, Y,和Z分量來指定平移速度。對于包括旋轉(zhuǎn)運動的問

22、題,在旋轉(zhuǎn)速度中指定旋轉(zhuǎn)速度。旋轉(zhuǎn)軸的定義請參閱指定旋轉(zhuǎn)軸一節(jié)。關(guān)于在移動參考系中模擬流動的詳細(xì)內(nèi)容請參閱移動區(qū)域的流動一節(jié)。定義輻射參數(shù)如果你使用DO輻射模型,你可以用參加輻射選項指定固體區(qū)域是否參加輻射的計算。詳情請參閱輻射邊界條件一節(jié)。多孔介質(zhì)條件多孔介質(zhì)模型可以應(yīng)用于很多問題,如通過充滿介質(zhì)的流動、通過過濾紙、穿孔圓盤、流量分配器以及管道堆的流動。當(dāng)你使用這一模型時,你就定義了一個具有多孔介質(zhì)的單元區(qū)域,而且流動的壓力損失由多孔介質(zhì)的動量方程中所輸入的內(nèi)容來決定。通過介質(zhì)的熱傳導(dǎo)問題也可以得到描述,它服從介質(zhì)和流體流動之間的熱平衡假設(shè),具體內(nèi)容可以參考多孔介質(zhì)中能量方程的處理一節(jié)。多孔

23、介質(zhì)的一維化簡模型,被稱為多孔跳躍,可用于模擬具有已知速度/壓降特征的薄膜。多孔跳躍模型應(yīng)用于表面區(qū)域而不是單元區(qū)域,并且在盡可能的情況下被使用(而不是完全的多孔介質(zhì)模型),這是因為它具有更好的魯棒性,并具有更好的收斂性。詳細(xì)內(nèi)容請參閱多孔跳躍邊界條件。多孔介質(zhì)模型的限制如下面各節(jié)所述,多孔介質(zhì)模型結(jié)合模型區(qū)域所具有的阻力的經(jīng)驗公式被定義為“多孔”。事實上多孔介質(zhì)不過是在動量方程中具有了附加的動量損失而已。因此,下面模型的限制就可以很容易的理解了。l 流體通過介質(zhì)時不會加速,因為事實上出現(xiàn)的體積的阻塞并沒有在模型中出現(xiàn)。這對于過渡流是有很大的影響的,因為它意味著FLUENT不會正確的描述通過介

24、質(zhì)的過渡時間。l 多孔介質(zhì)對于湍流的影響只是近似的。詳細(xì)內(nèi)容可以參閱湍流多孔介質(zhì)的處理一節(jié)。多孔介質(zhì)的動量方程多孔介質(zhì)的動量方程具有附加的動量源項。源項由兩部分組成,一部分是粘性損失項 (Darcy),另一個是內(nèi)部損失項:其中S_i是i向(x, y, or z)動量源項,D和C是規(guī)定的矩陣。在多孔介質(zhì)單元中,動量損失對于壓力梯度有貢獻,壓降和流體速度(或速度方陣)成比例。對于簡單的均勻多孔介質(zhì):其中a是滲透性,C_2時內(nèi)部阻力因子,簡單的指定D和C分別為對角陣1/a 和C_2其它項為零。FLUENT還允許模擬的源項為速度的冪率:其中C_0和C_1為自定義經(jīng)驗系數(shù)。注意:在冪律模型中,壓降是各向

25、同性的,C_0的單位為國際標(biāo)準(zhǔn)單位。多孔介質(zhì)的Darcy定律通過多孔介質(zhì)的層流流動中,壓降和速度成比例,常數(shù)C_2可以考慮為零。忽略對流加速以及擴散,多孔介質(zhì)模型簡化為Darcy定律:在多孔介質(zhì)區(qū)域三個坐標(biāo)方向的壓降為:其中1/a_ij為多孔介質(zhì)動量方程1中矩陣D的元素v_j為三個方向上的分速度,D n_x、 D n_y、以及D n_z為三個方向上的介質(zhì)厚度。在這里介質(zhì)厚度其實就是模型區(qū)域內(nèi)的多孔區(qū)域的厚度。因此如果模型的厚度和實際厚度不同,你必須調(diào)節(jié)1/a_ij的輸入。.多孔介質(zhì)的內(nèi)部損失在高速流動中,多孔介質(zhì)動量方程1中的常數(shù)C_2提供了多孔介質(zhì)內(nèi)部損失的矯正。這一常數(shù)可以看成沿著流動方向

26、每一單位長度的損失系數(shù),因此允許壓降指定為動壓頭的函數(shù)。如果你模擬的是穿孔板或者管道堆,有時你可以消除滲透項而只是用內(nèi)部損失項,從而得到下面的多孔介質(zhì)簡化方程: 寫成坐標(biāo)形式為: 多孔介質(zhì)中能量方程的處理對于多孔介質(zhì)流動,F(xiàn)LUENT仍然解標(biāo)準(zhǔn)能量輸運方程,只是修改了傳導(dǎo)流量和過度項。在多孔介質(zhì)中,傳導(dǎo)流量使用有效傳導(dǎo)系數(shù),過渡項包括了介質(zhì)固體區(qū)域的熱慣量:其中:h_f=流體的焓 h_s=固體介質(zhì)的焓 f=介質(zhì)的多孔性 k_eff=介質(zhì)的有效熱傳導(dǎo)系數(shù) Sh_f=流體焓的源項 Sh_s=固體焓的源項多孔介質(zhì)的有效傳導(dǎo)率多孔區(qū)域的有效熱傳導(dǎo)率k_eff是由流體的熱傳導(dǎo)率和固體的熱傳導(dǎo)率的體積平均

27、值計算得到:其中:f=介質(zhì)的多孔性 k_f=流體狀態(tài)熱傳導(dǎo)率(包括湍流的貢獻k_t)k_s=固體介質(zhì)熱傳導(dǎo)率如果得不到簡單的體積平均,可能是因為介質(zhì)幾何外形的影響。有效傳導(dǎo)率可以用自定義函數(shù)來計算。然而,在所有的算例中,有效傳導(dǎo)率被看成介質(zhì)的各向同性性質(zhì)。多孔介質(zhì)中的湍流處理在多孔介質(zhì)中,默認(rèn)的情況下FLUENT會解湍流量的標(biāo)準(zhǔn)守恒防城。因此,在這種默認(rèn)的方法中,介質(zhì)中的湍流被這樣處理:固體介質(zhì)對湍流的生成和耗散速度沒有影響。如果介質(zhì)的滲透性足夠大,而且介質(zhì)的幾何尺度和湍流渦的尺度沒有相互作用,這樣的假設(shè)是合情合理的。但是在其它的一些例子中,你會壓制了介質(zhì)中湍流的影響。如果你使用k-e模型或者

28、Spalart-Allmaras模型,你如果設(shè)定湍流對粘性的貢獻m_t為零,你可能會壓制了湍流對介質(zhì)的影響。當(dāng)你選擇這一選項時,F(xiàn)LUENT會將入口湍流的性質(zhì)傳輸?shù)浇橘|(zhì)中,但是它對流動混合和動量的影響被忽略了。除此之外,在介質(zhì)中湍流的生成也被設(shè)定為零。要實現(xiàn)這一解策略,請在流體面板中打開層流選項 。激活這個選項就意味著多孔介質(zhì)中的m_t為零,湍流的生成也為零。如果去掉該選項(默認(rèn))則意味著多孔介質(zhì)中的湍流會像大體積流體流動一樣被計算。概述模擬多孔介質(zhì)流動時,對于問題設(shè)定需要的附加輸入如下:1. 定義多孔區(qū)域2. 確定流過多孔區(qū)域的流體材料3. 設(shè)定粘性系數(shù)(多孔介質(zhì)動量方程3中的1/a_ij)

29、以及內(nèi)部阻力系數(shù)(多孔介質(zhì)動量方程3中的C_2_ij),并定義應(yīng)用它們的方向矢量。冪率模型的系數(shù)也可以選擇指定。4. 定義多孔介質(zhì)包含的材料屬性和多孔性5. 設(shè)定多孔區(qū)域的固體部分的體積熱生成速度(或任何其它源項,如質(zhì)量、動量)(此項可選)。6. 如果合適的話,限制多孔區(qū)域的湍流粘性。7. 如果相關(guān)的話,指定旋轉(zhuǎn)軸和/或區(qū)域運動。在定義粘性和內(nèi)部阻力系數(shù)中描述了決定阻力系數(shù)和/或滲透性的方法。如果你使用多孔動量源項的冪律近似,你需要輸入多孔介質(zhì)動量方程5中的C_0和C_1來取代阻力系數(shù)和流動方向。在流體面板中(下圖)你需要設(shè)定多孔介質(zhì)的所有參數(shù),該面板是從邊界條件菜單中打開的(詳細(xì)內(nèi)容請參閱邊

30、界條件的設(shè)定一節(jié))Figure 1:多孔區(qū)域的流體面板定義多孔區(qū)域正如定義邊界條件概述中所提到的,多孔區(qū)域是作為特定類型的流體區(qū)域來模擬的。亞表明流體區(qū)域是多孔區(qū)域,請在流體面板中激活多孔區(qū)域選項。面板會自動擴展到多孔介質(zhì)輸入狀態(tài)。定義穿越多孔介質(zhì)的流體在材料名字下拉菜單中選擇適當(dāng)?shù)牧黧w就可以定義通過多孔介質(zhì)的流體了。如果你模擬組分輸運或者多相流,流體面板中就不會出現(xiàn)材料名字下拉菜單了。對于組分計算,所有流體和/或多孔區(qū)域的混合材料就是你在組分模型面板中指定的材料。對于多相流模型,所有流體和/或多孔區(qū)域的混合材料就是你在多相流模型面板中指定的材料。定義粘性和內(nèi)部阻力系數(shù)粘性和內(nèi)部阻力系數(shù)以相同

31、的方式定義。使用笛卡爾坐標(biāo)系定義系數(shù)的基本方法是在二維問題中定義一個方向矢量,在三維問題中定義兩個方向矢量,然后在每個方向上指定粘性和/或阻力系數(shù)。在二維問題中第二個方向沒有明確定義,它是垂直于指定的方向矢量和z向矢量所在的平面的。在三維問題中,第三個方向矢量是垂直于所指定的兩個方向矢量所在平面的。對于三維問題,第二個方向矢量必須垂直于第一個方向矢量。如果第二個方向矢量指定失敗,解算器會確保它們垂直而忽略在第一個方向上的第二個矢量的任何分量。所以你應(yīng)該確保第一個方向指定正確。在三維問題中也可能會使用圓錐(或圓柱)坐標(biāo)系來定義系數(shù),具體如下:定義阻力系數(shù)的過程如下:1. 定義方向矢量。l 使用笛

32、卡爾坐標(biāo)系,簡單指定方向1矢量,如果是三維問題,指定方向2矢量。每一個方向都應(yīng)該是從(0,0)或者(0,0,0)到指定的(X,Y)或(X,Y,Z)矢量。(如果方向不正確請按上面的方法解決)l 對于有些問題,多孔介質(zhì)的主軸和區(qū)域的坐標(biāo)軸不在一條直線上,你不必知道多孔介質(zhì)先前的方向矢量。在這種情況下,三維中的平面工具或者二維中的線工具可以幫你確定這些方向矢量。1. 捕捉Snap平面工具(或者線工具)到多孔區(qū)域的邊界。(請遵循使用面工具和線工具中的說明,它在已存在的表面上為工具初始化了位置)。2. 適當(dāng)?shù)男D(zhuǎn)坐標(biāo)軸直到它們和多孔介質(zhì)區(qū)域成一條線。3. 當(dāng)成一條線之后,在流體面板中點擊從平面工具更新或

33、者從線工具更新按鈕。FLUENT會自動將方向1矢量指向為工具的紅(三維)或綠(二維)箭頭所指的方向。l 要使用圓錐坐標(biāo)系(比方說環(huán)狀、錐狀顧慮單元),請遵循下面步驟(這一選項只用于三維問題):1. 打開圓錐選項2. 指定圓錐軸矢量和在錐軸上的點。圓錐軸矢量的方向?qū)菑?0,0,0)到指定的(X,Y,Z)方向的矢量。FLUENT將會使用圓錐軸上的點將阻力轉(zhuǎn)換到笛卡爾坐標(biāo)系。3. 設(shè)定錐半角(錐軸和錐表面之間的角度,如下圖),使用柱坐標(biāo)系,錐半角為0.Figure 1:錐半角l 對于有些問題,錐形過濾單元的主軸和區(qū)域的坐標(biāo)軸不在一條直線上,你不必知道錐軸先前的方向矢量以及錐軸上的點。在這種情況下

34、,三維中的平面工具或者二維中的線工具可以幫你確定這些方向矢量。一種方法如下:1. 在點擊捕捉到區(qū)域按鈕之前,你可以在下拉菜單中選擇垂直于錐軸矢量的軸過濾單元的邊界區(qū)域。2. 點擊捕捉到區(qū)域按鈕,F(xiàn)LUENT會自動將平面工具捕捉到邊界。它也會設(shè)定錐軸矢量和錐軸上的點(需注意的是你還要自己設(shè)定錐半角)。l 另一種方法為: 1. 捕捉Snap平面工具到多孔區(qū)域的邊界。(請遵循使用面工具和線工具中的說明,它在已存在的表面上為工具初始化了位置)。2. 旋轉(zhuǎn)和平移工具坐標(biāo)軸,直到工具的紅箭頭指向錐的軸向。工具的起點在軸上。3. 當(dāng)軸和工具的起點成一條線時,在流體面板中點擊從平面工具更新按鈕。FLUENT會

35、自動設(shè)定軸向矢量以及在軸上的點(注意:你還是要自己設(shè)定錐的半角)。2. 在粘性阻力中指定每個方向的粘性阻力系數(shù)1/a,在內(nèi)部阻力中指定每一個方向上的內(nèi)部阻力系數(shù)C_2(你可能需要將滾動條向下滾動來查看這些輸入)。如果你使用錐指定方法,方向1為錐軸方向,方向2為垂直于錐表面(對于圓柱就是徑向)方向,方向3圓周(q)方向。在三維問題中可能有三種可能的系數(shù),在二維問題中有兩種:l 在各向同性算例中,所有方向上的阻力系數(shù)都是相等的(如海綿)。在各向同性算例中你必須將每個方向上的阻力系數(shù)設(shè)定為相等。l 在三維問題中只有兩個方向上的系數(shù)相等,第三個方向上的阻力系數(shù)和前兩個不等,或者在二維問題中兩個方向上的

36、系數(shù)不等,你必須準(zhǔn)確的指定每一個方向上的系數(shù)。例如,如果你得多孔區(qū)域是由具有小洞的細(xì)管組成,細(xì)管平行于流動方向,流動會很容易的通過細(xì)管,但是流動在其它兩個方向上(通過小洞)會很小。如果你有一個平的盤子垂直于流動方向,流動根本就不會穿過它而只在其它兩個方向上。l 在三維問題中還有一種可能就是三個系數(shù)各不相同。例如,如果多孔區(qū)域是由不規(guī)則間隔的物體(如針腳)組成的平面,那么阻礙物之間的流動在每個方向上都不同。此時你就需要在每個方向上指定不同的系數(shù)(請注意指定各向同性系數(shù)時,多孔介質(zhì)的解策略的注解)。推導(dǎo)粘性和內(nèi)部損失系數(shù)的方法在定義粘性和內(nèi)部阻力系數(shù)一節(jié)中介紹。當(dāng)你使用多孔介質(zhì)模型時,你必須記住F

37、LUENT中的多孔單元是100%打開的,而且你所指定1/a_ij和/或C_2_ij的值必須是基于這個假設(shè)的。然而,假如你知道通過真實裝置壓降和速度之間的的變化,它只是部分地對流動開放。下面的練習(xí)會告訴你如何對FLUENT模型計算適當(dāng)?shù)腃_2值。假定穿孔圓盤只有25%對流動開放。已知通過圓盤的壓降為0.5。在圓盤內(nèi)真實流體速度基礎(chǔ)上,即通過%25開放區(qū)域的的基礎(chǔ)上,損失系數(shù)由下式定義的損失系數(shù)K_L為0.5:要計算適當(dāng)?shù)腃_2值,請注意在FLUENT模型中:1. 通過穿孔圓盤的速度假定圓盤為100%開放的。2. 損失系數(shù)必須轉(zhuǎn)化為多孔區(qū)域每個單位長度的動壓頭損失。對于第一條,第一步是計算并調(diào)節(jié)損

38、失因子K_L,它應(yīng)該是在100%開放區(qū)域的速度基礎(chǔ)上的:或者注意對于相同的流速,v_25% open = 4 v_100% open,調(diào)節(jié)之后的損失系數(shù)為8。對于第二條,你必須將它轉(zhuǎn)換為穿孔圓盤每個單位厚度的損失系數(shù)。假定圓盤的厚度為1.0 mm。內(nèi)部損失系數(shù)為(國際標(biāo)準(zhǔn)單位):注意,對于各向異性介質(zhì),這些信息必須分別從每一個坐標(biāo)方向上計算。第二個例子,考慮模擬充滿介質(zhì)的流動。在湍流流動中,充滿介質(zhì)的流動用滲透性和內(nèi)部損失系數(shù)來模擬。推導(dǎo)適當(dāng)常數(shù)的方法包括了Ergun方程49的使用,對于在很大范圍雷諾數(shù)內(nèi)和許多類型的充滿形式,有一個半經(jīng)驗的關(guān)系式:當(dāng)模擬充滿介質(zhì)的層流流動時,上面方程中的第二項

39、可能是個小量,從而得到Blake-Kozeny方程49:在這些方程中,m是粘性,D_p是平均粒子直徑,e空間所占的分?jǐn)?shù)(即空間的體積除以總體積)。比較多孔介質(zhì)中Darcy定律的方程1和內(nèi)部損失系數(shù)為9的方程1,則每一方向上的滲透性和內(nèi)部損失系數(shù)定義為:第三個例子我們會考慮Van Winkle等人146,121的方程,并表明如何通過具有方孔圓盤的多孔介質(zhì)輸入來計算壓力損失。作者所聲明的應(yīng)用在通過在等邊三角形上的方洞圓盤的湍流中的表達式為:其中:m(dot)=通過圓盤的質(zhì)量流速A_f=剩下的面積或者洞的總面積A_p=圓盤的面積(固體和洞) C=對于不同D/t的不同雷諾數(shù)范圍被列成不同的表的系數(shù) D

40、/t=洞的直徑和圓盤厚度的比例對于t/D 1.6和Re 4000,系數(shù)C近似為0.98,其中雷諾數(shù)是基于洞的直徑與速度的使用下式整理方程17:除以圓盤的厚度D x = t有:其中v是表面速度而不是洞內(nèi)的速度。與多孔介質(zhì)內(nèi)部損失系數(shù)中的方程1比較可以看出,對于垂直于圓盤方向,常數(shù)C_2可由下式計算:考慮通過由隨機方向的纖維或者玻璃材料組成的墊子或者過濾器的層流。對于可以二選一的方程Blake-Kozeny(方程11),我們可能會選擇將實驗數(shù)據(jù)列成表。很多類型的纖維都由這一類相關(guān)的數(shù)據(jù)70。固體體積分?jǐn)?shù)f 玻璃絲織品的無量綱滲透性Q0.262 0.25 0.258 0.26 0.221 0.40

41、0.218 0.41 0.172 0.80 其中Q =,a為纖維直徑。使用多孔介質(zhì)的Darcy定律中的方程1可以很容易從給定的纖維直徑和體積分?jǐn)?shù)種計算出。使用冪律模型對于多孔介質(zhì)動量源項(多孔介質(zhì)動量方程中的方程5),如果你使用冪律模型近似,你只要在流體面板的冪律模型中輸入系數(shù)C_0和C_1就可以了。如果C_0或C_1為非零值,解算器會忽略面板中除了多孔介質(zhì)冪律模型之外的所有輸入。定義熱傳導(dǎo)如果你選擇在多孔介質(zhì)中模擬熱傳導(dǎo),你必須指定多孔介質(zhì)中的材料以及多孔性。要定義多孔介質(zhì)的材料,向下拉流體面板中阻力輸入下面的滾動條,然后在多孔熱傳導(dǎo)的固體材料下拉列表中選中適當(dāng)?shù)墓腆w。然后在多孔熱傳導(dǎo)下設(shè)定

42、多孔性。多孔性f是多孔介質(zhì)中流體的體積分?jǐn)?shù)(即介質(zhì)的開放體積分?jǐn)?shù))。多孔性用于介質(zhì)中的熱傳導(dǎo)預(yù)測,處理方法請參閱多孔介質(zhì)能量方程的處理一節(jié)。它還對介質(zhì)中的反應(yīng)源項和體力的計算有影響。這個源項和介質(zhì)中流體的體積成比例。如果你想要模擬完全開放的介質(zhì)(固體介質(zhì)沒有影響),你應(yīng)該設(shè)定多孔性為1.0。當(dāng)多孔性為1.0時,介質(zhì)的固體部分對于熱傳導(dǎo)和(或)熱源項/反應(yīng)源項沒有影響。注意:多孔性永遠不會影響介質(zhì)中的流體速度,這已經(jīng)在多孔介質(zhì)的動量方程一節(jié)中介紹了。不管你將多孔性設(shè)定為何值,F(xiàn)LUENT所預(yù)測的速度都是介質(zhì)中的表面速度。定義源項如果你想在多孔流動的能量方程中包括熱的影響,請激活源項選項并設(shè)定非零

43、的能量源項。FLUENT會計算多孔區(qū)域所生成的能量,該能量為能量源項值乘以組成多孔區(qū)域的單元所有體積值。你也可以定義質(zhì)量、動量、湍流、組分或者其它標(biāo)量的源項,詳細(xì)內(nèi)容請參閱、質(zhì)量、動量、能量和其它源項的定義。在多孔區(qū)域內(nèi)壓制湍流源項如多孔介質(zhì)的湍流處理中所討論的,湍流在多孔介質(zhì)中的計算和大量(bulk)流體流動是一樣的。如果你使用k-e模型或者Spalart-Allmaras模型,你想要壓制湍流在多孔區(qū)域的影響可以打開流體區(qū)域面板中的層流區(qū)域選項(從而使得多孔區(qū)域的湍流生成為零)。指定旋轉(zhuǎn)軸并定義區(qū)域運動旋轉(zhuǎn)軸和區(qū)域運動的輸入和標(biāo)準(zhǔn)流體區(qū)域的輸入是相同的,詳細(xì)情況可以參閱流體區(qū)域的輸入一節(jié)。多

44、孔介質(zhì)的解策略一般說來,在模擬多孔介質(zhì)時,你可以使用標(biāo)準(zhǔn)的解算步驟以及解參數(shù)的設(shè)置。然而你會發(fā)現(xiàn)如果多孔區(qū)域在流動方向上壓降相當(dāng)大(比如:滲透性a很低或者內(nèi)部因子C_2很大)的話,解的收斂速度就會變慢。這就表明由于動量源項中出現(xiàn)了多孔介質(zhì)的壓降(方程的矩陣不再是對角占優(yōu)了),收斂性問題就出現(xiàn)了。解決多孔介質(zhì)區(qū)域收斂性差最好的補救辦法就是對于通過介質(zhì)的流向壓降有一個很好初始預(yù)測。猜測的辦法之一就是,在介質(zhì)流體單元的上游或者下游補償一個壓力值,詳細(xì)內(nèi)容請參閱所選單元的補償值一節(jié)。必須記住的是,當(dāng)補償壓力時,你所輸入的壓力可以定義為解算器所使用的gauge壓力(即在操作條件面板中定義的相對于操作壓力

45、的壓力)。另一個處理收斂性差的方法是臨時取消多孔介質(zhì)模型(在流體面板中關(guān)閉多孔區(qū)域)然后獲取一個不受多孔區(qū)域影響的初始流場。取消多孔區(qū)域后,F(xiàn)LUENT會將多孔區(qū)域處理為流體區(qū)域并按相應(yīng)的流體區(qū)域來計算。一旦獲取了初始解,或者計算很容易收斂,你就可以激活多孔模型繼續(xù)計算包含多孔區(qū)域的流場(對于大阻力多孔介質(zhì)不推薦使用該方法)。對于高度各向異性的多孔介質(zhì),有時會造成收斂性的麻煩。對于這些問題你可以將多孔介質(zhì)的各向異性系數(shù)(1/a_ij和C_2_i,j)限制在二階或者三階的量級。即使在某一方向上介質(zhì)的阻力為無窮大,你也不需要將它設(shè)定超過初始流動方向上的1000倍。多孔介質(zhì)的后處理可以通過檢查速度分

46、量和壓力值來確定多孔區(qū)域?qū)τ诹鲌龅挠绊?。你可能對下列變量或函?shù)的圖形(XY圖,等值線圖或者矢量圖)或者文檔報告感興趣:l X,Y,Z速度(在速度類別中)l 靜壓(在壓力類別中)這些變量會在后處理面板的變量選擇下拉菜單制定類別中出現(xiàn)。需要注意的是多孔區(qū)域的熱報告不影響固體介質(zhì)的屬性。所報告的多孔區(qū)域內(nèi)的熱容、傳導(dǎo)率以及焓是流體的屬性不包括固體介質(zhì)的影響。排氣扇邊界條件排氣扇模型是集總模型,可用于確定具有已知特征的排氣扇對于大流域流場的影響。排氣扇邊界類型允許你輸入控制通過排氣扇單元頭部(壓升)和流動速率(速度)之間關(guān)系的經(jīng)驗曲線。你也可以制定排氣扇旋轉(zhuǎn)速度的徑向和切向分量。排氣扇模型并精確模擬經(jīng)

47、過排氣扇葉片的詳細(xì)流動。它所預(yù)測的是通過排氣扇的流量。排氣扇的使用可能和其它流動源項關(guān)聯(lián),或作為模擬中流動的唯一源項。在后面的算例中,系統(tǒng)的流動速度由系統(tǒng)的損失和排氣扇曲線之間的相互平衡決定。FLUENT還提供了與用戶自定義模型之間的連接,這個模型在計算時更新了壓力跳躍函數(shù)。該功能在自定義排氣扇模型一節(jié)介紹。排氣扇方程模擬通過排氣扇的壓升在FLUENT的排氣扇模型中,排氣扇被看成無限薄,通過排氣扇的不連續(xù)壓升被指定為通過排氣扇速度的函數(shù)。它們之間的關(guān)系可能是常數(shù),多項式、分段線性函數(shù)或者分段多項式函數(shù),也可以是自定義函數(shù)。多于多項式情況,關(guān)系式為:其中D p為壓力升高(單位:Pa),f_n為壓

48、力跳躍多項式系數(shù),v垂直于排氣扇的當(dāng)?shù)亓黧w速度。速度v既可以是正也可以是負(fù)。你必須正確的模擬排氣扇以保證從排氣扇流過之后流體有個壓力升高的現(xiàn)象。對于排氣扇區(qū)域內(nèi)所有表面,你可以選擇使用垂直于排氣扇的質(zhì)量平均速度來確定單獨的壓力跳躍值。模擬排氣扇漩渦速度對于三維問題,對流的切向何徑向速度值可以加到排氣扇表面來產(chǎn)生渦流。這些速度可以指定為到排氣扇中心的徑向距離的函數(shù)。它們之間的關(guān)系可以是常數(shù)、多項式函數(shù)或者自定義函數(shù)。注意:所有渦流速度輸入都使用國際單位。對于多項式函數(shù),切向何徑向速度公式為:其中U_q和U_r分別為排氣扇表面的切向和徑向速度,單位為。m/s,f_n和g_n是切向和徑向速度的多項式

49、系數(shù),r為到排氣扇中心的距離。排氣扇的用戶輸入概述一旦排氣扇區(qū)域被確定(在邊界條件面板),你需要在排氣扇面板(下圖)中設(shè)定所有的模型輸入。該面板是從邊界條件菜單中打開的,詳細(xì)內(nèi)容清參閱邊界條件的設(shè)定一節(jié)。Figure 1: 排氣扇面板對于排氣扇,需要輸入如下:1. 確定排氣扇區(qū)域2. 定義通過排氣扇的壓力跳躍3. 為排氣扇定義離散相邊界條件(對于離散相計算)4. 需要的話,定義漩渦速度(只用于三維)確定排氣扇區(qū)域因為排氣扇被定義為無限薄,所以它必須被模擬為單元之間的界面而不是單元區(qū)域。因此排氣扇區(qū)域是內(nèi)部表面區(qū)域類型(其中表面是是二維中的線段或者三維中的三角形/四邊形)。當(dāng)你將網(wǎng)格讀入到FLU

50、ENT中時,如果排氣扇區(qū)域被確定為內(nèi)部區(qū)域,請使用邊界條件(見改變邊界區(qū)域類型)將適當(dāng)?shù)膬?nèi)部區(qū)域改變?yōu)榕艢馍葏^(qū)域。菜單:Define/Boundary Conditions.。內(nèi)部區(qū)域改變?yōu)榕艢馍葏^(qū)域后,你可以打開排氣扇面板并指定壓力跳躍,以及(可選)漩渦速度。定義壓力跳躍要定義壓力跳躍,你需要指定速度的多項式函數(shù)、分段線性函數(shù)、分段多項式函數(shù)或者常數(shù),也可以是自定義函數(shù)。你還應(yīng)該檢查區(qū)域平均方向矢量,保證流過排氣扇有個壓力升高。由解算器計算的區(qū)域平均方向是排氣扇區(qū)域的表面平均方向矢量。如果這個方向指向和排氣扇吹的方向一致就不用選擇排氣扇翻轉(zhuǎn)方向了,否則選擇排氣扇翻轉(zhuǎn)方向。對于壓力跳躍,請遵循

51、下面的步驟定義多項式函數(shù)、分段線性函數(shù)、分段多項式函數(shù):1. 檢查排氣扇面板,其中的壓力跳躍輪廓指定選項是關(guān)閉的。2. 在壓力跳躍右邊的下拉菜單中選擇多項式、分段線性或者分段多項式(如果所所要選擇的類型已被選中,你就可以點擊編輯按鈕打開定義函數(shù)的面板了)。3. 在定義壓力跳躍函數(shù)的面板中(如下圖)輸入適當(dāng)?shù)臄?shù)值。這些輪廓輸入面板和溫度相關(guān)屬性的輪廓輸入面板用法相同。請參閱使用溫度相關(guān)函數(shù)定義屬性來查看如何使用它。Figure 1: 壓力跳躍定義的多項式輪廓面板4. 設(shè)定下面所述的任何可選參數(shù)(此步可選)。當(dāng)你用這些函數(shù)的任何類型來定義壓力跳躍時,你可以限定計算壓力跳躍的速度值的最大和最小極限。

52、打開多項式速度范圍極限選項就可以設(shè)定速度范圍的最大最小值了。如果計算的法向速度范圍超出了你所指定的最大/最小速度范圍,那么解算器就會用極限值來替換它。你也可以選用垂直于風(fēng)扇的質(zhì)量平均速度來確定風(fēng)扇區(qū)域內(nèi)所有表面的單一的壓力跳躍值。打開從平均條件計算壓力跳躍可以激活這個選項。要定義常數(shù)壓力跳躍,請遵循如下步驟:1. 在排氣扇面板中打開指定壓力跳躍輪廓選項。2. 在壓力跳躍右邊的下拉菜單中選擇常數(shù)。3. 輸入壓力跳躍場中的D p值。 如果更方便的話,你也可以使用如下步驟:1. 打開壓力跳躍的輪廓指定選項。2. 在壓力跳躍輪廓下面的下拉菜單中選擇常數(shù),然后輸入壓力跳躍輪廓場的D p值。對于自定義壓力

53、跳躍函數(shù)或者邊界輪廓中定義的函數(shù),請遵循如下步驟:1. 打開壓力跳躍的輪廓指定選項。2. 在壓力跳躍輪廓下面的下拉菜單中選擇適當(dāng)?shù)暮瘮?shù),然后輸入壓力跳躍輪廓場的D p值。關(guān)于自定義函數(shù)的信息請參閱自定義函數(shù)一節(jié),關(guān)于邊界輪廓文件的信息請參閱邊界輪廓一節(jié)。下面的例子告訴了我們?nèi)绾未_定壓力跳躍的函數(shù)。考慮簡單的二維管流(如圖2)。進入長2.0m寬0.4m的導(dǎo)管的常密度空氣的速度為15 m/s。管的中心是個排氣扇。Figure 2: 定位于二維導(dǎo)管的排氣扇當(dāng)風(fēng)扇的轉(zhuǎn)速是2000rpm時,假定風(fēng)扇的特征如下:Q (m3/s)D p (Pa) 250.0 20175 15350 10525 5700 0

54、875 其中Q時通過風(fēng)扇的流動,D p時通過風(fēng)扇的壓升。在本例中,風(fēng)扇的特征為壓力升高和速度呈線性關(guān)系。要將這些特征轉(zhuǎn)換為壓力和速度的關(guān)系,必須知道風(fēng)扇的截面信息。在本例中,假定導(dǎo)管是1.0米深,面積為0.4平方米。相應(yīng)的速度值如下:v (m/s)D p (Pa) 62.50.0 50.0175 37.5350 25.0525 12.5700 0875 下面的對于一條線的方程是關(guān)系的多項式形式:為風(fēng)扇定義離散相邊界條件如果你是模擬粒子的離散相問題,你可以設(shè)定粒子在風(fēng)扇處的軌跡。關(guān)于邊界條件的設(shè)定清參閱離散相邊界條件一節(jié)。定義排氣扇旋轉(zhuǎn)速度如果你想在風(fēng)扇表面設(shè)定切向和徑向速度來產(chǎn)生三維問題中的渦

55、流,步驟如下:1. 在排氣扇面板打開漩渦速度指定選項。2. 定義軸的起始點(風(fēng)扇的起始點)和方向矢量(風(fēng)扇的旋轉(zhuǎn)軸)來指定風(fēng)扇的旋轉(zhuǎn)軸。3. 設(shè)定風(fēng)扇旋轉(zhuǎn)軸的半徑值。默認(rèn)為以避免多項式中出現(xiàn)除零問題。4. 設(shè)定切向和徑向速度為半徑的多項式函數(shù),常數(shù)值或者自定義函數(shù)。注意:渦流的速度輸入必須是國際標(biāo)準(zhǔn)單位。要定義切向和徑向速度的多項式函數(shù),步驟如下:1. 在排氣扇面板中,檢查切向速度的輪廓指定或者徑向速度的輪廓指定選項是關(guān)閉的。2. 輸入模擬排氣扇渦流中的方程1的系數(shù)f_n,或者在切向或徑向速度多項式系數(shù)框中模擬風(fēng)扇漩渦速度的方程3的g_n。首先輸入f_-1然后是f_0等。記住用空格符將每一個系數(shù)分開,第一個系數(shù)是(1/r)。要定義常數(shù)切向或者徑向速度,步驟如下:1. 在排氣扇面板中打開切向速度或者徑向速度的輪廓指定選項。2. 在切向或者徑向速度輪廓中選擇常數(shù)。3. 在切向或者徑向速度輪廓中輸入相應(yīng)的U_q或者U_r值。如果更方便的話,你可以遵照如下步驟:1. 在排氣扇面板中打開切向速度或者徑向速度輪廓指定選項。2. 在切向或者徑向多項式速度系數(shù)中輸入U_q或者U_r的值。對于

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論