Conventional preytaxis systems assume that prey-tactic velocity is proportional to the prey density gradient. However, many experiments exploring the predator–prey interactions show that it is the predator’s acceleration instead of velocity that is proportional to the prey density gradient in the prey-tactic movement, which we call preytaxis with prey-induced acceleration. Mathematical models of preytaxis with prey-induced acceleration were proposed by Arditi et al. ((2001) Theor. Popul. Biol. 59(3), 207–221) and Sapoukhina et al. ((2003) Am. Nat. 162(1), 61–76) to interpret the spatial heterogeneity of predators and prey observed in experiments. This paper is devoted to exploring the qualitative behaviour of such preytaxis systems with prey-induced acceleration and establishing the global existence of classical solutions with uniform-in-time bounds in all spatial dimensions. Moreover, we prove the global stability of spatially homogeneous prey-only and coexistence steady states with decay rates under certain conditions on system parameters. For the parameters outside the stability regime, we perform linear stability analysis to find the possible patterning regimes and use numerical simulations to demonstrate that spatially inhomogeneous time-periodic patterns will typically arise from the preytaxis system with prey-induced acceleration. Noticing that conventional preytaxis systems are unable to produce spatial patterns, our results imply that the preytaxis with prey-induced acceleration is indeed more appropriate than conventional preytaxis to interpret the spatial heterogeneity resulting from predator–prey interactions.

Mathematical models of polyelectrolyte gels are often simplified by assuming the gel is electrically neutral. The rationale behind this assumption is that the thickness of the electric double layer (EDL) at the free surface of the gel is small compared to the size of the gel. Hence, the thin-EDL limit is taken, in which the thickness of the EDL is set to zero. Despite the widespread use of the thin-EDL limit, the solutions in the EDL are rarely computed and shown to match to the solutions for the electrically neutral bulk. The aims of this paper are to study the structure of the EDL and establish the validity of the thin-EDL limit. The model for the gel accounts for phase separation, which gives rise to diffuse interfaces with a thickness described by the Kuhn length. We show that the solutions in the EDL can only be asymptotically matched to the solutions for an electrically neutral bulk, in general, when the Debye length is much smaller than the Kuhn length. If the Debye length is similar to or larger than the Kuhn length, then phase separation can be initiated in the EDL. This phase separation spreads into the bulk of the gel and gives rise to electrically charged layers with different degrees of swelling. Thus, the thin-EDL limit and the assumption of electroneutrality only generally apply when the Debye length is much smaller than the Kuhn length.

We propose an extension of the anisotropic interaction model which allows for collision avoidance in pairwise interactions by a rotation of forces (Totzeck (2020) Kinet. Relat. Models 13(6), 1219–1242.) by including the agents’ body size. The influence of the body size on the self-organisation of the agents in channel and crossing scenarios as well as the fundamental diagram is studied. Since the model is stated as a coupled system of ordinary differential equations, we are able to give a rigorous well-posedness analysis. Then we state a parameter calibration problem that involves data from real experiments. We prove the existence of a minimiser and derive the corresponding first-order optimality conditions. With the help of these conditions, we propose a gradient descent algorithm based on mini-batches of the data set. We employ the proposed algorithm to fit the parameter of the collision avoidance and the strength parameters of the interaction forces to given real data from experiments. The results underpin the feasibility of the method.

Sonar systems are frequently used to classify objects at a distance by using the structure of the echoes of acoustic waves as a proxy for the object’s shape and composition. Traditional synthetic aperture processing is highly effective in solving classification problems when the conditions are favourable but relies on accurate knowledge of the sensor’s trajectory relative to the object being measured. This article provides several new theoretical tools that decouple object classification performance from trajectory estimation in synthetic aperture sonar processing. The key insight is that decoupling the trajectory from classification-relevant information involves factoring a function into the composition of two functions. The article presents several new general topological invariants for smooth functions based on their factorisations over function composition. These invariants specialise to the case when a sonar platform trajectory is deformed by a non-small perturbation. The mathematical results exhibited in this article apply well beyond sonar classification problems. This article is written in a way that supports full mathematical generality.

