有限体积法立体角离散误差的研讨.pdfVIP

中国工程热物理学会 传热传质学 86 学术会议论文 编号：0631 有限体积法立体角离散误差的研究木 张昊春 谈和平 甄灸 哈尔滨工业大学能源科学与工程学院，150001 1．8641 (电话：045 摘要：采用有限体积法求解辐射传输方程时需要同时对空间坐标和立体角进行近似处理，从而引起离 散误差，对于立体角离散误差的研究亟待深入。本文针对有限体积法立体角内热流密度的求解过程：l{ 发，运用傅立叶级数推导了近似误差的表达式，并分析影响离散误差的各种因素。包括离散射线数， 辐射强度分布，假散射及网格比。根据分析结果，采用辐射强度的连续分布模型，模拟了各种因素对 误差的影响。结果表明，辐射强度的分布函数对于有限体积法的离散误差有一定程度影响。当增加网 格数及减小网格尺寸比时，可以有效地减小离散误差。另外，假散射与离散误差问呈较强的非线性耦 合关系。 关键词；辐射换热有限体积法离散误差立体角假散射 1．．引言 Volume 近年来，求解辐射换热的有限体积法lI删(FiniteMethod，FVM)得到了快速 发展，广泛地应用在多种热辐射传输问题中p矧。FVM方法的基本思想是将4zr立体角 空间离散为不重叠的立体角，在控制体y和控制立体角国内对辐射传递方程积分，考虑 立体角的能量平衡关系，得到缈内辐射能量守恒方程的有限体积表达式。 作为一种数值计算方法，FVM在求解过程中，对空间坐标和立体角同时作了近似 处理，从而产生离散误差【71。近年来，FVM方法离散误差的研究成为计算热辐射界普遍 关注的问题【聃1。相对而言，对空间坐标离散引起假散射得到了较为深入的研究[1,3,6,8-9】， 但对空间立体角离散误差(在计算热辐射领域，由于空间立体角离散所引起的误差，通 常称之为射线效应)的深入研究则很少。FVM方法的创始人Q D．Raithby指出【81，由于 立体角离散所引起的误差的研究没有得到足够重视。除此之外，空间坐标离散引起的误 差同立体角离散误差间的交互影响也是FVM中出现的特殊误差现象，对此也需要进行 深入的分析。 本文针对FVM方法离散立体角内热流密度的求解过程，应用文献【10】中提出的辐 射强度连续分布函数模型，分析了影响有限体积法立体角离散误差的各种因素并做数值 模拟，并对假散射对立体角离散误差的影响做了讨论。 ‘’田末白然科擘基金重点资助项目lNo．503360l∞夏田容重点基础研宄{973专项经费{N02003CB214500)蚤助 631 2．空间立体角离散误差分析 本文主要研究立体角离散误差，因此本节对FVM的详细求解过程．-7以参考文献 度的精确解可以表达为： 时印啪躲僻．j--警仨簿邶捌cos“删删‘㈣ 该立体角的中心处的坐标为(oc，九)，在FVM中，采用如下方式进行近似计算： (2) 【g；Ⅲ一≈，(包。．，，仍√)cosec，』sint?c√sin(A0)A矽 利用Taylor级数，在(铭，允)点将，(口，≯)展开得到： ，(口，≯)=，(如，芘)『+(口一亿)·(V／o)c+(≯一九)·(V厶)c+… +志(p吲¨旷1械+志(≯一∥～矿劫c。’ 忽略高阶项j。仅保留二阶，这时，由于立体角奄散所造成的误差为钳： ‘ (4