Moving Meshes of One Dimensional Slab Symmetry Shock Waves
-
摘要: 移动网格方法可以有效提高冲击波阵面计算的分辨率,多年来一直受到数值计算研究领域的关注。针对冲击波在凝聚介质中传播的数值模拟,利用变分原理构建了基于移动网格的自适应算法。移动网格的生成依赖包含控制函数的欧拉方程的迭代计算,通过控制方程从物理域到计算域的映射来求解物理量,并研究了不同移动网格迭代方法对计算效率的影响。数值结果显示了该算法的有效性。Abstract: Recently, the moving mesh method has received plenty of attention in the numerical computation areas owing to its capability to improve effectively the calculation precision on the shock waveplane.This paper describes a moving mesh method that simulates the shock wave's propagation in condensed matter based on the variation principle.The generation of the moving meshes consists of the iterative computation of the Euler-Lagrange equation and the governing equation's mapping from the physical plane to the computational plane.We have studied the effect of computational efficiency using different iterative methods.Finally, the numerical results show the validity of the arithmetic.
-
Key words:
- moving mesh /
- iterative methods /
- numerical computation /
- governing functions
-
表 1 两种迭代方法的比较
Table 1. Comparison of two iterative methods
Iterative method D t/(ms) Brackbill 2 667 109.375 Jacobi 257 15.625 -
[1] TANG H, TANG T.Adaptive mesh methods for one- and two-dimensional hyperbolic conservation laws [J].Siam J Numer Anal, 2003, 41(2):487-515. doi: 10.1137/S003614290138437X [2] 安德森·J D.计算流体力学基础及其应用[M].吴颂平, 刘赵森, 译.北京: 机械工业出版社, 2007: 139-148.ANDERSON J D.Computational fluid mechanics [M].Translated by WU S P, LIU Z S.Beijing: China Machine PRESS, 2007: 139-148. [3] SEMPLICE M, COCO A, RUSSO G.Adaptive mesh refinement for hyperbolic systems based on third-order compact WENO reconstruction[J].J Sci Comput, 2016, 66(2):692-724. doi: 10.1007/s10915-015-0038-z [4] VALORANI M, GIACINTO M D.Performance assessment of an adaptive mesh refinement technique for detonation waves[C]//NAPOLITANO M, SABETTA F.Thirteenth International Conference on Numerical Methods in Fluid Dynamics.Berlin: Springer-Verlag Berlin Heidelberg, 1993: 300-304. [5] FRYXELL B, OLSON K, RICKER P, et al.FLASH:an adaptive mesh hydrodynamics code for modeling astrophysical thermonuclear flashes[J].Astrophys J Suppl, 2000, 131(1):273-334. doi: 10.1086/apjs.2000.131.issue-1 [6] 孙承纬.一维冲击波和爆轰波计算程序SSS[J].计算物理, 1986(2):18-30. http://www.cnki.com.cn/Article/CJFDTOTAL-JSWL198602001.htmSUN C W.Computational program SSS of one dimensional shock and detonation waves[J].Chinese Journal of Computational Physics, 1986(2):18-30. http://www.cnki.com.cn/Article/CJFDTOTAL-JSWL198602001.htm [7] 忻孝康, 朱士灿, 张慧生.一维Burgers方程的各种差分格式研究[J].上海力学, 1980(1):65-84. http://cdmd.cnki.com.cn/Article/CDMD-10459-1017129064.htmXIN X K, ZHU S C, ZHANG H S.A study of various finite difference schemes for Burgers' equation [J].Shanghai Mechanics, 1980(1):65-84. http://cdmd.cnki.com.cn/Article/CDMD-10459-1017129064.htm [8] WINSLOW A M.Numerical solution of the quasilinear Poisson equation [J].J Comput Phys, 1966, 1(2):149-172. doi: 10.1016/0021-9991(66)90001-5 [9] BRACKBILL J U, SALTZMAN J S.Adaptive zoning for singular problems in two dimensions[J].J Comput Phys, 1982, 46(3):342-368. doi: 10.1016/0021-9991(82)90020-1 [10] 陈国贤, 汤华中, 张平文.移动网格上的二阶Godunov型格式在一维爆轰波模拟中的应用[C]//刘凯欣.计算爆炸力学进展.2006: 73-79.CHEN G X, TANG H Z, ZHANG P W.The simulation of one dimensional detonation waves based on second-order Godunov-type scheme on moving meshes[C]//LIU K X.Advances in Computational Explosion Mechanics.2006: 73-79.