層流非預混擴散火焰課件_第1頁
層流非預混擴散火焰課件_第2頁
層流非預混擴散火焰課件_第3頁
層流非預混擴散火焰課件_第4頁
層流非預混擴散火焰課件_第5頁
已閱讀5頁,還剩329頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領

文檔簡介

九層流非預混擴散火焰

NonpremixedDiffusonal九層流非預混擴散火焰

Nonpremixed2022/11/192背景(Backgrounds)燃料的射流火焰;蠟燭的火焰;本生燈的外焰;碳煙(Soot)的形成。Jet

flameCandle

flameBunsenouterflame2022/11/102背景(Backgrounds)燃料的射2022/11/193預混火焰非預混火焰預混火焰和非預混燃燒的區(qū)別非預混火焰結構分三個區(qū):富燃料區(qū)、反應區(qū)和富氧區(qū).C.K.Law,Combustionphysics,2006火焰速度、火焰厚度火焰長度富燃料區(qū)反應區(qū)富氧區(qū)2022/11/103預混火焰非預混火焰預混火焰和非預混燃燒2022/11/194未燃常密度層流射流射流火焰的物理描述Burke-Schumann簡化理論描述圓口和狹縫燒嘴的火焰長度Soot形成及氧化分解對吹火焰的介紹內容2022/11/104未燃常密度層流射流內容2022/11/1951.未反應常密度層流射流NonreactingJet物理描述:燃料射入氧化劑中中心線上流速半徑方向的流速流體速度和密度等于出口值氣流核心:速度和燃料的質量分數(shù)不變且處處相等,等于噴嘴出口值。從該核到射流邊界,速度與燃料濃度均單調遞減到0。在氣流核心之外(x>xc),粘性剪切力和質量擴散在整個射流區(qū)域內都起作用。Vc為射流初始速度;Vx為射流方向x處徑向(r方向)某點的x方向速度2022/11/1051.未反應常密度層流射流Nonrea2022/11/19在整個流場里面,初始的射流動量是守恒的。當燃料流向周圍的空氣時,它的一部分動量就傳給了空氣,因此燃料的速度減小;同時在流場中越往下游,進入射流區(qū)域里的空氣量就越多。這可以用動量守恒的積分表達式來表示:速度場動量守恒2022/11/10在整個流場里面,初始的射流動量是守恒的。2022/11/197其中e

和ve

分別為噴管口燃料的密度和速度。9.1動量守恒射流出口初始動量9.1速度場2022/11/1079.1動量守恒射流出口初始動量9.1速2022/11/198影響速度場的是動量的對流和擴散,而影響燃料濃度場的則是物質的對流和擴散,這兩者具有一定的類似性,因此燃料的質量分數(shù)Yf(r,x)和無量綱速度vx(r,x)/ve也應該具有類似的分布規(guī)律(如圖9.1)。和初始射流動量類似,從噴嘴流入的燃料質量也是守恒的,即:9.2=1濃度場9.1質量守恒2022/11/108影響速度場的是動量的對流和擴散,而影響2022/11/199現(xiàn)在需要確定的就是方程9.1和9.2中的速度場和方程9.2中的燃料質量分數(shù)場的具體分布。

