Source author record

David J. W. Simpson

David J. W. Simpson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Border-collision bifurcations from stable fixed points to any number of coexisting chaotic attractors

In diverse physical systems stable oscillatory solutions devolve into more complicated dynamical behaviour through border-collision bifurcations. Mathematically these occur when a stable fixed point of a piecewise-smooth map collides with a switching manifold as parameters are varied. The purpose of this paper is to highlight the extreme complexity possible in the subsequent dynamics. We perturb instances of the border-collision normal form in $n \ge 2$ dimensions for which the $n^{\rm th}$ iterate is a direct product of identical skew tent maps that have chaotic attractors comprised of $k \ge 2$ disjoint intervals. The resulting maps have coexisting attractors and we use Burnside's lemma to count the number of mutually disjoint trapping regions produced by taking unions of Cartesian products of slight enlargements of the disjoint intervals. The attractors are shown to be chaotic by demonstrating that some iterate of the map is piecewise-expanding. The resulting transition from a stable fixed point to many coexisting chaotic attractors is shown to occur throughout open subsets of parameter space and not destroyed by adding higher order terms to the normal form, hence can be expected to arise generically in mathematical models.

preprint2022arXiv

Normal forms, differentiable conjugacies and elementary bifurcations of maps

We strengthen the standard bifurcation theorems for saddle-node, transcritical, pitchfork, and period-doubling bifurcations of maps. Our new formulation involves adding one or two extra terms to the standard truncated normal forms with coefficients determined by algebraic equations. These extended normal forms are differentiably conjugate to the original maps on basins of attraction and repulsion of fixed points or periodic orbits. This reflects common assumptions about the additional information in normal forms despite standard bifurcation theorems being formulated only in terms of topological equivalence.

preprint2021arXiv

On the stability of boundary equilibria in Filippov systems

The leading-order approximation to a Filippov system $f$ about a generic boundary equilibrium $x^*$ is a system $F$ that is affine one side of the boundary and constant on the other side. We prove $x^*$ is exponentially stable for $f$ if and only if it is exponentially stable for $F$ when the constant component of $F$ is not tangent to the boundary. We then show exponential stability and asymptotic stability are in fact equivalent for $F$. We also show exponential stability is preserved under small perturbations to the pieces of $F$. Such results are well known for homogeneous systems. To prove the results here additional techniques are required because the two components of $F$ have different degrees of homogeneity. The primary function of the results is to reduce the problem of the stability of $x^*$ from the general Filippov system $f$ to the simpler system $F$. Yet in general this problem remains difficult. We provide a four-dimensional example of $F$ for which orbits appear to converge to $x^*$ in a chaotic fashion. By utilising the presence of both homogeneity and sliding motion the dynamics of $F$ can in this case be reduced to the combination of a one-dimensional return map and a scalar function.

preprint2021arXiv

Robust Devaney chaos in the two-dimensional border-collision normal form

The collection of all non-degenerate, continuous, two-piece, piecewise-linear maps on $\mathbb{R}^2$ can be reduced to a four-parameter family known as the two-dimensional border-collision normal form. We prove that throughout an open region of parameter space this family has an attractor satisfying Devaney's definition of chaos. This strengthens existing results on the robustness of chaos in piecewise-linear maps. We further show that the stable manifold of a saddle fixed point, despite being a one-dimensional object, densely fills an open region containing the attractor. Finally we identify a heteroclinic bifurcation, not described previously, at which the attractor undergoes a crisis and may be destroyed.

preprint2020arXiv

Numerical Bifurcation Analysis of Pacemaker Dynamics in a Model of Smooth Muscle Cells

Evidence from experimental studies shows that oscillations due to electro-mechanical coupling can be generated spontaneously in smooth muscle cells. Such cellular dynamics are known as \textit{pacemaker dynamics}. In this article we address pacemaker dynamics associated with the interaction of $\text{Ca}^{2+}$ and $\text{K}^+$ fluxes in the cell membrane of a smooth muscle cell. First we reduce a pacemaker model to a two-dimensional system equivalent to the reduced Morris-Lecar model and then perform a detailed numerical bifurcation analysis of the reduced model. Existing bifurcation analyses of the Morris-Lecar model concentrate on external applied current whereas we focus on parameters that model the response of the cell to changes in transmural pressure. We reveal a transition between Type I and Type II excitabilities with no external current required. We also compute a two-parameter bifurcation diagram and show how the transition is explained by the bifurcation structure.

