五点隐式TVD格式地研究探索.pdfVIP

中国工程热物理学会 流体机械 2004年学术会议 编号：047052 五点隐式TVD格式的研究探索 陈宏冀 中国科学院工程热物理研究所，北京100080 电话：(010 E-mail：chj@mail．etp．ac．cn 摘要：本文探索研究在构造时间空间均为高精度的隐式TVD格式中采用较多的空间网格 基点，以减小差分格式高阶的色散和耗散。构造了一个五空间点隐式TVD基本格式，其 差分系数用拉格朗日插值多项式得出，利用限量机制和作者发展的推导论证方法，得出 限量因子的条件，使之在任意时间步长下是TVD的。实例计算验证了理论，显示了效果。 一、引 言 我们知道对于数值求解流动的定常问题时其收敛解的精度只取决于差分格式的空 间精度，而对非定常问题差分格式的时间精度和空间精度同样重要，另外人们为了不要 受到稳定性对于时间步长的过分限制而常常采用隐式格式。对于日渐广泛使用的高分辨 率TVD格式也是这样的。本作者在对隐式TVD格式的研究探索中曾提出了两个格式“11， 分析表明都是时间空间均为二阶精度的，且TVD性的保持不受时间步长限制，这超过了 过去己有的格式。它们仍是本质属于空间三点格式的．一般构造差分格式时使用较多的 网格基点可实现较高阶精度，而使用较高阶精度的格式可减少使用的网格数。也可见到 在TVD格式中使用较多网格基点的情形。然而鉴于TVD格式的机理，构造高于二阶精度 的格式是有困难的，但利用较多的网格基点仍可以达到减小高阶的色散和耗散。本文对 五空间点的隐式TVD格式进行探索研究，构造基本格式，采用一种简单的方法确定变量 的差分系数，利用限量因子机制和作者发展的推导论证方法，得出限量因子的条件，使 之在任意时间步长情况下保持TVD性。用计算实例对理论进行验证，观察效果。 二、五点隐格式 对于双曲型守恒律方程 皇竺+af(u__A：0或皇竺+口旦竺：o (1) at 敏 al 瓠 散项的大小正好消掉Euler格式的修正方程中的二阶导数项。而本文拟构造空间五点格 式、欲消掉更高阶导数项，用这种方法将很麻烦，因此本文采用另外一种可以简单地得 出差分方程各项的系数的方法。先看两个显式格式的情况，图1表示通过(j，n+1)节 用线性插值将“(P)以“■和“?表示，则得到一阶精度的迎风格式；若“(P)利用抛物 本工作得到国家重点基础研究项目(NoG1999022309)的资助 133 们将此具有普遍性的方法用于隐式五点格式。 0-1，11)P a，n) O+1，n) O,n) 图(1)显式格式 图(2)五点隐式格式 就得到差分方程 此差分方程的截断误差由拉格朗日插值多项式余项决定，为pb4)，即格式为4阶精 度。将其重新排列得到 ‘ ～ (3) ， 。 一三署k：一心。+鸭一q。+％：)I，“=彬 可以看出，方程相当于隐式一阶迎风格式2nT--个反各高阶误差的项。 三、构造TVD格式 我们不对上述格式进行一般的稳定性分析，而直接利用它构造具有TVD稳定性的格 式。现将式(3)进一步改写成守恒的形式 叶n“+协川∥“+半(△j+112u--△j-112uP“ 一兰半[(△j+3／2“--△j-l／2“)一¨∥一△j-312U)1I，“ (4) 这里使用了符号△，+1，2“=“J+1一“川以便于下面的步骤中限量机制的运用·各反高阶