vx(r,x)Yf(r,x)9.19.22022/11/109現(xiàn)在需要確定的就是方程9.1和9.2中2022/11/1910假設1、射流和周圍流體的摩爾質量相等。有了這個假設,加上理想氣體性質,并設流場內壓力和溫度都是常數(shù),那么整個流場內流體的密度也就是常數(shù)。2、物質之間的擴散為遵從費克定律的簡單二元擴散。3、動量和組分的擴散率都是常數(shù),且相等,即施密特數(shù)[Sc=ν(動量擴散率)/D

(質量擴散率)]等于1。4、只考慮物質的徑向動量擴散,忽略軸向擴散。因此,下面得出的結論只在距離噴嘴出口下游一定距離的地方,也就是隱核以外的地方才適用,因為噴嘴出口處隱核內部軸向擴散起著較大的作用,不能忽略。2022/11/1010假設1、射流和周圍流體的摩爾質量相等2022/11/1911守恒定律方程邊界層方程:(第七章)質量守恒9.3單位體積內軸向對流引起的質量流量增量單位體積內徑向對流引起的質量流量增量對比7.7式,消去了密度ρ(已假設密度為常數(shù))7.72022/11/1011守恒定律方程邊界層方程:(第七章)92022/11/1912守恒定律方程9.4兩邊除以密度ρ,動力粘度μ變?yōu)閯恿繑U散率ν單位體積內軸向對流引起的x方向動量增加速率單位體積內徑向對流引起的x方向動量增加速率單位體積內所受的粘性力單位體積內所受的浮力(可忽略)軸向動量守恒(徑向動量之和為0)7.482022/11/1012守恒定律方程9.4兩邊除以密度ρ,動2022/11/1913組分守恒:針對燃料射流(無反應)由于組分只有燃料和氧化劑兩種,所以兩者的質量分數(shù)相加應該為1。9.5單位體積內軸向對流引起的燃料的質量流量單位體積內徑向對流引起的燃料的質量流量單位體積內徑向擴散引起的燃料的質量流量7.202022/11/1013組分守恒:針對燃料射流(無反應)92022/11/1914邊界條件為了求解未知的vx(r,x),vr(r,x),YF(r,x)在軸線上:(r=0)r=0r=軸對稱軸對稱2022/11/1014邊界條件為了求解未知的vx(r,x2022/11/1915在徑向無窮遠處(r),流體速度為0,且沒有燃料:在射流出口(x=0處),假設噴嘴內部(r≤R)的軸向速度和燃料質量分數(shù)都相等,而在噴嘴外部其值均為零,即:2022/11/1015在徑向無窮遠處(r),流體速度為2022/11/1916x=0處2022/11/1016x=0處2022/11/1917速度場求解求解速度場可以通過相似理論來實現(xiàn),相似性就是指速度的內在分布在流場內的每個地方都是一樣的。對于現(xiàn)在討論的情況來說,這就意味著用中心線速度vx(0,x)將相同軸向距離處速度vx(r,x)無量綱化后,得到的無量綱速度僅僅是相似變量r/x的函數(shù)。2022/11/1017速度場求解求解速度場可以通過相似理論2022/11/1918其中Je是初始的動量流率(單位),0ξ里面包含著相似變量r/x,9.11速度場求解2022/11/100ξ里面包含著相2022/11/1919將(9.10)代入(9.8)并加以整理,即可得到軸向速度分布的無量綱形式:9.12此方程反映了任意x和r處的軸向速度與初始速度ve的比值速度場求解2022/11/10199.12此方程反映了任意x和r處的軸2022/11/1920再令r=0(ξ=0),即可得到中心線速度的衰減關系式:從(9.12)中我們可以看出,該無量綱速度vx,0/ve和軸向距離x成反比,和射流的雷諾數(shù)(Rej=eveR/

)成正比。另外對于(9.13),在靠近噴嘴的地方,因為x/R將變得很小,其倒數(shù)將變得很大,vx,0/ve會大于1,因此這個解不再適用。9.139.122022/11/1020再令r=0(ξ=0),即可得2022/11/1921由(9.13)作出的中心線速度變化曲線如圖9.2所示,從圖中可以明顯地看到,射流的雷諾數(shù)越低,中心線速度的衰減就越快。這是由于雷諾數(shù)減小時,射流的初始動量和粘性剪切力作用相比,其重要性變小。9.132022/11/1021由(9.13)作出的中心線速度變化曲2022/11/1922其它常用的射流參數(shù)有擴張率和噴射角α。我們先引入射流半寬r1/2概念:在射流的某一軸向距離處,當射流速度減小到該軸向距離處中心線速度一半時的徑向距離為此軸向距離處的射流半寬,如圖9.3所示。

射流半寬的表達式可以通過聯(lián)立(9.12)和(9.13),并令vx/vx,0

=1/2得到。9.139.122022/11/1022其它常用的射流參數(shù)有擴張率和噴射角α2022/11/1923有了射流半寬這一概念,我們就可以這樣來定義擴張率和噴射角:擴張率是射流半寬和軸向距離的比值,正切值等于擴張率的角度即為噴射角,其表達式分別如下所示:及9.14令上式左邊等于1.287,r變?yōu)閞1/2由結論:高雷諾數(shù)射流比低雷諾數(shù)射流窄。2022/11/1023有了射流半寬這一概念,我們就可以這樣2022/11/1924下面進行濃度場的求解。對比(9.4)和(9.5)可以看出,如果Sc=ν/D=1,那么燃料質量分數(shù)YF和無量綱速度vx/ve的表達式就完全相同。比較一下其中QF是燃料在噴嘴處的體積流率(QF=veR2).濃度場9.8濃度場求解2022/11/1024下面進行濃度場的求解。對比(9.4)2022/11/1925對于(9.16),因Sc=1,進一步簡化得到以射流雷諾數(shù)為參數(shù)的表達式:中心線(r=0)的質量分數(shù)衰減關系式:9.179.18濃度場求解2022/11/1025對于(9.16),因Sc=1,進一步2022/11/1926同前面一樣,這個解只在距離噴嘴一定距離以外地方才適用,因此這個范圍的無量綱軸向距離和雷諾數(shù)的關系為:由于(9.18)和(9.13)是一樣的,因此圖9.2也就是中心線質量分數(shù)的衰減曲線。因為YF,0須小于等于1,即9.189.13濃度場求解2022/11/1026同前面一樣,這個解只在距離噴嘴一定距2022/11/1927例題9.1乙烯射流從直徑10毫米的圓管射入300K,1atm的靜止空氣中。設初始射流速度為10cm/s和1.0cm/s,比較兩種速度下噴射角以及當射流中心線上的質量分數(shù)等于化學當量值時的高度。乙烯在300K,1atm下的粘度為102.3×10-7N-s/m2.略2022/11/1027例題9.1乙烯射流從直徑10毫米的2022/11/1928答案低速射流比較寬,擴張角是高速射流的9倍低速射流的燃料濃度衰減到與高速射流相同的值時,其軸向距離只是高速射流的1/10略2022/11/1028答案低速射流比較寬,擴張角是高速射流2022/11/1929例題9.2根據(jù)例題9.1中的情況2(ve=1.0cm/s,R=5mm),如果出口速度上升到10cm/s,試求為了保持流量,出口半徑需要變?yōu)槎啻蟆M瑫r,試求YF,0=YF,stoic的軸向位置并與情況2比較。答案A略2022/11/1029例題9.2根據(jù)例題9.1中的情況22022/11/1930燃料質量分數(shù)的空間分布取決于初始質量流率(體積流率),而非速度答案B略2022/11/1030燃料質量分數(shù)的空間分布取決于初始質量2022/11/19312.射流火焰的物理描述層流射流的燃燒情況在很大程度上和前面討論的等溫射流的一樣的。燃料一邊沿著軸向流動一邊快速向外擴散,同時氧化劑(如空氣)迅速向內擴散。概念:射流火焰的火焰面2022/11/10312.射流火焰的物理描述層流射流的燃2022/11/1932在流場中,燃料和氧化劑之比為化學當量比的點就構成了火焰表面。層流非預混火焰(一):火焰表面定義火焰表面2022/11/1032層流非預混火焰(一):火焰表面定義火2022/11/1933在流場中,燃料和氧化劑之比為化學當量比的點就構成了火焰表面。層流非預混火焰(一):火焰表面火焰表面2022/11/1033層流非預混火焰(一):火焰表面火焰表2022/11/1934在這里需要注意的是,雖然燃料和氧化劑在火焰處都消耗了,但是產(chǎn)物的組成成分只和Φ的取值有關,因此射流火焰初始的當量比仍然有意義。產(chǎn)物在火焰表面形成后,就向著內、外側快速擴散?;鹧姹砻?當量比Ф等于1的點的軌跡層流非預混火焰(一):火焰表面定義2022/11/1034在這里需要注意的是,雖然燃料和氧化劑2022/11/1935對于富氧燃燒,周圍存在著過量的氧化劑,火焰長度Lf可以這樣定義:層流非預混火焰(二):火焰長度2022/11/1035對于富氧燃燒,周圍存在著過量的氧化劑2022/11/1936在整個火焰中,發(fā)生化學反應的區(qū)域通常來說是很窄的,就像圖9.4中看到的一樣,在到達火焰頂部以前,高溫的反應區(qū)是一個環(huán)形的區(qū)域。這個區(qū)域可以通過一個簡單的實驗顯示出來,即在本生燈的火焰前面垂直于軸線放置一個金屬濾網(wǎng),在火焰區(qū)對應的地方濾網(wǎng)就會受熱而發(fā)光,就可以看到這種環(huán)形的結構。層流非預混火焰(二):火焰長度定義2022/11/1036在整個火焰中,發(fā)生化學反應的區(qū)域通常2022/11/1937在垂直火焰的上部,由于氣體較熱,浮力的作用就不能不考慮了。浮力的存在在加快氣體流動的同時,也使火焰變窄。由質量守恒定律我們知道,當速度變大時,流體的流線將變得彼此靠近,但也就導致了燃料的濃度梯度dYF/dr的增加,也增強了擴散作用。浮力作用和擴散作用對圓噴嘴火焰長度的影響相互抵消[6,7],因此本章后面推導的簡化理論雖然忽略了浮力,也能夠合理地計算出圓口或方口射流的火焰長度。層流非預混火焰(三):浮力的影響2022/11/1037在垂直火焰的上部,由于氣體較熱,浮力2022/11/1938在碳氫化合物的燃燒火焰中,由于經(jīng)常會有碳黑存在,火焰就可能呈現(xiàn)為橙色或黃色。層流非預混火焰(四):碳煙Soot的產(chǎn)生2022/11/1038在碳氫化合物的燃燒火焰中,由于經(jīng)常會2022/11/1939碳煙在反應區(qū)的燃料側形成,并在流向氧化區(qū)過程中不斷被氧化、消耗。形成條件:高溫、缺氧。圖9.5層流射流火焰中碳煙的形成區(qū)和消耗區(qū)2022/11/1039碳煙在反應區(qū)的燃料側形成,并在流向氧2022/11/1940由于燃料和火焰停留時間的不同,在燃料側形成的碳煙在向高溫氧化區(qū)移動的過程中可能無法被完全氧化。在這種情況下,soot就會沖出火焰而形成碳黑的“翼”,這部分從火焰中出來的碳黑就是我們通常說的說的煙。2022/11/1040由于燃料和火焰停留時間的不同,在燃料41圖9.6是一個乙烯火焰的照片,在圖中可以看到焰舌的右邊出現(xiàn)了碳黑翼。R.J.Santoro提供(PSU)41圖9.6是一個乙烯火焰的照片,在圖中可以看到焰舌的右邊出42R.J.Santoro實驗室,PSU42R.J.Santoro實驗室,PSU2022/11/1943對于圓口火焰來說,火焰長度和初始速度以及管徑都無關,但是和初始體積流量QF有關。由于QF=veπR2,不同ve和R的組合也可以得到相同的火焰長度。這一點的合理性可以從非反應層流射流中應用前面的分析結果(見例9.2)中得到驗證。層流非預混火焰(五):火焰長度-流量的關聯(lián)在層流射流非預混火焰中,還有一個突出的特點,就是火焰長度和初始條件之間的關系。2022/11/1043對于圓口火焰來說,火焰長度和初始速度2022/11/1944對前面得出的(9.16),如果忽略反應放熱,并將其中的YF改為YF,stoic,代入(9.16)就可以得到火焰表面軌跡的粗略描述方程。如果再令r=0,就可以得到火焰長度為:9.16Lf相當于9.16中的xr=0,則ξ=09.222022/11/1044對前面得出的(9.16),如果忽略反2022/11/1945由此可以看出,火焰長度確實是和體積流量成正比而且還和燃料的化學當量質量分數(shù)成反比。這就意味著如果燃料完全燃燒需要的空氣越少,那么燃燒的火焰也就越短。

