Abstract
We apply a Runge-Kutta-based waveform relaxation method to initial-value problems for implicit differential equations. In the implementation of such methods, a sequence of nonlinear systems has to be solved iteratively in each step of the integration process. The size of these systems increases linearly with the number of stages of the underlying Runge-Kutta method, resulting in high linear algebra costs in the iterative process for high-order Runge-Kutta methods. In our earlier investigations of iterative solvers for implicit initial-value problems, we designed an iteration method in which the linear algebra costs are almost independent of the number of stages when implemented on a parallel computer system. In this paper, we use this parallel iteration process in the Runge-Kutta waveform relaxation method. In particular, we analyse the convergence of the method. The theoretical results are illustrated by a few numerical examples.
Similar content being viewed by others
References
A. Bellen and F. Tagliaferro, A combined WR-parallel steps method for ordinary differential equations, in: Proc. Conf. Parallel Computing: Achievements, Problems and Prospects, Capri, Italy (1990).
K. Burrage, Parallel and Sequential Methods for Ordinary Differential Equations (Clarendon Press, Oxford, 1995).
W. M. Lioen, J. J. B. de Swart and W. A. van der Veen, Test set for IVP solvers, available via WWW at URL: http://www.cwi.nl/cwi/projects/IVPtestset.shtml (1996).
B. A. Gottwald, MISS: A simple simulation system for biological and chemical processes, EDV in Medizin und Biologie 3 (1977) 85–90 (in German).
E. Hairer and G. Wanner, Solving Ordinary Differential Equations II; Stiff and Differential-Algebraic Problems, Springer Series in Computational Mathematics 14 (Springer, Berlin, 1991).
K. J. in 't Hout, On the convergence of waveform relaxation methods for stiff nonlinear ordinary differential equations, Appl. Numer. Math. 18 (1995) 175–190.
E. Lelarasmee, The waveform relaxation method for the time domain analysis of large scale nonlinear dynamical systems, Ph.D. thesis, University of California, Berkeley (1982).
E. Lelarasmee, A. E. Ruehli and A. L. Sangiovanni-Vincentelli, The waveform relaxation method for time domain analysis of large scale integrated circuits, IEEE Trans. CAD IC Syst. 1 (1982) 131–145.
U. Miekkela and O. Nevanlinna, Convergence of dynamic iteration methods for initial value problems, SIAM J. Sci. Statist. Comput. 8 (1987) 459–482.
E. Schäfer, A new approach to explain the “High Irradiance Responses” of photomorphogenesis on the basis of phytochrome, J. Math. Biology 2 (1975) 41–56.
P. J. van der Houwen, B. P. Sommeijer and W. A. van der Veen, Parallelism across the steps in iterated Runge-Kutta methods for stiff initial value problems, Numer. Algorithms 8 (1994) 293–312.
P. J. van der Houwen and J. J. B. de Swart, Triangularly implicit iteration methods for ODE-IVP solvers, SIAM J. Sci. Comput. 18 (1997) 41–55.
P. J. van der Houwen and W. A. van der Veen, Solving implicit differential equations on parallel computers, CWI Report NM-R9526 (1995), submitted for publication.
J. K. White and A. L. Sangiovanni-Vincentelli, Relaxation Techniques for the Simulation of VLSI Circuits (Kluwer, Dordrecht, 1987).
Rights and permissions
About this article
Cite this article
van der Houwen, P.J., van der Veen, W.A. Waveform relaxation methods for implicit differential equations. Advances in Computational Mathematics 7, 183–197 (1997). https://doi.org/10.1023/A:1018942718589
Issue Date:
DOI: https://doi.org/10.1023/A:1018942718589