1、湍流燃燒數(shù)值模擬李 國(guó) 岫 的研究與進(jìn)展1一、基本概念與基本原理二、湍流燃燒數(shù)值模擬研究狀況三、湍流燃燒模擬的主要模型四、非預(yù)混湍流火焰的數(shù)值模擬五、其他湍流火焰數(shù)值模擬的實(shí)例六、湍流燃燒數(shù)值模擬的進(jìn)展內(nèi)容2 1、什么是湍流（turbulence） （1）現(xiàn)象：流體的運(yùn)動(dòng)狀態(tài)可以分為層流和湍 流，當(dāng)流動(dòng)的特征雷諾數(shù)（ ）超過(guò) 相應(yīng)的臨界值，流動(dòng)從層流轉(zhuǎn)捩到湍流。 （2）定義：無(wú)一、基本概念與基本原理3（3）基本特征 大雷諾數(shù)下的流體運(yùn)動(dòng) 湍流的渦團(tuán)結(jié)構(gòu)和小渦團(tuán)的隨機(jī)性 有旋的和三維的 極強(qiáng)的擴(kuò)散性和耗散性 是流動(dòng)的屬性，而非流體的屬性 大尺度問(wèn)題具有擬序性（Coherence）和間歇性 （in

2、termittence）一、基本概念與基本原理42、湍流燃燒（Turbulent Combustion） 又稱為：湍流反應(yīng)流（Turbulent Reactive Flows） （1）燃燒 一種帶劇烈放熱化學(xué)反應(yīng)的流動(dòng)現(xiàn)象。 （2）湍流燃燒 一種極其復(fù)雜的帶劇烈放熱化學(xué)反應(yīng)的湍流 流動(dòng)現(xiàn)象。 （湍流化學(xué)反應(yīng)傳熱傳質(zhì)）一、基本概念與基本原理5（a）特征： 強(qiáng)非線性 高度耦合和關(guān)聯(lián) 高度隨機(jī)性（b）復(fù)雜性： 湍流問(wèn)題 湍流與燃燒的相互作用（Interaction） 流動(dòng)參數(shù)與化學(xué)動(dòng)力學(xué)參數(shù)之間的耦合機(jī)理一、基本概念與基本原理63、湍流燃燒數(shù)值模擬（Numerical Simulation） 應(yīng)用計(jì)

3、算機(jī)為工具，將流體力學(xué)、傳熱學(xué)、化學(xué)反應(yīng) 動(dòng)力學(xué)和數(shù)值計(jì)算方法相結(jié)合所得到的求解化學(xué)流體 力學(xué)基本方程的理論和方法。計(jì)算機(jī) 流體力學(xué)傳熱傳質(zhì)學(xué)化學(xué)反應(yīng)動(dòng)力學(xué) 數(shù)值計(jì)算方法湍流燃燒數(shù)值模擬 一、基本概念與基本原理7 （1）通過(guò)數(shù)值模擬，可以求出各種模型的數(shù)值解， 檢驗(yàn)、發(fā)展和優(yōu)化理論模型。 （2）深入認(rèn)識(shí)現(xiàn)有湍流燃燒的特征，進(jìn)一步揭示 其共性。 （3）為工業(yè)裝置的設(shè)計(jì)和優(yōu)化提供理論分析依據(jù)。作用：一、基本概念與基本原理8二、湍流燃燒模擬研究狀況1、基本思想 分別獨(dú)立描述湍流流動(dòng)和化學(xué)反應(yīng)過(guò)程， 然后考慮湍流流動(dòng)和化學(xué)反應(yīng)的相互作用。9 在描述流動(dòng)的守恒方程組中的能量方程和組分方程中，由于存在化學(xué)

