圆环域中结构与非结构网格的热传导数值计算.PDFVIP

东 北 电 力 学 院 学 报 第 24 卷第 6 期 　 Journal Of Northeast China Vol. 24 ,No. 6 2004 年 12 月 Institute Of Electric Power Engineering Dec. ,2004 ( ) 文章编号 :1005 - 2992 2004 06 - 0015 - 05 圆环域中结构与非结构网格的热传导数值计算 1 2 1 张　敏 ,戴春林 ,范悦宏 ( 1. 南京理工大学 动力工程学院 ,江苏 南京 210094 ;2. 吉林市民营经济发展局 ,吉林 吉林 132012) 摘 　　　要 :用有限容积方法求解圆域和圆环域中的热传导问题。同时采用结构化网格和非结构化网 格进行数值计算 ,其中结构化网格采用极坐标 ,而非结构化网格采用笛卡尔直角坐标。所有数值计算结 果同精确解相比较 ,并得到相一致的结果。结构化网格求解问题因受求解区域几何形状的限制而需选 择不同的坐标系进行方程离散 ,而非结构化网格则没有这种限制。 关 　键 　词 :结构网格 ;非结构网格 ;热传导 ;有限容积法 中图分类号:O 552. 4 + 22 　　　　文献标识码 :A ( ) 数值网格生成方法 Numerical Grid Generation Method ,N GGM 着眼于网格节点和单元特征量的离 散数据。一般而言 ,数值网格生成方法大致可分为结构和非结构两种。所谓结构 ,就是网络数据的生成 是有数学表达式的 ,并且这种表达式和对应的离散方程依赖于所选择的坐标系而不同。它的特点是序 列性和数据结构上的顺序性质。非结构网格既没有数学上严格的递推关系。对于对应的离散方程 ,可 采用最简单的直角坐标系。它的数据结构具有一定的随机性和自适应特性。尽管结构化网格能采用某 些数学变换 ,如代数法、无限插值法、多面法等方法处理非规则边界问题 ,但它更多地表现为规则和有 序 ;而非结构化网格则反映出无规则和无序性[1 ,2 ] 。 本文通过基本导热方程在圆域和圆环域中的数值求解来论述结构化和非结构化网格在这方面的应 用。 1 　基本导热方程和离散方程 对一个标量物理变量 , 一般的稳态扩散方程可写成 1 Γ ·f ( x ) ( ) i + S = 0 1 f ( x i) x i x i ( ) Γ ( ) 其中, f x i 是几何形状因子, S 是单位体积中的净源项, 是对应于变量 的扩散系数。根据方程 1 , 在笛卡尔直角和圆柱坐标系中, 导热方程可分别为 ( )