一维各向同性散射系统内辐射传递分析的ES-DRESOR法

摘要

同其他传热方式相比,辐射传热包含了发射、吸收、散射、反射等复杂过程,其控制方程为微分积分方程。在辐射换热问题分析以及实际应用中,辐射传递方程的求解和辐射强度的计算至为关键。到目前为止,已经发展了多种方法用来求解辐射传递方程,如离散坐标法、蒙特卡洛法、热流法、区域法、球形谐波法等,这些方法各有优缺点,根据其特点被用于不同的辐射问题求解中。在一些辐射问题的研究中,方向辐射强度信息显得格外重要。基于蒙特卡洛法发展起来的DRESOR(Distributions of Ratio of Energy Scattered by the medium Or Reflected by the boundary surface,被介质散射或者被壁面反射的能量份额分布)法,因能够提供高方向分辨率的辐射强度,而被应用到一系列的辐射传递问题的研究中,如工业炉内燃烧温度场的可视化监控、平行脉冲入射辐射、各向异性散射、梯度折射率介质中等。然而因该方法基于蒙特卡洛法,计算耗时不可避免,同时存在着统计误差,限制了此方法的进一步发展。
　　本文致力于提高DRESOR法计算效率和计算精度的研究,在基本DRESOR法的基础上发展了ES-DRESOR(Equation-solving DRESOR)法。通过求解线性方程组直接得到DRESOR数(DRESOR法求解的核心),无需像先前那样由蒙特卡洛法的能束跟踪得到,以实现计算效率和计算精度上都能有大幅度的提升。本研究对于提高DRESOR法的分析求解能力,拓展其应用范围具有重要意义。具体来说,主要工作如下:
　　在基本DRESOR法的基础上提出了ES-DRESOR。给出了该方法的计算原理,分为透明壁面和漫壁面参与两种情况,对于前者,设定系统内某微元体具有单位黑体发射,其它部分没有发射,这样对DRESOR法给出的辐射强度简化后,可由入射辐射(对辐射强度的空间积分)的两种计算方法(即DRESOR法给出的计算公式和标准计算公式)所得结果相等,建立关于DRESOR数的方程。在漫壁面参与辐射过程后,对DRESOR法计算入射辐射表达式进行了修正,加入壁面的贡献部分。对壁面单元建立方程的原则是被壁面反射的总能量等于其它单元散射或反射到此壁面的总能量;对空间单元建立方程的原则依然是由两种方法计算的入射辐射相等。
　　在一维吸收发射各向同性散射的辐射系统内,对ES-DRESOR求解辐射传递方程进行了研究。数值部分,首先对ES-DRESOR法的可靠性进了验证,包括与文献中的精确解对比,与广泛使用的离散坐标法对比,以及将ES-DRESOR法计算的结果代入到辐射传递方程进行检验。而后将ES-DRESOR法计算出来的DRESOR数与蒙特卡洛法计算的结果进行对比,显示了很好的吻合性,且消除了蒙特卡洛法中因统计误差存在而导致的“抖动”现象。对比ES-DRESOR法同蒙特卡洛法计算出来的辐射强度相对误差,前者的精度高出一个数量级。更为重要的是ES-DRESOR法的计算效率在测试网格数下,与基于蒙特卡洛法的基本DRESOR法相比,有两个数量级的提升。

目录

封面

声明

中文摘要

英文摘要

主要符号表

目录

1 绪 论

1.1 课题研究背景

1.2 辐射传递方程求解方法研究现状

1.3 DRESOR法的发展历程及研究现状

1.4 本文的主要工作

2 ES-DRESOR法的求解原理

2.1 引言

2.2 基本DRESOR法

2.3 ES-DRESOR法的基本原理

2.4 本章小结

3 一维透明壁面下ES-DRESOR法研究

3.1 引言

3.2 网格离散及辐射强度公式

3.3 方程组的建立

3.4 数值结果及分析

3.5 本章小结

4 一维漫壁面条件下ES-DRESOR法求解

4.1 引言

4.2 辐射强度离散

4.3 组建方程组

4.4 数值结果及讨论

4.5 本章小结

5 全文总结及展望

5.1 全文总结

5.2 后续工作的展望

致谢

参考文献

附录A 攻读博士学位期间的科研成果

附录B 方程组AX=B的系数矩阵