Research on pseudo arc length algorithm for strong discontinuity based on coordinate transformation

DOI:	10.11809/bqzbgcxb2023.04.011
Keywords:	explosion shock wave; numerical simulation; pseudo arc length algorithm; coordinate transformation; accuracy verification
Abstract:	Aiming at strong discontinuity singularity often occurring in the numerical simulation of explosion shock wave, this paper introduces arc length parameters in one or multi dimensional space. By adding a constraint equation, the original grid can move and deform adaptively to achieve distribution of a great number of grid points in strong discontinuous areas so as to improve the resolution of large gradient physical quantity capture. Concerning the difficulty in constructing a scheme in non uniform deformation physical space, a coordinate transformation is proposed to map the control equation in physical space into the uniform orthogonal arc length computational space. In addition, the corresponding one and two dimensional programs of the pseudo arc length algorithm are compiled based on MUSCL, and the method of manufactured solutions is developed to verify this algorithm. The numerical results of the pseudo arc length algorithm are compared with the finite volume method. Through error analysis, it is verified that the pseudo arc length algorithm has high resolution near the discontinuity.
Issue:	Vol. 44 No. 4 (2023)
Published:	2023-04-28
PDF HTML