FLUENT中組分輸運(yùn)及化學(xué)反應(yīng)(燃燒)模擬

1、燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)Temperature in a gas furnaceFLUENT中中組分輸運(yùn)及化學(xué)反應(yīng)組分輸運(yùn)及化學(xué)反應(yīng)（燃燒）模擬燃燒）模擬熱科學(xué)與能源工程系2003年10月燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)概要u應(yīng)用u燃燒模擬簡(jiǎn)介 u化學(xué)動(dòng)力學(xué)u氣相燃燒模型u稀疏相燃燒模型 u污染物排放模擬u燃燒數(shù)值模擬步驟介紹燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)燃燒模擬u廣泛應(yīng)用與均相和非均相燃燒過(guò)程模擬l燃燒爐l鍋爐l加熱器l燃?xì)廨啓C(jī)l火箭發(fā)動(dòng)機(jī)u求解內(nèi)容l流場(chǎng)流動(dòng)特性及其混合特性l溫度場(chǎng)l組分濃度場(chǎng)l顆粒和污染物排

2、放Temperature in a gas furnaceCO2 mass fractionStream function燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)燃燒模型概要稀疏相模型液滴/顆粒動(dòng)力學(xué)非均相反應(yīng)液化蒸發(fā)輸運(yùn)控制方程質(zhì)量動(dòng)量 (湍流)能量化學(xué)組分燃燒模型預(yù)混局部預(yù)混非預(yù)混燃燒污染物模型輻射換熱模型燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)u氣相燃燒l有限速率模型 (Magnussen model)l守恒標(biāo)量的 PDF模型 (一個(gè)或兩個(gè)混合分?jǐn)?shù))l層流火焰面（小火焰）模型 (V5)lZimont model (V5)u稀疏相模型n湍流顆粒彌散s隨機(jī)軌道

3、模型（Stochastic tracking）s顆粒云團(tuán)模型（Particle cloud model） (V5)n粉煤與噴油燃燒子模型u 輻射模型: DTRM, P-1, Rosseland 和Discrete Ordinates (V5)u湍流模型: k-, RNG k-, RSM, Realizable k- (V5) and LES (V5)u污染物排放模型: NOx with reburn chemistry (V5) and sootFLUENT提供的燃燒模型燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)模擬燃燒過(guò)程的化學(xué)反應(yīng)動(dòng)力學(xué)u難點(diǎn)與挑戰(zhàn)l多數(shù)實(shí)際的燃燒過(guò)程是湍流l

4、化學(xué)反應(yīng)速率高度非線性; 湍流化學(xué)反應(yīng)高度耦合，相互作用很重要。l真實(shí)化學(xué)反應(yīng)機(jī)理包含數(shù)十個(gè)組分, 數(shù)百個(gè)基元反應(yīng)，并且方程組極具剛性 (基元化學(xué)反應(yīng)時(shí)間尺度相差大)u實(shí)際處理方法l簡(jiǎn)化化學(xué)反應(yīng)機(jī)理n有限速率燃燒模型l考慮湍流及其混合、弱化反應(yīng)化學(xué)n混合分?jǐn)?shù)模型s平衡化學(xué)的 PDF模型s層流火焰面模型n進(jìn)展變量模型sZimont 模型燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)有限速率模型u用總包機(jī)理反應(yīng)描述化學(xué)反應(yīng)過(guò)程.u求解化學(xué)組分輸運(yùn)方程.l求解當(dāng)?shù)貢r(shí)間平均的各個(gè)組分的質(zhì)量分?jǐn)?shù), mj.u組分 j的源項(xiàng) (產(chǎn)生或消耗）是機(jī)理中所有k個(gè)反應(yīng)的凈反應(yīng)速率 : uRjk (第k 個(gè)

5、化學(xué)反應(yīng)生成或消耗的j 組分)是根據(jù) Arrhenius速率公式、混合或渦旋破碎（EBU）速率的小值。. u混合速率與渦旋壽命相關(guān), k /.l物理意義是湍流渦旋是決定化學(xué)反應(yīng)的首要因素。對(duì)于非預(yù)混燃燒，湍流渦旋決定了組分混合；對(duì)于預(yù)混燃燒湍流決定了熱輸運(yùn)（高溫加熱低溫）。即：化學(xué)反應(yīng)決定于湍流混合組分（非預(yù)混燃燒）和熱量（預(yù)混燃燒）的速率。RRjjkk燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)Fluent燃燒模型u有限速率模型有限速率模型 求解組分的質(zhì)量分?jǐn)?shù)輸運(yùn)方程，化學(xué)反應(yīng)機(jī)理由用戶自己定義。u非預(yù)混燃燒模型非預(yù)混燃燒模型 該模型中并不求解單個(gè)組分的輸運(yùn)方程，而是求解一個(gè)或者

