CLEAR与SIMPLE系列的三种算法的收敛性和健壮性比较.pdfVIP

收敛标准为：船。：?痸?诼?局褂裥?蔓丛型丝立?詜≤?畂×?胡的收敛性和系列的三种算法健壮性比较????隨???张京兆屈治国何雅铃陶文铨???年得到了不断的发展，相继提出了???舻副和?????等算法。????惴ㄖ饕???荡?矸椒????算法自??年由?????蚐????岢鲆岳础?ǎ?诩扑懔魈宥?ρШ?数值传热学领域得到了广泛的应用。???算法引入了包括速度和压力初场独立假定、忽略邻点速度修正影响等简化假定【??】。正是对这些假定的改进使???算法在随后几改进了???方法中速度与压力初场独立假定的近似处理，提高了速度与压力初场的协调性；????惴ǜ慕?薙???椒ㄖ泻雎粤诘闼俣刃拚?跋斓募蚧??恚?欢ǔ?度上改善了算法的收敛特性。???惴ㄖ菔墙?晏岢龅囊恢中碌乃惴āＫ?肷鲜鯯???盗兴惴ú煌??挥薪?压力修正方程，而是通过直接改进压为来计算改进的速度，同时解决了上述???算法的两个简化假定。文献??詂??算法和????惴ㄗ髁吮冉舷晗傅氖樟残员冉希?得出与????啾菴??算法收敛性明显占优的结果；但从比较数据来看???惴?在一定程度上没有???八惴ń∽承院谩?对于???系列算法的比较【??】文献较多，结论为没有一种算法绝对占优，一般情况下???算法健壮性较差，密网格下????惴ㄊ樟残越虾谩?本文在文献??幕?∩?程序正确性考核同文献??，通过六个二维例子着重比较???惴ㄓ隨???????、????惴ǖ氖樟布敖∽程匦浴?为了合理的比较几种算法，本文所有算法的数值处理方法都保持一致。具体数值处?对流格式为乘方格式；?代数方程求解采用交替方向隐式算法辅以块修正技术；其中，“’，?为动量方程解出的临时速度值；?０凑瘴南?第??车亩ㄒ迳瓒ā???????瑉??????．??理方法如下：中国工程热物理学会??学术会议论文传热传质学编号：??西安交通大学能源与动力工程学院，???摘要：通过六个典型流动换热问题对???????????和????闹炙惴ń?惺樟残院徒?壮性比较。比较结果表明，对于开口系统????惴ㄗ酆闲Ч?虾茫??湛谙低砪??算法综合占优。关键词：收敛性健壮性奇异点?疚氖芄?易匀豢蒲Щ?鹬氐阆钅孔手??．?????钟¨ 采用时步倍率?豪甲魑1硎舅惴ū冉辖峁?暮嶙?辍Ｊ辈奖堵视胙撬沙谝蜃拥亩??闹炙惴ū冉???俣妊撬沙谝蜃觓，，?冢??口；?????惴ǖ难沽ρ撬沙谝蜃涌冢???；??阶匀欢粤魉憷?露茸餮撬沙诳冢???·????算法的第二亚松弛因子孱；屈，??????凇躱．?，属?伲?∥????凇???网格系统已经网格独立性考核。本文通过六个常物性、稳态、层流、不可压缩二维流动换热算例来对几种算法进行比较。这六个算例为：?圆柱坐标突扩流动；?直角坐标外掠后台阶流动；?方腔顶盖驱动流；?环腔顶盖驱动流；?环形空间自然对流；?方腔自然对流。其中前两个算例为开口系统，后四个为闭口系统算例。六个算例的物理模型完全同文献???韵卤冉辖峁?鹑≈捣段?狾．?～???，．收敛所用迭代次数和??奔?主要是??奔?多少用来判断算法的收敛特性；能得到收敛解的时步倍率?嶙??范围大小用来比较算法的健壮性，获得收敛解的时步倍率范围约宽则健壮性越优。??开口系统算例比较?圆柱坐标突扩流动图??中四种算法的收敛特性变化趋向一致，????惴ǖ慕∽承宰詈茫琒???次之，????虲??较差。???谛⊙撬沙谝蜃?速度?笔樟残宰詈茫?被∷沙??亚松弛因子的设定：应关系如表??尽?网格系统：????琑???不同时步倍率下的迭代次数及所用??奔浔冉辖峁?缤?，??尽?表?撬沙谝蜃印Ｓ胧辈奖堵蔈对应表主??图?仓??晖焕┝鞫?樟驳??问?冉?图?仓??晖焕┝鞫?樟彩奔浔冉??豢??????巴霉一’? ．??差?；芝墨因子大于??以后，?獬??惴ㄊ樟蔡匦宰詈茫琒?琇?算法最差。不同时步倍率下的迭代次数及所用??奔浔冉辖峁?缤?，??尽?从图??可知，对直角坐标外掠后台阶流动，????惴ǖ慕∽承宰詈茫琒?芁?和????次之，???畈睢???在小弧松弛因子时收敛性最好，当亚松弛因子大于??以后，?鮌??惴ㄊ樟蔡匦宰詈茫琒???算法最差。以上两个开口系统算例可以看出，?御??惴ǖ慕∽承宰詈茫?渌??炙惴ń∽?性相当：????惴ㄔ诔Ｓ醚撬沙谝蜃???～???段?谑樟残哉加牛?ㄒ榭?谙?统使用?押??惴ā１疚牡谋冉辖峁?胛南譡??慕崧凼且恢碌摹???闭口系统算例比较网格系统：?×?，???从图??可知，对方腔顶盖驱动流动，????蚐???算法的收敛范围最大，???和???惴ń喜睿坏?荢???蚐???算法出现Ⅱ松弛冈子奇异点?狭??直兔坐标外掠后台阶流动网格系统：????琑????方腔顼盖驱动流动咖图?苯亲?晖饴雍筇ń琢鞫?樟驳??问?冉?图?苯亲?晖饴雍筇ń琢鞫?樟彩奔浔冉?图?角欢ジ乔??鞫