preprint2019arXiv

Unfolding codimension-two subsumed homoclinic connections in two-dimensional piecewise-linear maps

For piecewise-linear maps, the phenomenon that a branch of a one-dimensional unstable manifold of a periodic solution is completely contained in its stable manifold is codimension-two. Unlike codimension-one homoclinic corners, such `subsumed' homoclinic connections can be associated with stable periodic solutions. The purpose of this paper is to determine the dynamics near a generic subsumed homoclinic connection in two dimensions. Assuming the eigenvalues associated with the periodic solution satisfy $0 < |λ| < 1 < σ< \frac{1}{|λ|}$, in a two-parameter unfolding there exists an infinite sequence of roughly triangular regions within which the map has a stable single-round periodic solution. The result applies to both discontinuous and continuous maps, although these cases admit different characterisations for the border-collision bifurcations that correspond to boundaries of the regions. The result is illustrated with a discontinuous map of Mira and the two-dimensional border-collision normal form.

preprint2016arXiv

Jitter in Dynamical Systems with Intersecting Discontinuity Surfaces

Mathematical models involving switches --- in the form of differential equations with discontinuities --- can accomodate real-world non-idealities through perturbations by hysteresis, time-delay, discretization, and noise. These are used to model the processes associated with switching in electronic control, mechanical contact, predator-prey preferences, and genetic or cellular regulation. The effect of such perturbations on rapid switching dynamics about a single switch are somewhat benign: in the limit that the size of the perturbation goes to zero the dynamics is given by Filippov's sliding solution. When multiple switches are involved, however, perturbations can have complicated effects, as shown in this paper. In the zero-perturbation limit, hysteresis, time-delay, and discretization cause erratic variation or `jitter' for stable sliding motion, whilst noise generates a relatively regular sliding solution similar to the canopy solution (an extension of Filippov's solution to multiple switches). We illustrate the results with a model of a switched power circuit, and showcase a variety of complex phenomena that perturbations can generate, including chaotic dynamics, exit selection, and coexisting sliding solutions.

preprint2016arXiv

The structure of mode-locking regions of piecewise-linear continuous maps: II. Skew sawtooth maps

In two-parameter bifurcation diagrams of piecewise-linear continuous maps on $\mathbb{R}^N$, mode-locking regions typically have points of zero width known as shrinking points. Near any shrinking point, but outside the associated mode-locking region, a significant proportion of parameter space can be usefully partitioned into a two-dimensional array of nearly-hyperbolic annular sectors. The purpose of this paper is to show that in these sectors the dynamics is well-approximated by a three-parameter family of skew sawtooth circle maps, where the relationship between the skew sawtooth maps and the $N$-dimensional map is fixed within each sector. The skew sawtooth maps are continuous, degree-one, and piecewise-linear, with two different slopes. They approximate the stable dynamics of the $N$-dimensional map with an error that goes to zero with the distance from the shrinking point. The results explain the complicated radial pattern of periodic, quasi-periodic, and chaotic dynamics that occurs near shrinking points.

preprint2016arXiv

Unfolding homoclinic connections formed by corner intersections in piecewise-smooth maps

The stable and unstable manifolds of an invariant set of a piecewise-smooth map are themselves piecewise-smooth. Consequently, as parameters of a piecewise-smooth map are varied, an invariant set can develop a homoclinic connection when its stable manifold intersects a non-differentiable point of its unstable manifold (or vice-versa). This is a codimension-one bifurcation analogous to a homoclinic tangency of a smooth map, referred to here as a homoclinic corner. This paper presents an unfolding of generic homoclinic corners for saddle fixed points of planar piecewise-smooth continuous maps. It is shown that a sequence of border-collision bifurcations limits to a homoclinic corner and that all nearby periodic solutions are unstable.

preprint2015arXiv

The influence of localised randomness on regular grazing bifurcations with applications to impacting dynamics

