Changes between Version 19 and Version 20 of doc/tec/advection
- Timestamp:
- Aug 28, 2012 4:18:39 PM (13 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/tec/advection
v19 v20 17 17 }}} 18 18 denotes the fluxes staggered half a grid length related to the advected quantity. \\ 19 Wicker and Skamarock (2002) di cretized the 6^th^ and 5^th^ order fluxes as follows:19 Wicker and Skamarock (2002) discretized the 6^th^ and 5^th^ order fluxes as follows: 20 20 {{{ 21 21 #!Latex … … 52 52 [[Image(prop.png, 700px, border=1)]] \\ 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 Piacsek and Williams (1970) (PW). 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 di cretization. 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).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 discretization. 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 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).
