一类新的高阶混合型加权紧致非线性格式的构造.pdfVIP

第十三届现代数学和力学学术会议暨钱伟长诞辰100周年纪念大会文集 一类新的高阶混合型加权紧致非线性格式的构造 邓小刚 刘化勇 张涵信 1 中国空气动力研究与发展中心，空气动力学国家重点实验室， 四川绵阳（621000 ） 2 国防科学技术大学，湖南长沙（410073 ） 摘要 本文基于 WCNS 格式，提出了节点与半节点混合构造思想，发展了一类新的加权紧致 非线性节点与半节点混合格式（HWCNS ）。HWCNS 格式的分辨率（耗散和色散特性）与 WCNS 格式相当，但 HWCNS-E-5 格式的模板更加紧凑，仅有 7 点，更有利于实现高阶边 界信息的传递。边界和边界附近，节点与半节点混合格式在相对较窄的模板上也能够达到高 阶精度。本文选用了多个算例对HWCNS-E-5 格式进行了数值实验，结果表明 HWCNS-E-5 格式具有良好的间断捕捉和粘性分辨能力，格式分辨率及模拟精度与 WCNS-E-5 格式相当。 HWCNS-E-5 格式也能满足几何守恒率，很容易应用于复杂流动的高精度数值模拟。 1 引言 高阶精度格式作为CFD 研究的一个重要方向，在过去的三十多年中得到了广泛、深入的研 究，国内外已经发展了一系列高阶精度格式，包括：ENO、WENO、FCT、ENN、加权紧致非 线性格式（WCNS ）、间断有限元方法（DG ）、高阶迎风偏置和高阶紧致格式等。尽管高阶精度 格式在简单外形的复杂流动机理研究方面取得了巨大的成功，但在复杂外形中的应用还远不够 成熟，工程上的应用还比较少。 1988 年张涵信构造了无波动、无自由参数的耗散差分格式 （NND ）[3, 6]。上世纪90 年 代初，张涵信提出了构造高阶精度格式的四条基本原则 [5, 7]，在该构造原则的指导下，贺国 宏、张涵信构造了 ENN 格式 [8]。1997 年邓小刚、张涵信 [9, 1] 又提出了加权紧致非线性格 式（WCNS ），目前已被广泛应用于复杂外形、复杂流动的高精度数值模拟 [11]。WCNS 格式 有许多优点，比如在半节点插值可采用多种变量插值方式，通量构造既可以采用矢通量分裂方 法，也可采用通量差分分裂方法，非常灵活。Nonomura 指出 [2]，与WENO 格式 [13] 相比， WCNS 主要的优势在于在一般坐标系下WCNS 格式可以满足几何守恒率（GCL），而WENO 格式却很难满足。这一特性对高阶精度格式在复杂外形、复杂流动中的应用非常重要，因为在 复杂网格条件下，高网格质量很难保证，几何守恒误差会比较大，严重时会导致计算终止 [12]。 WCNS 格式的构造过程可分为三个步骤：（1）从节点到半节点的加权非线性变量插 值，（2 ）半节点的通量计算，（3 ）半节点到节点的中心型通量差分。第三个步骤可通过紧致和 显式格式来实现。由于显式格式有着计算量小，实现简单，易于并行等优点，文献 [1] 推荐使用 显式的加权非线性插值和显式的通量差分格式。WCNS-E-5 是一个显式的五阶格式，该格式尽 管在采用了7 点五阶加权非线性插值，然而格式的整体模板宽度达到了11 个点。对复杂外形的 工程应用问题，减少WCNS 格式的模板宽度非常重要，尤其是对于多块对接边界处理。 本文基于 WCNS 格式，提出了节点与半节点混合构造思想，发展了一类新的混合型加权 紧致非线性格式（HWCNS ）。HWCNS 格式的分辨率（耗散和色散特性）与WCNS 格式相当， 但 HWCNS-E-5 格式的模板更加紧凑。HWCNS-E-5 格式的模板宽度仅有 7 点，更有利于实现 1 第十三届现代数学和力学学术会议暨钱伟长诞辰100周年纪念大会文集 高阶边界信息的传递。此外，HWCNS-E-5 格式也能满足几何守恒率，很容易应用于复杂流动 的计算。在边界和边界附近，节点与半节点混合格式在较窄的模板上也能够达到高阶精度。 2 五阶 HWCNS 格式的构造 在这一节中，首先回