越大9.222022/11/1045由此可以看出,火焰長度確實是和體積流2022/11/19463.簡化(數(shù)學)理論描述最早的關于層流噴射擴散火焰的理論描述是Burke和Schumann[4]在1928年發(fā)表的。雖然做了很多假設,例如認為速度場在每個地方的分布都是恒定,并且平行于火焰軸,他們的理論也能夠較為合理地對軸對稱(如圓口)火焰的火焰長度進行求解。在此之后,其他的研究者對其理論進行擴展提煉,但一直保留著恒速這一假設。2022/11/10463.簡化(數(shù)學)理論描述最早的關于2022/11/19471977年,Roper發(fā)表了他的理論學說,其中保留了Burke-Schumann理論中的基本簡化,但是去掉了恒速這一條件,他的結論給出了圓口和方口噴嘴燃燒的火焰長度的合理解。9.5節(jié)將將推導這個適用于工程計算的結論。2022/11/10471977年,Roper發(fā)表了他2022/11/1948首先,我們先來看看問題的數(shù)學描述,以便了解問題固有的難度,從而體會簡化理論的好處。先給出的是一些通用方程,其中包含了速度、質量分數(shù)和溫度這些大家熟悉的變量;然后根據(jù)合理的假設得到用守恒量來表示的關系式,再對這兩個偏微分方程求解就可以得到該問題的解了。層流射流擴散火焰的Burke和Schumann數(shù)學描述2022/11/1048首先,我們先來看看問題的數(shù)學描述,以2022/11/19493-1基本假設1、流動為穩(wěn)定的軸對稱層流,燃料由半徑為R的圓形噴嘴噴出,在靜止的、無限大的、充滿氧化劑的空間里燃燒。2、只有燃料、氧化劑和產(chǎn)物三種物質存在,火焰面內部只存在燃料和產(chǎn)物,火焰外部只存在氧化劑和產(chǎn)物。3、火焰表面,燃料和氧化劑按化學當量比進行反應?;瘜W動力學速度無限快,意味著火焰面只存在于一個無限薄的薄層里;這就是通常說的“火焰面近似”

。4、物質間的擴散為遵從費克定律的簡單二元擴散。5、熱擴散率和質量擴散率相等,即路易斯數(shù)(Le=α/D)等于1。6、忽略輻射換熱。7、只考慮徑向的動量、熱和物質擴散,而忽略軸向的各種擴散;8、火焰的軸線垂直向上。2022/11/10493-1基本假設1、流動為穩(wěn)定的軸對2022/11/1950“火焰面近似”(Flame-sheetapproximation)實際的火焰面火焰面近似2022/11/1050“火焰面近似”(Flame-shee2022/11/19513-2基本守恒方程質量守恒軸向動量守恒這個方程適用于整個空間,即火焰面內部和外部都適用,并且在火焰面處保持連續(xù)。9.24(7.7)(7.48)9.232022/11/10513-2基本守恒方程質量守恒9.242022/11/1952組分守恒:根據(jù)火焰面假設,火焰面內、外的產(chǎn)物質量產(chǎn)率()均為零(反應區(qū)體積為0)。所有的化學反應現(xiàn)象都體現(xiàn)在邊界條件里面,因此(7.20)可以寫作:產(chǎn)物的質量分數(shù)為:單位體積內化學反應引起的A組分的質量凈生成率9.25(7.20)2022/11/1052組分守恒:根據(jù)火焰面假設,火焰面2022/11/1953能量守恒(火焰面的前、后,可與組分守恒方程一樣簡化)這個方程在火焰面處不連續(xù)(有能量產(chǎn)生)?;鹧婷嫔戏磻懦龅臒崃浚ㄔ错棧┛梢赃吔鐥l件的形式給出。因此,反應放熱是邊界條件,其中邊界為火焰表面?;鹧婷媲昂?,源項等于0;上式僅適用于火焰面前后。單位體積內化學反應引起生成焓變化速率軸向對流引起的單位體積顯焓輸運率徑向對流引起的單位體積顯焓輸運率軸向擴散引起的單位體積顯焓輸運率2022/11/1053能量守恒(火焰面的前、后,可與組分守2022/11/1954附加關系和方程歸結狀態(tài)方程:5個方程:質量,動量,燃料,氧化劑和能量5個未知的函數(shù):聯(lián)立求解5個偏微分方程,任務艱巨,需要知道火焰面的邊界條件來求解燃料和氧化劑的組分守恒方程和能量守恒方程。但不知道火焰面的位置。引入守恒標量:只需要火焰軸線、遠離火焰和噴嘴出口的邊界條件即可。2022/11/1054附加關系和方程歸結狀態(tài)方程:5個方程2022/11/1955守恒標量的概念討論其中的兩個:一個是下面定義的混合物分數(shù)(f);一個是前面提到的混合物絕對焓(h)。略過2022/11/1055守恒標量的概念討論其中的兩個:略過2022/11/1956混合分數(shù):ff:在燃料中為1,在氧化劑中為0!!!推導過程Q-2略過2022/11/1056混合分數(shù):ff:在燃料中為1,在氧化2022/11/1957略過2022/11/1057略過2022/11/1958略過組分守恒方程變?yōu)?022/11/1058略過組分守恒方程變?yōu)?022/11/1959能量守恒方程變?yōu)樽钤际阶勇赃^絕對焓:h2022/11/1059能量守恒方程變?yōu)樽钤际阶勇赃^絕對焓2022/11/1960建立的4個標量守恒方程a.混合物分數(shù)守恒方程(New)b.絕對焓守恒方程(New)c.連續(xù)性方程d.動量守恒方程略過2022/11/1060建立的4個標量守恒方程a.混合物分數(shù)2022/11/1961(a)守恒標量的推導混合分數(shù):f的邊界條件因為在火焰處f=fstoic,所以只要知道了

f(r,x)的分布,也就能得到火焰的位置。略過2022/11/1061(a)守恒標量的推導混合分數(shù):因為2022/11/1962(b)絕對焓守恒方程能量方程可以寫為:和混合分數(shù)一樣,h在火焰面也連續(xù),其邊界條件如下:如何構造略過2022/11/1062(b)絕對焓守恒方程能量方程可以寫2022/11/1963(c-d)質量和動量守恒方程質量守恒軸向動量守恒這個方程適用于整個空間,即火焰面內部和外部都適用,并且在火焰片處保持連續(xù)。略過2022/11/1063(c-d)質量和動量守恒方程質量守恒2022/11/1964我們只需要用守恒量來代替物質組分方程和能量方程,這對連續(xù)方程和軸向動量方程沒有影響,因此不對前面得到的(9.23)和(9.24)做變換。速度的邊界條件和非反應射流的邊界條件一樣:

