We study the planar FitzHugh–Nagumo system with an attracting periodic orbit that surrounds a repelling focus equilibrium. When the associated oscillation of the system is perturbed, in a given direction and with a given amplitude, there will generally be a change in phase of the perturbed oscillation with respect to the unperturbed one. This is recorded by the phase transition curve (PTC), which relates the old phase (along the periodic orbit) to the new phase (after perturbation). We take a geometric point of view and consider the phase-resetting surface comprising all PTCs as a function of the perturbation amplitude. This surface has a singularity when the perturbation maps a point on the periodic orbit exactly onto the repelling focus, which is the only point that does not return to stable oscillation. We also consider the PTC as a function of the direction of the perturbation and present how the corresponding phase-resetting surface changes with increasing perturbation amplitude. In this way, we provide a complete geometric interpretation of how the PTC changes for any perturbation direction. Unlike other examples discussed in the literature so far, the FitzHugh–Nagumo system is a generic example and, hence, representative for planar vector fields.

Macroscopically, a Darcian unsaturated moisture flow in the top soil is usually represented by an one-dimensional volume scale of evaporation from a static water table. On the microscale, simple pore-level models posit bundles of small-radius capillary tubes of a constant circular cross-section, fully occupied by mobile water moving in the Hagen–Poiseuille (HP) regime, while large-diameter pores are occupied by stagnant air. In our paper, cross-sections of cylindrical pores are polygonal. Steady, laminar, fully developed two-dimensional flows of Newtonian water in prismatic conduits, driven by a constant pressure gradient along a pore gradient, are more complex than the HP formula; this is based on the fact that the pores are only partially occupied by water and immobile air. The Poisson equation in a circular tetragon, with no-slip or mixed (no-shear-stress) boundary conditions on the two adjacent pore walls and two menisci, is solved by the methods of complex analysis. The velocity distribution is obtained via the Keldysh–Sedov type of singular integrals, and the flow rate is evaluated for several sets of meniscus radii by integrating the velocity over the corresponding tetragons.

We consider a pair of identical theta neurons in the active regime, each coupled to the other via a delayed Dirac delta function. The network can support periodic solutions and we concentrate on solutions for which the neurons are half a period out of phase with one another, and also solutions for which the neurons are perfectly synchronous. The dynamics are analytically solvable, so we can derive explicit expressions for the existence and stability of both types of solutions. We find two branches of solutions, connected by symmetry-broken solutions which arise when the period of a solution as a function of delay is at a maximum or a minimum.

In this paper, we derive simple analytical bounds for solutions of x - \ln x = y -\ln y$x - \ln x = y -\ln y$, and use them for estimating trajectories following Lotka–Volterra-type integrals. We show how our results give estimates for the Lambert W function as well as for trajectories of general predator–prey systems, including, for example, Rosenzweig–MacArthur equations.

This study explores the dynamics of a simple mechanical oscillator involving a magnet on a spring constrained to an axis; this magnet is additionally subject to the attractive force from a second magnet, which is placed on a parallel offset axis. The moments of both magnets remain aligned. The dynamics of the first magnet is first analysed in isolation for an unforced situation in which the second magnet is static and its position is taken as a parameter. We find codimension-1 saddle-node bifurcations, as well as a codimension-2 cusp bifurcation. The system has a region of bistability which increases in size with increasing force ratio. Next, the parametrically forced situation is considered, in which the second magnet moves sinusoidally. A comprehensive analysis of the forced oscillator behaviour is presented from the dynamical-systems standpoint. The solutions are shown to include periodic, quasiperiodic and chaotic trajectories. Resonances are shown to exist and the effect of weak damping is explored. Layered stroboscopic maps are used to produce cross-sections of the chaotic attractor as the parametric forcing frequency is varied. The strange attractor is found to disappear for a narrow window of forcing frequencies near the natural frequency of the spring.

