一般坐标系下的二维euler方程时2空守恒格式.pdfVIP

第 16 卷　第 1 期 应 用 力 学 学 报 . 16　 . 1 V o l N o 1999 年 3 月 CH INESE JOURNAL OF APPL IED M ECHAN ICS M ar. 1999 一般坐标系下的二维Eu ler 方程时空守恒格式 张增产　沈孟育 (清华大学　北京　100084) 摘　　要 将文[ 13 ] 中的一维时空守恒格式推广到了二维情形, 得到了一般坐标系下的 二维 Eu ler 方程时空守恒格式, 并用几个典型算例进行了检验计算, 结果表明: 本文 得到的二维时空守恒格式保留了一维格式所有的优点, 格式简章, 通用性强, 而且对 激波等间断具有很高的分辨率。 关键词: Eu ler 方程; 时空守恒律; 数值方法; 激波分辩率 1　引　言 时空守恒元和解元方法( ) 是 . . 在文献[ 4 ] 中提出的一种新的数值方法。 CE SE S C Chang 该方法与传统的数值方法相比有许多不同之处, 例如, 它把时间与空间完全统一起来同等对 待, 并从守恒律积分型方程出发, 通过设立守恒元和解元, 使格式局部和全局都严格保证其物 理意义上的守恒律; 其次, 该方法把流场变量及其对空间偏导数都作为独立变量同时进行求 解; 另外, 它除利用了简单的台劳展开式外, 无需任何其它的数值逼近技术, 不需要任何单调性 限制或特征技术, 也不需要进行通量分裂, 因此该方法通用性好。但是, 该方法推广到高阶格式 ( ) 及多维情形 尤其是三维 时会遇到一些困难。为此, 文献[13 ]对其进行了改进。改变了格式的 构造方法, 由此得到的新格式既可保留原方法的优点, 同时格式又更简单, 通用性更强, 而且对 激波等间断具有很高的分辩率。本文将文[13 ] 中的一维时空守恒格式推广到了二维情形, 得 到了一般坐标系下的二维Eu ler 方程的时空守恒格式, 该格式具有一维格式全部优点。最后, 用几个二维典型算例进行了检验计算。 2　一般坐标系下的二维 Eu ler 方程时空守恒格式 完全气体二维非定常 Eu ler 方程组的无量纲化形式可写成如下: um t + f m x + g m y = 0, 　　m = 1, 2, 3, 4 ( 1) 国家自然科学基金资助项目 来稿日期: 修回日期: © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 第 1 期 一般坐标系下的二维 Euler 方程时空守恒格式 11 ( ) 对于物理区域不太规范的流动问题, 通常可以通过坐标变换 代数变换或数值变换 的方法把 ( ) 它变换到比较规则的区域 计算平面上 中来进行求解。现假设有如下坐标变换: = (x , y )