Stability in the Numerical Treatment of Volterra Integral and Integro-Differential Equations with emphasis on Finite Recurrence Relations.

In the last two decades the theory of Volterra integral equations and of integro-differential equations has developed extensively. New classes of methods for the numerical solution of such equations have been developed and at the same time there have been advances in the qualitative theory of these equations. More frequent use is being made of Volterra equations to model various physical and biological phenomenon as new insight has occurred into the asymptotic behaviour of solutions. In consequence, there has emerged a need for reliable and efficient methods for the numerical treatment of such equations. This thesis is concerned with an aspect of numerical solution of Volterra integral and integro-differential equations. In Chapters 1 and 2 we are concerned with background material. We provide results on the classical theory of Volterra equations in Chapter 1 and on numerical methods in Chapter 2. The original material is contained in Chapters 3, 4 and 5. Here, stability results which involve the construction and analysis of finite-term recurrence relations are presented. The techniques relate to the treatment of Volterra integral and integro-differential equations. They permit the analysis of classical and 7-modified numerical methods. The results presented should be viewed as a contribution towards an understanding of numerical stability for the methods considered. The area is one in which further work (subsequent to the present investigation and involving advanced techniques) has been performed and where open questions still remain. The techniques which are employed in this thesis are applicable in other areas of numerical analysis and therefore have intrinsic interest

A Rabies Model with Distributed Latent Period and Territorial and Diffusing Rabid Foxes

abstract: Rabies is an infectious viral disease. It is usually fatal if a victim reaches the rabid stage, which starts after the appearance of disease symptoms. The disease virus attacks the central nervous system, and then it migrates from peripheral nerves to the spinal cord and brain. At the time when the rabies virus reaches the brain, the incubation period is over and the symptoms of clinical disease appear on the victim. From the brain, the virus travels via nerves to the salivary glands and saliva. A mathematical model is developed for the spread of rabies in a spatially distributed fox population to model the spread of the rabies epizootic through middle Europe that occurred in the second half of the 20th century. The model considers both territorial and wandering rabid foxes and includes a latent period for the infection. Since the model assumes these two kinds of rabid foxes, it is a system of both partial differential and integral equations (with integration over space and, occasionally, also over time). To study the spreading speeds of the rabies epidemic, the model is reduced to a scalar Volterra-Hammerstein integral equation, and space-time Laplace transform of the integral equation is used to derive implicit formulas for the spreading speed. The spreading speeds are discussed and implicit formulas are given for latent periods of fixed length, exponentially distributed length, Gamma distributed length, and log-normally distributed length. A number of analytic and numerical results are shown pertaining to the spreading speeds. Further, a numerical algorithm is described for the simulation of the spread of rabies in a spatially distributed fox population on a bounded domain with Dirichlet boundary conditions. I propose the following methods for the numerical approximation of solutions. The partial differential and integral equations are discretized in the space variable by central differences of second order and by the composite trapezoidal rule. Next, the ordinary or delay differential equations that are obtained this way are discretized in time by explicit continuous Runge-Kutta methods of fourth order for ordinary and delay differential systems. My particular interest is in how the partition of rabid foxes into territorial and diffusing rabid foxes influences the spreading speed, a question that can be answered by purely analytic means only for small basic reproduction numbers. I will restrict the numerical analysis to latent periods of fixed length and to exponentially distributed latent periods. The results of the numerical calculations are compared for latent periods of fixed and exponentially distributed length and for various proportions of territorial and wandering rabid foxes. The speeds of spread observed in the simulations are compared to spreading speeds obtained by numerically solving the analytic formulas and to observed speeds of epizootic frontlines in the European rabies outbreak 1940 to 1980.Dissertation/ThesisDoctoral Dissertation Applied Mathematics 201

Symbolic-Numeric Methods for Nonlinear Integro-Differential Modeling

International audienceThis paper presents a proof of concept for symbolic and numeric methods dedicated to the parameter estimation problem for models formulated by means of nonlinear integro-differential equations (IDE). In particular, we address: the computation of the model input-output equation and the numerical integration of IDE systems

Effect of Random Parameter Switching on Commensurate Fractional Order Chaotic Systems

This is the author accepted manuscript. The final version is available from Elsevier via the DOI in this record.The paper explores the effect of random parameter switching in a fractional order (FO) unified chaotic system which captures the dynamics of three popular sub-classes of chaotic systems i.e. Lorenz, Lu and Chen's family of attractors. The disappearance of chaos in such systems which rapidly switch from one family to the other has been investigated here for the commensurate FO scenario. Our simulation study show that a noise-like random variation in the key parameter of the unified chaotic system along with a gradual decrease in the commensurate FO is capable of suppressing the chaotic fluctuations much earlier than that with the fixed parameter one. The chaotic time series produced by such random parameter switching in nonlinear dynamical systems have been characterized using the largest Lyapunov exponent (LLE) and Shannon entropy. The effect of choosing different simulation techniques for random parameter FO switched chaotic systems have also been explored through two frequency domain and three time domain methods. Such a noise-like random switching mechanism could be useful for stabilization and control of chaotic oscillation in many real-world applications

SciCADE 95: International conference on scientific computation and differential equations

This report consists of abstracts from the conference. Topics include algorithms, computer codes, and numerical solutions for differential equations. Linear and nonlinear as well as boundary-value and initial-value problems are covered. Various applications of these problems are also included

Resonances and resonant frequencies for a class of nonlinear systems

Resonant phenomena for a class of nonlinear systems, which can be described by a SDOF model with a polynomial type nonlinear stiffness, are investigated using Nonlinear Output Frequency Response Functions (NOFRFs). The concepts of resonance and resonant frequencies are proposed for the first time for a class of nonlinear systems. The effects of damping on the resonances and resonant frequencies are also analyzed. These results produce a novel interpretation of energy transfer phenomena in this class of nonlinear systems and show how the damping effect influences the system resonant frequencies and amplitudes. The results are important for the design and fault diagnosis of mechanical systems and structures which can be described by the SDOF nonlinear model