4、反應(yīng)而產(chǎn)生的源項(xiàng)，即稱為化學(xué)反應(yīng)源項(xiàng)。表示化學(xué)反應(yīng)中能量的釋放速率和組分的生成（消耗）率，它是反應(yīng)物濃度和反應(yīng)流體溫度的強(qiáng)非線性函數(shù)。如何準(zhǔn)確計(jì)算化學(xué)反應(yīng)速率便成為湍流燃燒數(shù)值模擬的中心任務(wù)。2、問(wèn)題描述二、湍流燃燒模擬研究狀況10 致力于非線性源項(xiàng)的封閉和湍流流動(dòng)與 化學(xué)反應(yīng)的耦合與關(guān)聯(lián)的研究3、根本目標(biāo)二、湍流燃燒模擬研究狀況114、目前研究所關(guān)注的主要問(wèn)題 （1） 湍流流動(dòng)與化學(xué)反應(yīng)之間的相互作用組分濃度及溫度脈動(dòng)而強(qiáng)化組分的混合與傳熱 迅速放熱而引起密度變化，同時(shí)使流體輸運(yùn)系數(shù)變化湍流流動(dòng)化學(xué)反應(yīng)二、湍流燃燒模擬研究狀況12（2）強(qiáng)非線性源項(xiàng) 湍流燃燒模擬的中心問(wèn)題 反應(yīng)率的時(shí)均值 時(shí)

5、均值的反應(yīng)率（3） 需考慮密度變化而引起的封閉問(wèn)題（4） 各標(biāo)量（如溫度、組分等）與速度的脈動(dòng) 耦合量及標(biāo)量的強(qiáng)非線性源項(xiàng)的封閉問(wèn)題二、湍流燃燒模擬研究狀況135、湍流燃燒數(shù)值模擬的核心問(wèn)題： 湍流反應(yīng)率的封閉問(wèn)題 湍流反應(yīng)率的影響因素 湍流混合 分子擴(kuò)散 化學(xué)動(dòng)力學(xué)二、湍流燃燒模擬研究狀況14（一） 湍流模擬（二） 湍流燃燒相互作用模型（三） 化學(xué)反應(yīng)模型三、湍流燃燒數(shù)值模擬的主要模型15直接數(shù)值模擬（DNS）大渦模擬（LES）湍流輸運(yùn)模型雙流體模型湍流粘性系數(shù)模型 雷諾應(yīng)力模型（RSM）代數(shù)RSM微分RSM零方程模型單方程模型雙方程模型湍流模擬（一） 湍流模擬三、湍流燃燒模擬的主要模型16

6、1、 直接數(shù)值模擬（DNS） 直接數(shù)值求解非定常（瞬態(tài)）三維 Navier Stokes（NS）方程，來(lái)確定和描述湍流流 動(dòng)過(guò)程。特點(diǎn)： 無(wú)需任何模型，計(jì)算結(jié)果足夠精確。特征： 空間尺度上，要求計(jì)算區(qū)域足夠大，以 包含和分辨最大尺度的湍流渦團(tuán)，同時(shí) 計(jì)算網(wǎng)格尺寸又要足夠小，以模擬湍流 的最小渦團(tuán)。三、湍流燃燒模擬的主要模型17特征： 時(shí)間尺度上，時(shí)間必須大于大渦的 時(shí)間尺度，同時(shí)計(jì)算步長(zhǎng)又要小于 小渦的時(shí)間尺度。缺點(diǎn)： 計(jì)算量是 的量級(jí)，計(jì)算成本極 高，目前只適用于低雷諾數(shù)和簡(jiǎn)單 幾何形狀的湍流流動(dòng)。三、湍流燃燒模擬的主要模型18基本思想：認(rèn)為大尺度渦是高度各向異性的，小尺 度渦是近似各向同性

7、的。采用濾波器對(duì) NS方程進(jìn)行濾波，保留大渦特征，消 除小渦特征，對(duì)大渦進(jìn)行直接模擬，對(duì) 小渦引入亞網(wǎng)格尺度模型進(jìn)行模擬。2、大渦模型（LES）特點(diǎn)： 對(duì)DNS的近似，保證一定精度的前提下， 可以降低計(jì)算成本。缺點(diǎn)： 計(jì)算成本仍較高。 三、湍流燃燒模擬的主要模型193、 湍流輸運(yùn)模型基本思想：將瞬態(tài)NS方程的瞬時(shí)量分解為時(shí)均值 和脈動(dòng)值之和（雷諾分解），再取時(shí) 間平均，得到雷諾時(shí)均方程。然后利用 某些模擬假設(shè)，將方程中的高階的未知 關(guān)聯(lián)項(xiàng)用低階項(xiàng)或時(shí)均量來(lái)表達(dá)，從而 使雷諾（Reynolds）時(shí)均方程封閉。 三、湍流燃燒模擬的主要模型20 雷諾時(shí)均方程 三、湍流燃燒模擬的主要模型21 湍流模型