6、兩個(gè)守恒標(biāo)量（混合分?jǐn)?shù)）的輸運(yùn)方程 u預(yù)混燃燒模型預(yù)混燃燒模型 模擬完全混合的燃燒問(wèn)題。充分混合的燃燒物和產(chǎn)物被火焰前鋒分隔，求解出的化學(xué)反應(yīng)進(jìn)展變量來(lái)描述該火焰前峰的位置 u部分預(yù)混燃燒模型部分預(yù)混燃燒模型 該模型用來(lái)處理系統(tǒng)中同時(shí)具有非預(yù)混和充分預(yù)混的情況。該方法同時(shí)求解了混合分?jǐn)?shù)和反應(yīng)進(jìn)展變量 燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)有限化學(xué)反應(yīng)速率模型設(shè)置u要求:l給出組分及其物性l給出化學(xué)反應(yīng)及其反應(yīng)速率在內(nèi)的化學(xué)反應(yīng)動(dòng)力學(xué)數(shù)據(jù)uFLUENT V5 在mixture material database里面提供了數(shù)據(jù)u對(duì)于常用的燃料，數(shù)據(jù)庫(kù)都會(huì)給定機(jī)理，組分物性等信息.

7、 u如果用戶需要給定個(gè)性化機(jī)理，則: l生成新的混合物.l改變已有混合物的物性/化學(xué)反應(yīng).燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)有限速率模型小節(jié)u優(yōu)點(diǎn):l可以應(yīng)用于nonpremixed, partially premixed和premixed combustionl簡(jiǎn)單、直觀l應(yīng)用廣泛u缺點(diǎn):l不適合混合速率與化學(xué)反應(yīng)動(dòng)力學(xué)時(shí)間尺度相當(dāng)時(shí)候的化學(xué)反應(yīng) (要求 Da 1).l沒(méi)有嚴(yán)格考慮湍流化學(xué)反應(yīng)之間的相互作用問(wèn)題l不能考慮中間產(chǎn)物或組分、不能考慮分裂影響.l模型常數(shù)不確定, 特別是用于計(jì)算多個(gè)化學(xué)反應(yīng)的時(shí)候尤為如此，模型常數(shù)通用性較差。.燃燒熱科學(xué)與能源工程系2003年1

8、0月計(jì)算流體與傳熱傳質(zhì)守恒標(biāo)量 (混合物分?jǐn)?shù)) 模型: PDF 模型u只適應(yīng)用于非預(yù)混 (擴(kuò)散) 火焰燃燒u假定化學(xué)反應(yīng)過(guò)程受混合速率控制l滿足局部化學(xué)平衡.l控制體（計(jì)算單元）組分、物性決定于燃料和氧化劑在該處的混合程度. u化學(xué)反應(yīng)機(jī)理不明確. l用化學(xué)平衡計(jì)算來(lái)處理化學(xué)反應(yīng) (prePDF).u只求解混合物分?jǐn)?shù)及其方差的輸運(yùn)方程, 無(wú)需求解組分的輸運(yùn)方程.u可以嚴(yán)格考慮湍流與化學(xué)反應(yīng)的相互作用燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)混合分?jǐn)?shù)定義u混合分?jǐn)?shù), f, 寫(xiě)成元素的質(zhì)量分?jǐn)?shù)形式:l其中， Zk 是元素k的質(zhì)量分?jǐn)?shù) ；下標(biāo) F 和O 表示燃料和氧化劑進(jìn)口流處的值。u

