Research Article  Open Access
Jürgen Geiser, "Iterative Splitting Methods for Integrodifferential Equations: Theory and Applications", Journal of Applied Mathematics, vol. 2014, Article ID 812137, 12 pages, 2014. https://doi.org/10.1155/2014/812137
Iterative Splitting Methods for Integrodifferential Equations: Theory and Applications
Abstract
We present novel iterative splitting methods to solve integrodifferential equations. Such integrodifferential equations are applied, for example, in scattering problems of plasma simulations. We concentrate on a linearised integral part and a reformulation to a system of first order differential equations. Such modifications allow for applying standard iterative splitting schemes and for extending the schemes, respecting the integral operator. A numerical analysis is presented of the system of semidiscretised differential equations as abstract Cauchy problems. In the applications, we present benchmark and initial realistic applications to transport problems with scattering terms. We also discuss the benefits of such iterative schemes as fast solver methods.
1. Introduction
The motivation is to apply fast splitting schemes that have been developed for systems of differential equations, to integrodifferential equations. The applications are related to transport problems with scattering terms, for example, plasma simulations with collision terms; see [1, 2]. Due to the physics of the situation, the numerical implementation allows decoupling the system into transport and scattering terms, so that the underlying methods can employ optimal splitting schemes.
We concentrate on a semidiscretised transportscattering equation, given as a Cauchy problem. In the following, we consider the abstract homogeneous Cauchy problem as a system of ordinary integrodifferential equations: where are bounded and nonlinear operators, which are derived from the semidiscretised transport and scattering operators. is a Banach space and is the corresponding norm in . The operator represents the transport processes, and the operator is a given scattering operator.
This paper is organised as follows.
In Section 2 we present the mathematical model and a possible reduced model for further approximations. The numerical schemes are presented in Section 3. The results of some numerical experiments are given in Section 4. The contents of Section 5 contain a summary of the results.
2. Analytical Approach to the Integrodifferential Equation
In the following, we discuss the transformation to a system of differential equations in order to apply iterative splitting schemes; see also the decomposition of matrices in [3].
For the linearisation, we make the following assumptions.
Assumption 1. We assume a weak nonlinearity and employ the following simplification.(1)The Jacobian of , given as , is linearised at , which means , where .(2)The operator is linearised at , which means , where .(3)The initial condition of the first derivative is given by .(4)For the linearised case, the following matrix operations can be applied: where and are simultaneous and diagonalisable, with the same matrix , to the diagonal matrices and . Furthermore, and have eigenvectors as a basis for . Further and commute.
Corollary 2. The integrodifferential equation (1) can be transformed, under Assumption 1, to the following differential equation:
where the analytical solution is
where .
While and commute, we also have the commutation of and and, for small , we have the solution
Proof. The analytical solution of the second order differential equation (4) is given as follows. We rewrite (4) into
The exponential matrix can be decoupled via the Cayley Hamilton theorem as
where is the identity matrix.
The characteristic polynomial of is given as
and the roots are given as
Then and ; see also [4]. Therefore the analytical solution is given by (7).
We apply this solution to our integrodifferential equation (1) and, with Assumption 1, the equation is satisfied.
A further simplification can be made to rewrite the integraldifferential equation as two first order differential equations. Later such a reduction will allow us to apply fast iterative splitting methods.
Corollary 3. The integrodifferential equation (1) can be transformed, with the help of Assumption 1, to two first order differential equations:
where and .
The analytical solution is
Proof. The analytical solution of the first order differential equations (14) and (16) is obtained using their characteristic polynomials:
hence the solution is, using the notations of and ,
and therefore the analytical solution is given by (18) with .
Therefore this is the solution of our integrodifferential equation (1) under Assumption 1.
3. Iterative Splitting Approach to the Integrodifferential Equation
Operator splitting methods are used to solve complex models in geophysical and environmental physics. They have been developed and applied in [5–7]. Such ideas are the basis for this paper: solving the simpler equations but using higher order discretisation methods for the remaining, more complicated equations. With this aim, we use the operator splitting method and decouple the system of equations as will be described in the following.
We will concentrate on the iterative splitting method. This algorithm is based on an iteration with fixed splitting discretisation stepsize ; namely, on the time interval we solve the following subproblems consecutively for (cf. [8, 9]): where is the known split approximation at the time level . The split approximation at the time level is defined by . (Clearly, the function depends on the interval , too, but, for the sake of simplicity, in our notation, we omit the dependence on .)
For the simplification, we apply two standard iterative splitting schemes to the first order differential equations.
The first iterative process is and the second iterative process is where is the known split approximation at the time level . The split approximation at the time level is defined by .
In the following, we will analyse the convergence and the rate of convergence of the methods (22) and (23) for tending to infinity for the linear operators using Assumption 1.
It suffices to discuss one case, since the other case is the same with only a change in the sign of the operator .
Theorem 4. Let us consider the abstract Cauchy problem in a Banach space where are given linear operators which are generators of the semigroup and is a given element. Then the iteration process (22) is convergent and the rate of convergence is of second order.
Proof. The proof is given in [10].
In the following Algorithm 5, we have combined the first and second iterative processes.
Algorithm 5. We apply and . The timesteps are , with , and where is a given number of iterative steps.
We have the following algorithm: ,
the starting values , results of the last iteration, , ,
,
we go to the next timestep , else .
4. Numerical Experiments
For the applications part, we will now treat test problems and reallife applications to transportscattering problems.
The first examples are to validate our novel splitting and solver methods; then, in the subsequent examples, we use more delicate coupled differential equations.
4.1. First Example: A Matrix Problem with Integral Term
We deal with simpler integrodifferential equations: where we take as a first order approximation of the integral and treat where and we have the analytical solution for the approximation
We split this into
We have the following solutions for the iterative scheme: where and while the timesteps are given by .
We employ the following recurrence relations, distinguishing between even and odd iterations.
For the odd iterations, , for ,
For the even iterations, , for , We have a blowup solution of exponential character in long time behavior; see Figure 1. Therefore, it is advisable that the integrator should be highly accurate for long time behavior.
(a)
(b)
Table 1 gives the numerical results of the iterative splitting scheme.