8、中心問(wèn)題 對(duì)于脈動(dòng)速度相關(guān)矩 （雷諾應(yīng)力），通過(guò)一定的假設(shè)，建立關(guān)于雷諾應(yīng)力的數(shù)學(xué)表達(dá)式或可以求解的輸運(yùn)方程。三、湍流燃燒模擬的主要模型22（1） 湍流粘性系數(shù)模型 Boussinesq于1887年提出 三、湍流燃燒模擬的主要模型23（1） 雙方程模型三、湍流燃燒模擬的主要模型24（1） 雙方程模型優(yōu)點(diǎn)：對(duì)于簡(jiǎn)單的流動(dòng)情況，（平壁邊界層、無(wú)浮 力平面射流、管流、尾跡流）能給出相當(dāng)滿 意的計(jì)算結(jié)果，計(jì)算工作量小。缺點(diǎn)：對(duì)于復(fù)雜的流動(dòng)（旋流、浮力流、曲壁邊界 層、圓射流） 無(wú)法準(zhǔn)確描述，甚至無(wú)法給出 正確的結(jié)果，缺乏普適性。 三、湍流燃燒模擬的主要模型25（2）雷諾應(yīng)力模型（Reynolds St

9、ress Model ）三、湍流燃燒模擬的主要模型26（2）雷諾應(yīng)力模型（Reynolds Stress Model ）基本思想：利用瞬時(shí)N-S方程和Reynolds時(shí)均方程， 直接推導(dǎo)出雷諾應(yīng)力的輸運(yùn)方程，然后 再對(duì)輸運(yùn)方程中的未知項(xiàng)進(jìn)行模擬。優(yōu)點(diǎn)： 可以較為準(zhǔn)確描述突擴(kuò)回流、鈍體回流、 旋流、浮力流等復(fù)雜流動(dòng)，具有更廣的 適用范圍、 更高的預(yù)測(cè)能力、更好的計(jì) 算精度。缺點(diǎn)： 計(jì)算工作量仍較大 。三、湍流燃燒模擬的主要模型27(a) 微分雷諾應(yīng)力模型（DRSM）(b) 代數(shù)雷諾應(yīng)力模型（ASM） 某些情況下（高剪切的流動(dòng)，局部平衡的湍流），可以消去微分雷諾應(yīng)力方程中的對(duì)流與擴(kuò)散項(xiàng)，方程轉(zhuǎn)化為

10、代數(shù)方程。使計(jì)算量減小。 三、湍流燃燒模擬的主要模型281 關(guān)聯(lián)矩模型2 Arrhenius 公式3、 EBU模型（EddyBreak-up Model）4快速反應(yīng)模型5特征時(shí)間模型6簡(jiǎn)化PDF模型7PDF輸運(yùn)方程模型（二）湍流燃燒相互作用模型三、湍流燃燒模擬的主要模型29 湍流燃燒數(shù)值模擬的核心問(wèn)題： 湍流反應(yīng)率的封閉問(wèn)題 湍流反應(yīng)率的影響因素 湍流混合 分子輸運(yùn) 化學(xué)動(dòng)力學(xué)三、湍流燃燒模擬的主要模型301、關(guān)聯(lián)矩模型 時(shí)均反應(yīng)率 相關(guān)矩 等是新的未知數(shù)，需要進(jìn)行模擬。三、湍流燃燒模擬的主要模型312、Arrhenius公式 簡(jiǎn)化化學(xué)反應(yīng)系統(tǒng)特征： 只考慮化學(xué)反應(yīng)動(dòng)力學(xué)的作用。缺點(diǎn)： 未考慮

11、湍流混合、分子輸運(yùn)兩方面 的因素。 三、湍流燃燒模擬的主要模型32假設(shè)：認(rèn)為化學(xué)反應(yīng)率取決于未燃?xì)馕F(tuán)在湍 流作用下破碎成更小微團(tuán)的速率公式：特征：突出湍流混合對(duì)燃燒速率的控制作用。缺點(diǎn)： 未能考慮分子輸運(yùn)和化學(xué)動(dòng)力學(xué)因素的 影響，過(guò)于粗糙。3、EBU模型（EddyBreak-up Model）三、湍流燃燒模擬的主要模型33 快速反應(yīng)假設(shè)：假設(shè)化學(xué)反應(yīng)速率與湍流混合（擴(kuò)散）速率相比無(wú)窮快。即湍流燃燒過(guò)程由燃料與氧化劑，或已燃?xì)怏w與未燃?xì)怏w的混合過(guò)程控制。4、 快速反應(yīng)模型 （1）反應(yīng)面模型 （2） 模型 （3） 模型 三、湍流燃燒模擬的主要模型34假設(shè)：認(rèn)為燃料與氧無(wú)論在同一時(shí)間還是同一空間內(nèi)