While constructing mathematical models, scientists usually consider biotic factors, but it is crystal-clear that abiotic factors, such as wind, are also important as biotic factors. From this point of view, this paper is devoted to the investigation of some bifurcation properties of a fractional-order prey–predator model under the effect of wind. Using fractional calculus is very popular in modelling, since it is more effective than classical calculus in predicting the system’s future state and also discretization is one of the most powerful tools to study the behaviour of the models. In this paper, first of all, the model is discretized by using a piecewise discretization approach. Then, the local stability of fixed points is considered. We show using the centre manifold theorem and bifurcation theory that the system experiences a flip bifurcation and a Neimark–Sacker bifurcation at a positive fixed point. Finally, numerical simulations are given to demonstrate our results.

Recently, we analysed spontaneous symmetry breaking (SSB) of solitons in linearly coupled dual-core waveguides with fractional diffraction and cubic nonlinearity. In a practical context, the system can serve as a model for optical waveguides with the fractional diffraction or Bose–Einstein condensate of particles with Lévy index \alpha <2$\alpha <2$. In an earlier study, the SSB in the fractional coupler was identified as the bifurcation of subcritical type, becoming extremely subcritical in the limit of \alpha \rightarrow 1$\alpha \rightarrow 1$. There, the moving solitons and collisions between them at low speeds were also explored. In the present paper, we present new numerical results for fast solitons demonstrating restoration of symmetry in post-collision dynamics.

The study applies a two-dimensional adaptive mesh refinement (AMR) method to estimate the coordinates of the locations of the centre of vortices in steady, incompressible flow around a square cylinder placed within a channel. The AMR method is robust and low cost, and can be applied to any incompressible fluid flow. The considered channel has a blockage ratio of 1/8$1/8$. The AMR is tested on eight cases, considering flows with different Reynolds numbers (5\le Re\le 50$5\le Re\le 50$), and the estimated coordinates of the location of the centres of vortices are reported. For all test cases, the initial coarse meshes are refined four times, and the results are in good agreement with the literature where a very fine mesh was used. Furthermore, this study shows that the AMR method can capture the location of the centre of vortices within the fourth refined cells, and further confirms an improvement in the estimation with more refinements.

We derive and analyse well-posed boundary conditions for the linear shallow water wave equation. The analysis is based on the energy method and it identifies the number, location and form of the boundary conditions so that the initial boundary value problem is well-posed. A finite-volume method is developed based on the summation-by-parts framework with the boundary conditions implemented weakly using penalties. Stability is proven by deriving a discrete energy estimate analogous to the continuous estimate. The continuous and discrete analysis covers all flow regimes. Numerical experiments are presented verifying the analysis.

The behaviour of an axisymmetric bubble in a pure liquid forced by an acoustic pressure field is analysed. The bubble is assumed to have a sharp deformable interface, which is subject both to surface tension and to Rayleigh viscosity damping. Two modelling regimes are considered. The first is a linearized solution, based on the assumption of small axisymmetric deformations to an otherwise spherical bubble. The second involves a semi-numerical solution of the fully nonlinear problem, using a novel spectral method of high accuracy. For large-amplitude nonspherical bubble oscillations, the fully nonlinear solutions show that a complicated resonance structure is possible and that curvature singularities may occur at the interface, even in the presence of surface tension. Rayleigh viscosity at the interface prevents singularity formation, but eventually causes the bubble to become purely spherical unless shape-mode resonances occur. An extended analysis is also presented for purely spherical bubbles, which allows for a more detailed study of the effects of resonance and the Rayleigh viscosity at the bubble surface.

