Impulsive perturbations to differential equations: stable/unstable pseudo-manifolds, heteroclinic connections, and flux

State-dependent time-impulsive perturbations to a two-dimensional autonomous flow with stable and unstable manifolds are analysed by posing in terms of an integral equation which is valid in both forwards- and backwards-time. The impulses destroy the smooth invariant manifolds, necessitating new definitions for stable and unstable pseudo-manifolds. Their time-evolution is characterised by solving a Volterra integral equation of the second kind with discontinuous inhomogeniety. A criteria for heteroclinic trajectory persistence in this impulsive context is developed, as is a quantification of an instantaneous flux across broken heteroclinic manifolds. Several examples, including a kicked Duffing oscillator and an underwater explosion in the vicinity of an eddy, are used to illustrate the theory

From regional pulse vaccination to global disease eradication: insights from a mathematical model of Poliomyelitis

Mass-vaccination campaigns are an important strategy in the global fight against poliomyelitis and measles. The large-scale logistics required for these mass immunisation campaigns magnifies the need for research into the effectiveness and optimal deployment of pulse vaccination. In order to better understand this control strategy, we propose a mathematical model accounting for the disease dynamics in connected regions, incorporating seasonality, environmental reservoirs and independent periodic pulse vaccination schedules in each region. The effective reproduction number, $R_e$, is defined and proved to be a global threshold for persistence of the disease. Analytical and numerical calculations show the importance of synchronising the pulse vaccinations in connected regions and the timing of the pulses with respect to the pathogen circulation seasonality. Our results indicate that it may be crucial for mass-vaccination programs, such as national immunisation days, to be synchronised across different regions. In addition, simulations show that a migration imbalance can increase $R_e$ and alter how pulse vaccination should be optimally distributed among the patches, similar to results found with constant-rate vaccination. Furthermore, contrary to the case of constant-rate vaccination, the fraction of environmental transmission affects the value of $R_e$ when pulse vaccination is present.Comment: Added section 6.1, made other revisions, changed titl

Three Lectures: Nemd, Spam, and Shockwaves

We discuss three related subjects well suited to graduate research. The first, Nonequilibrium molecular dynamics or "NEMD", makes possible the simulation of atomistic systems driven by external fields, subject to dynamic constraints, and thermostated so as to yield stationary nonequilibrium states. The second subject, Smooth Particle Applied Mechanics or "SPAM", provides a particle method, resembling molecular dynamics, but designed to solve continuum problems. The numerical work is simplified because the SPAM particles obey ordinary, rather than partial, differential equations. The interpolation method used with SPAM is a powerful interpretive tool converting point particle variables to twice-differentiable field variables. This interpolation method is vital to the study and understanding of the third research topic we discuss, strong shockwaves in dense fluids. Such shockwaves exhibit stationary far-from-equilibrium states obtained with purely reversible Hamiltonian mechanics. The SPAM interpolation method, applied to this molecular dynamics problem, clearly demonstrates both the tensor character of kinetic temperature and the time-delayed response of stress and heat flux to the strain rate and temperature gradients. The dynamic Lyapunov instability of the shockwave problem can be analyzed in a variety of ways, both with and without symmetry in time. These three subjects suggest many topics suitable for graduate research in nonlinear nonequilibrium problems.Comment: 40 pages, with 21 figures, as presented at the Granada Seminar on the Foundations of Nonequilibrium Statistical Physics, 13-17 September, as three lecture

Stabilization of systems with asynchronous sensors and controllers

We study the stabilization of networked control systems with asynchronous sensors and controllers. Offsets between the sensor and controller clocks are unknown and modeled as parametric uncertainty. First we consider multi-input linear systems and provide a sufficient condition for the existence of linear time-invariant controllers that are capable of stabilizing the closed-loop system for every clock offset in a given range of admissible values. For first-order systems, we next obtain the maximum length of the offset range for which the system can be stabilized by a single controller. Finally, this bound is compared with the offset bounds that would be allowed if we restricted our attention to static output feedback controllers.Comment: 32 pages, 6 figures. This paper was partially presented at the 2015 American Control Conference, July 1-3, 2015, the US