This paper concerns stochastic perturbations of piecewise-smooth ODE systems relevant for vibro-impacting dynamics, where impact events constitute the primary source of randomness. Such systems are characterised by the existence of switching manifolds that divide the phase space into regions where the system is smooth. The initiation of impacts is captured by a grazing bifurcation, at which a periodic orbit describing motion without impacts develops a tangential intersection with a switching manifold. Oscillatory dynamics near regular grazing bifurcations are described by piecewise-smooth maps involving a square-root singularity, known as Nordmark maps. We consider three scenarios where coloured noise only affects impacting dynamics, and derive three two-dimensional stochastic Nordmark maps with the noise appearing in different nonlinear or multiplicative ways, depending on the source of the noise. Consequently the stochastic dynamics differs between the three noise sources, and is fundamentally different to that of a Nordmark map with additive noise. This critical dependence on the nature of the noise is illustrated with a prototypical one-degree-of-freedom impact oscillator.

preprint2015arXiv

The structure of mode-locking regions of piecewise-linear continuous maps

The mode-locking regions of a dynamical system are the subsets of the parameter space of the system within which there exists an attracting periodic solution. For piecewise-linear continuous maps, these regions have a curious chain structure with points of zero width called shrinking points. In this paper we perform a local analysis about an arbitrary shrinking point. This is achieved by studying the symbolic itineraries of periodic solutions in nearby mode-locking regions and performing an asymptotic analysis on one-dimensional slow manifolds in order to build a comprehensive theoretical framework for the local dynamics. We obtain leading-order quantitative descriptions for the shape of nearby mode-locking regions, the location of nearby shrinking points, and the key properties of these shrinking points. We apply the results to the three-dimensional border-collision normal form, nonsmooth Neimark-Sacker-like bifurcations, and grazing-sliding bifurcations in a model of a dry friction oscillator.

preprint2014arXiv

On the relative coexistence of fixed points and period-two solutions near border-collision bifurcations

At a border-collision bifurcation a fixed point of a piecewise-smooth map intersects a surface where the functional form of the map changes. Near a generic border-collision bifurcation there are two fixed points, each of which exists on one side of the bifurcation. A simple eigenvalue condition indicates whether the fixed points exist on different sides of the bifurcation (this case can be interpreted as the persistence of a single fixed point), or on the same side of the bifurcation (in which case the bifurcation is akin to a saddle-node bifurcation). A similar eigenvalue condition indicates whether or not there exists a period-two solution on one side of the bifurcation. Previously these conditions have been combined to obtain five distinct scenarios for the existence and relative coexistence of fixed points and period-two solutions near border-collision bifurcations. In this Letter, it is shown that one of these scenarios, namely that two fixed points exist on one side of the bifurcation and a period-two solution exists on the other side of the bifurcation, cannot occur. The remaining four scenarios are feasible. Therefore there are exactly four distinct scenarios for fixed points and period-two solutions near border-collision bifurcations.

preprint2014arXiv

Scaling laws for large numbers of coexisting attracting periodic solutions in the border-collision normal form

A wide variety of intricate dynamics may be created at border-collision bifurcations of piecewise-smooth maps, where a fixed point collides with a surface at which the map is nonsmooth. For the border-collision normal form in two dimensions, a codimension-three scenario was described in previous work at which the map has a saddle-type periodic solution and an infinite sequence of stable periodic solutions that limit to a homoclinic orbit of the saddle-type solution. This paper introduces an alternate scenario of the same map at which there is an infinite sequence of stable periodic solutions due to the presence of a repeated unit eigenvalue in the linearization of some iterate of the map. It is shown that this scenario is codimension-four and that the sequence of periodic solutions is unbounded, aligning with eigenvectors corresponding to the unit eigenvalue. Arbitrarily many attracting periodic solutions coexist near either scenario. It is shown that if $K$ denotes the number of attracting periodic solutions, and $\varepsilon$ denotes the distance in parameter space from one of the two scenarios, then in the codimension-three case $\varepsilon$ scales with $λ^{-K}$, where $λ> 1$ denotes the unstable stability multiplier associated with the saddle-type periodic solution, and in the codimension-four case $\varepsilon$ scales with $K^{-2}$. Since $K^{-2}$ decays significantly slower than $λ^{-K}$, large numbers of attracting periodic solutions coexist in open regions of parameter space extending substantially further from the codimension-four scenarios than the codimension-three scenarios.

preprint2014arXiv

