非灰介質輻射特性計算方法的研究

非灰介質輻射特性計算方法的研究

非灰介質輻射特性在許多碳氫化合物的燃燒過程中，輻射熱流的計算對環(huán)境和燃燒性能產生了重大影響。同時，氣體輻射也是一些非接觸光學氣流診斷技術的基礎，在工程應用中正確計算氣體輻射具有重要意義。鍋爐燃燒室中不僅有燃燒氣體混合物(主要是二氧化碳和水蒸氣),還含有煤灰、炭黑粒子,這些粒子的吸收和散射顯著增強了爐內的輻射換熱。為了準確計算輻射,必須考慮混合物的非灰輻射特性?？紤]氣體的非灰輻射特性已使整個輻射計算量非常大,加上粒子使輻射問題變得更復雜。工程應用中的氣體輻射計算主要涉及三個不同層次的問題:(1)有效處理大量吸收氣體光譜線的輻射特性計算模型;(2)準確地吸收氣體光譜線數據;(3)準確有效求解多維空間輻射傳遞方程的方法。但是目前還不能對工程中的輻射熱傳遞進行精確和高效的計算,原因在于:(1)缺乏精確的氣體輻射光譜數據,由于高溫實驗難以實現,現有數據庫高溫區(qū)的光譜線數據大多是通過外推得到的,其正確性有待進一步研究;(2)缺乏有效、靈活的精確預測非灰介質輻射傳遞的計算模型。由于高溫實驗需要耗費大量的財力,而且容易受實驗設備的影響,使實驗的實施有一定的困難,所以目前的研究主要集中在理論模型方面。1譜帶模型的設置根據計算波數間距大小,可以把氣體輻射特性計算方法分成三類:(1)逐線計算,波數間隔一般為0.0002～0.02cm-1;(2)譜帶模型,窄譜帶模型的波數間隔一般為5～50cm-1,寬譜帶模型的波數間隔一般為100～1000cm-1;(3)總體模型,波數范圍包含整個光譜區(qū)間。計算氣體輻射特性最精確的方法是逐線計算(LBL),由于其巨大的計算量,一般只將逐線計算結果作為檢驗其他模型精確度和有效性的基準解。2采用k分步法和傳統(tǒng)帶模式法k分布法是近些年發(fā)展起來的一種先進的計算大氣分子吸收的處理方法,它把隨波長迅速變化的大氣吸收離散成有限項的贗單色吸收,通過對有限項的積分可以得到與逐線積分(LBL)非常接近的精度,大大提高了計算效率。同時,因為使用贗單色計算,k分布法或相關k分布法可以直接用于散射的處理,而傳統(tǒng)的帶模式法不能直接用于散射處理,必須作一定的近似。k分步法相對于傳統(tǒng)的帶模式法在計算精度和速度方面都有很大的優(yōu)勢,在氣體輻射學的研究中得到了廣泛的應用。在小的譜帶間隔內,普朗克函數基本不變,而光譜吸收系數劇烈變化,在波數間隔很小的情況下吸收系數數值多次相同,求解吸收系數相同的輻射傳遞方程所得到的輻射強度相同。吸收系數又稱氣體輻射減弱系數或衰減系數,它是由氣體輻射喝吸收再整個容積中進行所引起的。當氣體或容器壁發(fā)射出輻射能時,可以射入到氣體容積內的任何地方,但輻射能在射線行程過程中被有吸收能力的氣體分子所部分吸收而逐漸削弱。k分布思想就是將光譜吸收系數按數值大小重新排列成平滑單調上升的曲線,對數值相同的吸收系數就只用求解一次輻射傳遞方程,用很少的計算量就可以求解氣體輻射傳遞。3k分布法廣泛應用于氣體輻射特性的計算中3.1基于分離子的高溫輻射特性模型在非等溫介質中,如果不同溫度下累積k分布函數在相同位置的吸收系數是一一對應的,則說k分布是相互“關聯(lián)”的。對溫度介于200K到300K之間,壓力變化明顯的大氣科學研究,窄譜帶關聯(lián)k分布模型應用很成功。石廣玉等人論述了大氣輻射計算中的窄譜帶k分布模型,包括均勻路徑和非均勻路徑的關聯(lián)k分布模型。Taine等人利用EM2C數據庫得到適用于統(tǒng)計窄譜帶模型和窄譜帶關聯(lián)k模型的高溫氣體輻射特性參數。Dembele等人用窄譜帶關聯(lián)k分布方法研究了變濃度、非等溫燃燒產物中的輻射熱傳遞。Coelho等人在三維非灰氣體輻射熱傳遞的數值模擬中應用了窄譜帶關聯(lián)k分布模型。Goutiere等人在對二維氣體輻射進行研究時提出了統(tǒng)計窄譜帶關聯(lián)k分布模型(SNBCK),它是將統(tǒng)計窄譜帶模型和關聯(lián)k分布思想相結合的一種方法,其通過對統(tǒng)計窄譜帶明中的平均透射率反Laplace變換求得吸收系數的k分布函數。Liu等人考慮了積分格式的選取對SNBCK計算精度的影響,發(fā)現:7點Gauss-Lobatto積分結果和SNB吻合很好;4點Gauss-Lobatto積分和7點Gauss-Lobatto積分精度差不多,但明顯節(jié)約計算時間;對多維問題,綜合考慮計算精度和效率,推薦選用2點Gauss-Lobatto積分格式。3.2混合氣體重疊譜帶將k分布思想應用到整個振動—轉動譜帶,就得到寬譜帶k分布模型。Marin和Buckius等人提出了一種較有效準確求解寬帶k分布吸收系數的方法,模擬了CO2和H2O重要輻射譜帶的輻射。他們將累計k分布函數g(k)分四個區(qū)間,用關聯(lián)式來確定這些譜帶的k分布吸收系數,并討論了H2O和CO2混合氣體重疊譜帶吸收系數的計算。石廣玉等人利用AFGL(AirForceGeophysicalLaboratory)數據庫中吸收氣體的光譜特性參數,在Malkmus模型和Edwards寬譜帶假設基礎上得到寬帶k分布吸收系數。聶宇宏等人用Marin寬譜帶k分布模型和重新排列吸收系數模型,計算了H2O和CO2的有效帶寬,并將修正寬譜帶模型計算帶寬的方法引入到累積k分布函數的計算中。Cayan等人在計算氣體輻射特性時用了寬譜帶k分布方法。3.3多尺度全光譜標定及最優(yōu)分區(qū)的確定Modest和Zhang提出了一種新的全光譜k分布模型(FSK),為計算氣體輻射特性提供了一種較為準確有效的方法。其主要思想是引進普朗克函數加權吸收系數k的幾率分布函數,求解輻射傳遞時只需對少量的吸收系數進行數值積分,大大提高了輻射計算效率。在吸收系數k分布“關聯(lián)”假設成立的情況下,該方法較準確。全光譜k分布模型進一步發(fā)展到多尺度全光譜標定k分布模型——Multi-ScaleFull-SpectrumScaled-kDistributionModel(MSFSSK)和多組全光譜關聯(lián)k分布模型——MultiGroupFull-SpectrumCorrelated-kModel(MGFSCK)。MSFSSK模型是根據氣體種類的不同或一種氣體按分子躍遷的低能級能量E″(T)值分組,和假想氣體思想類似。對溫度場合摩爾分數的變化的一維或二維非均勻氣體混合物,只需幾組,該模型就能得到較準確的計算結果。但MSFSSK模型中關于如何選擇最優(yōu)參考狀態(tài)以及確定分區(qū)的E″(T)的數值,有待進一步研究。MGFSCK模型是根據氣體光譜吸收系數與溫度和局部壓力的關系將光譜位置分成M組,對M組分別計算全光譜關聯(lián)k分布吸收系數,當M≤32就可以獲得和LBL就按相當的計算精度。但對非均勻氣體,若考慮大的溫度范圍,需要建立大量的光譜組數據,使得計算量很大。MGFSCK模型中總體壓力變化對氣體光譜吸收系數的影響,以及該模型用于混合物、處理非灰邊界和散射問題都需要進一步研究。4特性計算方法的運用本文主要介