1、湍流邊界條件設置在流場的入口、出口和遠場邊界上，用戶需要定義流場的湍流參數。在 FLUENT中可以使用的湍流模型有很多種。在使用各種湍流模型時，哪些變量需 要設定，哪些不需要設定以及如何給定這些變量的具體數值，都是經常困擾用戶的問 題。本小節只討論在邊界上設置均勻湍流參數的方法，湍流參數在邊界上不是均勻分 布的情況可以用型函數和UDF(用戶自定義函數)來定義，具體方法請參見相尖章節 的敘述。在大多數情況下，湍流是在入口后面一段距離經過轉扳形成的，因此在邊界上 設置均勻湍流條件是一種可以接受的選擇。特別是在不知道湍流參量的分布規律時， 在邊界上采用均勻湍流條件可以簡化模型的設置。在設置邊界條件時

2、，首先應該定 性地對流動進行分析，以便邊界條件的設置不違背物理規律。違背物理規律的參數設 置往往導致錯誤的計算結果，甚至使計算發散而無法進行下去。在TurbulenceSpecificationMethod (湍流定義方法)下拉列表中，可以簡單地用 一個常數來定義湍流參數，即通過給定湍流強度、湍流粘度比、水力直徑或湍流特征 長在邊界上的值來定義流場邊界上的湍流。下面具體討論這些湍流參數的含義，以保 證在設置模型時不出現違背流動規律的錯誤設置：(1) 湍流強度(Turbulence Intensity)湍流強度I的定義為：l=Sqrt(u f*u '+vf*v ! +w *-w1)

3、9; )/u_avg(8上式中u；v，和w，是速度脈動量，u_avg是平均速度。湍流強度小于1%時，可以認為湍流強度是比較低的，而在湍流強度大于10% 時，則可以認為湍流強度是比較高的。在來流為層流時，湍流強度可以用 繞流物體的 幾何特征粗略地估算出來。比如在模擬風洞試驗的計算中，自由流的湍流強度可以用 風洞的特征長度估計出來。在現代的低湍流度風洞中，自由流的湍流強度通常低于0.05%。內流問題進口處的湍流強度取決于上游流動狀態。如果上游是沒有充分發 展的 未受擾流動，則進口處可以使用低湍流強度。如果上游是充分發展的湍 流，則進口處 湍流強度可以達到幾個百分點。如果管道中的流動是充分發展的 湍

4、流，則湍流強度可 以用公式(8-2)計算得到，這個公式是從管流經驗公式得到的：l=u ' /u_avg=0.16*Re_DHA-0.125(8-2)其中Re_DH是HydraulicDiameter (水力直徑)的意思，即式(8-2)中的雷諾數 是以水力直徑為特征長度求出的。(2) 湍流的長度尺度與水力直徑湍流能量主要集中在大渦結構中，而湍流長度尺度I則是與大渦結構相尖的物理量。在充分發展的管流中，因為漩渦尺度不可能大于管道直徑，所以 I是受 到管道尺寸制約的幾何量。湍流長度尺度I與管道物理尺寸L尖系可以表示為：0.07L( 8-3)式中的比例因子0.07是充分發展管流中混合長的最大值

5、，而L則是管道直徑。在管道截面不是圓形時，L可以取為管道的水力直徑。湍流的特征長取決于對湍流發展具有決定性影響的幾何尺度。在上面的討論中，管道直徑是決定湍流發展過程的唯一長度量。如果在流動中還存在其他對流動影 響更大的物體，比如在管道中存在一個障礙物，而障礙物對湍流的發生和發展過 程起 著重要的干擾作用。在這種情況下，湍流特征長就應該取為障礙物的特征長度。從上面的分析可知，雖然式（8-2）對于大多數管道流動是適用的，但并不是 普遍適用的，在某些情況下可以進行調整。在FLUENT中選擇特征長L或湍流長度尺度I的方法如下：1 ）對于充分發展的內流，可以用IntensityandHydraulicD

6、iameter （湍流強度 與水力直徑）方法定義湍流，其中湍流特征長度就是Hydraulic Diameter （水力直 徑）HD。2）對于導向葉片或分流板下游的流場，可以用IntensityandHydraulicDiameter （湍流強度與水力直徑）定義湍流，并在 HydraulicDiameter （水力直徑）中將導向葉片或分流板的開口部分的長度 L定義 為特征長度。3）如果進口處的流動為受到壁面限制且帶有湍流邊界層的流動，可以在Intensity and LengthScale面板中用邊界層厚度delta_99通過公式匸0.4*delta_99計算得到湍流長度尺度|。最后在Turbu

7、lenee Length Scale （湍 流 長度尺度）中輸入I的值。（3）湍流粘度比湍流粘度比mu_t/mu與湍流雷諾數Re_t成正比。湍流雷諾數的定義為：Re_t=k*k/ （Epsilon*nu） （8-4）在高雷諾數邊界層、剪切層和充分發展的管道流動中的數值較大，其量級大約 在100到1000之間。而在大多數外部流動的自由流邊界上，湍流粘度比的值很 小。在典型情況下，其值在1到10之間。4）推導湍流變量時采用的尖系式為了從前面講到的湍流強度I,湍流長度尺度L和湍流粘度比mu_t/mu求出其他 湍流變量'必須采用幾個經驗矢系式。在FLUENT中使用的經驗尖系式主要包括下面幾種：

8、1 )從湍流強度和長度尺度求出修正的湍流粘度在使用Spalart-Allmaras模型時，可以用湍流強度I和長度尺度I求出修正的 湍 流粘度，具體公式如下：nu=Sqrt(1.5)*u_avg*l*L(8-5)在使用FLUENT時，如果在Spalart-Allmaras模型中選擇Intensity andHydraulic Diameter (湍流強度與水力直徑)選項，則修正的湍流粘度就用這個公式求 出。其中的長度尺度I則用式(8-3)求出。2) 用湍流強度求出湍流動能湍流動能k與湍流強度I的尖系如下：k=1.5*(u_avg*l)A2(8-6)如果在使用FLUENT時沒有直接輸入湍流動能k和

