In many transport-chemistry models, a huge system ofODEs of the advection-diffusion-reaction type has to be integrated in time. Typically, this is done with the help of operator splitting. Rosenbrock schemes combined with approximate matrix factorization (ROS-AMF) are an alternative to operator splitting which does not suffer from splitting errors. However, implementation of ROS-AMF schemes often requires serious changes in the code. In this paper we test another classical second order splitting introduced by Strang in 1963, which, unlike the popular Strang splitting, seemed to be forgotten and rediscovered recently (partially due to its intrinsic parallellism). 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 favorably to the Strang splitting, but is less accurate than ROS-AMF.
展开▼