间断Galerkin法对前台阶绕流问题的研究.pdfVIP

LHD 2012 年度夏季学术研讨会 间断 Galerkin 法对前台阶绕流问题的研究 赵文赓，郑洪伟 中国科学院力学研究所，高温气体动力学国家重点实验室（筹），北京海淀区 100190 摘要 相对于有限差分法，间断 Galerkin 方法可以很容易的运用于非结构网格，并且在边界上通量的 计算只需要相邻两个单元的信息，因此在处理复杂的边界问题时具有很大优势。本文采用基于三角 形非结构网格的间断 Galerkin 方法计算前台阶绕流问题，用分片线性的多项式近似精确解，采用龙 格库塔法进行时间推进。从计算结果可以清晰地观察到流场建立的过程。 关键词 前台阶绕流，间断 Galerkin，龙格库塔法，非结构网格 U f g 引 言   0 t x y （1） 间断有限元方法的出现，最早可以追溯到 上式的矢量形式为： 1973 年Reed 和Hill 关于中子输运方程问题的论 U F 0 t 文。之后，在解决含有间断现象的问题中发挥 ， F f , g 着越来越大的作用，它被广泛地应用到了水动 其中 U 是状态变量，  ，具体为，       力学，气动力学，波传播等问题。    2  u u p DG 方法在很多方面有突出的优点，使得该 U   f       方法引起很多学者的关注。首先，DG 中，精度 v uv     E (  ) 的提高是通过适当选取基函数，即提高单元插   u E p  ， ， 值多项式的次数来实现，这克服了有限体积法    （   FVM ）中通过扩大基架点数目（比如k-exact,  uv  ENO, WENO ）来提高精度的不足；其次，由于 g v2 p  近似解的间断性假设，这