用k分布法計算氣體輻射特性的進(jìn)展

尹雪梅 閆立強 吳學(xué)紅 鄭州輕工業(yè)學(xué)院機(jī)電工程學(xué)院，河南省 鄭州市 450002

用k分布法計算氣體輻射特性的進(jìn)展

尹雪梅 閆立強 吳學(xué)紅 鄭州輕工業(yè)學(xué)院機(jī)電工程學(xué)院，河南省 鄭州市 450002

介紹了氣體輻射特性計算的分類，以及k分布法在各種氣體輻射特性計算中的應(yīng)用，重點闡述了國內(nèi)外用k分布法計算氣體輻射特性的研究現(xiàn)狀。

氣體輻射；k分布法

引言

在許多碳?xì)浠衔锏娜紵龍龊?，輻射熱流的計算結(jié)果對環(huán)境和燃燒性能影響很大，同時氣體輻射也是一些非接觸式光學(xué)氣流診斷技術(shù)的基礎(chǔ)，在工程應(yīng)用中對氣體輻射進(jìn)行準(zhǔn)確計算具有重要意義。鍋爐燃燒室中不僅有燃燒氣體混合物（主要是二氧化碳和水蒸氣），還含有煤灰、炭黑粒子，這些粒子的吸收和散射顯著增強了爐內(nèi)的輻射換熱。為了準(zhǔn)確計算輻射，必須考慮混合物的非灰輻射特性。考慮氣體的非灰輻射特性已使整個輻射計算量非常大，加上粒子使輻射問題變得更復(fù)雜。

工程應(yīng)用中的氣體輻射計算主要涉及三個不同層次的問題：(1)有效處理大量吸收氣體光譜線的輻射特性計算模型；(2)準(zhǔn)確地吸收氣體光譜線數(shù)據(jù)；(3)準(zhǔn)確有效求解多維空間輻射傳遞方程的方法。但是目前還不能對工程中的輻射熱傳遞進(jìn)行精確和高效的計算，原因在于：(1)缺乏精確的氣體輻射光譜數(shù)據(jù)，由于高溫實驗難以實現(xiàn)，現(xiàn)有數(shù)據(jù)庫高溫區(qū)的光譜線數(shù)據(jù)大多是通過外推得到的，其正確性有待進(jìn)一步研究；(2)缺乏有效、靈活的精確預(yù)測非灰介質(zhì)輻射傳遞的計算模型。由于高溫實驗需要耗費大量的財力，而且容易受實驗設(shè)備的影響，使實驗的實施有一定的困難，所以目前的研究主要集中在理論模型方面。

1 氣體輻射特性計算的分類

根據(jù)計算波數(shù)間距大小，可以把氣體輻射特性計算方法分成三類：(1)逐線計算，波數(shù)間隔一般為0.0002～0.02cm-1；(2)譜帶模型，窄譜帶模型的波數(shù)間隔一般為5～50cm-1，寬譜帶模型的波數(shù)間隔一般為100～1000cm-1；(3)總體模型，波數(shù)范圍包含整個光譜區(qū)間。計算氣體輻射特性最精確的方法是逐線計算(LBL)，由于其巨大的計算量，一般只將逐線計算結(jié)果作為檢驗其他模型精確度和有效性的基準(zhǔn)解。

2 k分布法簡介

k分布法是近些年發(fā)展起來的一種先進(jìn)的計算大氣分子吸收的處理方法，它把隨波長迅速變化的大氣吸收離散成有限項的贗單色吸收，通過對有限項的積分可以得到與逐線積分(LBL)非常接近的精度，大大提高了計算效率。同時，因為使用贗單色計算，k分布法或相關(guān)k分布法可以直接用于散射的處理，而傳統(tǒng)的帶模式法不能直接用于散射處理，必須作一定的近似。k分步法相對于傳統(tǒng)的帶模式法在計算精度和速度方面都有很大的優(yōu)勢，在氣體輻射學(xué)的研究中得到了廣泛的應(yīng)用。

