Changes between Version 12 and Version 13 of doc/tec/advection


Ignore:
Timestamp:
Jan 10, 2011 2:58:31 PM (14 years ago)
Author:
suehring
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • doc/tec/advection

    v12 v13  
    5353Fig. 1 shows the dispersion and dissipation error as a function of the dimensionless wavenumber κ Δx for WS3 (3^rd^ order scheme), WS4 (4^th^ order scheme), WS5, WS6 and the 2^nd^ order scheme of Piascek and Williams (1970) (PW).
    5454The dispersion error of the upwind schemes and the dispersion error of the next higher, even ordered scheme are identical. Generally the dispersion error decreases with increasing order of the dicretization. However, no of the depicted schemes is able to adequately resolve structures with wavelengths near 2-Δx (generally no scheme based on finite differences is capable to do this).
    55 The centered, even ordered schemes hold no dissipation errors. With increasing order the numerical dissipation is more local. So the maximum wavelength affected by the dissipation term is round about 8-Δx with WS5, whereas wavelength of 16-Δx are still affected with WS3. Accordingly to the maximum of the amplification factor at κ Δx = 1.69 (these waves become unstable at first) in conjunction with the used Runge-Kutta method (Baldauf, 2008), the 5^th^ order dissipation is sufficient to avoid instabilities.
     55The centered, even ordered schemes hold no dissipation errors. With increasing order the numerical dissipation is more local. So the maximum wavelength affected by the dissipation term is round about 8-Δx with WS5, whereas wavelength of 16-Δx are still affected with WS3. Accordingly to the maximum of the amplification factor at κ Δx = 1.69 (these waves become unstable at first) in conjunction with the used [../rk3 Runge-Kutta method] (Baldauf, 2008), the 5^th^ order dissipation is sufficient to avoid instabilities.
    5656The maximum stable Courant-number is C,,r,, = 1.4 (Baldauf, 2008).
    5757
    58 '''Note: A stable numerical solution can only be guaranteed with the 3 rd order Runge-Kutta method.'''
     58'''Note: A stable numerical solution can only be guaranteed with the 3 rd order [../rk3 Runge-Kutta method].'''
    5959
    6060=== Boundaries ===