We propose a novel time-asymptotically stable, implicit–explicit, adaptive, time integration method (denoted by the \theta $\theta$-method) for the solution of the fractional advection–diffusion-reaction (FADR) equations. The spectral analysis of the method (involving the group velocity and the phase speed) indicates a region of favourable dispersion for a limited range of Péclet number. The numerical inversion of the coefficient matrix is avoided by exploiting the sparse structure of the matrix in the iterative solver for the Poisson equation. The accuracy and the efficacy of the method is benchmarked using (a) the two-dimensional fractional diffusion equation, originally proposed by researchers earlier, and (b) the incompressible, subdiffusive dynamics of a planar viscoelastic channel flow of the Rouse chain melts (FADR equation with fractional time-derivative of order ) and the Zimm chain solution (). Numerical simulations of the viscoelastic channel flow effectively capture the nonhomogeneous regions of high viscosity at low fluid inertia (or the so-called “spatiotemporal macrostructures”), experimentally observed in the flow-instability transition of subdiffusive flows.

We study the influence of a low-frequency harmonic vibration on the formation of the two-dimensional rolling solitary waves in vertically co-flowing two-layer liquid films. The system consists of two adjacent layers of immiscible fluids with the first layer being sandwiched between a vertical solid plate and the second fluid layer. The solid plate oscillates harmonically in the horizontal direction inducing Faraday waves at the liquid–liquid and liquid–air interfaces. We use a reduced hydrodynamic model derived from the Navier–Stokes equations in the long-wave approximation. Linear stability of the base flow in a flat two-layer film is determined semi-analytically using Floquet theory. We consider sub-millimetre-thick films and focus on the competition between the long-wavelength gravity-driven and finite wavelength Faraday instabilities. In the linear regime, the range of unstable wave vectors associated with the gravity-driven instability broadens at low and shrinks at high vibration frequencies. In nonlinear regimes, we find multiple metastable states characterized by solitary-like travelling waves and short pulsating waves. In particular, we find the range of the vibration parameters at which the system is multistable. In this regime, depending on the initial conditions, the long-time dynamics is dominated either by the fully developed solitary-like waves or by the shorter pulsating Faraday waves.

We devise schemes for producing, in the least possible time, p identical objects with n agents that work at differing speeds. This involves halting the process to transfer production across agent types. For the case of two types of agent, we construct schemes based on the Euclidean algorithm that seeks to minimize the number of pauses in production.

For microscale heterogeneous partial differential equations (PDEs), this article further develops novel theory and methodology for their macroscale mathematical/asymptotic homogenization. This article specifically encompasses the case of quasi-periodic heterogeneity with finite scale separation: no scale separation limit is required. A key innovation herein is to analyse the ensemble of all phase-shifts of the heterogeneity. Dynamical systems theory then frames the homogenization as a slow manifold of the ensemble. Depending upon any perceived scale separation within the quasi-periodic heterogeneity, the homogenization may be done in either one step or two sequential steps: the results are equivalent. The theory not only assures us of the existence and emergence of an exact homogenization at finite scale separation, it also provides a practical systematic method to construct the homogenization to any specified order. For a class of heterogeneities, we show that the macroscale homogenization is potentially valid down to lengths which are just twice that of the microscale heterogeneity! This methodology complements existing well-established results by providing a new rigorous and flexible approach to homogenization that potentially also provides correct macroscale initial and boundary conditions, treatment of forcing and control, and analysis of uncertainty.

Clustering is a method of allocating data points in various groups, known as clusters, based on similarity. The notion of expressing similarity mathematically and then maximizing it (minimize dissimilarity) can be formulated as an optimization problem. Spectral clustering is an example of such an approach to clustering, and it has been successfully applied to visualization of clustering and mapping of points into clusters in two and three dimensions. Higher dimension problems remained untouched due to complexity and, most importantly, lack of understanding what “similarity” means in higher dimensions. In this paper, we apply spectral clustering to long timeseries EEG (electroencephalogram) data. We developed several models, based on different similarity functions and different approaches for spectral clustering itself. The results of the numerical experiment demonstrate that the created models are accurate and can be used for timeseries classification.

We find solutions that describe the levelling of a thin fluid film, comprising a non-Newtonian power-law fluid, that coats a substrate and evolves under the influence of surface tension. We consider the evolution from periodic and localized initial conditions as separate cases. The particular (similarity) solutions in each of these two cases exhibit the generic property that the profiles are weakly singular (that is, higher-order derivatives do not exist) at points where the pressure gradient vanishes. Numerical simulations of the thin film equation, with either periodic or localized initial condition, are shown to approach the appropriate particular solution.