最后,要確定密度:(r,x)?。。÷赃^2022/11/1064我們只需要用守恒量來代替物質組分方程2022/11/1965建立的無量綱標量守恒方程略過2022/11/1065建立的無量綱標量守恒方程略過2022/11/19663-4無量綱方程在這里,我們使用噴嘴半徑R為特征長度,噴嘴出口速度ve為特征速度來定義下面的無量綱空間坐標和速度:略過2022/11/10663-4無量綱方程在這里,我們使用噴2022/11/1967混合分數(shù)f(0≤f≤1)本身就是一個無量綱變量,可以直接使用。而混合物的絕對焓h具有量綱,對此則做如下定義:注意,在燒嘴出口,h=hF,e

且h*=1;而在環(huán)境中(r),h*=0略過2022/11/1067混合分數(shù)f(0≤f≤1)本身就是一個2022/11/1968再定義如下的無量綱密度,來使控制方程完全無量綱化:其中e是燒嘴出口的燃料密度略過2022/11/1068再定義如下的無量綱密度,來使控制方程2022/11/1969將這些無量綱變量和參數(shù)與其對應的物理量聯(lián)系起來,再代入前面的基本守恒方程里面,就可以得到下面的無量綱控制方程:連續(xù)性軸向動量雷諾數(shù)和Freud數(shù)略過2022/11/1069將這些無量綱變量和參數(shù)與其對應的物理2022/11/1970混合分數(shù)無量綱焓略過2022/11/1070混合分數(shù)無量綱焓略過2022/11/1971上述方程無量綱的邊界條件為:中心線R無窮遠中心線對稱出口略過2022/11/1071上述方程無量綱的邊界條件為:中心線R2022/11/1972通過觀察這些無量綱控制方程和邊界條件,我們可以得到一些有意思的結論。首先,混合分數(shù)和無量綱焓的方程以及邊界條件是一樣的形式,就是說f和h*在他們各自的控制方程里的所處的地位是一樣的。因此只需要對(9.38)和(9.39)中的一個進行求解就可以了,例如已經(jīng)解出了f(r*,x*),那么就可以認為h*(r*,x*)=f(r*,x*)。略過2022/11/1072通過觀察這些無量綱控制方程和邊界條件2022/11/1973附加假設如果忽略浮力的作用,那么軸向動量方程(9.37)的右邊項就為零,這個方程和混合分數(shù)以及無量綱焓方程相對比,就只有其中的μ和后者中的ρD是不一樣的了。如果再假設μ和ρD相等,那么問題就可以得到進一步的簡化。施密特數(shù)Sc定義如下:如果μ=ρD,那么即施密特數(shù)為1(Sc=1)。前面已經(jīng)做了熱擴散率和質量擴散率相等的假設(Le=1)。略過2022/11/1073附加假設如果忽略浮力的作用,那么軸向2022/11/1974忽略了浮力的作用,并做了Sc=1這個假設,前面的軸向動量、混合分數(shù)(物質)和焓(能量)方程(9.37~9.39)就可以用下面這個統(tǒng)一的式子來表達:式中的通用變量=vx*=f=h*,雷諾數(shù)Re=eveR/

。略過2022/11/1074忽略了浮力的作用,并做了Sc=12022/11/1975狀態(tài)關系式為了對上述射流火焰的問題進行求解,需要將無量綱密度ρ*(=ρ/ρe)和混合分數(shù)或其它任何一個守恒量聯(lián)系起來。為此將引入理想氣體狀態(tài)方程(9.28),但是這又必須知道各組分的質量分數(shù)和溫度,因此接下來要做的就是將各組分的Yi和T表示為混合分數(shù)的函數(shù),然后就可以得到所需要的函數(shù)關系式ρ=ρ(f)了。對于我們考慮的簡單系統(tǒng),在火焰內部只有燃料和產(chǎn)物,火焰外只有氧化劑和產(chǎn)物(參見假設2),需要確定的是下面的狀態(tài)關系式:略過2022/11/1075狀態(tài)關系式為了對上述射流火焰的問題進2022/11/1976略過2022/11/1076略過2022/11/1977由火焰面的假設(假設3),對于火焰內、火焰面處和火焰外的YF、YOx和YPr,都可以用混合分數(shù)的定義將其同f聯(lián)系起來,如圖9.7所示,其關系為:略過2022/11/1077由火焰面的假設(假設3),對于火焰內2022/11/1978火焰內部(fstoic<f1)從定義式出發(fā)略過2022/11/1078火焰內部(fstoic<f1)從2022/11/1979在火焰面上(f=fstoic)略過2022/11/1079在火焰面上(f=fstoic)略過2022/11/1980火焰外部(0f<fstoic)略過2022/11/1080火焰外部(0f<fstoic)略過2022/11/1981由上面的結論可知所有組分的質量分數(shù)都與混合分數(shù)成線性關系,如圖9.8所示。略過2022/11/1081由上面的結論可知所有組分的質量分數(shù)都2022/11/1982為了將混合物的溫度表示為混合物分數(shù)f的函數(shù),還需要引入熱量狀態(tài)方程(2.4)。像前面幾章那樣,這里也應用Spalding理論,并做下面的假設:1、所有組分的比熱均為常數(shù),并且彼此相等,即cp,F=cp,Ox=cp,

Pr

≡cp。2、氧化劑和產(chǎn)物的生成焓均為零,即hof,

Ox=hof,Pr=0。這就使得燃料的生成焓和燃燒熱相等。狀態(tài)關系式-溫度分布略過2022/11/1082為了將混合物的溫度表示為混合物分數(shù)f2022/11/1983將(9.48)代入無量綱焓h*的定義式(9.35e)中,并利用控制方程中h*=f這一相似性,可以得到:略過2022/11/1083略過2022/11/1984式中用到了和這兩個焓定義式。通過求解(9.49),再注意到YF也是混合分數(shù)f的函數(shù),就可以得到關于T的狀態(tài)關系式:略過2022/11/1084式中用到了2022/11/1985將YF在火焰內、火焰面處和火焰外的表達式(9.44a、9.45a和9.46a)代入(9.50),就可以得到下面的結論:火焰內部(fstoic<f1)略過2022/11/1085將YF在火焰內、火焰面處和火焰外的表2022/11/1986在火焰面上(f=fstoic)火焰外部(0f<fstoic)略過2022/11/1086在火焰面上(f=fstoic)略過2022/11/1987溫度在火焰面內部和外部均為線性分布,并在火焰面處達到最大值略過2022/11/1087溫度在火焰面內部和外部均為線性分布,2022/11/1988得到無量綱方程組連續(xù)性軸向動量混合分數(shù)無量綱焓2022/11/1088得到無量綱方程組連續(xù)性軸向動量2022/11/19894各種不同的解法4.1

Burke-Schumann解;4.2常密度近似解;4.3變密度近似解;4.4

Ropper的近似解;4.5數(shù)值解(不做要求)4.6

P.Sutherland的最新工作P.B.Sunderland,et.al,Combust.Flame116(1999)376–386;Combust.Flame152(2008)60–68針對無量綱方程組的不同解法2022/11/10894各種不同的解法4.1Burke2022/11/1990在忽略浮力作用的基礎上,再加上速度不變這個假設,就可以去掉軸向動量方程(9.24)了。雖然那個時候還沒有正式地提出守恒標量這個概念,但他們對物質方程的處理后也得到了和守恒標量方程一樣的形式。4.1Burke-Schumann的簡化解Burke和Schumann最早得到了層流圓管射流火焰問題的近似解[4]。對于軸對稱和二維模型,他們都使用了火焰面近似,并簡單地認為流體的速度不變,即vx=v,vr

