Authors: Mike Botchev, Istvan Farago, Agnes Havasi
Addresses: Department of Applied Mathematics, University of Twente, The Netherlands. ' Eotvos Lorand University, Pazmany P. setany 1/C, Budapest, Hungary. ' Eotvos Lorand University, Pazmany P. setany 1/C, Budapest, Hungary
Abstract: In many transport-chemistry models, a huge system of one-dimensional equations of the advection-diffusion-reaction type has to be integrated in time. Typically, this is done with the help of operator splitting. Operator splitting is attractive for complex large-scale transport-chemistry models because it can handle different processes separately in different parts of the computer program. Rosenbrock schemes combined with approximate matrix factorisation (ROS-AMF) are an alternative to operator splitting which does not suffer from splitting errors. However, since the ROS2-AMF schemes are not based on operator splitting, implementation of these methods often requires major changes in the code. In this paper we test another second-order splitting introduced by Strang in 1963, which seemed to be forgotten and rediscovered recently (partially due to its intrinsic parallelism). This splitting, called symmetrically weighted sequential (SWS) splitting, is simple and straightforward to apply, independent of the order of the operators and has an operator-level parallelism. In the experiments, the SWS scheme compares favourably with Strang splitting, but is less accurate than ROS-AMF.
Keywords: numerical comparison; parallelisation; ROS-AMF method; Strang splitting; chemistry transport model; weighted splitting.
International Journal of Environment and Pollution, 2004 Vol.22 No.1/2, pp.3 - 16
Available online: 12 Oct 2004 *Full-text access for editors Access for subscribers Purchase this article Comment on this article