一种新的宽带k分布模型.docxVIP

一种新的宽带k分布模型 0 理论内固定cg 气体辐射计算在工程应用中起着重要作用。气体分子通过一定数量的振动-旋转光谱带之间吸收和发射辐射。光谱带中包含数万条光谱线，因此整个光谱吸收体系的传输非常快。rbl的计算是最精确的。然而，对于每个光谱带的详细光谱特性参数，如光谱位置、光谱强度、光谱半宽、光谱梯度能量，必须是光谱系列的详细光谱特征。以智能出版社（2000）为例，有效载荷为1032009条，ho有1283468条光谱。由于其巨大的计算量,一般只将LBL结果作为检验其它模型准确度和有效性的基准解.工程应用中,综合考虑计算精度和速度,宽带模型是比较合理的选择.传统宽谱带模型(EWB)是将小波数间隔内的平均吸收系数重新排列后得到的近似关系式,误差很大,且因其用谱带有效带宽表示辐射特性,和路径长度有关,只能用射线踪迹法求解辐射传递方程,适用于黑体壁面内非散射介质.现有的宽谱带k分布模型基本都是建立在指数宽谱带模型基础上,并作了一些假设,使其精度和适用性受到限制. 气体辐射特性计算只占整个辐射计算的小部分,应尽可能地提高其计算速度及与辐射传递方程求解方法的兼容性,本文建立了一种新的宽带k分布模型,从高温气体的数据库得到模型参数,准确快速并能与任意辐射传递方程求解方法相容. 1 谱带密度的影响 在小的谱带间隔内普朗克函数基本不变,而吸收系数在很窄的谱带内剧烈变化,吸收系数值在波数间隔很小的情况下多次相同,每次都得到相同的辐射强度和热流。k分布模型就是将吸收系数重新排列成平滑单调上升的函数,对相同的吸收系数只进行一次强度计算。窄带k分布模型,波数间隔小,比较准确,但其计算量大,不太适合工程应用。考虑精度和计算速度,不少学者将k分布方法应用于宽谱带,但基本都建立在指数宽谱带模型基础上,使其精度和适用性受到限制。 表1是H2O重要谱带的谱带参数,ηu、ηl和η0分别为谱带的上下限波数和中心波数.表2是不同温度下,最大普朗克函数Ibηmax和最小普朗克函数Iηmin差值与谱带中心波数η0的普朗克函数Ibη0的比值Idiff。从表2可以看出:由于宽谱带的波数间距较大,假设普朗克函数不变会导致很大的误差。 Modest等人提出了全光谱k分布模型(FSK),使k分布模型能在整个光谱应用,但由于谱带范围太大,温度对不同波数区间影响不同,导致其对非等温介质误差很大.本文考虑从高温气体数据库HITEMP得到宽带k分布模型的模型参数,提高其精度和适用性.由于宽谱带的波数间隔较大,引入普朗克函数加权f(T,k): 式中,IbΔη和Ibη分别表示谱带区间总的黑体强度和光谱黑体强度,δ(k-κη)是Dirac-delta函数。 对式(1)积分可得累积分布函数g(T,k), 对于任意的只依赖于吸收系数的辐射变量φη,它在波数Δη上的积分可以变换为: 常用Gauss型的数值积分方法来计算,其中7点Gauss-Lobatto积分应用最广泛,考虑到波数间隔较大,本文采用7点和12点Gauss-Lobatto积分: 式中,N是积分点数,wi是求积的权. 以H2O的2.7μm谱带为例,图1所示的光谱吸收系数随光谱剧烈变化.对应的累积分布函数如图2所示,吸收系数随累积分布函数变化是条平滑曲线,可以方便地进行曲线拟合.图1和图2的计算条件为:总压力p=1.013 MPa,H2O分压力pH2O=0.18p,温度T=1700 K。 2 同分法的比较 有效带宽是宽带模型的重要参数,定义为 式中,L是路径长度.对于宽带k分布模型,有效带宽可按下式计算: 图3～图5给出了宽带k分布模型计算的水蒸气3个重要谱带的有效带宽,并与逐线计算等计算方法进行比较.计算条件与上节相同。从图中可以看出,对于宽带k分布模型,12点Gauss-Lobatto求积结果和7点Gauss-Lobatto求积结果基本没有差别。与逐线计算比较,除6.3μm谱带最大误差大约为10%外,其它谱带本文提出的宽带k分布模型结果和逐线计算结果吻合很好.指数宽谱带模型在所计算的谱带上误差均很大,相对误差高达50%。文献应用Marin等人提出的宽带关联k分布模型计算了相同条件下H2O的2.7μm谱带的有效带宽.图4对本文提出的宽带k分布模型结果和文献中的关联k分布模型结果作了比较,结果显示,本文提出的宽带k分布模型结果更准确. 3 谱带有效带宽 本文建立了一种新的宽带k分布模型,从高温气体数据库得到模型参数,容易实施,且提高了计算精度.用该模型计算了H2O的谱带有效带宽。本文提出的宽带k分布模型结果和逐线计算结果吻合很好,比以前的宽带k分布模型更准确.对于宽带k分布模型,12点Gauss-Lobatto求积结果和7点Gauss-Lobatto求积结果差别不大。