12、均不 共存，亦即存在二者濃度均為零的反應(yīng)面或火焰面 （1） 模型混合物分?jǐn)?shù)f 的時(shí)均值方程 與k， 方程聯(lián)立，組成 模型缺點(diǎn)：過(guò)于簡(jiǎn)化，假設(shè)成立得太勉強(qiáng)。 三、湍流燃燒模擬的主要模型35物質(zhì)m的質(zhì)量分量質(zhì)量分量的瞬時(shí)的熱力平衡狀態(tài)質(zhì)到達(dá)平衡度的特征時(shí)間層流時(shí)間尺度湍流時(shí)間尺度延遲系數(shù)5、特征時(shí)間模型 由化學(xué)反應(yīng)引起的某種質(zhì)量分?jǐn)?shù)的變化率： 三、湍流燃燒模擬的主要模型36 （1）PDF模型6、 簡(jiǎn)化PDF模型基本思想： 化學(xué)反應(yīng)率是熱力學(xué)狀態(tài)量 ，T和各組分質(zhì)量分?jǐn)?shù)的非線性函數(shù)，而這些量的隨機(jī)脈動(dòng)對(duì)平均反應(yīng)率有強(qiáng)烈的影響，因此考慮采用概率統(tǒng)計(jì)的方法來(lái)描述。采用概率密度函數(shù)（PDF）可以得到平均反

13、應(yīng)率的表達(dá)式。三、湍流燃燒模擬的主要模型37簡(jiǎn)化PDF（設(shè)定PDF ） PDF輸運(yùn)方程（2）確定PDF的方法三、湍流燃燒模擬的主要模型38 分布 函數(shù)截尾Gauss分布 即先假定輸運(yùn)變量脈動(dòng)的概率密度函數(shù)的具體形式，通過(guò)確定其中的一些待定參數(shù)獲得概率分布。 設(shè)定PDF常采用以下幾種分布（3）簡(jiǎn)化PDF模型三、湍流燃燒模擬的主要模型39 由于湍流脈動(dòng)的存在，使瞬時(shí)值不共存，平均值共存。 在 模型的基礎(chǔ)上，進(jìn)一步求解混合物脈 動(dòng)的均方值 的方程。 a、 模型三、湍流燃燒模擬的主要模型（4）幾個(gè)典型的簡(jiǎn)化PDF模型40對(duì)于任意瞬時(shí)標(biāo)量是f的函數(shù) 的時(shí)均值即為 三、湍流燃燒模擬的主要模型41特征：基于

14、快速反應(yīng)假設(shè)，突出了湍流混合過(guò)程在湍流 燃燒中的控制作用。優(yōu)點(diǎn)：某些應(yīng)用實(shí)例與實(shí)驗(yàn)結(jié)果符合較好，優(yōu)于 模型。缺點(diǎn)：完全忽略燃燒與氧化劑分子的相互擴(kuò)散作用。 三、湍流燃燒模擬的主要模型42 基本原理： 把湍流擴(kuò)散火焰看成是由無(wú)數(shù)個(gè)層流擴(kuò)散小火 焰組成的渦團(tuán)。湍流與化學(xué)反應(yīng)相互作用的問(wèn) 題分成兩部分：（1）按照守恒標(biāo)量定義進(jìn)行準(zhǔn)穩(wěn)定一維層流擴(kuò) 散小火焰結(jié)構(gòu)的計(jì)算；（2）湍流火焰中這種結(jié)構(gòu)出現(xiàn)的概率分布。b、層流小火焰模型（Flamelet Model ）三、湍流燃燒模擬的主要模型43 基本思想： 把整個(gè)燃燒場(chǎng)看成是隨機(jī)分布的小火焰的集合，從組分和焓的守恒方程入手，引入非平衡參數(shù)瞬時(shí)標(biāo)量耗散率（考慮

