高精度迎风紧致差分格式二维全离散Fourier分析研究.pdfVIP

高精度迎风紧致差分格式二维全离散Fourier分析· 王 3910.2) 阵德薰 劝‘ 马延文(2) (t北京空气:3r力研究所.北京 100074) 2〔中国科学院力学研究所.LNM，北京 100080) 摘要 以Fourier分析方法导出了一类A李削优迎风紧致格式二维全离散色散关系式。分析了格式 的耗散、色散及各向异性等数值行为特征。折出这些特征主要体现于Fourier空间的高波数区域 内，且各向异性区域小于阳近梢反的对称紧致洛式。提高时间精度往往可减小这些误差。以可压 平面混合层算例检验了所研究洛式对线性扰动波及夏杂流动尺度较好的分辨能力. 关titi. Euler方程 迎风紧致格式 Fourier分析 非定常多尺度流动现象在实际应用中经常可以遇到。对这类问题的模拟不仅需要具有较高 时空精度的数值算法，还要求所用数值算法的耗散、色散等数值行为特性对所关心的流动现象 (特别是流场中的高波数小尺度物理量)能较好地反映。因此、计算中对数值算法的数值行为特性 的充分认识是必要的，并且这还有助于对所得数值结果进行正确地解释和说明。 近年，紧致格式因类谱性质而庄复杂流动直r}R值模拟中受到广泛关注。文I【3中、Lele通川 归纳总结，得到了一类七点对称型紧致格式、并给出了详细的空间半离散行为分析结果，随后 Y。等t=t对其中部分格式分别研究了它们在一维和二维情形下的全离散厅为特性。为了克服对称 格式无耗散所引起的数值混淆误差，又不损害紧致格式的类谱性质，Fu和Ma3[)基于矢通量分裂 思想先后提出了一类高精度迎风型紧致洛式.井研究了其在空间半离散情况下的数值特性。最近、 作者在文h[4]中又进一步对该类格式进行了一维全离散数值行为分析。 为了认识文 3〔3高精度迎风紧致筹分洛式公抓份清形下的数值特性，本文给出一个相应的 维全离散行为分析结果。分析表明.泣维隋况下，所研究格式构成的数值色散介质中有三组性质 不同的数值波存在.并在高波数区均受到明显的耗散和色散作用，其中以传播方向同流动方向 相反的声波最为明显。格式的各向异性误差也主要集中于高波数区，与阶数相当的对称紧致格SI 相比，该误差可忽略区域更大。最后以二维Igliil和空间发展可压平面混合层算例显示了高精度迎 风紧致格式对线性扰动波及流动复杂尺度的分辨能力。 数值方法 Cartesian坐标系下，二维Euler方程非守恒型特征形式为 阅 .*aQ.。阅 一 宁 气」 — 宁 厄J - 二 ()) a ax 即 式中Q是基本流动元未知向量.A和B是,f:粘通if!系数矩阵。为了使方程(I)线化，以适于线性分 析，A和B肢设作常矩阵，同时，假L,C它们的特征rt均非负 对空间一阶导数，文[3]的一类迎 。国家自然科学基金资助课题(编号:NNSrc-19391)10) 298 风紧致差分格式中最常用的两个3阶和5阶格式为 (对特征值非负情形) UCD3 3 二 ， 8}u+上S;u, 2 ( a ) 合··-; 6h 3 口 蜘 UCD5 5 一 普··-;60h(-6夏u+:+118,u,,,+476,u,+35,u,-,) 其中u为任一基本流动元分量，u为‘该分6r 阶睁数，上标+，一分别表示一阶前差和后差。本又 主要对此两种格式加以分析。方程(I)的m级显式Runge-Kutta时间推进全离散形式为 Ql=Q+4t芝a,sR(Q-) 倒=L二，m) (3) 这里s对应n至n+l时间层内顺序中间步.ct，是确定时间精度的对流非齐次项加权平均系数。 2Fourier分析 2.1全离散Fourier变换 设数值解为无限区