Dispersion and dissipation error in high-order Runge-Kutta ...
Dispersion and dissipation error in high-order Runge-Kutta ...
Dispersion and dissipation error in high-order Runge-Kutta ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
DoF / λ<br />
DoF / λ<br />
1<br />
0.9<br />
−3e−05<br />
14<br />
−0.0001<br />
12<br />
10<br />
−0.0003<br />
8<br />
−0.001<br />
6<br />
−0.01<br />
4<br />
1<br />
0.9<br />
1<br />
0.9<br />
−3e−07<br />
−1e−06<br />
15<br />
−3e−06<br />
−1e−05<br />
10<br />
5<br />
1<br />
0.9<br />
0.8<br />
0.8<br />
0.8<br />
0.8<br />
∆ t / ∆ t max<br />
0.7<br />
0.6<br />
0.5<br />
−1e−05<br />
−3e−05<br />
−0.0001<br />
−0.0003<br />
−0.001<br />
−0.01<br />
0.7<br />
0.6<br />
0.5<br />
∆ t / ∆ t max<br />
∆ t / ∆ t max<br />
0.7<br />
0.6<br />
0.5<br />
−3e−07<br />
−1e−06<br />
−3e−06<br />
−1e−05<br />
−0.0001<br />
0.7<br />
0.6<br />
0.5<br />
∆ t / ∆ t max<br />
0.4<br />
0.3<br />
0.2<br />
25<br />
−3e−05<br />
20<br />
−0.0001<br />
15<br />
λ / h<br />
10<br />
−0.001<br />
−0.01<br />
5<br />
0.4<br />
0.3<br />
0.2<br />
0.4<br />
0.3<br />
0.2<br />
−1e−07<br />
25<br />
−3e−07<br />
20<br />
−1e−06<br />
15<br />
λ / h<br />
10<br />
−0.0001<br />
5<br />
0.4<br />
0.3<br />
0.2<br />
(a) p = 2<br />
(b) p = 3<br />
DoF / λ<br />
DoF / λ<br />
1<br />
0.9<br />
20<br />
−1e−08<br />
−1e−07<br />
15<br />
−1e−06<br />
10<br />
−1e−05<br />
−0.0001<br />
5<br />
1<br />
−0.001<br />
0.9<br />
1<br />
0.9<br />
20<br />
−1e−09<br />
−1e−08<br />
15<br />
−1e−07<br />
10<br />
1<br />
0.9<br />
∆ t / ∆ t max<br />
0.8<br />
0.7<br />
0.6<br />
0.5<br />
0.4<br />
0.3<br />
−1e−09<br />
0.2<br />
35<br />
30<br />
−1e−08<br />
25<br />
−1e−07<br />
20<br />
λ / h<br />
−1e−06<br />
15<br />
−1e−05<br />
10<br />
−0.0001<br />
5<br />
0.8<br />
0.7<br />
0.6<br />
0.5<br />
0.4<br />
0.3<br />
0.2<br />
∆ t / ∆ t max<br />
∆ t / ∆ t max<br />
0.8<br />
0.7<br />
0.6<br />
0.5<br />
−1e−10<br />
0.4<br />
0.3<br />
20<br />
−1e−09<br />
15<br />
−1e−08<br />
−1e−08<br />
λ / h<br />
−1e−07<br />
10<br />
−1e−06<br />
−1e−05<br />
−1e−06<br />
−1e−05<br />
5<br />
−0.0001<br />
−0.0001<br />
0.8<br />
0.7<br />
0.6<br />
0.5<br />
0.4<br />
0.3<br />
∆ t / ∆ t max<br />
(c) p = 4<br />
(d) p = 5<br />
Figure 8: Dissipation <strong>error</strong> as a function of wave length per mesh size (λ/h), degrees of<br />
freedom per wavelength (DoF/λ) <strong>and</strong> relative time step (∆t/∆t max ) for polynomial <strong>order</strong>s<br />
p = 2, 3, 4, 5 on a non-uniform r<strong>and</strong>om mesh<br />
28