Changes between Version 12 and Version 13 of doc/tec/advection
- Timestamp:
- Jan 10, 2011 2:58:31 PM (14 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/tec/advection
v12 v13 53 53 Fig. 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). 54 54 The 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.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 [../rk3 Runge-Kutta method] (Baldauf, 2008), the 5^th^ order dissipation is sufficient to avoid instabilities. 56 56 The maximum stable Courant-number is C,,r,, = 1.4 (Baldauf, 2008). 57 57 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].''' 59 59 60 60 === Boundaries ===