Many industrial design problems are characterized by a lack of an analytical expression defining the relationship between design variables and chosen quality metrics. Evaluating the quality of new designs is therefore restricted to running a predetermined process such as physical testing of prototypes. When these processes carry a high cost, choosing how to gather further data can be very challenging, whether the end goal is to accurately predict the quality of future designs or to find an optimal design. In the multi-fidelity setting, one or more approximations of a design’s performance are available at varying costs and accuracies. Surrogate modelling methods have long been applied to problems of this type, combining data from multiple sources into a model which guides further sampling. Many challenges still exist; however, the foremost among them is choosing when and how to rely on available low-fidelity sources. This tutorial-style paper presents an introduction to the field of surrogate modelling for multi-fidelity expensive black-box problems, including classical approaches and open questions in the field. An illustrative example using Australian elevation data is provided to show the potential downfalls in blindly trusting or ignoring low-fidelity sources, a question that has recently gained much interest in the community.

Singularly perturbed ordinary differential equations often exhibit Stokes’ phenomenon, which describes the appearance and disappearance of oscillating exponentially small terms across curves in the complex plane known as Stokes lines. These curves originate at singular points in the leading-order solution to the differential equation. In many important problems, it is impossible to obtain a closed-form expression for these leading-order solutions, and it is therefore challenging to locate these singular points. We present evidence that the analytic leading-order solution of a linear differential equation can be replaced with a numerical rational approximation using the adaptive Antoulas–Anderson (AAA) method. Despite such an approximation having completely different singularity types and locations, we show that the subsequent exponential asymptotic analysis accurately predicts the exponentially small behaviour present in the solution. For sufficiently small values of the asymptotic parameter, this approach breaks down; however, the range of validity may be extended by increasing the number of poles in the rational approximation. We present a related nonlinear problem and discuss the challenges that arise due to nonlinear effects. Overall, our approach allows for the study of exponentially small asymptotic effects without requiring an exact analytic form for the leading-order solution; this permits exponential asymptotic methods to be used in a much wider range of applications.

We consider planar flow involving two viscous fluids in a porous medium. One fluid is injected through a line source at the origin and moves radially outwards, pushing the second, ambient fluid outwards. There is an interface between the two fluids and if the inner injected fluid is of lower viscosity, the interface is unstable to small disturbances and radially directed unstable Saffman–Taylor fingers are produced. A linearized theory is presented and is compared with nonlinear results obtained using a numerical spectral method. An additional theory is also discussed, in which the sharp interface is replaced with a narrow diffuse interfacial region. We show that the nonlinear results are in close agreement with the linearized theory for small-amplitude disturbances at early times, but that large-amplitude fingers develop at later times and can even detach completely from the initial injection region.

Mathematical modelling of microwaves travelling through bauxite ore provides a way to compute moisture content in the free space transmission method given data on signal attenuation, phase shift and variable bauxite depth. We extend a recently developed four-layer model that uses coupled ordinary differential wave equations for the electric field together with continuity boundary conditions at interfaces between ore, air and antenna to find a solution that incorporates multiple internal reflections in ore and air. The model provides good fits to data, depending on ore permittivity and conductivity.

Our extensions are to use effective medium models to obtain electromagnetic properties of the ore mixture from moisture content and to incorporate the damping effects of scattering from the ore surface. Our model leads to a formula for the received signal showing how signal strengths SS and phase shifts depend on the moisture content of the bauxite ore, through the effects of moisture on permittivity and conductivity. We show that SS may be noninvertible, indicating that attenuation data alone cannot be used to infer moisture content. Combining with phase data typically corrects the noninvertibility. Reducing the operating frequency dramatically improves the usefulness of signal strength data for inferring moisture content.