矩形通道充分发展层流的精确与数值解.pdfVIP

矩形通道充分发展层流的精确和数值解 1 2 1 1,2 1 张敏 ，John C. Chai ，张钧波 ，许彬 ，刘少培 1 南京理工大学动力工程学院，南京 (210094) 2 南洋理工大学机械与航天工程学院，新加坡 (639815) E-mail ：zhangm@ 摘 要：用有限容积法，在结构/非结构化网格中，求解矩形通道内充分发展层流流动的速 度分布和流阻情况，并将计算结果与精确解进行比较，从而证明本方法和所编写程序的正确 性和实用性。 关键词：结构/非结构网格，矩形通道，充分发展层流，有限容积法 矩形通道内的充分发展层流流动在工业界有着广泛的应用，并且对它的研究由来已久。 如美国学者Patankar[1-2] [3-4] [5-6] 和Manglik 等，中国学者陶文铨 等。由于这种流动自身具有的 特性，使得其动量方程可以得到很大的简化，因而寻找其精确解变为可能。这些正好可以为 [7-8] 数值计算的正确性提供参考依据 。 1. 基本概念和方程 所谓充分发展的层流流动，是指通道内截面上流体的无量纲速度分布与流动方向上的坐 标无关的这样一种流动状况。这种状况大致可以分为简单的和复杂的两大类。所谓简单的充 分发展层流流动是指在垂直于主流方向的截面上，速度分量为零的这一类情况。此时求解的 纳维-斯托克斯方程可简化为导热型方程[7]，即z 方向的动量方程可以化简为如下形式， ⎛ 2 2 ⎞ µ⎜∂ w +∂ w ⎟−dp 0 (1.1) ⎜⎝∂x 2 ∂y 2 ⎠⎟ dz 将动量方程(1.1)化为无量纲直角坐标形式如下， ∂ ∂W ∂ ∂W ⎛ ⎞ ⎛ ⎞ ⎜ ⎟+ ⎜ ⎟+1 0 (1.2) ∂X ∂X ∂Y ∂Y ⎝ ⎠ ⎝ ⎠ 其中，无量纲坐标X ，Y ，以及无量纲速度W 定义如下， x y µw X ,Y ,W 2 (1.3) D D − dp / dz D ( ) 其中D f 为通道截面的某一特性尺寸，摩擦因子 ，以及雷诺数Re 定义如下， dp −D h dz ρwDh f ，Re (1.4) 1 2 µ ρw 2 Dh 为当量直径, w 为横截面上