一种可压缩流动的高阶加权紧致非线性格式(WCNS)的加速收敛方法.pdfVIP

物 理 计算 .244.6四月 VoI . 2A .同。ι CHINESE JOURNAL OF COMPU时TIONAL PHYSl臼 笃IJ7年 1111 N~. 犁赔? [:It..~￥1 1∞1-246X(2007)06-伪98-07 一种可压缩流动的高阶加权紧致非线性格式 (WCNS) 的加速收敛方法 张毅锋， 邓小刚， 毛枚良， 陈坚强 [画面 要] 以高趋向速精位统流的数值模拟为例，研究ω-SCS、高斯- J睡德尔点松地、线始他以及 CMR岛电F .式求解方法在空闲项采用离阶#幅度格式 WCNS~国散时的收iUt，并对GMR臼(伊附世-1 IDÌTÙ mal reaidual) 方法 中的子迭代辱电作了对比计算.结果褒明.采用准确的第析ft亮比矩僻的杰、辈革松殊的收敛速度优于 LU-SCS， 以线松弛为预处理的 CMRES算Il具有良好的收敛将位. [关键姐 1 满阶精度格式; A是松也z 线松也; CM皿S 算法币收敛性 [中图分.-If) V21 1.3 [文献标识确1 A o 51 富 在过去的几十年中.为获得高阶精度的计算结果，发展了多种高阶格式，包括四O， WENO， B 佯条，FCf (通囊修正输运〉、离阶迎风偏慧和高阶紧敖格式等[，并成功用于各种流动计算.紧致格式以极少的网格点 数和较高精度受到众多学者的关注和研究.邓小刚等以单元中心裂紧致格式为基础，构造了一类基于单元边 界紧致非线位插值的高阶精度格式 CN$1231 ，之后又应用加权技术发展了 WCNS 格式[川1 该格式具有良好 的耗散、包徽特性.在光滑区与高阶迎风偏重显式线性格式的精度很接近，同时又具有很好的间断捕捉能力， 且变量插值独立于通量分裂格式，可以根据流动特点选用各种j!萤分裂格式，在对激波、涡波干扰、边界层等 典型流动的大量数值模拟[6.11 中，表现出良好的佳能，特别是在含有强激泼、大梯度变量的计算中有很强的优 势，如对于然流初摩阻均有准确的模拟，因而在高超声速绕流的数值计算中有很大的吸引力. 然而在应用过程中发现，WCNS 格式收敛很慢，尽管采用在二阶格式计算中具有很好稳定性和很高收敛 效率的 LU-民slll式方法，但依然存在收敛慢的困难，时间步长稍大就出现残差不降的现象，甚至计算失败， 而较小的时间步长往往使-个简单的绕流计算也需要几万步才能收敛.初步分析其主要原因是，在隐式方法 中，左端空间项的离散采用-阶格式，与右端项的高阶离散不匹配.另外在 LU-民S 方法中，左编项因采用後 量是大特征值分裂的近似雅克比矩阵而产生对角化的标tt矩阵，这虽然大大简化了计算，却引人较大的数值挺 徽，降低了左蝙项精度，而且近似因式分解法也引入了分解误楚，这些都会影响收敛性[1.) 因此，为了提高高 阶格式的收敛性，有必要研究适于篇阶格式的隐式求解方法. fl\~离散方程的求解方法主要有:直接解法、近似因式分解法、迭代法和Krylov 子空间法.尽管目前多采 用近似因式分解法和迭代法，但沁ylov 子空间法以其求解大型线性方程组的有效性吸引着众多学者的目光， 其中广义极小残差(GMR臼)算法，，-u) 的汁算量少、稳定性离.每步具有最优佐.是簸为有效可靠的一种方 法，近年来得到了广泛的研究和应用.各秒隐式方法的收敛位是不同的，而且离散方程左端项的精度也会影 响收敛性，据此，我们在左瑞项采用按Sceger-W由四ng 分裂的准确解析雅克比矩阵来提高摘自使.应用寓斯-赛 德尔点松弛、线松弛选代以及 GMR臼算法求解离徽方程，希望能够改善离阶精度格式WCNS 的收敛性能. 本文以二维、三维高超声速粘住绕流为例.空间项采用 WCNS-E习格式离徽，时间一阶后差隐式离散协S 方程，应用以上几种隐式方法求解.计算结果表明，采用