15、流動(dòng)的影響），利用混合百分?jǐn)?shù)的守恒方程得出小火焰結(jié)構(gòu)，再?gòu)恼w上考慮，綜合概率密度函數(shù)，從而得到湍流燃燒過(guò)程中各瞬時(shí)值在整個(gè)燃燒區(qū)域內(nèi)的統(tǒng)計(jì)行為。三、湍流燃燒模擬的主要模型44 模型實(shí)現(xiàn)： Flamelet模型從組分和焓的守恒方程出發(fā)，將焓的守恒方程轉(zhuǎn)化為溫度的偏微分方程。通過(guò)定義混合百分?jǐn)?shù)，并將其作為獨(dú)立變量，按照坐標(biāo)轉(zhuǎn)換的原則得到一維準(zhǔn)穩(wěn)定結(jié)構(gòu)的組分和溫度的偏微分方程，同時(shí)引入非平衡參數(shù)瞬時(shí)標(biāo)量耗散率，利用混合百分?jǐn)?shù)和耗散率的概率密度函數(shù)，得到湍流燃燒過(guò)程中各瞬時(shí)值在整個(gè)燃燒區(qū)域內(nèi)的統(tǒng)計(jì)行為。三、湍流燃燒模擬的主要模型45基本思想：直接求解關(guān)于概率密度函數(shù)PDF的輸運(yùn) 方程，求出所有有關(guān)流

16、動(dòng)與燃燒的參量。7、 PDF 輸運(yùn)方程模型取消了其他模型的假設(shè)前提。對(duì)于守恒方程中的對(duì)流項(xiàng)、非線形化學(xué) 反應(yīng) 項(xiàng)、平均壓力項(xiàng)可以精確處理，不 需模擬?？梢蕴峁┝鲌?chǎng)的完整信息?？梢阅M著火、熄火、湍流燃燒和排放污染物生成過(guò)程。優(yōu)點(diǎn)：三、湍流燃燒模擬的主要模型46PDF模型分類：（1）速度PDF方程（2）標(biāo)量PDF方程（3）速度標(biāo)量聯(lián)合PDF方程（4）速度、耗散率和標(biāo)量聯(lián)合PDF方程三、湍流燃燒模擬的主要模型47（1）速度PDF方程 對(duì)流項(xiàng)是封閉的。 避免梯度擴(kuò)散模型的模擬。三、湍流燃燒模擬的主要模型48（2）標(biāo)量PDF方程 對(duì)于任意復(fù)雜的化學(xué)反應(yīng)機(jī)理都可以 精確計(jì)算，無(wú)需模擬。 無(wú)法計(jì)算湍流流場(chǎng)

17、，仍需采用湍流模 型計(jì)算流場(chǎng)。三、湍流燃燒模擬的主要模型49Micro-mixingReactionSource of fa at boundary of composition spaceConservation of probability三、湍流燃燒模擬的主要模型50（3）速度標(biāo)量聯(lián)合PDF方程 將速度PDF模型與標(biāo)量PDF模型結(jié)合起 來(lái)，具備兩者的優(yōu)勢(shì)。 對(duì)流項(xiàng)是封閉的 化學(xué)反應(yīng)項(xiàng)封閉 無(wú)需湍流模型三、湍流燃燒模擬的主要模型51（3）速度標(biāo)量聯(lián)合PDF方程 缺點(diǎn)：沒(méi)有包含有關(guān)湍流長(zhǎng)度尺度和時(shí)間尺度的信息，需要應(yīng)用耗散率方程，而耗散率方程的建立幾乎沒(méi)有任何客觀依據(jù)。從而造成在解決一些復(fù)雜

18、的情況時(shí)，速度標(biāo)量聯(lián)合PDF方程效果不理想三、湍流燃燒模擬的主要模型52三、湍流燃燒模擬的主要模型53（4）速度、耗散率和標(biāo)量聯(lián)合PDF方程 以提供湍流時(shí)間尺度和長(zhǎng)度的信息。 適應(yīng)于各向異性湍流。三、湍流燃燒模擬的主要模型541 單步化學(xué)反應(yīng)模型2 簡(jiǎn)化多步化學(xué)反應(yīng)模型3詳細(xì)化學(xué)反應(yīng)模型（三）化學(xué)反應(yīng)模型三、湍流燃燒模擬的主要模型55非預(yù)混湍流火焰（Non-premixed Flame)： 燃料和氧化劑沒(méi)有預(yù)先混合好，而是邊混合邊燃燒。四、非預(yù)混湍流火焰的數(shù)值模擬56非預(yù)混Bluff-body Flame四、非預(yù)混湍流火焰的數(shù)值模擬57簡(jiǎn)介 1. Burner: bluff-body burn