Figure 2 presents the onesided and twosided iterative results.
(a)
(b)
Remark 6. In the experiments, we obtain improved results with each additional step. By the way, the solution blows up and we have to use also very fine timesteps to control the errors. Optimal results are obtained by using the integral part (stiff part) as the implicit part in the iteration (onesided over ).
4.2. Second Example: Integrodifferential Equation
We treat the system of integrodifferential equations where we have the analytical solution where and, for example, we take .
We have employed the integral approximation
The matrices are given by where the parameters are , and .
The splitting into the two operators and is made in order to accelerate the computation of the exponential matrices and .
(1) The operator is computed as where .
We have the following solutions for the iterative scheme: where and while the timesteps are given by .
We employ the following recurrence relations, distinguishing between even and odd iterations.
For the odd iterations, , for ,
For the even iterations, , for ,
For the integral computations, we apply higher order schemes, the NewtonCotes formula (e.g., Simpson’s rule and Boole’s rule); see [11].
We acchieve for the second example (vectorial benchmark problem) the same accuracy for its first component as in the first example (scalar benchmark problem).
the same accuracy for its first component as in the first example (scalar benchmark problem).
Figure 3 presents the onesided and twosided iterative results.
(a)
(b)
The computational results are given in Figure 4, which presents the onesided and twosided iterative results.
(a)
(b)
Remark 7. In the experiments, we obtain improved results with each additional step. By the way, the solution blows up and so we have to use very fine timesteps to control the errors. The optimal results are obtained by using the integral part (stiff part) as the implicit part in the iteration (onesided over ).
4.2.1. Exact Integral Part Based on the Splitting Approach
We treat a system of integrodifferential equations where we assume that the primitive of is .
We have the analytical solution where and, for example, we take .
The treatment of the square root for matrices is given in [12, 13].
We have the following matrices: where the parameters are , and .
The exponentials of these operators are computed using Padé approximants; see the discussion in [11].
We have the following solutions for the two iterative schemes.
Algorithm 8. Iterative process for : For the odd iterations, , we have, for , For the even iterations, , we have, for , where , , and the timesteps are given by .
Algorithm 9. Iterative process for : For the odd iterations, , we have, for , For the even iterations, , we have, for , where , , and the timesteps are given by .
The solution after each iterative step with the two single solutions is
For the integral computations, we apply higher order schemes, NewtonCotes formulas (e.g., Simpson’s rule and Boole’s rule); see [11].
Table 2 shows the numerical results of the iterative splitting scheme. We also acchieve similar results for the first component of the vectorial benchmark problem as in its scalar form in the first experiment.

Figure 5 presents the onesided and twosided iterative results for an optimised approach.
(a)
(b)
The computational results are given in Figure 6; we present the onesided and twosided iterative results for an optimised approach.
(a)
(b)
Remark 10. In the experiments, we obtain improved results with each additional step. By the way, the solution blows up and we have to use very fine timesteps to control the errors. The optimal results are obtained by using the integral part (stiff part) as the implicit part in the iteration (onesided over ).
Remark 11. To optimise the algorithms, we combine both of the iterative cycles with the idea of and in the Padé approximants we apply ; see Algorithm 5.
Remark 12. Table 3 presents the computational times of the three experiments. While the exact approach based on the splitting idea of higher order differential equations obtains more accurate solutions, its drawback is its higher computational costs. We reduce such costs by using the symmetry of the underlying iterative schemes and reach nearly the same computational time as in the numerical approach scheme.

