|
Article information
2018 , Volume 23, ¹ 2, p.37-54
Kovyrkina O.A., Ostapenko V.V.
On monotonicity and accuracy of CABARET scheme for calculation of weak solutions with shocks
We studied the monotonicity and the accuracy of the modified CABARET scheme approximating the quasilinear hyperbolic system of conservation laws. We obtained conditions under which this scheme preserves the monotonicity of the finite-difference solution with respect to invariants of the linear approximation of the considering system. As a specific example, we considered the approximation for the system of conservation laws for shallow water theory. To determine the accuracy of the scheme in the regions influenced by shocks that propagate with variable velocity and behind the fronts of which a non-constant solution is formed, series of test calculations on the sequence of contracting grids were carried out. These calculations enabled us to apply the Runge rule to define the order of convergence of the finite-difference solution. To estimate the accuracy of the modified CABARET scheme, the orders of its integral and local convergence are calculated, as well as local disbalances in calculating of the absolute value of averaged basis functions vector and local disbalances in calculating of the absolute value of the components of the averaged vector of invariants. The order of the integral convergence makes it possible to estimate the accuracy of finite-difference schemes in translation of the RankineHugoniot conditions across a non-stationary shock wave front. In this case, the order of integral convergence is calculated for the finite-difference solution itself (rather than for its absolute value, as in the L1 norm) on spatial intervals crossing the shock wave. It is shown that, like TVD-schemes of high order approximation on smooth solutions, the CABARET scheme being formally of the second order, in spite of high accuracy in the localization of shocks, reduces the order of integral convergence to the first order on the integration intervals, one of the boundaries of which is in the region influenced by a shock. The main reason of this decrease in accuracy is that the flux correction used in the CABARET scheme to its monotonization leads to the decrease in the smoothness of the finite-difference fluxes, which in turn leads to the decrease in the order for approximation of Rankine-Hugoniot conditions on the shocks. As a result, the local accuracy of the CABARET scheme in the regions of shocks influence also decreases to about the first order.
[full text]
Keywords: hyperbolic system of conservation laws, monotonicity and accuracy of CABARET scheme, shallow water flows, shocks
doi: 10.25743/ICT.2018.23.12757
Author(s):
Kovyrkina Olyana Aleksandrovna
PhD.
Position: Senior Research Scientist
Office: Lavrentyev Institute of Hydrodynamics SB RAS
Address: 630090, Russia, Novosibirsk, 15 Lavrentyev pr.
Phone Office: (383)333-22-01
E-mail: olyana@ngs.ru
Ostapenko Vladimir Viktorovich
Dr. , Professor
Position: General Scientist
Office: Lavrentyev Institute of Hydrodynamics of the Siberian Branch of the RAS
Address: 630090, Russia, Novosibirsk, Ac. Lavrentyev Ave., 15
Phone Office: (383)333-22-01
E-mail: Ostapenko_vv@ngs.ru
SPIN-code: 1676-5882 References:
[1] Iserles, A. Generalized leapfrog methods. IMA Journal of Numerical Analysis. 1986; 6(3):381–392.
[2] Goloviznin, V.M., Samarskii, A.A. Finite difference approximation of convective transport equation with space splitting time derivative. Matematicheskoe Modelirovanie. 1998; 10(1):86–100. (In Russ.)
[3] Goloviznin, V.M., Samarskii, A.A. Some characteristics of finite difference scheme “cabaret”. Matematicheskoe Modelirovanie. 1998; 10(1):101–116. (In Russ.)
[4] Rozhdestvenskiy, B.L., Yanenko, N.N. Sistemy kvazilineynykh uravneniy [Systems of Quasilinear Equations and Their Applications to Gas Dynamics]. Moscow: Nauka; 1978: 687. (In Russ.)
[5] Kulikovskiy, A.G., Pogorelov, N.V., Semenov, F.Yu. Matematicheskie voprosy chislennogo resheniya giperbolicheskikh sistem uravneniy [Mathematical Aspects of Numerical Solution of Hyperbolic Systems]. Moscow: Fizmatlit; 2001: 607. ( In Russ.)
[6] Goloviznin, V.M. Balanced characteristic method for systems of hyperbolic conservation laws. Doklady Mathematics. 2005; 72(1):619–623.
[7] Woodward, P., Colella, P. The numerical simulation of two-dimensional fluid flow with strong shocks. Journal of Computational Physics. 1984; 54(1):115–173.
[8] Karabasov, S.A., Goloviznin, V.M. New efficient high-resolution method for nonlinear problems in aeroacoustics. American Institute of Aeronautics and Astronautics journal. 2007; 45(12):2861–2871.
[9] Semiletov, V., Karabasov, S.A. CABARET scheme with conservation-flux asynchronous time-stepping for nonlinear aeroacoustics problems. Journal of Computational Physics. 2013; (253):157–165.
[10] Faranosov, G., Goloviznin, V., Karabasov, S., Kondakov, V., Kopiev, V., Zaitsev, M. CABARET method on unstructured hexahedral grids for jet noise computation. Comput. Fluids. 2013. Vol. 88. P. 165–179.
[11] Karabasov, S.A., Berloff, P.S., Goloviznin, V.M. Cabaret in the ocean gyres. Ocean Modelling. 2009; 30(2):155–168.
[12] Goloviznin, V.M., Zaytsev, M.A., Korotkin, I.A. Novye algoritmy vychislitel'noy gidrodinamiki dlya mnogoprotsessornykh vychislitel'nykh kompleksov [New CFD Algorithms for Multiprocessor Computer Systems]. Moscow: Izd-vo Moskovskogo universiteta; 2013: 467. (In Russ.)
[13] Kovyrkina, O.A., Ostapenko, V.V. On monotony of two layer in time cabaret scheme. Mathematical Models and Computer Simulations. 2013; 5(2):180–189.
[14] Kovyrkina, O.A., Ostapenko, V.V. On the monotonicity of the CABARET scheme in the multidimensional case. Doklady Mathematics. 2015; 91(3):323–328.
[15] Zuzina, N.A., Ostapenko, V.V. On the monotonicity of the CABARET scheme approximating a scalar conservation law with a convex flux. Doklady Mathematics. 2016; 93(1):69–73.
[16] Ostapenko, V.V. Convergence of finite-difference schemes behind a shock front. Computational Mathematics and Mathematical Physics. 1997; 37(10):1161–1172.
[17] Casper, J., Carpenter, M.H. Computational consideration for the simulation of shockinduced sound. SIAM Journal on Scientific Computing. 1998; 19(1):813–828.
[18] Engquist, B., Sjogreen, B. The convergence rate of finite difference schemes in the presence of shocks. SIAM Journal on Numerical Analysis. 1998; (35):2464–2485.
[19] Godunov, S.K. A difference method for numerical calculation of discontinuous solutions of the equations of hydrodynamics. Matematicheskii Sbornik. 1959; 47(3):271–306.
[20] Van Leer, B. Toward the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method. Journal of Computational Physics. 1979; 32(1):101–136.
[21] Harten, A. High resolution schemes for hyperbolic conservation laws. Journal of Computational Physics. 1983; (49):357–393.
[22] Nessyahu, H., Tadmor, E. Non-oscillatory Central Differencing for Hyperbolic Conservation Laws. Journal of Computational Physics. 1990; 87(2):408–463.
[23] Jiang, G.S., Shu, C.W. Efficient implementation of weighted ENO schemes. Journal of Computational Physics. 1996; (126):202–228.
[24] Kovyrkina, O.A., Ostapenko, V.V. On the convergence of shock-capturing difference schemes. Doklady Mathematics. 2010; 82(1):599–603.
[25] Kovyrkina, O.A., Ostapenko, V.V. On the practical accuracy of shock-capturing schemes. Mathematical Models and Computer Simulations. 2014; 6(2):183–191.
[26] Mikhailov, N.A. On convergence rate of WENO schemes behind a shock front. Mathematical Models and Computer Simulations. 2015; 7(5):467–474.
[27] Samarskiy, A.A. Teoriya raznostnykh skhem [The theory of difference schemes]. Moscow: Nauka; 1983: 616. (In Russ.)
[28] Ostapenko, V.V. Finite-difference approximation of the Hugoniot conditions on a shock front propagating with variable velocity. Computational Mathematics and Mathematical Physics. 1998; 38(8):1299–1311.
[29] Ostapenko, V.V. Construction of high-order accurate shock-capturing finite difference schemes for unsteady shock waves. Computational Mathematics and Mathematical Physics. 2000; 40(12):1784–1800.
[30] Stoker, J.J. Water Waves: The Mathematical Theory with Applications. New York: Wiley; 1957: 600.
[31] Karabasov, S.A. On the powerof second-order accurate numerical methods for model problems of gas- and hydrodynamics. Mathematical Models and Computer Simulations. 2011; 3(1):92–112.
[32] Zyuzina, N.A., Ostapenko, V.V. Modification of the CABARET scheme ensuring its strong monotonicity and high accuracy on local extrema. Doklady Mathematics. 2014; 90(1):453–457.
[33] Zuzina, N.A., Ostapenko, V.V. Modification of the Cabaret scheme ensuring its high accuracy on local extrema. Mathematical Models and Computer Simulations. 2016; 8(3):231–237.
[34] Lax, P., Wendroff, B. Systems of conservation laws. Communications on Pure and Applied Mathematics. 1960; (13):217–237.
[35] Rusanov, V.V. Difference schemes of third order accuracy for the through calculation of discontinuous solutions. Doklady Akademii Nauk SSSR. 1968; 180(6):1303–1305.
[36] MacCormack, R.W. The effect of viscosity in hypervelosity impact cratering. AIAA Paper. 1969: 69–354.
[37] Ostapenko, V.V. On the strong monotonicity of nonlinear difference schemes. Computational Mathematics and Mathematical Physics. 1998; 38(7):1119–1133.
Bibliography link:
Kovyrkina O.A., Ostapenko V.V. On monotonicity and accuracy of CABARET scheme for calculation of weak solutions with shocks // Computational technologies. 2018. V. 23. ¹ 2. P. 37-54
|