Stochastic Perturbations of Periodic Orbits with Sliding

Vector fields that are discontinuous on codimension-one surfaces are known as Filippov systems and can have attracting periodic orbits involving segments that are contained on a discontinuity surface of the vector field. In this paper we consider the addition of small noise to a general Filippov system and study the resulting stochastic dynamics near such a periodic orbit. Since a straight-forward asymptotic expansion in terms of the noise amplitude is not possible due to the presence of discontinuity surfaces, in order to quantitatively determine the basic statistical properties of the dynamics, we treat different parts of the periodic orbit separately. Dynamics distant from discontinuity surfaces is analyzed by the use of a series expansion of the transitional probability density function. Stochastically perturbed sliding motion is analyzed through stochastic averaging methods. The influence of noise on points at which the periodic orbit escapes a discontinuity surface is determined by zooming into the transition point. We combine the results to quantitatively determine the effect of noise on the oscillation time for a three-dimensional canonical model of relay control. For some parameter values of this model, small noise induces a significantly large reduction in the average oscillation time. By interpreting our results geometrically, we are able to identify four features of the relay control system that contribute to this phenomenon.

preprint2013arXiv

Non-Filippov dynamics arising from the smoothing of nonsmooth systems, and its robustness to noise

Switch-like behaviour in dynamical systems may be modelled by highly nonlinear functions, such as Hill functions or sigmoid functions, or alternatively by piecewise-smooth functions, such as step functions. Consistent modelling requires that piecewise-smooth and smooth dynamical systems have similar dynamics, but the conditions for such similarity are not well understood. Here we show that by smoothing out a piecewise-smooth system one may obtain dynamics that is inconsistent with the accepted wisdom --- so-called Filippov dynamics --- at a discontinuity, even in the piecewise-smooth limit. By subjecting the system to white noise, we show that these discrepancies can be understood in terms of potential wells that allow solutions to dwell at the discontinuity for long times. Moreover we show that spurious dynamics will revert to Filippov dynamics, with a small degree of stochasticity, when the noise magnitude is sufficiently large compared to the order of smoothing. We apply the results to a model of a dry-friction oscillator, where spurious dynamics (inconsistent with Filippov's convention or with Coulomb's model of friction) can account for different coefficients of static and kinetic friction, but under sufficient noise the system reverts to dynamics consistent with Filippov's convention (and with Coulomb-like friction).

preprint2013arXiv

On resolving singularities of piecewise-smooth discontinuous vector fields via small perturbations

A two-fold singularity is a point on a discontinuity surface of a piecewise-smooth vector field at which the vector field is tangent to the surface on both sides. Due to the double tangency, forward evolution from a two-fold is typically ambiguous. This is an especially serious issue for two-folds that are reached by the forward orbits of a non-zero measure set of initial points. The purpose of this paper is to explore the concept of perturbing the vector field so that forward evolution is well-defined, and characterising the perturbed dynamics in the limit that the size of the perturbation tends to zero. This concept is applied to a two-fold in two dimensions. Three forms of perturbation: hysteresis, time-delay, and noise, are analysed individually. In each case, the limit leads to a novel probabilistic notion of forward evolution from the two-fold.

preprint2013arXiv

Sequences of Periodic Solutions and Infinitely Many Coexisting Attractors in the Border-Collision Normal Form

The border-collision normal form is a piecewise-linear continuous map on $\mathbb{R}^N$ that describes dynamics near border-collision bifurcations of nonsmooth maps. This paper studies a codimension-three scenario at which the border-collision normal form with $N=2$ exhibits infinitely many attracting periodic solutions. In this scenario there is a saddle-type periodic solution with branches of stable and unstable manifolds that are coincident, and an infinite sequence of attracting periodic solutions that converges to an orbit homoclinic to the saddle-type solution. Several important features of the scenario are shown to be universal, and three examples are given. For one of these examples infinite coexistence is proved directly by explicitly computing periodic solutions in the infinite sequence.

preprint2013arXiv

The Positive Occupation Time of Brownian Motion with Two-Valued Drift and Asymptotic Dynamics of Sliding Motion with Noise

