Application of operator splitting to the Maxwell equations including a source term

Motivated by numerical solution of the lime-dependent Maxwell equations, we consider splitting methods for a linear system of differential equations omega'(t) = A omega(t) + f(t), A is an element of R-nxn split into two subproblems omega(1)'(t) = A(1)omega(1)(t) + f(1)(t) and omega(2)' (t) = A(2)omega(2) + f(2)(t), A = A(1) + A(2), f = f(1) + f(2). First, expressions for the leading term ofthe local errorare derived for the Strang-Marchuk and the symmetrically weighted sequential splitting methods. The analysis, done in assumption that the subproblems are solved exactly, confirms the expected second order global accuracy of both schemes. Second, several relevant numerical tests are performed for the Maxwell equations discretized in space either by finite differences or by finite elements. An interesting case is the splitting into the Subproblems omega(1)', = A omega(1) and omega(2)' = f (with the split-off source term f). For the central finite difference staggered discretization. we consider second order splitting schemes and compare them to the classical Yee scheme on a test problem with loss and source terms. For the vector Nedelec finite element discretizations, we test the Gautschi-Krylov time integration scheme. Applied in combination with the split-off source term, it leads to splitting schemes that are exact per split step. Thus. the time integration error of the schemes consists solely of the splitting error. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.