在小的譜帶間隔內(nèi)，普朗克函數(shù)基本不變，而光譜吸收系數(shù)劇烈變化，在波數(shù)間隔很小的情況下吸收系數(shù)數(shù)值多次相同，求解吸收系數(shù)相同的輻射傳遞方程所得到的輻射強度相同。吸收系數(shù)又稱氣體輻射減弱系數(shù)或衰減系數(shù)，它是由氣體輻射喝吸收再整個容積中進(jìn)行所引起的。當(dāng)氣體或容器壁發(fā)射出輻射能時，可以射入到氣體容積內(nèi)的任何地方，但輻射能在射線行程過程中被有吸收能力的氣體分子所部分吸收而逐漸削弱。

k分布思想就是將光譜吸收系數(shù)按數(shù)值大小重新排列成平滑單調(diào)上升的曲線，對數(shù)值相同的吸收系數(shù)就只用求解一次輻射傳遞方程，用很少的計算量就可以求解氣體輻射傳遞。

3 k分布法在氣體輻射特性計算中的應(yīng)用

3.1 窄譜帶k分布模型

在非等溫介質(zhì)中，如果不同溫度下累積k分布函數(shù)在相同位置的吸收系數(shù)是一一對應(yīng)的，則說k分布是相互“關(guān)聯(lián)”的。對溫度介于200K到300K之間，壓力變化明顯的大氣科學(xué)研究，窄譜帶關(guān)聯(lián)k分布模型應(yīng)用很成功。石廣玉[1]等人論述了大氣輻射計算中的窄譜帶k分布模型，包括均勻路徑和非均勻路徑的關(guān)聯(lián)k分布模型。Taine等人利用EM2C數(shù)據(jù)庫得到適用于統(tǒng)計窄譜帶模型和窄譜帶關(guān)聯(lián)k模型的高溫氣體輻射特性參數(shù)。Dembele等人用窄譜帶關(guān)聯(lián)k分布方法研究了變濃度、非等溫燃燒產(chǎn)物中的輻射熱傳遞。Coelho等人在三維非灰氣體輻射熱傳遞的數(shù)值模擬中應(yīng)用了窄譜帶關(guān)聯(lián)k分布模型。

Goutiere等人[2]在對二維氣體輻射進(jìn)行研究時提出了統(tǒng)計窄譜帶關(guān)聯(lián)k分布模型（SNBCK），它是將統(tǒng)計窄譜帶模型和關(guān)聯(lián)k分布思想相結(jié)合的一種方法，其通過對統(tǒng)計窄譜帶明中的平均透射率反Laplace變換求得吸收系數(shù)的k分布函數(shù)。Liu等人[3]考慮了積分格式的選取對SNBCK計算精度的影響，發(fā)現(xiàn)：7點Gauss-Lobatto積分結(jié)果和SNB吻合很好；4點Gauss-Lobatto積分和7點Gauss-Lobatto積分精度差不多，但明顯節(jié)約計算時間；對多維問題，綜合考慮計算精度和效率，推薦選用2點Gauss-Lobatto積分格式。

3.2 寬譜帶k分布模型

將k分布思想應(yīng)用到整個振動—轉(zhuǎn)動譜帶，就得到寬譜帶k分布模型。Marin和Buckius等人[4]提出了一種較有效準(zhǔn)確求解寬帶k分布吸收系數(shù)的方法，模擬了CO2和H2O重要輻射譜帶的輻射。他們將累計k分布函數(shù)g(k)分四個區(qū)間，用關(guān)聯(lián)式來確定這些譜帶的k分布吸收系數(shù)，并討論了H2O和CO2混合氣體重疊譜帶吸收系數(shù)的計算。石廣玉等人利用AFGL(Air Force Geophysical Laboratory)數(shù)據(jù)庫中吸收氣體的光譜特性參數(shù)，在Malkmus模型和Edwards寬譜帶假設(shè)基礎(chǔ)上得到寬帶k分布吸收系數(shù)。聶宇宏[5]等人用Marin寬譜帶k分布模型和重新排列吸收系數(shù)模型，計算了H2O和CO2的有效帶寬，并將修正寬譜帶模型計算帶寬的方法引入到累積k分布函數(shù)的計算中[6]。Cayan等人在計算氣體輻射特性時用了寬譜帶k分布方法。

3.3 全光譜k分布模型