4.3. RealLife Example: OnePhase Problem
The next example is a simplified reallife problem with a neutron transport equation, which includes the gain and loss of a neutron.
We concentrate on the computational benefits of a fast computation of the iterative scheme, given with matrix exponential functions.
The equation is
In the following, we deal with the semidiscretised equation given by the matrices where is the solution of the species in the mobile phase in each spatial discretisation point .
We have the following two operators for the splitting method: The reaction and scattering matrices are where is the number of spatial points.
In the first algorithm, we use the following approximation of the integral term:
We use the following recurrence relations, distinguishing between even and odd iterations.
For the odd iterations, , with , we have
For the even iterations, , with , we have where .
For the reference solution, we apply a fine timespatial scale without decoupling the equations.
For the algorithms, we employ novel Algorithm 5 and we obtain and .
The matrices for the novel algorithm are
Figure 7 presents the onesided and twosided iterative results.
(a)
(b)
The computational results are given in Figure 8, which presents the onesided and twosided iterative results.
(a)
(b)
Figure 9 presents the onesided and twosided iterative results for an optimised approach.
(a)
(b)
The computational results are given in Figure 10, which presents the onesided and twosided iterative results for an optimised approach.
(a)
(b)
Table 4 presents the computational times of the three experiments.

Remark 13. For all iterative schemes, we can reach faster results using the optimised approach. Such iterative schemes with fast computations of the exponential matrices accelerate standard splitting schemes. With four or five iterative steps we obtain more accurate results than we did using the expensive standard schemes. Onesided iterative schemes yield the best convergence results.
5. Conclusions and Discussion
We presented a coupled model for a transport and kinetic model for deposition species in a plasma environment. We proposed numerical schemes to decompose the flow field and the reactioncollision term. A generalised iterative scheme was discussed and their analysis and benefits were presented.
Such combined splitting schemes help to reduce the computational time required for delicate plasma problems. In the future, we will deal with nonlinear schemes of collision problems in plasma reactors.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
References
 J. Geiser, “An iterative splitting approach for linear integrodifferential equations,” Applied Mathematics Letters, vol. 26, no. 11, pp. 1048–1052, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 J. Geiser, “Multiscale splitting method for the BoltzmannPoisson equation: application to the dynamics of electrons,” International Journal of Differential Equations, vol. 2014, Article ID 178625, 8 pages, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 S. Blanes, F. Casas, and A. Murua, “Error analysis of splitting methods for the time dependent Schrödinger equation,” SIAM Journal on Scientific Computing, vol. 33, no. 4, pp. 1525–1548, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 R. Bhatia, Matrix Analysis, vol. 169 of Graduate Texts in Mathematics, Springer, 1997. View at: Publisher Site  MathSciNet
 G. Strang, “On the construction and comparison of difference schemes,” SIAM Journal on Numerical Analysis, vol. 5, pp. 506–517, 1968. View at: Publisher Site  Google Scholar  MathSciNet
 J. Verwer and B. Sportisse, “A note on operator splitting in a stiff linear case,” MAS Report R9830, 1998. View at: Google Scholar
 Z. Zlatev, Computer Treatment of Large Air Pollution Models, Kluwer Academic, New York, NY, USA, 1995.
 J. F. Kanney, C. T. Miller, and C. T. Kelley, “Convergence of iterative splitoperator approaches for approximating nonlinear reactive problems,” Advances in Water Resources, vol. 26, no. 3, pp. 247–261, 2003. View at: Publisher Site  Google Scholar
 E. J. Dean, R. Glowinski, and J. L. Lions, “An operator splitting approach to multilevel methods,” Applied Mathematics Letters, vol. 15, no. 4, pp. 505–511, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 J. Geiser, Iterative Splitting Methods for Differential Equations, CRC Press, Chapman & Hall, 2011. View at: Publisher Site  MathSciNet
 J. Geiser, “Computing exponential for iterative splitting methods: algorithms and applications,” Journal of Applied Mathematics, vol. 2011, Article ID 193781, 27 pages, 2011. View at: Publisher Site  Google Scholar
 S. H. Cheng, N. J. Higham, C. S. Kenney, and A. J. Laub, “Approximating the logarithm of a matrix to specified accuracy,” SIAM Journal on Matrix Analysis & Applications, vol. 22, no. 4, pp. 1112–1125, 2001. View at: Publisher Site  Google Scholar  MathSciNet
 E. D. Denman and J. Beavers, “The matrix sign function and computations in systems,” Applied Mathematics and Computation, vol. 2, no. 1, pp. 63–94, 1976. View at: Google Scholar  MathSciNet
Copyright
Copyright © 2014 Jürgen Geiser. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.