=0。9.242022/11/1090在忽略浮力作用的基礎上,再加上速度不2022/11/1991由于vr=0,質量守恒方程(9.23)決定了ρvx為常數(shù),4.1Burke-Schumann的簡化解9.230因此變密度的組分守恒方程(9.25)就可以寫作:2022/11/1091由于vr=0,質量守恒方程(9.2022/11/1992這個方程里面不包含組分生成源項,求解還需要知道火焰邊界。為了回避火焰邊界的問題,Burke和Schumann將燃料質量分數(shù)的定義擴展到整個流場里面,認為燃料質量分數(shù)在燃料流中內取1,火焰面處取0,純氧化劑中取-1/v,或–fstoic/(1–fstoic)。這樣,在火焰外面就出現(xiàn)了為負值的燃料濃度。9.25ν為化學當量系數(shù)9.52CxHy+vOx=(v+1)Pr0火焰面上當量燃燒:2022/11/1092這個方程里面不包含組分生成源項,求解2022/11/1993定義:混合分數(shù)f(無量綱數(shù))f:規(guī)定在燃料中為1,在氧化劑中為0!!!燃料、產(chǎn)物、氧化劑中源于燃料的各種元素質量之和混合物=燃料+產(chǎn)物+氧化劑2022/11/1093定義:混合分數(shù)f(無量綱數(shù))f:規(guī)定2022/11/1994將這個定義式代入(9.52)中,就可以得到前面的f守恒標量方程(9.30)。9.53火焰面處氧化劑的質量分數(shù)Yox=0火焰面處:9.302022/11/10949.53火焰面處氧化劑的質量分數(shù)Yo2022/11/1995用ρrefDref替換(9.52)中的ρD,并將其從微分符號里面提出,再注意到ρvx=constant=ρref

vx,ref,就可以從中消去ρref從而得到下面這個式子:進一步假設密度和質量擴散系數(shù)的乘積為常數(shù),即ρD=constant=ρrefDref。9.549.522022/11/1095用ρrefDref替換(9.52)中2022/11/1996這個式子里面,J0和J1分別為第零階和第一階貝塞耳函數(shù),可參見相應的數(shù)學參考書(如參考書19);λm為方程J1(λmR0)=0的所有正根[19];R和R0分別為燃料流動和外部流動的半徑;S在噴嘴和外部分別為燃料和氧化劑的化學當量摩爾比。這個偏微分方程的解YF(x,r)的表達式較復雜,并含貝塞耳函數(shù)?;鹧骈L度Lf并不能直接解出,需求解下面的超越方程:2022/11/1096這個偏微分方程的解YF(x,r)的2022/11/1997由于做了一些近似的假設,Burke和Schumann理論得出的火焰長度和圓管燃燒器理論得出的結果是基本一致的。Roper[6]提到過浮力的存在會使火焰變窄,使擴散作用增強。其實Burke和Schumann也注意到了這種現(xiàn)象。Kee和Miller[10]對有無浮力作用的情況做了數(shù)值分析并進行了對比,也精確地給出了浮力所帶來的影響。Burke-Schumann的簡化解總結Sunderland,etal(1999)2022/11/1097由于做了一些近似的假設,Burke和2022/11/19984.2常密度解法如果認為燃燒射流流體的密度為常數(shù),那么方程組(9.23)、(9.24)和(9.30)的解就和未反應射流(方程組9.3-9.5)的解一樣。此時,火焰長度由(9.22)給出:2022/11/10984.2常密度解法如果認為燃燒射流流體2022/11/19994.3變密度近似解Fay[5]對變密度的層流噴射火焰問題進行了求解。在他的求解中,忽略了浮力,因此軸向動量方程得到了化簡。對于熱物性參數(shù),假設施密特數(shù)和路易斯數(shù)都為1,這和建立控制方程過程中的假設是一致的,并設絕對粘度μ和溫度成正比,即:略2022/11/10994.3變密度近似解Fay[5]對變2022/11/19100Fay給出的火焰長度的解為:其中是噴嘴的質量流率,是遠離火焰處環(huán)境流體的密度,I(/f)是Fay解里面通過數(shù)值積分得到的函數(shù)。略2022/11/10100略2022/11/19101對于不同的環(huán)境與火焰密度比/f

下的/ref和I(/f)的值都列與表9.2注意到和則(9.57)可以寫成和常密度的解(9.56)相似的形式:對比略2022/11/10101對于不同的環(huán)境與火焰密度比2022/11/19102變密度層流射流火焰動量積分近似Ambient-to-flame變密度理論的求解結果比常密度理論的解大,二者之間相差倍數(shù)為:對于碳氫化合物在空氣中燃燒火焰,可取/f=5和F=

,此時根據(jù)變密度理論求得的火焰長度約為常密度解的2.4倍。略2022/11/10102變密度層流射流火焰動量積分近似Am2022/11/19103

不管哪種計算(常密度、變密度)的結果更接近實際數(shù)值,這兩種理論都指出了火焰長度和噴嘴的體積流量成正比,和噴射流體的燃料當量質量分數(shù)成反比,而與噴嘴直徑的大小無關。火焰長度的特征:FromUCBerkeley的學生實驗:J.Y.Chen2022/11/10103不管哪種計算(常密度、變密度2022/11/191044.4Roper的簡化解Roper[6]沿用了Burke-Schumann分析的主要思想,并加以擴展,考慮了浮力的作用,認為特征速度會隨著軸向距離改變,這和連續(xù)性也是一致的。除了圓口以外,Roper還分析了長方形口和弧形口[6,8]燃燒器。他給出的分析解和經(jīng)過實驗修正后的解將在下面單獨列出。略2022/11/101044.4Roper的簡化解Rope2022/11/191054.5數(shù)值解法(僅了解)火焰面假設不再必要熱輻射擴散比不考慮輻射時低150K略2022/11/101054.5數(shù)值解法(僅了解)火焰面2022/11/191065.圓口和扁口燃燒器的火焰長度Roper關聯(lián)式對于各種幾何形狀(圓形、方形、槽形和弧形)的燃燒器,Roper都做了研究[6,8],得出了對應于不同控制情況(動量控制、浮力控制和過渡區(qū))下的層流射流火焰長度,并且通過實驗進行了校核[7,8]。表9.3中簡單列舉了Roper得到的結果,在下面還將對這些結果進行較為詳細的討論。2022/11/101065.圓口和扁口燃燒器的火焰長度R2022/11/19107噴口幾何尺寸條件可用的計算方程圓口動量或浮力控制圓口方程:式(9.59)和式(9.60)方口動量或浮力控制方口方程:式(9.61)和式(9.62)槽型口動量控制浮力控制混合動量-浮力控制式(9.63)和式(9.64)式(9.65)和式(9.66)式(9.70)表9.3垂直射流火焰長度計算的經(jīng)驗與理論公式對于圓口和方口燃燒器,可用用表9.3中的公式計算火焰長度。這些結果適用于氧氣過量的情況(富氧燃燒),不管燃燒中動量和浮力哪個占主導,也不管燃料噴射進入的是靜止的氧化劑空間還是和燃料和氧化劑同軸射流。2022/11/10107噴口幾何尺寸條件可用的計算方程圓口2022/11/191085.1圓口燃燒器其中S是化學當量氧化劑-燃料摩爾比,D是氧化劑在溫度T下的平均擴散系數(shù)(T

是遠離火焰處環(huán)境流體溫度),TF和Tf分別是燃料流溫度和火焰平均溫度。9.599.60經(jīng)驗公式中去掉了D對比2022/11/101085.1圓口燃燒器其中S是化學當量2022/11/191095.2方口燃燒器式中inverf是反誤差函數(shù)。表9.4中列出了誤差函數(shù)erf的部分取值。從誤差函數(shù)表中查取反誤差函數(shù)和查反三角函數(shù)的方法一樣,即有ω=inverf(erfω)。和上面一樣,式中的所有量均采用國際單位制單位。9.619.62經(jīng)驗公式中去掉了D2022/11/101095.2方口燃燒器9.619.622022/11/191105.3槽形口-動量控制式中的b為槽的寬度,h為槽的長度,β由下面的函數(shù)給出bh9.639.64經(jīng)驗公式中去掉了D方程9.63和9.64僅適用于靜止氧化劑的情形,不適用于同軸射流。2022/11/101105.3槽形口-動量控制bh9.62022/11/19111I為實際流動時槽流出的初始動量流率與均勻流動時值的比值,即:如果流動是均勻的,則有I=1;若h>>b,則在流動充分發(fā)展段內有I=1.5。(9.63)和(9.64)僅適用于氧化劑靜止的情況。對于流動氧化劑的情況,讀者可以參考文獻[6]和[7]。2022/11/10111I為實際流動時槽流出的初始動量流率2022/11/19112槽形口-浮力控制式中a為平均浮力加速度,由下面的式子來計算g為重力加速度;由9.65,9.66,火焰長度隨著a的變化不大,僅和a的-1/3次方成正比,且與槽型口尺寸b無關。9.659.669.67經(jīng)驗公式中去掉了Dbh2022/11/10112槽形口-浮力控制g為重力加速度;92022/11/19113槽形口-過渡區(qū)控制為了判斷火焰是受動量控制還是受浮力控制,需要計算火焰的弗勞德(Froude)數(shù)。從物理意義上來說,F(xiàn)roude數(shù)為射流初始動量流率和火焰處的浮力作用之比。對于噴入靜止介質中的層流噴射火焰,有Frf,而流動受何種控制則可以由下面的標準來判斷:注意:在判斷槽形口流動控制的時候,需要用到Lf,因此需要對Lf進行預估,然后計算弗勞德數(shù),判斷適用何種控制并進行計算,最后對計算結果進行檢驗,看選取的控制是否正確。2022/11/10113槽形口-過渡區(qū)控制為了判斷火焰是受2022/11/19114對過渡區(qū),動量和浮力都起著比較重要的作用,Roper給出了下面的處理方法:其中角標M,B和T分別代表動量控制、浮力控制和過渡區(qū)(先計算出動量控制和浮力控制的火焰長度,再根據(jù)9.70式計算過渡控制的火焰長度)。9.702022/11/10114對過渡區(qū),動量和浮力都起著比較重要2022/11/19115例題9.3實驗室一個方口燃燒器,預得到一個500毫米長的丙烷擴散火焰,求所需的體積流率,同時試確定火焰的放熱量。如果是甲烷火焰又如何?略2022/11/10115例題9.3實驗室一個方口燃燒器,2022/11/19116解運用公式如果假設:T=TF=300K,那么在求QF

之前,我們需求空氣-燃料的化學當量摩爾比,S:化學當量燃燒時氧化劑(空氣和燃料的摩爾比)。a=x+y/4因此,獲得1mol氧氣需要4.76mol空氣,獲得amol需要4.76amol空氣空氣中O2和N2的摩爾比為1:3.76,4.76mol空氣中有1mol氧氣,3.76mol氮氣(相同狀態(tài)下混合氣體總摩爾數(shù)等于各組分氣體摩爾數(shù)之和)9.62略2022/11/10116解運用公式a=x+y/2022/11/19117解對于丙烷:S=(3+8/4)×4.76/1=23.8kmol/kmol.則有:inverf[(1+23.8)-0.5]=inverf(0.2008)=0.18解方程9.62求出QF(體積流量)9.62略2022/11/10117解對于丙烷:S=(3+8/4)×42022/11/19118根據(jù)理想氣體狀態(tài)方程求丙烷的密度(求出質量流量)(P=1atm,T=300K),并根據(jù)附錄B查出燃燒熱值:則反應放熱為:略2022/11/10118根據(jù)理想氣體狀態(tài)方程求丙烷的密度(2022/11/19119對于甲烷,可以求出S=9.52,F(xiàn)=0.65,hC=50,016,000J/kg;則:和評述:從結果中可以看出,雖然甲烷的體積流率是丙烷的2.4倍,但二者的反應放熱卻很接近。略2022/11/10119對于甲烷,可以求出S=9.52,2022/11/19120作業(yè)02022/11/10120作業(yè)9.12022/11/19121影響火焰高度的因素a)流量和幾何形狀的影響圖9.9中是圓口燃燒器和不同長寬比的槽形口燃燒器的火焰長度的對比.(各種情況下的噴口面積都相等,因此平均出口速度也都相等)從圖中可以看到:1.圓口燃燒器的火焰長度和燃料的體積流率成線性關系,而槽形口燃燒器的火焰長度隨燃料體積流量的變化率呈比線性更強的上升趨勢。圖中火焰的弗勞德數(shù)都很?。髁啃。?,即火焰是受浮力控制的2.當h/b變大時,槽型口火焰會明顯地變短。bh2022/11/10121影響火焰高度的因素a)流量和幾何2022/11/19122b)影響化學當量的因素前面得出式子里面,用到了化學當量摩爾比S這個概念,它是用噴射流體和環(huán)境流體來定義的:可以看出,噴射流體和環(huán)境流體的化學成分都會影響到S的值。例如,對于純燃料和用氮氣稀釋后的燃料在空氣中燃燒這兩種情況,它們的S的取值就不同。類似的,氧氣在環(huán)境流體中的摩爾分數(shù)也會影響到S。在大多數(shù)的應用中,我們關心的主要是下面的幾個參數(shù)對S的影響。9.712022/11/10122b)影響化學當量的因素前面得出式2022/11/19123b-1)燃料類型對于碳氫化合物CxHy