Modest和Zhang[7]提出了一種新的全光譜k分布模型(FSK)，為計算氣體輻射特性提供了一種較為準(zhǔn)確有效的方法。其主要思想是引進(jìn)普朗克函數(shù)加權(quán)吸收系數(shù)k的幾率分布函數(shù)，求解輻射傳遞時只需對少量的吸收系數(shù)進(jìn)行數(shù)值積分，大大提高了輻射計算效率。在吸收系數(shù)k分布“關(guān)聯(lián)”假設(shè)成立的情況下，該方法較準(zhǔn)確。

全光譜k分布模型進(jìn)一步發(fā)展到多尺度全光譜標(biāo)定k分布模型——Multi-Scale Full-Spectrum Scaled-k Distribution Model(MSFSSK)和多組全光譜關(guān)聯(lián)k分布模型——Multi-Group Full-Spectrum Correlated-k Model(MGFSCK)。MSFSSK模型是根據(jù)氣體種類的不同或一種氣體按分子躍遷的低能級能量E″（T）值分組，和假想氣體思想類似。對溫度場合摩爾分?jǐn)?shù)的變化的一維或二維非均勻氣體混合物，只需幾組，該模型就能得到較準(zhǔn)確的計算結(jié)果。但MSFSSK模型中關(guān)于如何選擇最優(yōu)參考狀態(tài)以及確定分區(qū)的E″（T）的數(shù)值，有待進(jìn)一步研究。

MGFSCK模型是根據(jù)氣體光譜吸收系數(shù)與溫度和局部壓力的關(guān)系將光譜位置分成M組，對M組分別計算全光譜關(guān)聯(lián)k分布吸收系數(shù)，當(dāng)M≤32就可以獲得和LBL就按相當(dāng)?shù)挠嬎憔?。但對非均勻氣體，若考慮大的溫度范圍，需要建立大量的光譜組數(shù)據(jù)，使得計算量很大。MGFSCK模型中總體壓力變化對氣體光譜吸收系數(shù)的影響，以及該模型用于混合物、處理非灰邊界和散射問題都需要進(jìn)一步研究。

4 小結(jié)

本文主要介紹了用k分步法計算氣體輻射特性的研究進(jìn)展，并討論了k分布法在各種輻射特性計算中的具體運用。氣體輻射特性計算只占整個輻射計算的一部分，應(yīng)盡可能地提高其與輻射傳遞方程求解方法的兼容性及計算效率，所以建立一種能與任意輻射傳遞方程求解方法相容并且準(zhǔn)確快速的計算模型具有重要作用，k分布法在將來會有更廣泛的應(yīng)用。

[1]HUA Zhang, SHI Guang-yu. A new approach to solve correlated k-distribution function[J]. Journal of Quantitative Spectroscopy and Radiative,2005,96(2):311～324

[2]GOUTIERE V, LIU Feng-shan, CHARETTE A. An assessment of real-gas modelling in2D enclosures[J]. Journal of Quantitative Spectroscopy and Radiative,2000,64(3):299～326

[3]LIU Feng-shan, GUO Hong-sheng, SMALLWOOD J G. Effects of radiation model on the modeling of a laminar coflow methane/air diffusion flame[J]. Combustion and Flame,2004,138(1-2):136～154

[4]JING He, CHENG W L, BUCKIUS O R. Wide band cumulative absorption coefficient distribution model for overlapping absorption in H2O and CO2 mixtures[J]. International Journal of Heat and Mass Transfer,2008,51(5-6):1115～1129

[5]聶宇宏，陳海耿.非灰氣體的寬帶吸收系數(shù)模型[J].東北大學(xué)學(xué)報(自然科學(xué)版)，2002，23（4）：375～378

[6]聶宇宏，陳海耿.氣體輻射寬帶關(guān)聯(lián)k模型的修正算法[J].計算物理，2002，19（5）：439～442

[7]WANG L, YANG J, MODEST M F, et al. Application of the full-spectrum -distribution method to photon Monte Carlo solvers[J]. Journal of Quantitative Spectroscopy and Radiative,2007,104(2):297～304

10.3969/j.issn.1001-8972.2012.07.020

國家自然科學(xué)基金(21006099)；鄭州市科技創(chuàng)新團(tuán)隊項目(10CXTD151)；鄭州輕工業(yè)學(xué)院博士基金資助項目(2008BSJJ017)

尹雪梅(1979-)，女，四川資中人，機(jī)電工程學(xué)院，博士。