浅谈CFD中离散方程的求解方法.pdfVIP

浅谈CFD中离散方程的求解方法 总装备部工程设计研究总院 满孝新4 哈尔滨建筑大学 高甫生o。 捅要：本文介绍了CFD中离散方程的求解方法，把复杂的三维大型矩阵简化为三对角带 状稀疏矩阵r采用TDMA法和逆向ID／vIA洁相结合+百『大太加快离散方程的收敛速度。 ， 、 关键字： ／cFD、)离散方程 (TDM射一护’ 、 7 ^ ^ 一、概述 土巴年来，随着工业、民用建筑复杂化和计算机计算性能的飞速更新，CFD已 经在热能动力、航空航天、机械、土木水利、环境化工等备个领域方面得到了广 泛的应用，日本调查[1]显示59％的企业已经利用CFD，还有26．4％的企业准备 使用CFD，我们国家对CFD使用率也逐步增加．CFD已经成为r稃研究与设计 中备受期待和不可缺少的技术手段。 目前，工业、民用建筑的复杂化、现代化和智能化，常常推陈出新，超出了 现有的规范，使得建筑改计、设备设计上缺乏经验和技术，迫切要求利用CFD对 工程进行预测、分行和评价，目前的设汁思路开始由单一工况瑷计向过程设计、 由宏观设计向微观设计转化．舒适、先进和节能己成为不可扭转的潮流。通过CFD 对气流的形象化图形演示，可以絮助工程技术人员进行优化设计，提高工程设计 的质量。 STREAM等，总体上说，这些商业软件对于一般的问题，基本上可咀得剑较为满 意的模拟结果。但是，这些商业软件普遍存在价格昂贵．计算周期太长，使用过 于复杂．通用性较差等问题，某些软件往往还需要对其进行改编以后才能使用． 所有这些都制约着的CFD发展。针对这些问题，还需耍对CFD的计算方法进行 优化．寻找一些便于求解，加快离散方程的收敛速度，防止迭代发散的求解方法， 以利于CFD的发展手u推广。 二、离散方程的模型分析及其简化求解 描述流体运动的基本方程组，主要是根据质量守衡定律、动量定理、能量守 衡定律、粘性规律推导出的连续性方程、运动方程和能量方程。事实上，由f 流体运动特别是紊流运动的复杂性和脉动性，使得方程求解非常困难，几乎是 不可能直接求解。于是，国内外的许多学者对方程进行了许多科学假设和简化， 产生了不同的计算模型．主要有K_-6、K·L、K·W、K．f、DSM、ASM等模型， 在室内气流模拟中应用较为广泛的是K．．￡模型，其对流场的描述主要有运动方 程、能量方程、K方程、s方程和压力修正方程等构成。 cn ☆通信地址： 100028北京4702信箱2号(010· E-mail：manxx@ihw．com 50008 ☆☆通信地址 哈尔滨建筑大学新匡644信箱 (0451．)6282504 CFD眄iq较为I泛的领域主要是人空间的气流组鲥模拟中，盯r一“大笠i 闭={勺模拟对象，其闷格的划分常常有儿万甚至几皇万个网格倚成．也裁意味着 对于每一个方程组是由儿万个(甚至儿百万个)方程均成，也可以说是由，lL万 【甚至儿|白万)相乘的大型矩阵构成。对丁一个实际肥摸拟对象，常常需要六、 七个这样的凡型方珲组相互耦台．要使这几个锅台方程组最终达到整体收敛是 非常困难的，稍育不慎戟可能引起迭代发散，或者收敛不到～个物理量上真实 的解。因此，选择一种简便有效的求解方法．显碍尤为重要。 对丁一个三维的模拟对象，面对如此巨最而复杂j0同格17点．可能会感判 有些茫然干¨无从r手，甚至感纠简育不可求钾。此时．我¨肯定台鹰j!如果理 是一维的问题，就简量多了。按照这个思想，我刺下面介7if～i’方法——交擎 方向的逐行迭代洼．可巧妙地将二维问题简化为一维问题进fi求解， 对离散方程求解的选代方注很多，在计算流体力学和I计算传热学甲常用构 方法大致可分为竞选代法、块迭代法L逐行法)、交替方向慧式逐}i迭代法及 强隐选代法等。透}i迭代法