方形通道充发流动中的精确解和数值解.pdfVIP

方形通道充发流动中的精确解和数值解 1 2 1 张敏 ，John C. Chai ，许彬 1 南京理工大学动力工程学院，南京 (210094) 2 南洋理工大学机械与航天工程学院，新加坡 (639798) E-mail ：mz2455@163.com 摘 要：用有限容积法，在结构/非结构化网格中，求解正方形通道内充分发展层流流动的 情况，并将计算结果与精确解进行比较，从而证明本方法和所给程序的正确性和实用性。 关键词：结构/非结构网格，正方形通道，充分发展层流，有限容积法 正方形状通道内的充分发展层流流动在工业界有着广泛的应用，并且对它的研究由来 已久。如美国学者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 为横截面上的平均流速，定义如下，