计算辐射传热学发展.pptVIP

计算热辐射学的发展 一、概述 描述辐射能传输过程的基本方程为辐射传递方程。在三维半透明介质中,辐射强度是七维变量的函数，多数情况下只能通过数值计算的途径进行近似求解。数值模拟方法已成为半透明介质内热辐射传递理论研究和工程应用的重要手段。 均匀折射率灰介质内的辐射传递方程 若将s看作速度矢量,则 相当于对流项,因而辐射传递方程可以看作一个特殊的对流扩散方程. 二、辐射传递方程数值求解方法 1. 基于射线跟踪的方法 2. 基于微分形式辐射传递方程全局离散的方法 基于射线跟踪的方法 此方法一般需要通过跟踪光束传播轨迹来进行求解，如区域法、蒙特卡洛法和离散传递法。由于需要进行繁重的光束传播轨迹计算，这类方法计算速度极慢，不适合用来与热流体问题中的其他过程进行耦合计算。 基于微分形式辐射传递方程全局离散的方法 此方法的求解过程类似于一般偏微分方程的离散和求解，如离散坐标法、有限体积法、有限元法、谱元法、无网格法等。这类方法计算速度较快，适合与对流和导热过程耦合计算，是目前计算热辐射学研究的重点。 三、有限元法 基本原理： 将复杂求解域划分为若干规则形状的单元，然后分别在每个单元上构造插值函数对待求函数进行近似，最后结合加权余量法对方程进行离散和求解。 有限元法的优点 1. 对复杂形状求解域和边界条件的适应能力强 2. 与有限体积法和有限差分法相比,更易于在单元上实现高精度近似 3. 易于发展与其他模拟兼容的耦合分析 有限元应用于热辐射的发展历史 有限元法应用于辐射换热问题始于20世纪70年代。Reddy和Murty首次使用有限单元法计算了一维辐射传递的积分方程。 Razzaque等人采用有限元法研究了充满非散射性介质的正方形灰体空腔内的辐射换热。 有限元应用于热辐射的发展历史 传统的辐射换热求解的有限元法基于区域法的思想，其中辐射强度和辐射源项表示为一个与所有面元和体元相关的形式积分。这类有限元法在计算中需要进行十分复杂的几何积分和可视性的判断，尤其对于多维不规则几何问题，计算极为繁琐。 有限元应用于热辐射的发展现状 由于传统的有限元法存在的不足之处，目前的趋势是发展基于微分形式辐射传递方程离散的有限元法。微分形式的辐射传递方程可以看作一个特殊的对流扩散方程且具有强对流特性。 有限元应用于热辐射的发展现状 Fiveland在离散坐标法的基础上发展了一个基于偶宇称形式辐射传递方程的有限元法，首先将离散坐标法的思想与有限元法结合起来。该方法虽然有较好的稳定性，但由于偶宇称方程自身对各向异性散射处理的困难，因而很难用于处理各向异性散射问题。 有限元应用于热辐射的发展现状 Liu发展了可以模拟各向异性散射问题的有限元法。该方法建立在离散坐标和标准伽辽金有限元法的基础上，然而由于辐射传递方程的强对流特性，该方法存在一定的数值稳定性问题。 有限元应用于热辐射的发展现状 为了彻底解决原始辐射传递方程强对流特性引起的数值稳定性问题，最近，Zhao等人基于标准形式的辐射传递方程，使用一种与偶宇称分解不同的处理方式，推导了一个均匀折射率介质内以辐射强度为原始变量的二阶辐射传递方程。 有限元应用于热辐射的发展现状 该二阶方程中不存在对流项，取而代之的是一个二阶扩散项，因而完全消除了原始辐射传递方程中对流项带来的数值不稳定性。研究表明基于该二阶方程的有限元法具有很好的数值稳定性，并可以方便地用于求解吸收、发射及各向异性散射介质内的辐射传递问题。 四、其他热辐射计算方法 谱方法 谱元法 间断元法 无网格法 谢谢 * * r为空间位置矢量，s为辐射传递方向单位矢量，c0为真空中的光速，n为折射率，I为辐射强度，Ib为黑体辐射强度，κa为吸收系数，κs为散射系数，Φ为散射相函数，t为时间 由于辐射强度是波长、时间、空间位置坐标、角度方向的函数，一般难以获得理论解,多数情况下只能通过数值计算的途径进行近似求解。随着高速计算机的普及,目前数值模拟方法已成为热辐射传递理论研究和工程应用的重要手段。 * * *