High-Order Pseudo Arc-Length Method for Strongly Discontinuity of Multi-Medium Flows

被引：0
|
作者：
Wang C. [1 ]
Ma T. [1 ]
Li K. [1 ]
机构：
[1] State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing
来源：
Beijing Ligong Daxue Xuebao/Transaction of Beijing Institute of Technology | 2022年 / 42卷 / 09期
关键词：
coordinate transformation; high-order; pseudo arc-length method (PALM); real ghost fluid method (RGFM); weighted essentially non-oscillatory (WENO);
D O I：
10.15918/j.tbit1001-0645.2021.205
中图分类号：
学科分类号：
摘要：
Considering the strong discontinuities appeared in the solutions of hyperbolic conservation law equations with the development of time and the ability of pseudo arc-length method (PALM) in weakening the singularity of the equation, a high-order PALM was developed based on the combination of the PALM and the high-order weighted essentially non-oscillatory (WENO) scheme to effectively improve the accuracy and resolution of the calculation. To overcome the complexity of a high-order scheme construction in the deformed mesh, a coordinate transformation method was adopted to make the governing equation mapped to the orthogonal uniform arc-length space, completing the calculation in the arc-length space. Combined with Level Set technology and real ghost fluid method (RGFM) interface processing, the algorithm was extended to the calculation of multi-medium flows. For the interpolation of Level Set distance function after grid movement, a third-order non-conservative interpolation scheme was proposed. The results show that the high-order PALM has high convergence order, which can effectively reduce the numerical oscillation at the discontinuity and improve the discontinuity resolution. © 2022 Beijing Institute of Technology. All rights reserved.
引用
收藏
页码：900 / 908
页数：8
相关论文
共 18 条
  • [1] JIANG G S, SHU C W., Efficient implementation of weighted ENO schemes[J], Journal of Computational Physics, 126, 1, pp. 202-228, (1996)
  • [2] LL Peng, WANG Cheng, Improved WENO schemes for one-dimensional detonation simulations[J], Transactions of Beijing Institute of Technology, 37, 12, pp. 1211-1216, (2017)
  • [3] RIKS E., An incremental approach to the solution of snapping and buckling problems[J], International Journal of Solid Structures, 15, 7, pp. 529-551, (1979)
  • [4] CHAN T F., Newton-like pseudo-arclength methods for computing simple turning points[J], Siam Journal on Scientific & Statistical Computing, 5, 1, pp. 135-148, (1984)
  • [5] WU Jike, XU Weihou, DING Hongli, Arc-length method for differential equations[J], Applied Mathematics and Mechanics, 8, pp. 108-113, (1999)
  • [6] TANG H Z, TANG T., Adaptive mesh methods for one- and two-dimensional hyperbolic conservation laws[J], Siam Journal on Numerical Analysis, 41, 2, pp. 487-515, (2003)
  • [7] YUAN X P, NING J G, MA T B, Et al., Stability of Newton TVD Runge-Kutta scheme for one-dimensional Euler equations with adaptive mesh[J], Applied Mathematics & Computation, 282, (2016)
  • [8] LIN Xiankun, LI Jian, RONG Jili, Et al., Application of fifth-order WENO finite difference method in the moving-interface tracking[J], Transactions of Beijing Institute of Technology, 35, 4, pp. 341-346, (2015)
  • [9] ZHANG Q H, FOLGELSON A., Fourth-order interface tracking in two dimensions via an improved polygonal area mapping method[J], SIAM Journal on Scientific Computing, 36, 5, pp. A2369-A2400, (2014)
  • [10] OSHER S, FEDKIW R P., Level set methods[J], Geometric Level Set Methods in Imaging Vision & Graphics, 32, 11, pp. 3-20, (2003)
12