PISO算法的实现及与SIMPLESIMPLERSIMPLEC算法收敛性的比较.pdfVIP

第28卷第4期 华北水利水 电学 院学报 V01．28No．4 ChinaInstjtuteofWa￡er and Power 2007年8月 Jounlal《Nor￡h ConservancyHydmelectric Aug．2007 文章编号：1002—5634(2007)04一0033—04 SIMPLEC算法收敛性的比较 王为术1…，徐维晖1，翟肇江1，孟庆东2，陈听宽3，罗毓珊3 (1．华北水利水电学院，河南郑州4500ll；2．山东省阳谷县供电公司，山东阳谷252300 3．西安交通大学，动力工程多相流国家重点实验室，陕西西安710049) 摘 要：针对PISO算法的实现程序，通过2个典型流动换热例子，验证了其正确性，并将PIS0与sIMPLE， sIMPLER，SIMPLEC算法的收敛性和健壮性进行了比较．结果表明，不同算法的收敛性和健壮性与具体的问题 PLEc耗用的CPU时间少；PIs0算法的健壮性因不同的问题而异． 关键词：PIsO算法；收敛性；健壮性 中图分类号：TV543+．6文献标识码：A 知的速度场通过求解压力Poisson方程获得． 自SIMPLE算法提出以来¨o，又相继有sIM— PLER，SIMPLEc，sIMPLEx，PIs0等改进算法提出．3．求解动量离散方程以获得M’，”+． sIMPLE系列算法不仅在不可压缩流体Ⅳ一s方程4．求解压力修正值p’方程，要求与p’相应的 数值求解中广泛应用，而且已成功推广到可压缩流 “’，秽’能使(“++“’)，(秽++秽’)满足连续性方程． 场的数值计算中，可计算任何流速的流动．因此，了 解SIMPLE系列算法原理及特性有很重要的实际意 d。(p：一p：)，d。=A。／口。，邻点速度的影响忽略不计； 义．笔者在sIMPLE教学程序基础上实现PISO算 法，用2个典型二维流动传热问题检验该算法，并比 邻点速度修正值的影响已基本考虑在内． 较PIsO和SIMPLE，sIMPLER，SIMPLEc等算法的收 6．采用下列速度、压力值作为开始下一层次的 敛性和健壮性． 迭代计算(设速度的亚松驰已经组合在u，勘求解过 1 SIMPLE系列算法 程中)： 在Ⅳ一5方程的求解方法中，SIMPLE系列是以 压力为基本变量的原始变量法中的压力修正法． 7) SIMPLER的改进在于压力由假定速度通过求解压 SIMPLER：(“++u’)，(口’+秽’) 力Poisson方程获得．sIMPLEC改进了速度修正式， 解决了速度修正不协调一致问题，其压力不再需亚 松驰∽J．SIMPLE系列算法的主要计算步骤： 1．由假定的或上一次计算得到的速度u∞’，移‘0’ 确定动量离散方程的系数及常数项． 修正，可认为具有两步修正的算法．但其第二步修正 2．确定压力场．对于sIMPLE，sIMPLEC采用假不是通过求解方程，而是显示地进行，不是典型的2 步校正算法． 定的或用上次计算的值p+；对于sIMPLER则由已 收稿日期：2006一11—15；修订日期：2007—01一10 基金项目：华北水利水电学院青年基金项目(HsQJ2006003)． 作者简介：王为术(1972一)，男，重庆开县人，华北水利水电学院讲师，博士，主要从事多相流动和传热方面的研究 华北水利水电学院学报