Some fast variants of TRAP scheme for solving advection equation - comparison with other schemes. (English) Zbl 1142.76441
Summary: The advection scheme TRAP was elaborated for the Bulgarian three-dimensional PC-oriented Eulerian air pollution model. The TRAP scheme uses polynomial fit of concentration profile as the BOTT scheme but differs from it in calculating fluxes. Instead of integrating the polynomial fit over the neighboring grid values as in the BOTT scheme, the flux area is supposed to be trapezoidal and is determined as a product of the Courant number and a single value of the approximation polynomial referring the middle of the passed distance. The obtained scheme is explicit, positively definite and conservative with limited numerical dispersion and good transport ability. Displaying the same properties as the BOTT scheme, the TRAP scheme turns out to be several times faster. Following the TRAP approach some schemes are developed and presented. Some of them are obtained by optimization of the so-called “normalization” procedure. Second- and third-order polynomials are used for fitting concentration profile around each grid point. The performance quality of these as well as of some other well-known schemes is determined by exploiting the rotation and deformation tests. Recommendations for the proper usage of different schemes are made.Finally, an approach for treating non-homogeneous grids with advection schemes developed for the homogeneous ones is described and demonstrated on the TRAP scheme.
MSC:
| 76M25 | Other numerical methods (fluid mechanics) (MSC2010) |
References:
| [1] | WMO-TCSU, (Numerical Methods Used in Atmospheric Models, Vol. I and II. Numerical Methods Used in Atmospheric Models, Vol. I and II, GARP Publication series, vol. 17 (1979)) |
| [2] | Rood, R. B., Numerical advection algorithms and their role in atmospheric transport and chemistry models, Reviews on Geophysics, 25, 71-100 (1987) |
| [3] | E. Berge, L. Tarrason, An evaluation of eulerian advection methods for the modelling of long-range transport of air pollution, EMEP/MSC-W Note 2/1992, 1992; E. Berge, L. Tarrason, An evaluation of eulerian advection methods for the modelling of long-range transport of air pollution, EMEP/MSC-W Note 2/1992, 1992 |
| [4] | Peters, L. K.; Berkowitz, C. M.; Carmichael, G. R.; Easter, R. C.; Fairweather, G.; Ghan, S. J.; Hales, J. M.; Leung, L. R.; Pennell, W. R.; Potra, F. A.; Saylor, R. D.; Tsang, T. T., The current state and future direction of Eulerian models in simulation the tropospheric chemistry and transport of trace species: A review, Atmospheric Environment, 29, 2, 189-222 (1995) |
| [5] | Riehtmyer, R. D., Difference Methods for Initial Value Problems (1957), Interscience Publishers: Interscience Publishers NY, 238 pp · Zbl 0079.33702 |
| [6] | R.D. Riehtmyer, A survey of difference methods for non-steady fluid dynamics, NCAR Technical Note 63-2, National Center for Atmospheric Research, Boulder, CO, 1962, 325 pp; R.D. Riehtmyer, A survey of difference methods for non-steady fluid dynamics, NCAR Technical Note 63-2, National Center for Atmospheric Research, Boulder, CO, 1962, 325 pp |
| [7] | Lax, P.; Wendroff, B., Difference schemes for hyperbolic equations with high order of accuracy, Communications on Pure and Applied Mathematics, 381-398 (1964) · Zbl 0233.65050 |
| [8] | Roach, P., Computational hydrodynamics (1980), Mir: Mir Moscow, Russian translation |
| [9] | Holmgren, P., An advection algorithm and an atmospheric airflow application, Journal of Computational Physics, 115, 27-42 (1994) · Zbl 0830.76065 |
| [10] | Russell, G. L.; Lerner, J. A., A new finite-differencing scheme for the tracer transport equation, Journal of Applied Meteorology, 20, 1483-1498 (1981) |
| [11] | Leith, C. E., (Numerical Simulation of the Earth’s Atmosphere. Numerical Simulation of the Earth’s Atmosphere, Methods in Computational Physics, vol. 4 (1965), Academic Press: Academic Press NY), 1-28 |
| [12] | Crowley, W. P., Second order numerical advection, Journal of Computational Physics, 1, 4, 471-484 (1967) · Zbl 0149.44403 |
| [13] | Crowley, W. P., Numerical advection experiments, Monthly Weather Review, 96, 1, 1-11 (1968) · Zbl 0177.56504 |
| [14] | Eliassen, A., The OECD study of long range transport of air pollutants: Long range transport modeling, Atmospheric Environment, 12, 4, 479-487 (1978) |
| [15] | Egan, A. B.; Mahoney, J. R., Numerical modeling of advection and diffusion of urban area source pollutants, Journal of Applied Meteorology, 11, 312-322 (1972) |
| [16] | Pedersen, L. B.; Prahm, L. P., A method for numerical solution of the advection equation, Tellus, XXVI, 5, 594-602 (1974) |
| [17] | Pepper, D. W.; Long, P. E., A comparison of results using second-order moments with and without width correction to solve the advection equation, Journal of Applied Meteorology, 17, 228-233 (1978) |
| [18] | Prather, M. J., Numerical advection by conservation of second-order moment, Journal of Geophysical Research, 91, D6, 6671-6681 (1986) |
| [19] | Kreiss, H. O.; Oliger, J., Comparison of accurate methods for integration of hyperbolic equation, Tellus, XXIV, 2, 199-215 (1972) |
| [20] | Prahm, L. P.; Christensen, O., Long range transmission of pollutants simulated by a two-dimensional pseudospectral dispersion model, Journal of Applied Meteorology, 16, 896-910 (1977) |
| [21] | Zlatev, Z., Mathematical model for studying the sulphur pollution in Europe, Journal of Computational and Applied Mathematics, 12, 651-666 (1985) · Zbl 0545.92018 |
| [22] | Zlatev, Z.; Berkowicz, R., Numerical treatment of large-scale air pollution model, Journal of Computational and Applied Mathematics, 16, 93-109 (1988) · Zbl 0647.76025 |
| [23] | Bott, A., A positive definite advection scheme obtained by non-linear renormalization of the advection fluxes, Monthly Weather Review, 117, 1006-1012 (1989) |
| [24] | Bott, A., Monotone flux limitation in the area-preserving flux form advection algorithm, Monthly Weather Review, 120, 2592-2602 (1992) |
| [25] | Bott, A., The monotone area-preserving flux-form advection algorithm: Reducing the time-splitting error in two-dimensional flow fields, Monthly Weather Review, 121, 2637-2641 (1993) |
| [26] | Tremback, C. J.; Powell, J.; Cotton, W. R.; Pielke, R. A., The forward-in-time upstream advection scheme: Extension to higher orders, Monthly Weather Review, 115, 540-555 (1987) |
| [27] | Syrakov, D., On a PC-oriented eulerian multi-level model for long-term calculations of the regional sulphur deposition, (Gryning, S. E.; Schiermeier, F. A., Air Pollution Modelling and its Application XI (1995), Plenum Press: Plenum Press NY, London), 645-646 |
| [28] | Syrakov, D., On the TRAP advection scheme — Description, tests and applications, (Geernaert, G.; Walloe-Hansen, A.; Zlatev, Z., Regional Modelling of Air Pollution in Europe (1996), National Environmental Research Institute: National Environmental Research Institute Denmark), 141-152 |
| [29] | Syrakov, D.; Galperin, M., On a new BOTT-type advection scheme and its further improvement, (Hass, H.; Ackermann, I. J., Proceedings of the first GLOREAM Workshop (1997), Ford Forschungszentrum Aachen: Ford Forschungszentrum Aachen Aachen, Germany), 103-109 |
| [30] | Syrakov, D.; Galperin, M., On some explicit advection schemes for dispersion modelling applications, International Journal of Environment and Pollution, 14, 267-277 (2000) |
| [31] | Syrakov, D., Once more on the advection schemes: Description of trap-schemes, (Borrego, C.; Incecik, S., Air Pollution Modeling and Its Application XVI (2004), Kluwer: Kluwer NY), 291-304 |
| [32] | H. Kirova-Galabova, S. Petrova, D. Syrakov, M. Prodanova, On some new advection schemes for air pollution modeling application, in: Proc. of the Third Intern. Symposium on Air Quality Management at Urban, Regional and Global Scales, 26-30 September 2005, Istanbul-Turkey, v.I, 2005 pp. 598-607; H. Kirova-Galabova, S. Petrova, D. Syrakov, M. Prodanova, On some new advection schemes for air pollution modeling application, in: Proc. of the Third Intern. Symposium on Air Quality Management at Urban, Regional and Global Scales, 26-30 September 2005, Istanbul-Turkey, v.I, 2005 pp. 598-607 |
| [33] | D. Syrakov, Hr. Kirova, S. Petrova, M. Prodanova, Creation and testing of flux-type advection schemes for air pollution modeling application, in: Proc. of the 10th Int. Conf. on “Harmonization within Atmospheric Dispersion Modelling for Regulatory Purposes”, Sissi (Malia), Crete, Greece, 17-20 October, 2005, pp. 278-282; D. Syrakov, Hr. Kirova, S. Petrova, M. Prodanova, Creation and testing of flux-type advection schemes for air pollution modeling application, in: Proc. of the 10th Int. Conf. on “Harmonization within Atmospheric Dispersion Modelling for Regulatory Purposes”, Sissi (Malia), Crete, Greece, 17-20 October, 2005, pp. 278-282 |
| [34] | MacCormack, R. W., (Proceedings of the Second International Conference on Numerical Methods in Fluid Dynamics. Proceedings of the Second International Conference on Numerical Methods in Fluid Dynamics, Lecture Notes in Physics, vol. 8 (1971), Springer-Verlag: Springer-Verlag NY) |
| [35] | Galperin, M. V., Two numerical advection schemes for atmospheric models: Development and comparison testing, (Syrakov, D.; Batchvarova, E.; Wiman, B., Proceedings of Swedish-Bulgarian Workshop. Proceedings of Swedish-Bulgarian Workshop, Sozopol, Bulgaria (1997), Pensoft Publishers: Pensoft Publishers Sofia), 49-61 |
| [36] | Galperin, M. V., (The Approaches to Correct Computation of Airborne Pollution Advection. The Approaches to Correct Computation of Airborne Pollution Advection, Problems of Ecological Monitoring and Ecosystem Modelling, vol. XVII (2000), Gidrometeoizdat: Gidrometeoizdat St. Petersburg), 54-68 |
| [37] | Smolarkiewiecz, P. K., The multidimensional Crowley advection scheme, Monthly Weather Review, 113, 1109-1130 (1982) |
| [38] | Staniforth, A.; Côté, J.; Pudikiewicz, J., Comments on “Smolarkiewicz’s deformational flow”, Monthly Weather Review, 115, 894-900 (1987) |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