9、對(duì)于簡(jiǎn)單的 fuel/oxidizer系統(tǒng), 混合物分?jǐn)?shù)代表計(jì)算控制體里的燃料質(zhì)量分?jǐn)?shù).u混合物分?jǐn)?shù)是守恒標(biāo)量（conserved scalar）: l組分輸運(yùn)方程中沒(méi)有化學(xué)反應(yīng)源項(xiàng).OkFkOkkZZZZf,燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)可以用單個(gè)混合物分?jǐn)?shù)模擬的燃燒系統(tǒng)uFuel/air 擴(kuò)散火焰:u多氧化劑入口的擴(kuò)散火焰:u多燃料進(jìn)口的擴(kuò)散火焰:60% CH4 40% CO21% O2 79% N2f = 1f = 035% O2 65% N260% CH4 40% CO35% O2 65% N2f = 1f = 0f = 060% CH4 20% CO 10%

10、 C3H8 10% CO221% O2 79% N2f = 1f = 0f = 160% CH4 20% CO 10% C3H8 10% CO2燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)系統(tǒng)化學(xué)平衡假設(shè)u化學(xué)反應(yīng)很快到達(dá)平衡. u可以考慮中間組分.燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)PDF 模擬Turbulence-Chemistry相互作用uFluctuating mixture fraction is completely defined by its probability density function (PDF).up(V), the PDF, r

11、epresents fraction of sampling time when variable, V, takes a value between V and V + V.up(f) can be used to compute time-averaged values of variables that depend on the mixture fraction, f:lSpecies mole fractionslTemperature, densityp VVTTii( )lim1iipffdf()()01燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)PDF Model

12、Flexibilityu非絕熱系統(tǒng):lIn real problems, with heat loss or gain, local thermo-chemical state must be related to mixture fraction, f, and enthalpy, h.lAverage quantities now evaluated as a function of mixture fraction, enthalpy (normalized heat loss/gain), and the PDF, p(f).u第二守恒標(biāo)量:lFLUENT用第二守恒標(biāo)量可以模擬:nTw

13、o fuel streams with different compositions and single oxidizer stream (visa versa)nNonreacting stream in addition to a fuel and an oxidizernCo-firing a gaseous fuel with another gaseous, liquid, or coal fuelnFiring single coal with two off-gases (volatiles and char burnout products) tracked separate

14、ly燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)混合分?jǐn)?shù)/PDF模型小節(jié)u優(yōu)點(diǎn):l可以計(jì)算中間組分.l考慮分裂影響.l考慮湍流化學(xué)反應(yīng)之間作用.l無(wú)需求解組分輸運(yùn)方程（特別是多組分），簡(jiǎn)化計(jì)算量l性能好，經(jīng)濟(jì)u缺點(diǎn):l系統(tǒng)必須滿足（靠近）局部平衡.l不能用于可壓速或非湍流流動(dòng).l不能用于預(yù)混燃燒.燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)層流火焰面模型u用混合分?jǐn)?shù)和標(biāo)量耗散率來(lái)求解（絕熱）溫度、密度和組分等量。l對(duì)于混合分?jǐn)?shù) PDF 模型 (絕熱), 熱化學(xué)狀態(tài)只是 f 的函數(shù)lc 與當(dāng)?shù)貞?yīng)變率有關(guān)u把混合分?jǐn)?shù) PDF擴(kuò)展到模擬中度化學(xué)非平衡燃燒模擬中u用層流拉伸火焰

