Author:
Hansen A. B.,Brandt J.,Christensen J. H.,Kaas E.
Abstract
Abstract. Various semi-Lagrangian methods are tested with respect to advection in air pollution modeling. The aim is to find a method fulfilling as many of the desirable properties by Rasch and Williamson (1990) and Machenhauer et al. (2008) as possible. The focus in this study is on accuracy and local mass conservation. The methods tested are, first, classical semi-Lagrangian cubic interpolation, see e.g. Durran (1999), second, semi-Lagrangian cubic cascade interpolation, by Nair et al. (2002), third, semi-Lagrangian cubic interpolation with the modified interpolation weights, Locally Mass Conserving Semi-Lagrangian (LMCSL), by Kaas (2008), and last, semi-Lagrangian cubic interpolation with a locally mass conserving monotonic filter byKaas and Nielsen (2009). Semi-Lagrangian (SL) interpolation is a classical method for atmospheric modeling, cascade interpolation is more efficient computationally, modified interpolation weights assure mass conservation and the locally mass conserving monotonic filter imposes monotonicity. All schemes are tested with advection alone or with advection and chemistry together under both typical rural and urban conditions using different temporal and spatial resolution. The methods are compared with a current state-of-the-art scheme, Accurate Space Derivatives (ASD), see Frohn et al. (2002), presently used at the National Environmental Research Institute (NERI) in Denmark. To enable a consistent comparison only non-divergent flow configurations are tested. The test cases are based either on the traditional slotted cylinder or the rotating cone, where the schemes' ability to model both steep gradients and slopes are challenged. The tests showed that the locally mass conserving monotonic filter improved the results significantly for some of the test cases, however, not for all. It was found that the semi-Lagrangian schemes, in almost every case, were not able to outperform the current ASD scheme used in DEHM with respect to accuracy.
Reference30 articles.
1. Bartnicki, J.: A {S}imple {F}iltering {P}rocedure for {R}emoving {N}egative {V}alues from {N}umerical {S}olutions of the {A}dvection {E}quation, Environ. Softw., 4, 187–201, 1989.
2. Bott, A.: A {P}ositive {D}efinite {A}dvection {S}cheme {O}btained by {N}onlinear {R}enormalization of the {A}dvective {F}luxes, Mon. Weather Rev., 117, 1006–1015, 1989.
3. Brandt, J.: Modelling {T}ransport, {D}ispersion and {D}eposition of {P}assive {T}racers from {A}ccidental {R}eleases, PhD thesis, Ministry of Environment and Energy National Environmental Research Institute and Ministry of Research and Information Technology Risø National Laboratory, pp. 307, 1998.
4. Brandt, J., Dimov, I., Goergiev, K., Wasniewski, J., and Zlatev, Z.: Coupling the {A}dvection and the {C}hemical {P}arts of {L}arge {A}ir {P}ollution {M}odels, Lecture Notes in Computer Science, Applied Parallel Computing, Industrial Computation and Optimization, 1184, 65–76, 1996a.
5. Brandt, J., Wasniewski, J., and Zlatev, Z.: Handling the {C}hemical {P}art in {L}arge {A}ir {P}ollution {M}odels, Appl. Math. Comp. Sci., 6, 331–351, 1996b.