Abstract
We investigate the use of splitting methods for the numerical integration of three-dimensional transport-chemistry models. In particular, we investigate various possibilities for the time discretization that can take advantage of the parallelization and vectorization facilities offered by multi-processor vector computers. To suppress wiggles in the numerical solution, we use third-order, upwind-biased discretization of the advection terms, resulting in a five-point coupling in each direction. As an alternative to the usual splitting functions, such as co-ordinate splitting or operator splitting, we consider a splitting function that is based on a three-coloured hopscotch-type splitting in the horizontal direction, whereas full coupling is retained in the vertical direction. Advantages of this splitting function are the easy application of domain decomposition techniques and unconditional stability in the vertical, which is an important property for transport in shallow water. The splitting method is obtained by combining the hopscotch-type splitting function with various second-order splitting formulae from the literature. Although some of the resulting methods are highly accurate, their stability behaviour (due to horizontal advection) is quite poor. Therefore we also discuss several new splitting formulae with the aim to improve the stability characteristics. It turns out that this is possible indeed, but the price to pay is a reduction of the accuracy. Therefore, such methods are to be preferred if accuracy is less crucial than stability; such a situation is frequently encountered in solving transport problems. As part of the project TRUST (Transport and Reactions Unified by Splitting Techniques), preliminary versions of the schemes are implemented on the Cray C98 4256 computer and are available for benchmarking.
Similar content being viewed by others
REFERENCES
Douglas. J. jr. (1962). Alternating direction methods for three space variables. Numer. Math. 4, 41–63.
Dunsbergen, D. (1994). Particle models for transport in three-dimensional shallow water flow. Ph.D. Thesis. Delft Technical University.
Gourlay, A. R. (1970). Hopscotch: a fast second order partial differential equation solver. J. Inst. Math. Appl. 6, 375–390.
van der Houwen, P. J., and Verwer, J. G. (1979). One-step splitting methods for semidiscrete parabolic equations. Computing 22, 291–309.
Hundsdorfer, W. H. Trapezoidal and midpoint splittings for initial-boundary value problems. CWI Report NM-R9605 (to appear in Math. Comp. 1998).
van Leer, B. (1983). Upwind-difference methods for aerodynamic problems governed by the Euler equations. Proceedings of the 15th AMS-SIAM Summer Seminar on Applied Mathematics. Scripps Institution of Oceanography (B. E. Engquist, S.Osher and R. C. J. Sommerville, eds.). Lectures in Applied Mathematics 22-Part 2 (1985). 327–336, AMS. Providence, RI.
Lu, T., Neittaanmäki, P., and Tai, X.-C. (1992). A parallel splitting-up method for partial differential equations and its applications to Navier Stokes equations. Math. Model. Numer. Anal. 26, 673–708.
Sommeijer, B. P., van der Houwen, P. J., and Kok, J. (1994). Time integration of three-dimensional numerical transport models. Appl. Numer. Math. 16, 201–225.
Sommeijer, B. P., and Kok, J. (1995). Implementation and performance of the time integration of a 3D numerical transport model. Intern. J. Numer. Meth. Fluids 21, 349–367.
Sommeijer, B. P., and Kok, J. A vector parallel method for a three-dimensional transport model coupled with bio-chemical terms. CWI Report NM-R9503. CWI. Amsterdam (1995).
Strang, G. (1968). On the construction and comparison of difference schemes. SIAM J. Numer. Anal. 5, 506–517.
Verwer, J. G., and Sommeijer, B. P. (1997). Stability analysis of an odd-even line hopscotch method for three-dimensional advection-diffusion problems. SIAM J. Numer. Anal. 34, 376–388.
Yanenko, N. N. (1971). The method of fractional steps. Springer Verlag, Berlin-Heidelberg-New York. 1971.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
van der Houwen, P.J., Sommeijer, B.P. Splitting Methods for Three-Dimensional Transport Models with Interaction Terms. Journal of Scientific Computing 12, 215–231 (1997). https://doi.org/10.1023/A:1025645326705
Issue Date:
DOI: https://doi.org/10.1023/A:1025645326705