9、湍流耗散率Epsilon的值5則可以使用Intensity and Hydraulic Diameter (湍流強度與水力直徑)、Inten sity and Len gth Scale (湍流強度與長度尺度)或 Inten sitya ndViscosityRatio(湍流強度與粘度比)等方法確定湍流動能，而確定的辦法就是使用上面的公式(8- 6) °3) 用長度尺度求出湍流耗散率長度尺度I與湍流耗散率之間的尖系為：epsilon=C_muA0.75*kA1.5/l(8-7)式中C_mu為湍流模型中的一個經驗常數，其值約等于0.09。在沒有直接輸入湍流動能k和湍流耗散率epsilo

10、n的情況下，可以用IntensityandHydraulicDiameter （:帶流強度與水力直徑）或 Intensity and LengthScale （湍流強度與長度尺度）等辦法，利用上述公式確定湍流耗散率epsilon。4）用湍流粘度比求出湍流耗散率湍流耗散率epsilon與湍流粘度比mu_t/nriu和湍流動能k的尖系如下：epsilo n=rho*C_mu*"2/mu*（mu_t/mu）7Vl （8-8）式中C_mu為湍流模型中的一個經驗常數，其值約等于0.09。在沒有直接輸入湍流動能k和湍流耗散率epsilon的情況下，可以用Intensity and Viscosi

11、tyRatio（湍流強度與粘度比）定義湍流變量，實際上就是利用上述公式算 出湍流耗散率epsilon。5）湍流衰減過程中湍流耗散率的計算如果計算風洞阻尼網下游試驗段中的流場，可以用下式求出湍流耗散率Epsilon :epsilo n=delta_k*Ui_farfield/L_farfield（8-9）式中delta_k是湍流動能k的衰減量，比如可以設為入口處k值的10%, UJarfield是自由流速度，LJarfield是自由流區域的長度。（8-9）式是對高雷諾數各 向同性湍流衰減指數律的線性近似，其理論基礎是衰減湍流中湍流動能k的方程：U*（partial derivative of U

12、 with respect to x）= -epsilon（8-10）如果用這種方法計算epsilon,還需要用（8-8）式檢驗計算結果，以保證 湍流粘 度比不過大。雖然這種方法在FLUENT中沒有使用，但是可以用這種方法 估算出自由流中的湍流耗散率epsilon，然后再用（8-6）式確定k，最后在Turbulenee Specification Method （湍流定義方法）下拉列表中選擇 KandEpsilon （k和Epsilon）并k和Epsilon的計算結果輸入到相應的欄目中。6）用長度尺度計算比耗散率如果知道湍流長度尺度I，可以用下式確定omega:omega=kA0.5/ （C_

13、muA0.25*1）（8-11）式中C_mu和長度尺度I的取法與前面段落中所述相同。在使用IntensityandHydraulicDiameter （；需流弓雖度與水力直徑）或 Intensity and Length Scale （湍流強度與長度尺度）定義湍流時，FLUENT用的就是這種方法。7）用湍流粘度比計算比耗散率omega的值還可以用mu_t/mu和k通過下式計算得出：omega=rho*k/mu* （mu_t/mu） A-1 （8-12）在使用Intensity and Viscosity Ratio （湍流強度與粘度比）方法定義湍流時， FLUENT就是使用上述尖系式對湍流進行定義的。8）用湍流動能定義雷諾應力分量在使用RSM （雷諾應力模型）時，如果用戶沒有在Reynolds-StessSpecificationMethod （雷諾應力定義方法）的 Reynolds-Stress Components （雷諾應力分量）選項中直接定義雷諾應力的值，則雷諾應力的值將由給定的k值 計算得出。假定湍流是各向同性的，即：Average （uu !-_1j3） =） 0 （8且：Average（u ' _aphla* u ' _aphla）-=124k） /3（8如果用戶在Reynolds-Stress Specification Method （