谱元法求解辐射传递方程的各种稳定格式比较研究.pdfVIP

中困工程热物理学会’ 传热传质学 学术会议论文 编吁：063260 谱元法求解辐射传递方程的各种稳定 格式比较 赵军明刘林华 (哈尔滨工业人学能源科学与丁程学院，黑龙江哈尔滨150001) Tel．：045 1Email：Ihliu@hit．edu．cn 摘要：研究比较了最小二乘稳定格式和三种流向迎风型稳定格式，即流向迎风格式、彼得罗大伽辽会 流向迎风格式和伽辽金最j,--乘格式，及一种各向删性格式对于讲冗泫求解辐射传递方程的低阶多项 式稳定特性和商阶多项式稳定特性．结果表明，当多项式阶数较低时，最,j,--乘格式及流向迎风型稳 定格式显示f{j了一定的稳定效果。对于流向迎风型稳定格式，迎风【大J予的增加其稳定性能有所增强。 但当多项式阶数较商时，最／j、乘格式及流向迎风型稳定格式均不能起到自．效的稳定作用。在多项式 高和低的两种情形，各向同性稳定格式均显示出了较好的稳定性能，优于其他格式。随者迎风网了的 增加，扩散效应增加，由于不一致性，各向同性稳定格式的结果有所变差。 关键词：潜元法，辐射传递，迎风格式 1 引言 ． 由于求解半透明介质内的辐射传递具有很广泛的工程应用背景，发展求解二卜透明介 质内辐射传递的高效数值方法一直是一项重要的研究内容。近几十年来，对予半透明介 质内的辐射换热，各国学者做了大量的研究，已发展了许多辐射换热的计算方法，主要 有蒙特卡洛法、区域法、离散坐标法、有限体积法、有限元法等。综合来看，离散坐标 法，有限体积法具有计算效率高，易于求解多维辐射传递，易处理复杂边界条件及可以 方便地处理散射项等优点因而得到了广泛的关注与应用。但这些方法都是低阶方法(一 般仅具有l阶或2阶的精度)，同时仅具有h收敛性，即为了获得更高的精度只能通过 加密网格来实现，这往往需要繁杂的网格重新划分的操作。 谱元法首先由PateraIl】提出并在之后被成功的应用于计算力学和传热等领域。该方 法具有细收敛特性。为了克服低阶方法的不足，谱元法受到广泛重视弗开始应用于半 透明介质内辐射传递的求解。由于辐射传递方程为可以看作是一种特殊的对流扩散方 程，且具有强对流特性，这类方程的解往往会出现人梯度甚至间断。对于该类方程的求 解需要特别的稳定技术，比如设计迎风格式，添加人：’r扩散项等，否则结果可能会出现 折射率介质内的辐射传递。赵军明和刘林华pJ发展了一种最小二乘Chebyshev谱元法并 应用·丁-求解多维?卜透明介质内的辐射传递。这些+l：作显示出谱元法求解一维及多维辐射 传递问题具有很女f的前景，对于具有光滑解的辐射传递问题，谱元法的P收敛特性显示 nl了指数阶的收敛速率。 对丁．谱元法应川较为J“泛的稳定格式包括最：j,--乘法(LeastSquare，LS)，主要还 仃流向迎风格式(Streamline 552 Petrov UpwindGalerkin，SUPG)14】，伽辽金最小二乘格式(Galerkin 除了这些流向扩散的稳定格式之外，还有各向同性扩散稳定格式。对，丁．各向同性扩散， 一般为住方程Zi端添加一正枯性的各向同性扩散项。就我们所知，目前还没有对谱元法 求解辐射传递方程的各种稳定格式进行比较研究的J：作。 为了便于实际应用中对各种稳定方案进行择优选择，本文对最小二乘格式，上述三 种常JH流向扩散型格式及一种各向同性稳定格式应用‘r谱元法求解辐射强度存在人梯 度的辐射传递问题的性能进行比较。 2 离散坐标形式的辐射传递方程 考虑吸收、发射和散射灰介质内的辐射传递，直角坐标系下离散坐标形式的辐射传 递方程可以写为 戗 鲫 q万孑j 矿等叫7等∥娶Oz=砭厶一∥+番薹，∥卅w∥ 对于不透明漫射边界给定边界条件为