Spatial and temporal discretization has a big impact to the overall accuracy of the numerical scheme. The spatial discretization affects the calculation of physical flux or other space-derived variables, which often require high order scheme in the non-linear problem to achieve the desired accuracy. While temporal discretization also has a significant impact on the stability of the numerical scheme and total computational cost. Strong stability preserving Runge Kutta scheme often used as temporal discretization scheme for high order numerical scheme. In this work, we are motivated to study the significance of temporal discretization order (up to fourth order) to the overall numerical scheme accuracy and total computational cost while high order weighted essentially non-oscillatory scheme used as spatial discretization technique in hyperbolic conservation law. We found that the third order SSP Runge Kutta scheme is an effective scheme with respect to accuracy and computational cost.

1.
J. A.
Ekaterinaris
,
Progress in Aerospace Sciences
41
(
3
),
192
300
(
2005
).
2.
G. A.
Sod
,
Journal of Computational Physics
27
(
1
),
1
31
(
1978
).
3.
S.
Gottlieb
and
D. I.
Ketcheson
, in
Handbook of Numerical Analysis
, edited by
R.
Abgrall
and
C.-W.
Shu
(
Elsevier
,
2016
), Vol.
17
, pp.
549
583
.
4.
C.-W.
Shu
, in
Advanced numerical approximation of nonlinear hyperbolic equations
(
Springer
,
1998
), pp.
325
432
.
5.
C.
Huang
and
L. L.
Chen
,
Journal of Computational Physics
352
,
498
515
(
2018
).
6.
E. F.
Toro
,
M.
Spruce
and
W.
Speares
,
Shock waves
4
(
1
),
25
34
(
1994
).
7.
S.
Gottlieb
,
Journal of Scientific Computing
25
(
1
),
105
128
(
2005
).