,這個比值根據(jù)下面的式子來計算,其中xO2是空氣中的氧氣摩爾分數(shù)(一般取1/4.76)。燃料分子越大,S值越大9.724.76摩爾空氣中有1摩爾氧氣,因此:2022/11/10123b-1)燃料類型對于碳氫化合物C2022/11/19124在圖9.10給出了由圓口表達式(9.60)計算得出的氫氣、一氧化碳以及含1-4個碳的烷烴的火焰相對長度(其中的每種情況均具有相同的燃料流率,并都以甲烷的火焰長度為標準)。在這個式子中,認為各種混合物都具有相同的平均擴散系數(shù),這只是一個近似的假設;而對于氫氣來說,這個假設根本就不合理。2022/11/10124在圖9.10給出了由圓口表達式(92022/11/19125當燃料的氫碳比(摩爾比和質量比)減少時(碳氫燃料分子變大),火焰長度增加,如丙烷火焰是甲烷火焰長度的2.5倍。對于同一類的(比如烷烴類或烯烴類)高碳氫化合物,當碳原子數(shù)增加時,對火焰長度的影響相對低碳氫化合物碳原子變化對火焰產(chǎn)生的影響?。ㄒ驗闅涮急茸兓。?。一氧化碳和氫氣的火焰和碳氫化合物的相比要短得多。2022/11/10125當燃料的氫碳比(摩爾比和質量比)減2022/11/19126b-2)一次風對一個應用層流射流擴散火焰的燃氣設備,通常在氣體燃料在燃燒以前都要和空氣進行部分預混,這部分預混的空氣就是一次風。一次風量一般為完全燃燒所需空氣量的40-60℅,它使得燃燒的火焰變短,防止了碳黑的生成,通常會產(chǎn)生藍色的火焰。引入一次風量的最大值受到了安全性的限制,如果加入的量過大,就可能超過可燃上限(richflammabilty),即此時混合物具備形成預混火焰條件(一次風氧氣量足夠多,且與燃料混合),此時,燃燒中可能有回火現(xiàn)象(流速較低時);當氣流速度足夠大就可能產(chǎn)生類似Bunsen燈內焰的預混火焰。2022/11/10126b-2)一次風對一個應用層流射流2022/11/19127一次風量會影響到火焰長度,圖9.11中給出了一次風對圓口燃燒器中甲烷火焰長度的影響情況。一次風量為40-60℅時,Lf和不加一次風相比,減小了85-90℅。在加入一次風情況下,可以將噴射流體當作是純燃料和空氣的混合物,來計算(9.71)定義的化學當量摩爾比S:其中ψpri是一次風量占所需空氣的百分比,Spure是純燃料對應的化學當量摩爾比。推導一下9.719.732022/11/10127一次風量會影響到火焰長度,圖9.12022/11/19128b-3)氧化劑的含氧量氧化劑中的含氧量對火焰長度的影響很大,由圖9.12中,空氣中的含氧量為21℅,如果氧氣含量在此基礎上的減少一點,所產(chǎn)生的火焰長度就會大大增加。以甲烷在純氧中和空氣中燃燒為例,前者的火焰長度只有后者的四分之一左右。通過(9.72)來計算碳氫化合物在不同含氧量下化學當量摩爾比,觀察含氧量的影響。2022/11/10128b-3)氧化劑的含氧量氧化劑中的2022/11/19129b-4)惰性氣體的稀釋用惰性氣體來對燃料進行稀釋,也會影響到化學當量比,從而影響到火焰長度。對碳氫化合物來說,其中χdil是燃料流中稀釋劑(惰性氣體)的摩爾分數(shù)。9.74稀釋2022/11/10129b-4)惰性氣體的稀釋用惰性氣體2022/11/19130例題9.4設計一個餐飲業(yè)用的天然氣的燃燒器,多個圓口燒嘴呈圓形分布。圓管直徑為160mm(6.3in.)。在滿負荷時,燃燒器功率必須達到2.2KW,且混入40%的一次風。為了穩(wěn)定燃燒,每一個圓口燒嘴的功率不得超過10W/mm2。(設計時參見圖8.23)同時,滿負荷時火焰高度不得超過20毫米。試求需要布置多少個圓形燒嘴,每一個直徑是多少?略2022/11/10130例題9.4設計一個餐飲業(yè)用的天然2022/11/19131解假設燃料為甲烷,當然,實際設計中應該用天然氣的性質。接下來要做的就是將燒嘴數(shù)量N和它們的直徑D與題目給出的條件建立聯(lián)系;選擇合適的N和D;然后檢驗火焰長度是否超標,最后,我們要看所做設計是否符合實際的物理意義。Step1根據(jù)符合的限制,總的燒嘴橫截面積為:略2022/11/10131解假設燃料為甲烷,當然,實際設計中2022/11/19132Step1其中,要求即:現(xiàn)在,我們任意設定一個N(D)來計算D(N)。如設:N=36求得D=2.79mm.略2022/11/10132Step1其中,要求現(xiàn)在,我們2022/11/19133Step2.求體積流率。根據(jù)設計的放熱率求燃料的體積流率:用一次風確定預混入燃料的空氣流量:略2022/11/10133Step2.求體積流率。根據(jù)設2022/11/19134總的體積流率為:求解,我們運用理想氣體方程,其中,平均摩爾質量根據(jù)空氣-燃料混和物的比例確定:其中Z是一次空-燃摩爾比:上下除以NF略2022/11/10134總的體積流率為:上下除以NF略2022/11/19135狀態(tài)方程:略2022/11/10135狀態(tài)方程:略2022/11/19136Step3.檢驗火焰長度。每個燒嘴的流率為:環(huán)境空氣和燒嘴流體的當量摩爾比S由公式9.73給出:略2022/11/10136Step3.檢驗火焰長度。每個2022/11/19137利用公式9.60我們可以求出火焰長度:19.6mm的火焰長度符合要求Lf<20mm。略2022/11/10137利用公式9.60我們可以求出火焰長2022/11/19138Step4檢驗設計的可行性。如果在直徑160毫米的圓周上均勻布置36個燒嘴,則燒嘴間距為:略2022/11/10138Step4檢驗設計的可行性。如果1396.碳黑(Soot)的形成與破壞本章開頭對層流射流火焰的概述中提到過碳黑的形成和氧化,這是碳氫化合物非預混火焰的一個重要特點?;鹧胬锩娴奶己谑軣岷髸l(fā)光,這是擴散火焰發(fā)光體的主要來源,歷史上最早應用在油燈燃燒。碳黑發(fā)出的波長大多位于光譜的紅外區(qū)域,因此還能導致火焰的輻射熱損失。在煤氣爐、柴油機、燃煤爐等實際應用中,應盡量避免碳黑的形成,但是在對燃燒中碳黑形成的基礎研究中,因為層流擴散火焰和碳黑的形成有很大的關系,因此層流非預混火焰通常是一種研究soot工具。1396.碳黑(Soot)的形成與破壞本章開頭對層流射流火2022/11/19140在碳氫化合物的燃燒火焰中,由于經(jīng)常會有碳黑存在,火焰就可能呈現(xiàn)為橙色或黃色。Q-2:乙烯、甲烷、航空汽油柴油機煤燃燒層流非預混火焰:碳煙Soot的產(chǎn)生2022/11/10140在碳氫化合物的燃燒火焰中,由于經(jīng)常2022/11/19141研究普遍認為擴散火焰中的碳黑是在一定的溫度范圍形成的,這個范圍大致為1300K<T<1600K;圖9.13以乙烯為例描述了這種現(xiàn)象。圖中給出了燃燒器和焰舌之間的兩條變化曲線(離噴口軸向距離5mm處),一條是溫度沿徑向的變化曲線,另一條是被碳黑顆粒散射的光線強度曲線,而后一條曲線的峰值所處區(qū)域里面,存在著大量的碳黑。碳黑含量峰值對應的溫度大約為1600K,并且位于溫度峰值的徑向的內側。含有碳黑的區(qū)域很窄,并且只存在于一定的溫度范圍之內。圖9.13層流乙烯射流火焰的溫度和散射光的徑向分布,碳煙存在于散光強度高的區(qū)域柴油機低溫燃燒2022/11/10141研究普遍認為擴散火焰中的碳黑是在一2022/11/19142Santoro教授最近的研究2022/11/10142Santoro教授最近的研究2022/11/19143