We study the turnpike phenomenon for optimal control problems with mean-field dynamics that are obtained as the limit $N\rightarrow \infty$ of systems governed by a large number $N$ of ordinary differential equations. We show that the optimal control problems with large time horizons give rise to a turnpike structure of the optimal state and the optimal control. For the proof, we use the fact that the turnpike structure for the problems on the level of ordinary differential equations is preserved under the corresponding mean-field limit.

In this paper, we study consensus-based optimisation (CBO), a versatile, flexible and customisable optimisation method suitable for performing nonconvex and nonsmooth global optimisations in high dimensions. CBO is a multi-particle metaheuristic, which is effective in various applications and at the same time amenable to theoretical analysis thanks to its minimalistic design. The underlying dynamics, however, is flexible enough to incorporate different mechanisms widely used in evolutionary computation and machine learning, as we show by analysing a variant of CBO which makes use of memory effects and gradient information. We rigorously prove that this dynamics converges to a global minimiser of the objective function in mean-field law for a vast class of functions under minimal assumptions on the initialisation of the method. The proof in particular reveals how to leverage further, in some applications advantageous, forces in the dynamics without loosing provable global convergence. To demonstrate the benefit of the herein investigated memory effects and gradient information in certain applications, we present numerical evidence for the superiority of this CBO variant in applications such as machine learning and compressed sensing, which en passant widen the scope of applications of CBO.

The one-dimensional transient solidification of a binary alloy undergoing shrinkage is well-known as an invaluable benchmark for the testing of numerical codes that model macrosegregation. Here, recent work that considered the small-time behaviour of this problem is extended until complete solidification, thereby determining the solute profile across the entire solidified domain. The small-time solution is used as the initial condition for the numerical integration of a problem having three moving boundaries. Of particular significance is the so-called inverse segregation that is observed at the start of solidification, and the extreme segregation that is observed at the end; in the case of the example presented, which is for the often-cited Al–Cu system, the macrosegregation is found to be positive or negative, depending on whether Scheil’s equation or the lever rule is assumed at the microscale, respectively. The relevance of these results for the modelling of steady-state continuous casting processes – in particular, the phenomenon of centreline segregation – is also discussed.

In this paper, we investigate an initial-boundary value problem of a reaction–diffusion equation in a bounded domain with a Robin boundary condition and introduce some particular parameters to consider the non-zero flux on the boundary. This problem arises in the study of mosquito populations under the intervention of the population replacement method, where the boundary condition takes into account the inflow and outflow of individuals through the boundary. Using phase plane analysis, the present paper studies the existence and properties of non-constant steady-state solutions depending on several parameters. Then, we prove some sufficient conditions for their stability. We show that the long-time efficiency of this control method depends strongly on the size of the treated zone and the migration rate. To illustrate these theoretical results, we provide some numerical simulations in the framework of mosquito population control.

In this paper, we consider reaction-diffusion epidemic models with mass action or standard incidence mechanism and study the impact of limiting population movement on disease transmissions. We set either the dispersal rate of the susceptible or infected people to zero and study the corresponding degenerate reaction-diffusion model. Our main approach to study the global dynamics of these models is to construct delicate Lyapunov functions. Our results show that the consequences of limiting the movement of susceptible or infected people depend on transmission mechanisms, model parameters and population size.

In Chen and Liang previous work, a model, together with its well-posedness, was established for credit rating migrations with different upgrade and downgrade thresholds (i.e. a buffer zone, also called dead band in engineering). When positive dividends are introduced, the model in Chen and Liang (SIAM Financ. Math. 12, 941–966, 2021) may not be well-posed. Here, in this paper, a new model is proposed to include the realistic nonzero dividend scenarios. The key feature of the new model is that partial differential equations in Chen and Liang (SIAM Financ. Math. 12, 941–966, 2021) are replaced by variational inequalities, thereby creating a new free boundary, besides the original upgrading and downgrading free boundaries. Well-posedness of the new model, together with a few financially meaningful properties, is established. In particular, it is shown that when time to debt paying deadline is long enough, the underlying dividend paying company is always in high grade rating, that is, only when time to debt paying deadline is within a certain range, there can be seen the phenomenon of credit rating migration.

