非结构网格中非牛顿流体在偏心圆环管中数值传热计算.pdfVIP

非结构网格中非牛顿流体在偏心圆环管中 数值传热计算 马倩，张敏，许彬 南京理工大学动力工程学院，南京（210094 ） 摘 要：在非结构化网格中，用基元中心有限容积法进行非牛顿流体在偏心圆环管中的数值 传热计算。在简单充分发展层流和不同的边界条件下，对圆环管不同的偏心率和半径比，寻 找温度场的分布规律，所得结果对石油化工等工业十分有益。 关键词：非结构网格，有限容积法，非牛顿流体 1. 引言 非牛顿流体在偏心圆环管中的流动与传热现象，在石油、化工和食品加工等工业中十分 普遍。由于制造加工和振动等原因，使得原本同心的环形通道发生了偏心（如图 1 所示）， 这也使得流动区域具有不规则性。 在以前的数值传热计算中，往往采用结构化网格和有限差分方法，为了保证网格的正交 性，必须进行坐标变换，这使得控制方程、边界条件以及准则数等都要转化成新坐标系下的 形式[1-3]，这些都使得计算工作复杂和困难。非结构化网格对不规则区域有很强的适应性， 它对方程的离散在直角坐标系下进行，既可适用于一切其它形式的坐标系，因此相对于结构 [3-7] 化网格来说，非结构化网格有其独特的优势 。 图1 偏心圆环管截面几何形状和坐标系 2. 基本控制方程和离散 在直角坐标系中，非稳态的扩散方程或导热方程，对一个标量物理变量φ 可写成： ∂φ ∂ ∂⎛⎜ φ ⎞⎟ ρcp ∂ ∂ ⎜Γ∂ ⎟+Sφ (1.1) t x x i ⎝ i ⎠ 其中，Sφ 是单位体积中的净源项，Γ 是对应于变量φ 的扩散系数。在此仅对非结构化 -1- 网格的离散方程进行讨论，控制方程（1.1）可以写成： ∂φ ρcp D +S (1.2) ∂t 在一个控制容积P 中，积分上式有： ∂ nb ρc φdV ∑D + ( S + S φ ) ∆V (1.3) ∫∫V p i c p p p ∂t p i 1 其中，Di 是扩散项，它可以表示为基本扩散项Dp ,i 和二次扩散项Ds ,i 之和： Di Dp ,i +Ds ,i (1.4) 它们的物理意义分别是某一交界面上的一次扩散项（基本扩散项）和二次扩散项（切向 扩散项）。对于正交网格二次扩散项为零，对某交界面i 它的计算方法如下： → → φ φ ( − ) A i •A i Dp ,i Γi Eds P →