4、顆粒被氧化3、顆粒的長大和團聚2、開始生成初始顆粒1、顆粒前驅體(PAHs)形成碳煙形成過程分析2022/11/10143碳煙形成過程分析144第一個步驟是前驅物的形成。在這一步中多環(huán)芳烴(PAH)是燃料分子向初始碳黑顆粒轉變過程中的一種重要的中間產(chǎn)物[20],其中化學動力學起著很重要的作用。確切的PAH種類還有待研究。(1)形成碳煙顆粒的前驅物1-形成多環(huán)芳烴2-多環(huán)芳烴長大H.WangandM.Frenklach,1997144第一個步驟是前驅物的形成。在這一步中多環(huán)芳烴(PAH)2022/11/19145在顆粒的生成中,PAHs通過化學和凝結成核作用,形成了臨界尺寸(3000-10000原子質量單位)的小顆粒。通過這一步,大分子轉變成為了顆粒。(2)PAHs轉變成小顆粒2022/11/10145在顆粒的生成中,PAHs通過化學和2022/11/19146

接下來的第三步中,小的初始碳黑顆粒被置于燃料流中,在隨著燃料流向火焰運動的過程中,不斷長大和聚合,并且通過火焰的某個氧化區(qū)(如射流火焰的焰舌)。(3)顆粒聚合Coagulation:collision-coalescenceTime2022/11/10146接下來的第三步中,小的初2022/11/19147(4)碳煙氧化Warnatz,Maas&Dibble,combustion,4thedition包含了Soot的4個過程。