We study a system of nonlocal aggregation cross-diffusion PDEs that describe the evolution of opinion densities on a network. The PDEs are coupled with a system of ODEs that describe the time evolution of the agents on the network. Firstly, we apply the Deterministic Particle Approximation (DPA) method to the aforementioned system in order to prove the existence of solutions under suitable assumptions on the interactions between agents. Later on, we present an explicit model for opinion formation on an evolving network. The opinions evolve based on both the distance between the agents on the network and the ’attitude areas’, which depend on the distance between the agents’ opinions. The position of the agents on the network evolves based on the distance between the agents’ opinions. The goal is to study radicalisation, polarisation and fragmentation of the population while changing its open-mindedness and the radius of interaction.

We analyse a nonlinear partial differential equation system describing the motion of a microswimmer in a nematic liquid crystal environment. For the microswimmer’s motility, the squirmer model is used in which self-propulsion enters the model through the slip velocity on the microswimmer’s surface. The liquid crystal is described using the well-established Beris–Edwards formulation. In previous computational studies, it was shown that the squirmer, regardless of its initial configuration, eventually orients itself either parallel or perpendicular to the preferred orientation dictated by the liquid crystal. Furthermore, the corresponding solution of the coupled nonlinear system converges to a steady state. In this work, we rigorously establish the existence of steady state and also the finite-time existence for the time-dependent problem in a periodic domain. Finally, we will use a two-scale asymptotic expansion to derive a homogenised model for the collective swimming of squirmers as they reach their steady-state orientation and speed.

In this paper, we derive the effective model describing a thin-domain flow with permeable boundary through which the fluid is injected into the domain. We start with incompressible Stokes system and perform the rigorous asymptotic analysis. Choosing the appropriate scaling for the injection leads to a compressible effective model. In this paper, we derive the effective model describing a thin-domain flow with permeable boundary through which the fluid is injected into the domain. We start with incompressible Stokes system and perform the rigorous asymptotic analysis. Choosing the appropriate scaling for the injection leads to a compressible effective model.

This paper is devoted to the study of the propagation dynamics of a mutualistic model of mistletoes and birds with nonlocal dispersal. By applying the theory of asymptotic speeds of spread and travelling waves for monotone semiflows, we establish the existence of the asymptotic spreading speed $c^*$, the existence of travelling wavefronts with the wave speed $c\ge c^*$ and the nonexistence of travelling wavefronts with $c\lt c^*$. It turns out that the spreading speed coincides with the minimal wave speed of travelling wavefronts. Moreover, some lower and upper bound estimates of the spreading speed $c^*$ are provided.

In numerical linear stability investigations, the rates of change of the kinetic and thermal energy of the perturbation flow are often used to identify the dominant mechanisms by which kinetic or thermal energy is exchanged between the basic and the perturbation flow. Extending the conventional energy analysis for a single-phase Boussinesq fluid, the energy budgets of arbitrary infinitesimal perturbations to the basic two-phase liquid–gas flow are derived for an axisymmetric thermocapillary bridge when the material parameters in both phases depend on the temperature. This allows identifying individual transport terms and assessing their contributions to the instability if the basic flow and the critical mode are evaluated at criticality. The full closed-form energy budgets of linear modes have been derived for thermocapillary two-phase flow taking into account the temperature dependence of all thermophysical parameters. The influence of different approximations to the temperature dependence on the linear stability boundary of the axisymmetric flow in thermocapillary liquid bridges is tested regarding their accuracy. The general mechanism of symmetry breaking turns out to be very robust.

We consider a parabolic-parabolic chemotaxis system with singular chemotactic sensitivity and source functions, which is originally introduced by Short et al to model the spatio-temporal behaviour of urban criminal activities with the particular value of the chemotactic sensitivity parameter $\chi =2$. The available analytical findings for this urban crime model including $\chi =2$ are restricted either to one-dimensional setting, or to initial data and source functions with appropriate smallness, or to initial data and source functions with some radial symmetry. In the present work, our first result asserts that for any $\chi \gt 0$ the initial-boundary value problem of this urban crime model possesses a global generalised solution in the two-dimensional setting, without imposing any small or radial conditions on initial data and source functions. Our second result presents the asymptotic behaviour of such solution, under some additional assumptions on source functions.

