离散坐标SRAPN法及其与TN法的比较.pdfVIP

第 15 卷 第 6 期 　 计 　　算 　　物 　　理 　　 Vol . 15 ,No . 6 1998 年 11 月 　 CH IN ESE J OU RNAL OF COM PU TA T IONAL PH YSICS 　　 　Nov . , 1998 　 离散坐标 SRAPN 法及其与 TN 法的比较 李本文 魏小林 　 (浙江大学热能工程研究所 ,杭州 　3 10027) ( 中国科学院力学研究所 ,北京 　100080) 摘 　要 　提出了一种新的坐标离散方法 —球带等差数列微元等分法 。并经双动量理论检验 ,与较 精确的区域法数值解进行了比较计算 。 关键词 　辐射换热 　离散坐标法 　数值计算 中图分类号 　T K12 1 0 　引　言 十多年来 ,关于离散坐标法[ 1 ] 的研究主要集中在以下几个方面 : ①从一维到多维 、方腔式 到圆筒式的应用 ; ②介质从单一吸收 —发射到非灰或各向异性散射 ; ③对有遮蔽效应等复杂几 何系统的应用 ; ④差分格式的构造及对计算精度的影响; ⑤射线效应和假散射的分析 ; ⑥优化 离散坐标系列及新的坐标离散方法的研究 ; ⑦离散坐标法在曲线坐标中的应用 。 本文提出了一种新的 SRA PN ( successive sp herical rings arit hmetic p rogression) 坐标离散方 ( ) 法 —球带等差数列微元等分法 。经双动量 Double Moment s 理论检验以及典型算例的辐射换 热计算表明 ,球带等差数列微元等分法具有精度高 、易于实施等优点 。 1 　球带等差数列微元等分离散坐标法 球带等差数列微元等分法的基本思想是将单位球体表面划分成许多面积相等的球面微 元 ,然后确定各表面所对应立体角的重心 。连结单位球体的中心与各立体角之重心 ,形成一个 离散坐标系列 。由于离散坐标的反射对称性 ,其构造可以在单位球体的一个卦限内进行 。 ( 首先 ,将单位球体的半球从上到下划分成许多球带 最顶端为球冠 , 图 1 中只画出一个卦 ) 限 ,然后从球顶端开始 ,各球带等分为更多的球面微元 。各球带上的微元数 目从上到下按等 ( ) 差数列分布 。根据作者的数值实验 ,第一个球带 实为球冠 分成 2 等份 , 以下球带按 3 ,4 ,5 … 等差数列划分为宜 。 离散坐标方向的计算属于纯几何问题 。求出每一个球面微元所对应立体角的重心 。坐标 原点与各重心的连线即为对应的离散坐标方向。 各离散坐标的权值即为对应的球面微元面积 。例如 ,若将球面从上到下分成 N = 5 个球 ( ) 带 ,则一个卦限的微元总数为 2 + 3 + 4 + 5 + 6 = 20 ,那么这 20 个离散坐标具有相同的权值 ( π ) / 2 / 20 。 　收稿 日期 ;修回日期 :1998042 1 　李本文 　男　33 　副教授 　博士 © 1994-2009 China Academic Journal Electronic Publishing House. All rights reserved. 736 计 　　算 　　物 　　理 第 15 卷 　 2 　离散坐标的双动量检验及与 TN 法的比较 　　最近 Thurgood 等[2 ] 提出了一种比传统方法