Suspended fibres significantly alter fluid rheology, as exhibited in for example solutions of DNA, RNA and synthetic biological nanofibres. It is of interest to determine how this altered rheology affects flow stability. Motivated by the fact thermal gradients may occur in biomolecular analytic devices, and recent stability results, we examine the problem of Rayleigh–Bénard convection of the transversely isotropic fluid of Ericksen. A transversely isotropic fluid treats these suspensions as a continuum with an evolving preferred direction, through a modified stress tensor incorporating four viscosity-like parameters. We consider the linear stability of a stationary, passive, transversely isotropic fluid contained between two parallel boundaries, with the lower boundary at a higher temperature than the upper. To determine the marginal stability curves the Chebyshev collocation method is applied, and we consider a range of initially uniform preferred directions, from horizontal to vertical, and three orders of magnitude in the viscosity-like anisotropic parameters. Determining the critical wave and Rayleigh numbers, we find that transversely isotropic effects delay the onset of instability; this effect is felt most strongly through the incorporation of the anisotropic shear viscosity, although the anisotropic extensional viscosity also contributes. Our analysis confirms the importance of anisotropic rheology in the setting of convection.

Increasing evidence suggests that the presence of mobile ions in perovskite solar cells (PSCs) can cause a current–voltage curve hysteresis. Steady state and transient current–voltage characteristics of a planar metal halide CH3NH3PbI3 PSC are analysed with a drift-diffusion model that accounts for both charge transport and ion vacancy motion. The high ion vacancy density within the perovskite layer gives rise to narrow Debye layers (typical width ~2 nm), adjacent to the interfaces with the transport layers, over which large drops in the electric potential occur and in which significant charge is stored. Large disparities between (I) the width of the Debye layers and that of the perovskite layer (~600 nm) and (II) the ion vacancy density and the charge carrier densities motivate an asymptotic approach to solving the model, while the stiffness of the equations renders standard solution methods unreliable. We derive a simplified surface polarisation model in which the slow ion dynamics are replaced by interfacial (non-linear) capacitances at the perovskite interfaces. Favourable comparison is made between the results of the asymptotic approach and numerical solutions for a realistic cell over a wide range of operating conditions of practical interest.

Neural mass models are ubiquitous in large-scale brain modelling. At the node level, they are written in terms of a set of ordinary differential equations with a non-linearity that is typically a sigmoidal shape. Using structural data from brain atlases, they may be connected into a network to investigate the emergence of functional dynamic states, such as synchrony. With the simple restriction of the classic sigmoidal non-linearity to a piecewise linear caricature, we show that the famous Wilson–Cowan neural mass model can be explicitly analysed at both the node and network level. The construction of periodic orbits at the node level is achieved by patching together matrix exponential solutions, and stability is determined using Floquet theory. For networks with interactions described by circulant matrices, we show that the stability of the synchronous state can be determined in terms of a low-dimensional Floquet problem parameterised by the eigenvalues of the interaction matrix. Moreover, this network Floquet problem is readily solved using linear algebra to predict the onset of spatio-temporal network patterns arising from a synchronous instability. We further consider the case of a discontinuous choice for the node non-linearity, namely the replacement of the sigmoid by a Heaviside non-linearity. This gives rise to a continuous-time switching network. At the node level, this allows for the existence of unstable sliding periodic orbits, which we explicitly construct. The stability of a periodic orbit is now treated with a modification of Floquet theory to treat the evolution of small perturbations through switching manifolds via the use of saltation matrices. At the network level, the stability analysis of the synchronous state is considerably more challenging. Here, we report on the use of ideas originally developed for the study of Glass networks to treat the stability of periodic network states in neural mass models with discontinuous interactions.

The displacement λ-convexity of a non-standard entropy with respect to a non-local transportation metric in finite state spaces is shown using a gradient flow approach. The constant λ is computed explicitly in terms of a priori estimates of the solution to a finite-difference approximation of a non-linear Fokker–Planck equation. The key idea is to employ a new mean function, which defines the Onsager operator in the gradient flow formulation.

We consider the Gierer–Meinhardt system with precursor inhomogeneity and two small diffusivities in an interval