19、er2. Fuel: CH4/H2( 1/1 by volume)3. Inlet: Vjet = 118 m/s, Vair = 40 m/s4. 湍流模型：微分雷諾應(yīng)力模型（DRSM)5. 湍流燃燒模型：簡(jiǎn)化PDF模型 PDF輸運(yùn)方程模型6. 化學(xué)反應(yīng)模型：簡(jiǎn)化化學(xué)反應(yīng)模型四、非預(yù)混湍流火焰的數(shù)值模擬58 湍流模型 （1） 微分雷諾應(yīng)力模型（DRSM)四、非預(yù)混湍流火焰的數(shù)值模擬59（2）計(jì)算中所應(yīng)用的雷諾應(yīng)力模型1. LRR-IP 模型2. LRR-IP（ BM-M1) 模型3. JM 模型4. SSG 模型5. BM-M2 模型四、非預(yù)混湍流火焰的數(shù)值模擬60Bluff-body st

20、abilised flame CH4/H2AirFuel96 cells72 cells0.1 m0.1 m0.2 m61Bluff-body stabilised flame CH4/H2Mean axial velocity:Full second moment closure / Fast chemistrym/s62Stream functionFull second moment closure / Fast chemistryBluff-body stabilised flame CH4/H263Turbulent kinetic energyFull second moment

21、closure / Fast chemistryBluff-body stabilised flame CH4/H2m2/s264Turbulent dissipationFull second moment closure / Fast chemistryBluff-body stabilised flame CH4/H2m2/s365Mean mixture fraction:Full second moment closure / Fast chemistryBluff-body stabilised flame CH4/H266Mixture fraction variance:Ful

22、l second moment closure / Fast chemistryBluff-body stabilised flame CH4/H267Mean temperatureFull second moment closure / Fast chemistryBluff-body stabilised flame CH4/H2K68Reynolds stress:Full second moment closure / Fast chemistryBluff-body stabilised flame CH4/H2m2/s269Reynolds stress:Full second

23、moment closure / Fast chemistryBluff-body stabilised flame CH4/H2m2/s270Reynolds stress:Full second moment closure / Fast chemistryBluff-body stabilised flame CH4/H2m2/s271Reynolds stress:Full second moment closure / Fast chemistryBluff-body stabilised flame CH4/H2m2/s272Reynolds flux:Full second mo

24、ment closure / Fast chemistryBluff-body stabilised flame CH4/H2m/s73Reynolds flux:Full second moment closure / Fast chemistryBluff-body stabilised flame CH4/H2m/s74Reynolds flux:k-e model / Fast chemistryBluff-body stabilised flame CH4/H2m/s75k-e model / Fast chemistryReynolds flux:Bluff-body stabil

25、ised flame CH4/H2m/s76Turbulent viscosityFull second moment closure / Fast chemistryBluff-body stabilised flame CH4/H2kg/m.s77Mean mass fraction of OHFull second moment closure / Fast chemistryBluff-body stabilised flame CH4/H278（2）計(jì)算結(jié)果 軸向速度U（x/D = 0.2）四、非預(yù)混湍流火焰的數(shù)值模擬79（2）計(jì)算結(jié)果 軸向速度U（x/D = 1.0）四、非預(yù)混湍流火焰的數(shù)值模擬80（2）計(jì)算結(jié)果 軸向速度U（x/D = 2.4）四、非預(yù)混湍流火焰的數(shù)值模擬81（2）計(jì)算結(jié)果 軸向雷諾應(yīng)力（x/D = 0.2）四、非預(yù)混湍流火焰的數(shù)值模擬82（2）計(jì)算結(jié)果 軸向雷諾應(yīng)力（x/D = 1.0）四、非預(yù)混湍流火焰的數(shù)值模擬83（2）計(jì)算結(jié)果 軸向雷諾應(yīng)力（x/D = 2.4）四、非預(yù)混湍流火焰的數(shù)值模擬84（2）計(jì)算結(jié)果 平均溫度T