fluent61的二次開發(fā)在湍流燃燒rans模擬和大渦模擬les）中的應用

1、2005 Fluent 中國用戶大會論文集 156 FLUENT 軟件的二次開發(fā)用于 湍流燃燒的 RANS 模擬和大渦模擬（LES）研究 周力行 1，王方 1，胡礫元 1，2，王文麗 1 1． 清華大學工程力學系，北京 100084 2． 上海交通大學機械和能源學院，上海 200024 摘 要: 摘 要: 本文作者們利用 FLUENT 軟件的 UDF/UDS 功能， 對該軟件進行了二次開發(fā)， 將我們自己提出的二階矩燃燒模型分

2、別和 FLUENT 軟件對流動的雷諾應力方程 RANS 模擬和流動的大渦模擬（LES）相結合，加入到該軟件的平臺上，先后分別對甲 烷-空氣旋流擴散燃燒、 甲烷-空氣射流擴散燃燒以及丙烷-空氣旋流擴散燃燒進行了 RANS 模擬和大渦模擬， 并且和相應的實驗數(shù)據(jù)進行了對照。研究結果表明，二階矩燃燒模型的 RANS 模擬結果優(yōu)于通常廣泛使用的 EBU-Arrhenius 湍流燃燒模型和 簡化 PDF 湍流反應模型的結果。 作為大渦模擬的亞網(wǎng)格

3、燃燒模型， 二階矩燃燒模型也優(yōu)于 EBU-Arrhenius 燃燒模型。 大渦模擬瞬態(tài)結果揭示了射流和旋流流動和擴散燃燒的瞬態(tài)湍流結構和火焰結構的特點。 大渦模擬和標量 PDF 方程模型對關聯(lián)量的 統(tǒng)計結果證實了二階矩燃燒模型的封閉模型的合理性。 1．引 ．引 言 尋求經(jīng)濟而合理的湍流燃燒模型，是工業(yè)裝置內湍流燃燒和污染物生成 CFD 數(shù)值模擬的關鍵。近年來發(fā)展的直接模擬（DNS）、大渦模擬（LES）和概率密度函數(shù)（PDF）輸運方程模擬

4、比較嚴格，合理，但是計算量太大，難以直接應用于工程中大尺寸復雜流動中，只適合作為基礎研究手段。層流小火焰模型和條件矩封閉（CMC）模型是比較有希望的模型，正處于發(fā)展過程中，有待實驗檢驗和進一步完善。約 30 年前Spalding 等人發(fā)展的 EBU（EDC），EBU-Arrhenius（E-A）模型及其 Magnusen 修正和簡化（設定）PDF 模型等 [1]至今仍然被廣泛用于商業(yè)軟件中， 用來預報工程裝置中復雜流動的湍流燃燒和 NO

5、 生成。 但是二者給出的預報結果有時和實驗結果出入較大。另一方面，多年前人們就想到了非條件二階矩湍流燃燒模型。它采取類似于湍流模型的封閉方法。早期的二階矩模型 [1]對反應率表達式中的非線性指數(shù)項采取級數(shù)展開近似，假設 E/RT> ，同時 T //T 也并非總是遠遠小于 1，因此這種級數(shù)展開近似法有嚴重的誤差。原因是，在流場的大部分地區(qū) E-A 模型事實上主要考慮了湍流的作用，難以計入化學反應動力學，或者有限反應率的作用。簡化 P

6、DF 模型用溫度的單變量 PDF 和濃度的單變量 PDF 的乘積對溫度-濃度聯(lián)合 PDF 取近似，低估了湍流反應率。近年來我們系統(tǒng)地對二階矩湍流-反應模型進行了研究，提出了統(tǒng)一二階矩（USM）湍流燃燒/反應模型[3]。USM 模型在求解時均反應率時，對濃度脈動關聯(lián)以及反應率系數(shù) k 的脈動和濃度脈動關聯(lián)都用統(tǒng)一形式的輸運方程加以封閉，同時考慮了化學反應率對關聯(lián)量耗散的影響。該模型用于模擬 Sandia 實驗室測量的甲烷-空氣湍流射流燃燒

7、和 NO 生成，取得了很好的效果。為了進一步在復雜流動的湍流燃燒中應用和檢驗 USM 模型，我們應用 FLUENT 軟件，對其進行二次開發(fā)。首先將 USM 湍流燃燒模型通過 UDS/UDF 功能和流動的 Reynolds 應力方程模型的 RANS 模擬相結合，加入到 FLUENT 的平臺上，分別對本實驗室測量的甲烷-空氣旋流擴散燃燒和 NO 生成以及文獻中測量的丙烷-空氣擴散燃燒進行數(shù)值模擬， 將模擬結果和實驗結果以及 E-A 模型和簡

8、化 PDF 模型的結果讀對照， 驗證 USM2005 Fluent 中國用戶大會論文集 158 ) / 10 027 . 2 exp( 10 119 . 2 8 3 . 1 2 2 . 0 4 2 11 RT Y Y W O CH fu × ? × = ρ（5） 對熱 NO 生成的層流反應率采用擴展的 Zeldovich 機理[5]，層流反應率的表達式為 ? ?? ? ??? × = ? + T T Y Y

9、 W O N Th NO 69160 exp 10 35 . 1 1 5 . 0 16 , 2 2 ρ（6） ? ?? ? ?? ? = ? ? ? T Y Y T W O N Th NO 47355 exp 6 . 22 5 . 0 2 1 , 2 2 ρ（7） 該式也可以簡化為如下的一步反應 ? ?? ? ??? × = T T Y Y W O N Z NO 38440 exp 10 9 3 . 0 2 12 , 2 2

10、ρ（8） 對 NH3轉化成的 NO，采用以下反應機理[6]： O H NO O NHO H N NO NH2 2 32 2 3 6 4 5 46 5 6 4+ ? ++ ? +（9） 對應的燃料 NO 層流反應率表達式為： ? ?? ? ??? × =? ?? ? ??? × =+?T Y Y WT Y Y WO NH fuel NONO NH fuel NO50000 exp 10 23 . 525000 exp

11、10 05 . 12 332 22 ,2 16 ,ρρ（10） 4. EBU-Arrhenius 燃燒模型和簡化 燃燒模型和簡化 PDF 的 NO 湍流 湍流-反應模型 反應模型 為了對照不同的湍流-反應模型， 同時采用廣泛應用的 E-A 湍流燃燒模型和簡化 PDF 的 NO 湍流反應率模型 [1]。E-A 湍流燃燒模型的表達式是 ] , min[ sT sA s w w w =（11） 其中 ) / exp( 2 122 12 T

12、R E Y Y B Y Y k wsA ? = = ρ ρ（12） ? ? ?? ? ?+ ? = β βε ρ 1 , , min 2 4 P EBU O CH EBU sT Y B Y Y k A w（13） 此處，β是化學反應的化學當量比，一般 AEBU取為 4.0、BEBU取為 0.5。 簡化概率密度函數(shù)模型 NO 生成的時均反應率為： ∫∫ = 2 2 2 ) ( ) ( ) , ( 2 1 O O O NO NO dTdY