传热学第4章 导热问题的数值解法.pptVIP

第四章　导热问题的数值解法 上海交通大学 （3） CFDNHT的区别与联系 计算传热学或数值传热学（Numerical Heat Transfer, NHT）与计算流体动力学（Computational Fluid Dynamics, CFD）的主要研究内容是一致的，不少文献把热流问题的数值计算一概称为CFD 。 就它们的发展史来说，计算传热的发展史在很大程度上也就是计算流体动力学的发展史 。 CFD与NHT之间也存在不少区别 ：计算传热学不讨论无粘流动及跨、超音速流动数值计算中的一些特殊问题（如激波的捕获等），而后者却是CFD中的一个重要研究内容。 5、CFDNHT发展简史 早在1933年，英国科学家Thom应用手摇计算机完成了对一个外掠圆柱流动的数值计算，但应用计算机和数值方法求解流动及传热问题在全世界范围内逐渐形成规模而且得出有益的结果，大致始于60年代。 从60年代开始至今，可以把CFDNHT的发展过程分为三阶段，每一个阶段中都有主要的发展事件。 ◆ 第一阶段，萌芽初创时期(1965－1974)； ◆ 第二阶段，开始走向工业应用阶段(1975－1984)； ◆ 第三阶段，蓬勃发展与大量工业运用阶段(1985－) （6）1972年SIMPLE算法问世 SIMPLE算法采用分离式的求解方法，解决了速度与压力的耦合问题。其基本思想是：先求解有关一个速度分量，而把其他作为常数，随后再逐一求解其它变量；在流场迭代求解的任何一个层次上，速度场都必须满足质量守恒方程，这是保证流场迭代计算收敛的一个十分重要的原则。 （7）1974年生成贴体坐标的方法(TTM)的提出 1974年美国学者Thompson、Thames和Mastin提出了采用微分方程来生成贴体坐标的方法(TTM)。为有限差分法与有限容积法处理不规则边界问题提供了一条崭新的道路——通过变换把物理平面上的不规则区域（二维问题）变换到计算平面上的规则区域，从而在计算平面上完成计算，再将结果传递到物理平面上。 TTM方法提出以后，逐渐地在CFDNHT领域中形成了“网格生成技术”这一分支，并成为目前世界上很活跃的研究方向。 （1）1979年Leonard发表了著名的QUICK格式。 这是一个具有三阶精度的（从界函面数插值而言）的对流项离散格式，其稳 定性优于中心差分。目前QUICK已在CFD/NHT研究与应用中得到广泛的应用。 （2）1979年美国Illinois大学Minkowycz教授主编的学术期刊“Numerical Heat Transfer”创刊。 （3）1979年Spalding发布的PHOENICS 1.0问世。 PHOENICS是英语Parabolic, Hyperbolic or Elliptic Numerical Integration Code Series的缩写（意为对抛物型、双曲型、椭圆型方程进行数值积分的系列程序）。 （4）1980年Patankar教授的名著“Numerical Heat Transfer and Fluid Flow”出版。 这本书内容精炼，说理透彻，注重物理概念的阐述，深受全世界数值传热的研究者与使用者的欢迎。出版后不久，被相继译成俄文、日文、波兰文及中文等，成为数值传热学领域中的一本经典著作。 （5）1981年英国的CHAM公司把PHOENICS软件正式投入市场，开创了CFDNHT商用软件市场的先河。 （6）1982年Rhie与Chou提出了同位网格方法。 随着计算机工业的进一步发展，CFD/NHT的计算逐步由二维向三维，由规则区域向不规则区域，由正交坐标系向非正交坐标系发展。于是，为克服棋盘形压力场而引入的交错网格的一些弱点，1982年Rhie与Chou提出了同位网格方法[50]。这种方法吸取了交错网格成功的经验而又把所有的求解变量布置在同一套网格上，目前在非正交曲线坐标系的计算中得到广泛的应用。 （7） SIMPLER与SIMPLEC算法 关于处理不可压缩流场计算中流速与压力的耦合关系的算法，在这一段时期内也有进一步的发展，先后提出了SIMPLER，SIMPLEC算法。 1. 边界节点离散方程的建立： qw x y qw (1) 平直边界上的节点 * - * - 第4章 导热问题的数值解法———§4-2边界节点离散方程的建立 (2) 外部角点 x y qw * - * - 第4章 导热问题的数值解法———§4-2边界节点离散方程的建立 (3) 内部角点 * - * - 第4章 导热问题的数值解法———§4-2边界节点离