Numerical dissipation control of a hybrid large-particle method in vortex instability problems

Current trends in the development of numerical schemes are associated with a decrease in dissipative and dispersion errors, as well as an improvement in the grid convergence of the solution. Achieving the computational properties is not an easy problem since a decrease in scheme viscosity is often associated with an increase in the oscillations of gas dynamic parameters. The paper presents a study of the issues of numerical dissipation control in gas dynamics problems in order to increase the resolution in numerical reproduction of vortex instability at contact boundaries. To solve this problem, a hybrid large-particle method of the second order of approximation in space and time on smooth solutions is used. The method is constructed with splitting by physical processes in two stages: gradient acceleration and deformation of the finite volume of the medium; convective transfer of the medium through its facets. An increase in the order of approximation in time is achieved by a time correction step. The regularization of the numerical solution of problems at the first stage of the method consists in the nonlinear correction of artificial viscosity which, regardless of the grid resolution, tends to zero in the areas of smoothness of the solution. At the stage of convective transport, the reconstruction of fluxes was carried out by an additive combination of central and upwind approximations. A mechanism for regulating the numerical dissipation of the method based on a new parametric limiter of artificial viscosity is proposed. The optimal adjustment of the method by the ratio of dissipative and dispersive properties of the numerical solution is achieved by setting the parameter of the limiting function. The efficiency of the method was tested on two-dimensional demonstrative problems. In one of them, the contact surfaces are twisted into a spiral on which the Kelvin-Helmholtz vortex instability develops. Another task is the classic problem with a double Mach reflection of a strong shock wave. Comparison with modern numerical schemes has shown that the proposed variant of the hybrid large-particle method has a high competitiveness. For example, in the problem with double Mach reflection, the considered version of the method surpasses in terms of vortex resolution the popular WENO (Weighted Essentially Non-Oscillatory) scheme of the fifth order and is comparable to the numerical solution of WENO of the ninth order of approximation. The proposed method can be the basis of a convective block of a numerical scheme when constructing a computational technology for modeling turbulence.

1. Godunov S.K. A finite difference method for the computation of discontinuous solutions of the equations of fluid dynamics. Mathematics of the USSR – Sbornik, 1959, vol. 47, no. 8-9, pp. 357–393.

2. Roe P.L. Approximate Riemann solvers, parameter vectors and difference schemes. Journal of Computational Physics, 1981, vol. 43, no. 2, pp. 357–372. https://doi.org/10.1016/0021-9991(81)90128-5

3. Belotserkovskii O.M., Davydov Yu.M. A non-stationary “Coarse particle” method for gas-dynamical computations. USSR Computational Mathematics and Mathematical Physics, 1971, vol. 11, no. 1, pp. 241–271. https://doi.org/10.1016/0041-5553(71)90112-1

4. Kolgan V.P. An application of the minimal derivative values concept to construction of the finite-difference schemes for calculating the discontinuous solutions of gas dynamics. Uchenyye Zapiski TsAGI, 1972, vol. 3, no. 6, pp. 68–77. (in Russian)

5. Harten A. High resolution schemes for hyperbolic conservation laws. Journal of Computational Physics, 1983, vol. 49, no. 3, pp. 357–393. https://doi.org/10.1016/0021-9991(83)90136-5

6. Liu X.-D., Osher S., Chan T. Weighted essentially non-oscillatory schemes. Journal of Computational Physics, 1994, vol. 115, no. 1, pp. 200–212. https://doi.org/10.1006/JCPH.1994.1187

7. Jiang G.-S., Shu C.-W. Efficient implementation of weighted ENO schemes. Journal of Computational Physics, 1996, vol. 126, no. 1, pp. 202–228. https://doi.org/10.1006/jcph.1996.0130

8. Zhang S., Zhu J., Shu C.-W. A brief review on the convergence to steady state solutions of Euler equations with high-order WENO schemes. Advances in Aerodynamics, 2019, vol. 1, no. 1, pp. 16. https://doi.org/10.1186/s42774-019-0019-2

9. Minoshima T., Miyoshi T. A low-dissipation HLLD approximate Riemann solver for a very wide range of Mach numbers. Journal of Computational Physics, 2021, vol. 446, pp. 110639. https://doi.org/10.1016/j.jcp.2021.110639

10. Li J., Shu C.-W., Qiu J. Multi-resolution HWENO schemes for hyperbolic conservation laws. Journal of Computational Physics, 2021, vol. 446, pp. 110653. https://doi.org/10.1016/j.jcp.2021.110653

11. Shi J., Zhang Y.-T., Shu C.-W. Resolution of high order WENO schemes for complicated flow structures. Journal of Computational Physics, 2003, vol. 186, no. 2, pp. 690–696. https://doi.org/10.1016/S0021-9991(03)00094-9

12. Deng X., Xie B., Loubère R., Shimizu Y., Xiao F. Limiter-free discontinuity-capturing scheme for compressible gas dynamics with reactive fronts. Computers & Fluids, 2018, vol. 171, pp. 1–14. https://doi.org/10.1016/j.compfluid.2018.05.015

13. Sadin D.V. A modification of the large-particle method to a scheme having the second order of accuracy in space and time for shockwave flows in a gas suspension. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming & Computer Software, 2019, vol. 12, no. 2, pp. 112–122. (in Russian). https://doi.org/10.14529/mmp190209

14. Sadin D.V., Golikov I.O., Shirokova E.N. Testing of the hybrid large-particle method using two-dimensional Riemann problems. St. Petersburg State Polytechnical University Journal: Physics and Mathematics, 2021, vol. 14, no. 1, pp. 58–71. (in Russian). https://doi.org/10.18721/JPM.14104

15. Zhang D., Jiang C., Liang D., Cheng L. A review on TVD schemes and a refined flux-limiter for steady-state calculations. Journal of Computational Physics, 2015, vol. 302, pp. 114–154. https://doi.org/10.1016/j.jcp.2015.08.042

16. Sadin D.V. Simulation of physical instability on contact boundaries in multicomponent compressible gas flows by a hybrid large-particle method. Numerical Methods and Programming, 2020, vol. 21, no. 2, pp. 129–137. (in Russian). https://doi.org/10.26089/NumMet.v21r211

17. Semenov A.N., Berezkina M.K., Krasovskaya I.V. Classification of shock wave reflections from a wedge. Part 2: Experimental and numerical simulations of different types of Mach reflections. Technical Physics, 2009, vol. 54, no. 4, pp. 497–503. https://doi.org/10.1134/S1063784209040094

18. Sadin D.V. Analysis of dissipative properties of a hybrid large-particle method for structurally complicated gas flows. Computer Research and Modeling, 2020, vol. 12, no. 4, pp. 757–772. (in Russian). https://doi.org/10.20537/2076-7633-2020-12-4-757-772

19. Sadin D.V. Efficient implementation of the hybrid large-particle method. Matematicheskoe Modelirovanie, 2022, vol. 34, no. 4, pp. 113–127. (in Russian). https://doi.org/10.20948/mm-2022-04-08