保持对称性的新积分梯度格式研究.pdfVIP

第四届全国计算爆炸力学会议论文集 * 保持对称性的新积分梯度格式 勇珩 袁国兴 （北京应用物理与计算数学研究所，北京 100088 ） 摘 要：在数值模拟内爆物理过程中，一维对称性是物理问题的基本要求。因此在二维柱坐标系下， 一个好的计算格式向一维问题退化时应该保证一维对称性。在现有积分梯度格式 IGA Integral （ Gradient Average IGT Integral Gradient Total IGTSP ）和 （ ）的基础上，提出了一种新的积分梯度格式 Integral Gradient Total Symmetry -Preserving IGTSP （ ）。理论上证明了：新积分梯度格式 不仅能保持 O h 一维球对称性，而且能较好地保持系统的总动量守恒，系统总动量守恒误差为 （）。数值试验进 一步显示了该格式的优点。 关键词：球对称性，动量守恒性，积分梯度法，柱坐标系 A new integrated gradient scheme for the preservation of symmetry YONG Heng YUAN guo-xing Institute of Applied Physics and Computational Mathematics, Beijing 100088,China （ ） Abstract: A new integrated gradient scheme IGTSP, which combines the advantages of the IGA and IGT schemes, is presented. The IGTSP scheme gives not only the spherical motion in a spherical problem with equal angular spacing in the cylindrical coordinate system, but also the reasonable fluid motion even if the initial cell mass are widely disparate, and the accuracy of the total momentum conservation is of order one. Numerical results demonstrate the advantages and disadvantages of the new scheme. Key word: spherical symmetry, momentum conservation, integrated gradients method, cylindrical coordinate system 1 引 言 研究二维柱坐标系下可压缩多介质流体力学问题，采用交错网格（速度定义在网格结点，内能、 质量定义在网格中心），可以准确确定物质界面位置。在任意拉格朗日-欧拉（Arbitrary Lagrangian- Eulerian ，简称ALE ）方法的Lagrange 计算[1，2]过程中，网格跟随流体一起运动，初始规