高阶精度格式WCNS-E-5在复杂网格中的应用的研究.pdfVIP

第十四届全国计算流体力学会议论文集 高阶精度格式WONS-E-5在复杂 网格中的应用研究 邓小刚，毛枚良，涂国华，张毅锋 (中国空气动力研究与发展中心，四川绵阳621000) 摘 要：高阶格式常被用来计算流动复杂但是网格简单的问题，在解决复杂网格流动问题方面，国 内外都还处于探索阶段。本文总结了研究集体在采用高阶WCNS-E一5格式解决复杂网格问题的相关基础研 究方面的最新进展。主要包括WCNS—E-5的内点格式和边界格式，该格式对网格的总体要求，网格导数求 解方法，和分块特征对接技术等。通过对多种复杂网格问题的计算，表明WCNS-E一5格式已经初步具备了 计算复杂网格和复杂流动问题的能力。 关键词：高阶格式；复杂流动；对接网格；几何守恒律；特征边界 0 引 言 众所周知，在湍流的直接数值模拟、大涡模拟、计算声学和计算电磁学等需要精细计算的地方， 采用高阶格式是非常必要的。但是，由于高阶格式的数值耗散小，数值稳定性比低阶格式弱，对网 格质量的要求高，很少见有人在复杂网格体系下采用高阶格式。迄今为止，高阶格式的应用范围主 要集中在流动复杂但网格相对简单的问题上。 我们从上世纪90年代中期就开始一直致力于研究高阶精度格式的相关问题。首先，通过在中心 型紧致格式中加入耗散项，于1996年构造了单参数线性耗散紧致格式(DissipativeCompact 色散特性。同时，通过结合自适应非线性插值技术，构造了一类高阶紧致非线性格式(Compact Nonlinear 1997年，通过采用非线性加权技术，构造了一系列加权紧致非线性格式WCNS(WeightedCompact Nonlinear Scheme)∽1，包括隐式和显式两种类型。高阶精度的内点格式同时需要高阶精度的边界 格式才能达到全局的高阶精度，如P阶的内点格式需要(p-1)阶的边界格式，为此我们还建立了一套 适用于WCNS的边界格式，并且对边界格式做了渐进稳定性分析，结果表明在均匀网格或网格拉伸不 biasedfifth-order Pad6格式和EUW-5(ExplicitUpwind scheme)的耗散误差与色散误差特性， 分析了格式在分辨效率方面的特性，结果表明这种格式具有良好的耗散、色散特性。文献[7]还分析 了WCNs格式的多维相位误差的各向异性特性，研究表明WCNS格式在多维方向上的各向异性相位误 差主要体现在较窄的短波范围内。文献[8]通过GMRES方法来提高WCNS的收敛速度，并采用高阶格 式模拟了高超声速边界层对声波的感受性。 虽然WCNS格式具有较高的精度和优越的分辨率等性能，但是，与其他高阶格式类似，wCNs的 数值耗散小，在解决复杂外形复杂网格流动时的鲁棒性较低阶格式差。为了提高高阶格式对复杂外 形的适应能力，国内外主要有两种处理方法。一种方法是采用重叠网格(oversetgrid)，比如Sherer 区域内的物理信息。但是重叠网格技术在解决复杂外形上需要许多经验，而且在三维情况下的复杂 程度将大大增加，另外，插值技术也可能造成新的数值误差，且高阶插值易造成数值不稳定。另外 一种方法是采用对接网格或拼接网格，通过算术对接或特征特征对接等技术解决物理信息在不同网 17 第十四届全国计算流体力学会议论文集 格子块之间的传播。但是，文献所中给的与高阶格式相关的对接或拼接网格算例仍然比较简单，离 实际问题仍有不少差距。总之，虽然在高阶格式的应用方面取得了一定的进步，但是，目前国际上 采用高阶格式解决复杂网格复杂流动问题方面仍然还处于探索阶段。 文本围绕显式高阶WCNs格式(即WcNS—E一5)探讨解决复杂网格系统问题的方法。主要从高阶 边界格式入手，探讨高阶格式的网格策略，网格导数的求解方法，以及分块特征对接技术等。本文 分五个部分：第二部分介绍高阶精度的计算方法，包括无粘项和粘性项的计算格式和边界格式；第 三部分介绍了WcNs—E一5对网格的总体要求，网格守恒算法和特征对接技术；第四部分为一些典型算 例的计算结果；最后为总结。 1 计算方法 控制万栏为任葸坐杯糸卜册Navier-Stokes或