21.11.2014 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!