JürgenWarnatz在氧化區(qū)(如射流火焰的焰舌),碳煙被氧化,如果完全被氧化,就不產(chǎn)生碳煙。2022/11/10147(4)碳煙氧化Warnatz,2022/11/19148在圖9.14中,對丙烯和丁烯的層流非預混火焰,分別對有煙和無煙的不同情況做了分析。在焰舌外面(x/xstoic≥1.1),如果碳黑的體積流量不為零,就意味著SootingFlame,即含煙火焰。從圖9.14中可以看出,擴散火焰中是否會有煙的生成和燃料類型有很大的關系。圖9.14在成煙和不成煙條件下,燃燒丁烯和丙烯時測量到的碳煙體積分數(shù)隨燃燒器高度的變化1.12022/11/10148在圖9.14中,對丙烯和丁烯的層流2022/11/19149如果將表9.5中的燃料進行分類,就可以看到不同種類燃料的發(fā)煙趨勢從小到大依次為烷烴、烯烴、炔類和芳香烴。(4)發(fā)煙點Smokepoint定義發(fā)煙點:逐漸增大擴散火焰中燃料的流量,直到焰舌處開始出現(xiàn)碳煙,此時燃料的流量稱為發(fā)煙點。發(fā)煙點越大,表示這種燃料越不容易形成碳煙。2022/11/10149如果將表9.5中的燃料進行分類,就2022/11/19150表9.6中給出了這些燃料的分類,顯然燃料的分子結構對其發(fā)煙趨勢起著很重要的作用。根據(jù)煙形成的化學機理,環(huán)狀化合物及其通過和乙炔反應而長大是其中的重要特征,在這點上面和這里得出的分類也是一致的。(4)發(fā)煙點Smokepoint定義2022/11/10150表9.6中給出了這些燃料的分類,顯2022/11/191516-1碳煙顆粒的形貌分析——TEM燃燒器出口沿程變化Zhao,Uchikawa&HaiWang;Proc.Combust.Inst.31(2007)851–860HaiWang2022/11/101516-1碳煙顆粒的形貌分析——TEM2022/11/197.對吹火焰過去幾十年里,人們對圖9.15中所示的對吹火焰,即燃料和氧化劑對射產(chǎn)生的火焰,做了大量的理論和實驗研究。這種火焰可以近似看作是一維的,并且其火焰區(qū)的停留時間較容易調節(jié),因此被當作是一種基礎理論研究而被人們所重視。前面討論二維(軸對稱)射流火焰是很復雜的,而一維對吹火焰在試驗和計算等方面都容易實現(xiàn)。例如在試驗中,只需要測量一條線上面的溫度和物質質量分數(shù);理論研究中,即使是使用復雜的化學動力學(見表5.3)來計算,也不需要太多的時間。NonpremixedPremixed2022/11/107.對吹火焰過去幾十年里,人們對圖9.2022/11/19153

在對對吹火焰進行數(shù)學描述以前,有必要知道它的一些基本特征。典型的對吹火焰實驗如圖9.15所示,燃料和氧化劑相對噴射,在兩個噴嘴之間形成了一個滯止面(vx=0),其位置由燃料和氧化劑的初始噴射動量通量的相對大小來決定。如果二者的動量通量相等,即那么這個滯止面就位于兩個噴嘴的中點處,否則滯止面就會向低動量通量的流體側移動。由于火焰位于化學當量比下混合物分數(shù)fstoic的地方,就可以根據(jù)這些條件來確定火焰在噴嘴之間的位置。2022/11/10153在對對吹火焰進行數(shù)學描述以2022/11/19154對于大多數(shù)碳氫化合物在空氣中的燃燒來說,其化學當量比fstoic≈0.06(燃料摩爾數(shù)與燃料+空氣摩爾數(shù)的比值),因此需要空氣量比燃料量多很多才能滿足這個條件。在這種情況下的靜止面處,空氣的量達不到當量燃燒所需的空氣的量,而只有在靜止面的空氣側才能滿足這個條件,此處燃料已經(jīng)擴散并穿越靜止面,并發(fā)生燃燒(這樣導致火焰面在靜止面的空氣側),如圖9.15所示。ν為化學當量系數(shù)相反的,如果按化學當量比混合時所需的燃料量大于氧化劑量,即fstoic>0.5,那么火焰面就位于靜止面的燃料側。2022/11/10154對于大多數(shù)碳氫化合物在空氣中的燃燒2022/11/191559-129-16課后作業(yè)預習第10章2022/11/101559-12課后作業(yè)預習第10章2022/11/191567.1數(shù)學描述對對吹火焰的建模由兩種不同的處理方法。第一種將無限遠處點源產(chǎn)生的駐點勢流(stagnation-point)和邊界層分析結合起來[31],在這種分析里面,兩個噴嘴之間的有限距離無法得以考慮。第二種方法[32,33]明確指出了流體是從噴嘴流出的,而不是由無限遠處的點源產(chǎn)生的。本節(jié)重點介紹這種方法。這個模型最初是針對預混火焰建立的[32],后來才擴展到了非預混火焰[33]。下面將對第二種方法做簡單的概述,并使用這個模型來檢驗數(shù)值解給出的甲烷-空氣擴散火焰的具體結構。對該方法的詳細說明可以參見參考文獻32和33。略過2022/11/101567.1數(shù)學描述對對吹火焰的建模由2022/11/19157SAND96-8243UnlimitedReleasePrintedMay1997SAND85-8240,UC-4UnlimitedReleasePrintedApril1995對吹火焰的算法目前國際上最公認的是Sandia的PREMIX和OPPDIX.略過2022/11/10157SAND96-8243SAND852022/11/19158這種分析方法總的目標是將用于軸對稱的偏微分控制方程組轉換為常微分方程組,并把它看作邊界值問題來求解。首先將對用連續(xù)性和動量守恒方程的軸對稱形式進行處理,這些方程均在第七章中給出,其中連續(xù)性方程為(7.7),軸向和徑向動量守恒方程分別為(7.43)和(7.44)。為了實現(xiàn)方程組的轉換,引入下面的流函數(shù):對吹火焰的算法目前國際上最公認的是Sandia的PREMIX和OPPDIF.略過2022/11/10158這種分析方法總的目標是將用于軸對稱2022/11/19159VonKarman:可壓縮流體-了解應用流體力學的本質連續(xù)性方程和r無關(Eq.9.76)(Eq.9.77)略過2022/11/10159VonKarman:可壓縮流體-2022/11/19160將(9.76a)、(9.76b)和(9.77)代入動量守恒方程(7.43)和(7.44)中,再忽略浮力的作用,就可以得到下面的式子:略過2022/11/10160將(9.76a)、(9.76b)和2022/11/19161通過這兩個式子可以獲得關于徑向壓力梯度的特征方程,其左邊的表達式之間的關系為:而該問題為一維的,P/x和(1/r)(P/r)都只是x的函數(shù),因此可得略過2022/11/10161通過這兩個式子可以獲得關于徑向壓力2022/11/19162其中H為徑向壓力梯度的特征值,將在常微分方程組里面以下面的形式出現(xiàn):當流體的馬赫數(shù)較低時,可以近似認為流場里各處

溫馨提示

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

評論

0/150

提交評論