$$\begin{equation*} \left\{ \begin{array}{ll} A_t=\epsilon^2 A''- \mu(x) A+\frac{A^2}{H}, &x\in(-1, 1),\,t>0,\\[3mm] \tau H_t=D H'' -H+ A^2, & x\in (-1, 1),\,t>0,\\[3mm] A' (-1)= A' (1)= H' (-1) = H' (1) =0, \end{array} \right. \end{equation*}$$

$$\begin{equation*}\mbox{where } \quad 0<\epsilon \ll\sqrt{D}\ll 1, \quad \end{equation*}$$

$$\begin{equation*} \tau\geq 0 \mbox{ and $\tau$ is independent of $\epsilon$. } \end{equation*}$$

$O(\sqrt{D})$

Entropy maximising spatial interaction models have been widely exploited in a range of disciplines and applications: from trade and migration flows to the spread of riots and the understanding of spatial patterns in archaeological sites of interest. When embedded into a dynamic system and framed in the context of a retail model, the dynamics of centre growth poses an interesting mathematical problem, with bifurcations and phase changes, which may be addressed analytically. In this paper, we present some analysis of the continuous retail model and the corresponding discrete version, which yields insights into the effect of space on the evolving system, and an understanding of why certain retail centres are more successful than others. The slowly developing growths and the fast explosive growths that are of particular concern are explained in detail.

Using mathematical methods to understand and model crime is a recent idea that has drawn considerable attention from researchers during the last five years. From the plethora of models that have been proposed, perhaps the most successful one has been a diffusion-type differential equations model that describes how the number of criminals evolves in a specific area. We propose a more detailed form of this model that allows for two distinct criminal types associated with major and minor crime. Additionally, we examine a stochastic variant of the model that represents more realistically the ‘generation’ of new criminals. Numerical solutions from both models are presented and compared with actual crime data for the Greater Manchester area. Agreement between simulations and actual data is satisfactory. A preliminary statistical analysis of the data also supports the model's potential to describe crime.

The master stability function is a powerful tool for determining synchrony in high-dimensional networks of coupled limit cycle oscillators. In part, this approach relies on the analysis of a low-dimensional variational equation around a periodic orbit. For smooth dynamical systems, this orbit is not generically available in closed form. However, many models in physics, engineering and biology admit to non-smooth piece-wise linear caricatures, for which it is possible to construct periodic orbits without recourse to numerical evolution of trajectories. A classic example is the McKean model of an excitable system that has been extensively studied in the mathematical neuroscience community. Understandably, the master stability function cannot be immediately applied to networks of such non-smooth elements. Here, we show how to extend the master stability function to non-smooth planar piece-wise linear systems, and in the process demonstrate that considerable insight into network dynamics can be obtained. In illustration, we highlight an inverse period-doubling route to synchrony, under variation in coupling strength, in globally linearly coupled networks for which the node dynamics is poised near a homoclinic bifurcation. Moreover, for a star graph, we establish a mechanism for achieving so-called remote synchronisation (where the hub oscillator does not synchronise with the rest of the network), even when all the oscillators are identical. We contrast this with node dynamics close to a non-smooth Andronov–Hopf bifurcation and also a saddle node bifurcation of limit cycles, for which no such bifurcation of synchrony occurs.

We consider interacting particle systems and their mean-field limits, which are frequently used to model collective aggregation and are known to demonstrate a rich variety of pattern formations. The interaction is based on a pairwise potential combining short-range repulsion and long-range attraction. We study particular solutions, which are referred to as flocks in the second-order models, for the specific choice of the Quasi-Morse interaction potential. Our main result is a rigorous analysis of continuous, compactly supported flock profiles for the biologically relevant parameter regime. Existence and uniqueness are proven for three space dimensions, while existence is shown for the two-dimensional case. Furthermore, we numerically investigate additional Morse-like interactions to complete the understanding of this class of potentials.

We study of a new type of multi-bump blowup solutions of the Ginzburg–Landau equation. Multi-bump blowup solutions have previously been found in numeric simulations, asymptotic analysis and were proved to exist via geometric construction. In the geometric construction of the solutions, the existence of two types of multi-bump solutions was shown. One type is exponentially small at ξ=0, and the other type of solutions is algebraically small at ξ=0. So far, the first type of solutions were studied asymptotically. Here, we analyse the solutions which are algebraically small at ξ=0 by using asymptotic methods. This construction is essentially different from the existing one, and ideas are obtained from the geometric construction. Hence, this is a good example of where both asymptotic analysis and geometric methods are needed for the overall picture.