曲型守恒律的一个新的通量限制器.pdfVIP

ｄｏｉ：１０３９６９／ｊｉｓｓｎ１００１９２３５２０１４０１００５　　　　　　　　　　　　 　 　 　２０１４年第１期 ·ＰＥＡＲＬＲＩＶＥＲ　人民珠江 双曲型守恒律的一个新的通量限制器 １ ２ ３ 朱 嵩 ，万五一 ，毛根海 （１．珠江水利科学研究院，广东 广州 ５１０６１１；２．浙江大学水利与海洋工程学系，浙江 杭州 ３１００２７； ３．浙江大学玉泉校区水利实验室，浙江 杭州 ３１００２７） 摘 要：对于双曲型守恒律，当间断或激波存在时，线性的二阶差分格式往往会导致非物理的振荡而一阶迎风格式 耗散性则过大。为了抑制虚假振荡使得数值解锐利而逼真，离散格式往往需要引入非线性的限制器。限制器的性 能对计算结果仍有着较大的影响，构造了一个新的限制器Ｍ７，它是一种耗散型限制器，在满足Ｃ．Ｆ．Ｌ．条件时具有 ＴＶＤ的性质。数值计算结果表明，该限制器能够有效抑制虚假的振荡，减轻数值群速度效应，同时该限制器和一些 传统限制器相比提高了显格式计算的稳定性条件。 关键词：通量限制器；双曲型守恒律；Ｃ．Ｆ．Ｌ．条件；ＴＶＤ；有限体积法 中图分类号：ＴＶ１３１．６５　　文献标识码：Ａ　　文章编号：１００１９２３５（２０１４）０１００１９０５ ０　引言 限制器在一定的条件下，能够取得较优的计算结果。 双曲型偏微分方程系统可以用来描述波的运动或物质 １　数学模型 的对流输运等一大类现象，各向同性双曲型方程称为守恒 １１　模型方程 ［１］ 律 。虽然形式简单，但完全可以用它来检验数值方法的效 一维双曲型守恒律方程可以表示如下： 果和功能，因而在发展离散方法、数值格式方面具有广泛的 Ｕ＋Ｆ（Ｕ） ＝０ （１） ｔ ｘ 应用。可压缩流体动力学的Ｅｕｌｅｒ方程和无粘 Ｂｕｒｇｅｒｓ方程 式中　Ｕ———守恒变量；Ｆ———相应的通量函数。 都是它的特殊形式，其他类型的方程如描述浅水运动的Ｓａｉｎｔ 假设网格均匀，定义单元ｉ为 则在单元ｉ （ｘ ，ｘ ）， ｉ－１／２ ｉ＋１／２ ［２］ －Ｖｅｎａｎｔ方程也是该方程的派生形式 。 上对方程进行积分可得： 现代计算流体力学为了满足计算精度的要求，往往需要 ｎ＋１ ｎ Δｔ ２） Ｕ ＝Ｕ － （Ｆ －Ｆ ） （ ｉ ｉ ｉ＋１／２ ｉ－１／２ 其离散格式具有二阶或二阶以上的计算精度。但一般线性 Δｘ 的二阶精度的离散格式在间断附近产生非物理的振荡，理论 式中　 、 ———为时间步长和空间步长； ———变量