Researcher profile

David J. W. Simpson

David J. W. Simpson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 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&#39; 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.

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.

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.

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&#39;s convention or with Coulomb&#39;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&#39;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

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.