Abstract Systems of activator–inhibitor reaction–diffusion equations posed on an infinite line are studied using a variety of analytical and numerical methods. A canonical form is considered, which contains all known models with simple cubic autocatalytic nonlinearity and arbitrary constant and linear kinetics. Restricting attention to models that have a unique homogeneous equilibrium, this class includes the classical Schnakenberg and Brusselator models, as well as other systems proposed in the literature to model morphogenesis. Such models are known to feature Turing instability, when activator diffuses more slowly than inhibitor, leading to stable spatially periodic patterns. Conversely in the limit of small feed rates, semi-strong interaction asymptotic analysis shows existence of isolated spike-like patterns. This paper describes the broad bifurcation structures that connect these two regimes. A certain universal two-parameter state diagram is revealed in which the Turing bifurcation becomes sub-critical, leading to the onset of homoclinic snaking. This regime then morphs into the spike regime, with the outer-fold being predicted by the semi-strong asymptotics. A rescaling of parameters and field concentrations shows how this state diagram can be studied independently of the diffusion rates. Temporal dynamics is found to strongly depend on the diffusion ratio though. A Hopf bifurcation occurs along the branch of stable spikes, which is subcritical for small diffusion ratio, leading to collapse to the homogeneous state. As the diffusion ratio increases, this bifurcation typically becomes supercritical and interacts with the homoclinic snaking and also with a supercritical homogeneous Hopf bifurcation, leading to complex spatio-temporal dynamics. The details are worked out for a number of different models that fit the theory using a mixture of weakly nonlinear analysis, semi-strong asymptotics and different numerical continuation algorithms.
more »
« less
Curvature effects and radial homoclinic snaking
Abstract In this work, we revisit some general results on the dynamics of circular fronts between homogeneous states and the formation of localized structures in two dimensions (2D). We show how the bifurcation diagram of axisymmetric structures localized in radius fits within the framework of collapsed homoclinic snaking. In 2D, owing to curvature effects, the collapse of the snaking structure follows a different scaling that is determined by the so-called nucleation radius. Moreover, in the case of fronts between two symmetry-related states, the precise point in parameter space to which radial snaking collapses is not a ‘Maxwell’ point but is determined by the curvature-driven dynamics only. In this case, the snaking collapses to a ‘zero surface tension’ point. Near this point, the breaking of symmetry between the homogeneous states tilts the snaking diagram. A different scaling law is found for the collapse of the snaking curve in each case. Curvature effects on axisymmetric localized states with internal structure are also discussed, as are cellular structures separated from a homogeneous state by a circular front. While some of these results are well understood in terms of curvature-driven dynamics and front interactions, a proper mathematical description in terms of homoclinic trajectories in a radial spatial dynamics description is lacking.
more »
« less
- Award ID(s):
- 1908891
- NSF-PAR ID:
- 10333626
- Date Published:
- Journal Name:
- IMA Journal of Applied Mathematics
- Volume:
- 86
- Issue:
- 5
- ISSN:
- 0272-4960
- Page Range / eLocation ID:
- 1094 to 1111
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract The passive conserved Swift–Hohenberg equation (or phase-field-crystal [PFC] model) describes gradient dynamics of a single-order parameter field related to density. It provides a simple microscopic description of the thermodynamic transition between liquid and crystalline states. In addition to spatially extended periodic structures, the model describes a large variety of steady spatially localized structures. In appropriate bifurcation diagrams the corresponding solution branches exhibit characteristic slanted homoclinic snaking. In an active PFC model, encoding for instance the active motion of self-propelled colloidal particles, the gradient dynamics structure is broken by a coupling between density and an additional polarization field. Then, resting and traveling localized states are found with transitions characterized by parity-breaking drift bifurcations. Here, we briefly review the snaking behavior of localized states in passive and active PFC models before discussing the bifurcation behavior of localized states in systems of (i) two coupled passive PFC models with common gradient dynamics, (ii) two coupled passive PFC models where the coupling breaks the gradient dynamics structure and (iii) a passive PFC model coupled to an active PFC model.more » « less
-
Abstract Localized coherent structures can form in externally driven dispersive optical cavities with a Kerr-type non-linearity. Such systems are described by the Lugiato–Lefever (LL) equation, which supports a large variety of dynamical states. Here, we review our current knowledge of the formation, stability and bifurcation structure of localized structures in the one-dimensional LL equation. We do so by focusing on two main regimes of operation: anomalous and normal second-order dispersion. In the anomalous regime, localized patterns are organized in a homoclinic snaking scenario, which is eventually destroyed, leading to a foliated snaking bifurcation structure. In the normal regime, localized structures undergo a different type of bifurcation structure, known as collapsed snaking. The effects of third-order dispersion and various dynamical regimes are also described.more » « less
-
null (Ed.)Abstract An activator–inhibitor–substrate model of side branching used in the context of pulmonary vascular and lung development is considered on the supposition that spatially localized concentrations of the activator trigger local side branching. The model consists of four coupled reaction–diffusion equations, and its steady localized solutions therefore obey an eight-dimensional spatial dynamical system in one spatial dimension (1D). Stationary localized structures within the model are found to be associated with a subcritical Turing instability and organized within a distinct type of foliated snaking bifurcation structure. This behavior is in turn associated with the presence of an exchange point in parameter space at which the complex leading spatial eigenvalues of the uniform concentration state are overtaken by a pair of real eigenvalues; this point plays the role of a Belyakov–Devaney point in this system. The primary foliated snaking structure consists of periodic spike or peak trains with $N$ identical equidistant peaks, $N=1,2,\dots \,$, together with cross-links consisting of nonidentical, nonequidistant peaks. The structure is complicated by a multitude of multipulse states, some of which are also computed, and spans the parameter range from the primary Turing bifurcation all the way to the fold of the $N=1$ state. These states form a complex template from which localized physical structures develop in the transverse direction in 2D.more » « less
-
Dense suspensions of particles in viscous liquid often demonstrate the striking phenomenon of abrupt shear thickening, where their viscosity increases strongly with increase of the imposed stress or shear rate. In this work, discrete-particle simulations accounting for short-range hydrodynamic, repulsive, and contact forces are performed to simulate flow of shear thickening bidisperse suspensions, with the packing parameters of large-to-small particle radius ratio δ = 3 and large particle fraction ζ = 0.15, 0.50, and 0.85. The simulations are carried out for volume fractions 0.54 ≤ ϕ ≤ 0.60 and a wide range of shear stresses. The repulsive forces, of magnitude F R , model the effects of surface charge and electric double-layer overlap, and result in shear thinning at small stress, with shear thickening beginning at stresses σ ∼ F R a −2 . A crossover scaling analysis used to describe systems with more than one thermodynamic critical point has recently been shown to successfully describe the experimentally-observed shear thickening behavior in suspensions. The scaling theory is tested here on simulated shear thickening data of the bidisperse mixtures, and also on nearly monodisperse suspensions with δ = 1.4 and ζ = 0.50. Presenting the viscosity in terms of a universal crossover scaling function between the frictionless and frictional maximum packing fractions collapses the viscosity for most of the suspensions studied. Two scaling regimes having different exponents are observed. The scaling analysis shows that the second normal stress difference N 2 and the particle pressure Π also collapse on their respective curves, with the latter featuring a different exponent from the viscosity and normal stress difference. The influence of the fraction of frictional contacts, one of the parameters of the scaling analysis, and its dependence on the packing parameters are also presented.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/imamat/hxab028