Understanding the generation of mechanical stress in drying, particle-laden films is important for a wide range of industrial processes. One way to study these stresses is through the cantilever experiment, whereby a thin film is deposited onto the surface of a thin plate that is clamped at one end to a wall. The stresses that are generated in the film during drying are transmitted to the plate and drive bending. Mathematical modelling enables the film stress to be inferred from measurements of the plate deflection. The aim of this paper is to present simplified models of the cantilever experiment that have been derived from the time-dependent equations of continuum mechanics using asymptotic methods. The film is described using nonlinear poroelasticity and the plate using nonlinear elasticity. In contrast to Stoney-like formulae, the simplified models account for films with non-uniform thickness and stress. The film model reduces to a single differential equation that can be solved independently of the plate equations. The plate model reduces to an extended form of the Föppl-von Kármán (FvK) equations that accounts for gradients in the longitudinal traction acting on the plate surface. Consistent boundary conditions for the FvK equations are derived by resolving the Saint-Venant boundary layers at the free edges of the plate. The asymptotically reduced models are in excellent agreement with finite element solutions of the full governing equations. As the Péclet number increases, the time evolution of the plate deflection changes from $t$ to $t^{1/2}$, in agreement with experiments.

We examine several extensions to the basic Susceptible-Infected-Recovered model, which are able to induce recurrent outbreaks (the basic Susceptible-Infected-Recovered model by itself does not exhibit recurrent outbreaks). We first analyse how slow seasonal variations can destabilise the endemic equilibrium, leading to recurrent outbreaks. In the limit of slow immunity loss, we derive asymptotic thresholds that characterise this transition. In the outbreak regime, we use asymptotic matching to obtain a two-dimensional discrete map which describes outbreak times and strength. We then analyse the resulting map using linear stability and numerics. As the frequency of forcing is increased, the map exhibits a period-doubling route to chaos which alternates with periodic outbreaks of increasing frequency. Other extensions that can lead to recurrent outbreaks include the addition of noise, state-dependent variation and fine-graining of model classes.

Flowering plants depend on some animals for pollination and contribute to nourish the animals in natural environments. We call these animals pollinators and build a plants-pollinators cooperative model with impulsive effect on a periodically evolving domain. Next, we define the ecological reproduction index for single plant model and plants-pollinators system, respectively, whose threshold dynamics, including the extinction, persistence and coexistence, is established by the method of upper and lower solutions. Theoretical analysis shows that a large domain evolution rate has a positive influence on the survival of pollinators whether or not the impulsive effect occurs, and the pulse eliminates the pollinators even when the evolution rate is high. Moreover, some selective numerical simulations are still performed to explain our theoretical results.

The concurrency of edges, quantified by the number of edges that share a common node at a given time point, may be an important determinant of epidemic processes in temporal networks. We propose theoretically tractable Markovian temporal network models in which each edge flips between the active and inactive states in continuous time. The different models have different amounts of concurrency while we can tune the models to share the same statistics of edge activation and deactivation (and hence the fraction of time for which each edge is active) and the structure of the aggregate (i.e., static) network. We analytically calculate the amount of concurrency of edges sharing a node for each model. We then numerically study effects of concurrency on epidemic spreading in the stochastic susceptible-infectious-susceptible and susceptible-infectious-recovered dynamics on the proposed temporal network models. We find that the concurrency enhances epidemic spreading near the epidemic threshold, while this effect is small in many cases. Furthermore, when the infection rate is substantially larger than the epidemic threshold, the concurrency suppresses epidemic spreading in a majority of cases. In sum, our numerical simulations suggest that the impact of concurrency on enhancing epidemic spreading within our model is consistently present near the epidemic threshold but modest. The proposed temporal network models are expected to be useful for investigating effects of concurrency on various collective dynamics on networks including both infectious and other dynamics.