15、系綜來(lái)模擬湍流火焰, 對(duì)撞擴(kuò)散火焰2)/(xf c),(cfii燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)層流火焰面模型(2)u用指定概率密度函數(shù)（PDF） P(f,c) 的方法來(lái)決定層流火焰面系綜。假定的概率包括計(jì)算為：Pf (f) Pc (c), 其中， Pf (f) 用Beta 函數(shù)， Pc (c) 用delta函數(shù)u只適合絕熱系統(tǒng)（FLUENT V5）uImport strained flame calculationslprePDF or Sandias OPPDIF codeuSingle or multiple flameletslSingle:user speci

16、fied strain, alMultiple: strained flamelet library, 0 a aextinctionna=0 equilibriumna= aextinction is the maximum strain rate before flame extinguishesuPossible to model local extinction pockets (e.g. lifted flames)10 0)()(),(dfdPfPffiicccc燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)預(yù)混燃燒的Zimont模型u用單個(gè)過(guò)程變量來(lái)模擬熱化學(xué)過(guò)程, u

17、平均反應(yīng)速率, u湍流火焰?zhèn)鞑ニ俣? Ut, 根據(jù)貧燃預(yù)混燃燒推導(dǎo)得到，并考慮l預(yù)混燃料當(dāng)量比l湍流引起的火焰前鋒皺折和增厚l湍流拉伸引起的火焰前鋒淬熄l分子擴(kuò)散u適合絕熱燃燒,u非絕熱燃燒必須求解焓方程tcxu cxSccxRciiittic01RUccunburnttpadpppYYc/adunburntTcTcT)1(燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)稀疏相模型u在拉格朗日坐標(biāo)系下，求解顆粒、液滴、氣泡的軌道. l并與連續(xù)相（氣相）進(jìn)行熱、質(zhì)量和動(dòng)量的耦合求解。u稀疏相體積分?jǐn)?shù)必須 10%l質(zhì)量分?jǐn)?shù)可以比較高l不考慮顆粒之間作用、顆粒破碎u模擬湍流彌撒l隨機(jī)軌道（S

18、tochastic tracking）l顆粒云團(tuán)模型（Particle cloud (V5)） uRosin-Rammler or linear size distributionu非定常流動(dòng)中的顆粒軌道 (V5)u模擬顆粒分離, 噴霧烘干, 液體燃料或煤粉的燃燒過(guò)程。.計(jì)算連續(xù)相流動(dòng)場(chǎng)計(jì)算顆粒軌道更新連續(xù)相源項(xiàng)燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)uMonte-Carlo方法模擬湍流顆粒彌散 (discrete random walks)u顆粒運(yùn)動(dòng)計(jì)算中考慮氣體的平均速度及隨機(jī)湍流脈動(dòng)速度的影響。u每個(gè)軌道包含了一群具有相同特性的顆粒，如相同的初始直徑，密度等. u考慮湍流

19、彌散尤為重要l更復(fù)合真實(shí)物理過(guò)程，但計(jì)算量更大。l可以通過(guò)光滑源項(xiàng)、消除與氣相的耦合來(lái)強(qiáng)化計(jì)算穩(wěn)定性；顆粒彌散: 隨機(jī)軌道模型Coal particle tracks in an industrial boiler燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)顆粒彌散: 顆粒云團(tuán)模型u 用平均速度決定顆粒的平均軌道u假定該平均軌道為 3D多變量的 Gaussian 分布，計(jì)算顆粒偏離的范圍（3倍標(biāo)準(zhǔn)偏差）uRigorously accounts for inertial and drift velocitiesuA particle cloud is required for eac

20、h particle type (e.g. initial d, etc.)uParticles can escape, reflect or trap (release volatiles) at wallsuEliminates (single cloud) or reduces (few clouds) stochastic trackinglDecreased computational expenselIncreased stability since distributed source terms in gas phaseBUT decreased accuracy sincel

21、Gas phase properties (e.g. temperature) are averaged within cloud lPoor prediction of large recirculation zones燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)非定常流動(dòng)的顆粒軌道模擬uEach particle advanced in time along with the flowuFor coupled flows using implicit time stepping, sub-iterations for the particle tracking are perf