We derive the probability density function of the positive occupation time of one-dimensional Brownian motion with two-valued drift. Long time asymptotics of the density are also computed. We use the result to describe the transitional probability density function of a general N-dimensional system of stochastic differential equations representing stochastically perturbed sliding motion of a discontinuous, piecewise-smooth vector field on short time frames. A description of the density at larger times is obtained via an asymptotic expansion of the Fokker-Planck equation.

preprint2012arXiv

Stochastic Regular Grazing Bifurcations

A grazing bifurcation corresponds to the collision of a periodic orbit with a switching manifold in a piecewise-smooth ODE system and often generates complicated dynamics. The lowest order terms of the induced Poincare map expanded about a regular grazing bifurcation constitute a Nordmark map. In this paper we study a normal form of the Nordmark map in two dimensions with additive Gaussian noise of amplitude, epsilson [e]. We show that this particular noise formulation arises in a general setting and consider a harmonically forced linear oscillator subject to compliant impacts to illustrate the accuracy of the map. Numerically computed invariant densities of the stochastic Nordmark map can take highly irregular forms, or, if there exists an attracting period-n solution when e = 0, be well approximated by the sum of n Gaussian densities centred about each point of the deterministic solution, and scaled by 1/n, for sufficiently small e > 0. We explain the irregular forms and calculate the covariance matrices associated with the Gaussian approximations in terms of the parameters of the map. Close to the grazing bifurcation the size of the invariant density may be proportional to the square-root of e, as a consequence of a square-root singularity in the map. Sequences of transitions between different dynamical regimes that occur as the primary bifurcation parameter is varied have not been described previously.

preprint2012arXiv

Stochastically Perturbed Sliding Motion in Piecewise-Smooth Systems

Sliding motion is evolution on a switching manifold of a discontinuous, piecewise-smooth system of ordinary differential equations. In this paper we quantitatively study the effects of small-amplitude, additive, white Gaussian noise on stable sliding motion. For equations that are static in directions parallel to the switching manifold, the distance of orbits from the switching manifold approaches a quasi-steady-state density. From this density we calculate the mean and variance for the near sliding solution. Numerical results of a relay control system reveal that the noise may significantly affect the period and amplitude of periodic solutions with sliding segments.

preprint2011arXiv

Dynamics of Simple Balancing Models with State Dependent Switching Control

Time-delayed control in a balancing problem may be a nonsmooth function for a variety of reasons. In this paper we study a simple model of the control of an inverted pendulum by either a connected movable cart or an applied torque for which the control is turned off when the pendulum is located within certain regions of phase space. Without applying a small angle approximation for deviations about the vertical position, we see structurally stable periodic orbits which may be attracting or repelling. Due to the nonsmooth nature of the control, these periodic orbits are born in various discontinuity-induced bifurcations. Also we show that a coincidence of switching events can produce complicated periodic and aperiodic solutions.

preprint2010arXiv

Mixed-Mode Oscillations in a Stochastic, Piecewise-Linear System

We analyze a piecewise-linear FitzHugh-Nagumo model. The system exhibits a canard near which both small amplitude and large amplitude periodic orbits exist. The addition of small noise induces mixed-mode oscillations (MMOs) in the vicinity of the canard point. We determine the effect of each model parameter on the stochastically driven MMOs. In particular we show that any parameter variation (such as a modification of the piecewise-linear function in the model) that leaves the ratio of noise amplitude to time-scale separation unchanged typically has little effect on the width of the interval of the primary bifurcation parameter over which MMOs occur. In that sense, the MMOs are robust. Furthermore we show that the piecewise-linear model exhibits MMOs more readily than the classical FitzHugh-Nagumo model for which a cubic polynomial is the only nonlinearity. By studying a piecewise-linear model we are able to explain results using analytical expressions and compare these with numerical investigations.

preprint2009arXiv

Simultaneous Border-Collision and Period-Doubling Bifurcations

We unfold the codimension-two simultaneous occurrence of a border-collision bifurcation and a period-doubling bifurcation for a general piecewise-smooth, continuous map. We find that, with sufficient non-degeneracy conditions, a locus of period-doubling bifurcations emanates non-tangentially from a locus of border-collision bifurcations. The corresponding period-doubled solution undergoes a border-collision bifurcation along a curve emanating from the codimension-two point and tangent to the period-doubling locus here. In the case that the map is one-dimensional local dynamics are completely classified; in particular, we give conditions that ensure chaos.