22、ormed within each time stepuFor non-coupled flows or coupled flows with explicit time stepping, particles are advanced at the end of each time step燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)煤/油燃燒模型u變化模型可以模擬煤或油燃燒llDroplet 油燃燒模擬lCombusting particle 煤燃燒模擬uSeveral devolatilization and char burnout models provided.lNot

23、e: These models control the rate of evolution of the fuel off-gas from coal/oil particles. Reactions in the gas (continuous) phase are modeled with the PDF or finite rate combustion model.Particle TypeDescriptionInertinert/heating or coolingDroplet (oil)heating/evaporation/boilingCombusting (coal)he

24、ating;evolution of volatiles/swelling;heterogeneous surface reaction燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)NOx ModelsuNOx 中含量最多的是一氧化氮 (NO). l引起煙霧l引起酸雨l導(dǎo)致臭氧損耗uNOx 排放有三個(gè)機(jī)理:l熱機(jī)理（Thermal Nox） - Zeldovich mechanism (氧化空氣中的 N)n高溫燃燒最明顯l快速機(jī)理（Prompt NOx ）- De Soete, Williams提出的經(jīng)驗(yàn)機(jī)理. n通常該機(jī)理對(duì)Nox排放量影響較小n富燃燃燒區(qū)域該機(jī)理有明顯作用l燃料

25、NOx 機(jī)理- De Soete, Williams等人提出的經(jīng)驗(yàn)機(jī)理.n在煤燃燒中， 燃料中含N高，并且燃燒溫度不是很高時(shí)起作用。uNOx 再燃 (V5)l在燃料富足的燃燒區(qū)域，NO會(huì)與碳?xì)浠衔锓磻?yīng)，使得NO 減少燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)煙（Soot） 模型u兩種煙模型:lOne-step 模型 (Khan and Greeves)n1個(gè)輸運(yùn)方程，煙質(zhì)量分?jǐn)?shù)lTwo-Step 模型(Tesner)nTransport equations for radical nuclei and soot mass fraction concentrationsuSoot

26、 生成采用經(jīng)驗(yàn)公式其中, C, pf, 和 F 是模型常數(shù), 燃料分壓力和當(dāng)量比。u采用Magnussen模型計(jì)算 Soot 燃燒 (消耗)uSoot 對(duì)輻射的影響l激活Soot面板上的Soot-Radiation 選項(xiàng)RTEnfformationepCR/F燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)燃燒模擬步驟及求解策略u(píng)打開(kāi)求解器（打開(kāi)求解器（2D或或3D）l檢查物理模型適用性l網(wǎng)格分辨率需求 (resolve shear layers)l求解參數(shù)、收斂條件設(shè)置u邊界條件邊界條件l燃燒經(jīng)常對(duì)進(jìn)口邊界條件十分敏感n正確的進(jìn)口標(biāo)量和速度的分布l很難估計(jì)壁面熱條件；盡可能優(yōu)先給定邊界溫度，很難給定對(duì)流換熱或者輻射換熱量 u初始條件初始條件l穩(wěn)態(tài)解跟初始值無(wú)關(guān)，但給的初始值不好，由于各個(gè)方程之間不和諧和輸運(yùn)方程的非線性，導(dǎo)致解分歧或發(fā)散不收斂。先求解冷態(tài)流動(dòng)、接著氣體燃燒、顆粒燃燒、再考慮輻射；l對(duì)于強(qiáng)旋流動(dòng)，漸漸增加旋流速度；燃燒熱科學(xué)與能源工程系2003年10月計(jì)算流體與傳熱傳質(zhì)燃燒模擬步驟及求解策略（2）u低松弛因子低松弛因子l松弛因子改變對(duì)結(jié)果影響是高度非線性的n采用混合分?jǐn)?shù)PDF模型時(shí)，密度采用低松弛 (0.5)n考慮浮力流動(dòng)中，速度采用低松弛因子n高速流動(dòng)中壓力采用低松弛因子l一旦解較為穩(wěn)定, 盡可能讓松弛因子提高到系統(tǒng)默