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.

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 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.

Two-scale models pose a promising approach in simulating reactive flow and transport in evolving porous media. Classically, homogenised flow and transport equations are solved on the macroscopic scale, while effective parameters are obtained from auxiliary cell problems on possibly evolving reference geometries (micro-scale). Despite their perspective success in rendering lab/field-scale simulations computationally feasible, analytic results regarding the arising two-scale bilaterally coupled system often restrict to simplified models. In this paper, we first derive smooth dependence results concerning the partial coupling from the underlying geometry to macroscopic quantities. Therefore, alterations of the representative fluid domain are described by smooth paths of diffeomorphisms. Exploiting the gained regularity of the effective space- and time-dependent macroscopic coefficients, we present local-in-time existence results for strong solutions to the partially coupled micro–macro system using fixed-point arguments. What is more, we extend our results to the bilaterally coupled diffusive transport model including a level-set description of the evolving geometry.

The dynamics of interfaces in slow diffusion equations with strong absorption are studied. Asymptotic methods are used to give descriptions of the behaviour local to a comprehensive range of possible singular events that can occur in any evolution. These events are: when an interface changes its direction of propagation (reversing and anti-reversing), when an interface detaches from an absorbing obstacle (detaching), when two interfaces are formed by film rupture (touchdown) and when the solution undergoes extinction. Our account of extinction and self-similar reversing and anti-reversing is built upon previous work; results on non-self-similar reversing and anti-reversing and on the various types of detachment and touchdown are developed from scratch. In all cases, verification of the asymptotic results against numerical solutions to the full PDE is provided. Self-similar solutions, both of the full equation and of its asymptotic limits, play a central role in the analysis.

Stability is among the most important concepts in dynamical systems. Local stability is well-studied, whereas determining the ‘global stability’ of a nonlinear system is very challenging. Over the last few decades, many different ideas have been developed to address this issue, primarily driven by concrete applications. In particular, several disciplines suggested a web of concepts under the headline ‘resilience’. Unfortunately, there are many different variants and explanations of resilience, and often, the definitions are left relatively vague, sometimes even deliberately. Yet, to allow for a structural development of a mathematical theory of resilience that can be used across different areas, one has to ensure precise starting definitions and provide a mathematical comparison of different resilience measures. In this work, we provide a systematic review of the most relevant indicators of resilience in the context of continuous dynamical systems, grouped according to their mathematical features. The indicators are also generalised to be applicable to any attractor. These steps are important to ensure a more reliable, quantitatively comparable and reproducible study of resilience in dynamical systems. Furthermore, we also develop a new concept of resilience against certain nonautonomous perturbations to demonstrate how one can naturally extend our framework. All the indicators are finally compared via the analysis of a classic scalar model from population dynamics to show that direct quantitative application-based comparisons are an immediate consequence of a detailed mathematical analysis.

In this paper, we consider the time decay of the solutions to some problems arising in strain gradient thermoelasticity. We restrict to the two-dimensional case, and we assume that two dissipative mechanisms are introduced, the temperature and the mass dissipation. First, we show that this problem is well-posed proving that the operator defining it generates a contractive semigroup of linear operators. Then, assuming that the function involving the coupling terms is elliptic, the exponential decay of the solutions is concluded as well as the analyticity of the solutions. Finally, we describe how to obtain the exponential stability in the case of hyperbolic dissipation.

Motivated by applications in data science, we study partial differential equations on graphs. By a classical fixed-point argument, we show existence and uniqueness of solutions to a class of nonlocal continuity equations on graphs. We consider general interpolation functions, which give rise to a variety of different dynamics, for example, the nonlocal interaction dynamics coming from a solution-dependent velocity field. Our analysis reveals structural differences with the more standard Euclidean space, as some analogous properties rely on the interpolation chosen.

The time-global unique classical positive solutions to the reaction–diffusion equations for prey–predator models with dormancy of predators are constructed. The feature appears on the nonlinear terms of Holling type $\rm I\!I$ functional response. The crucial step is to establish time-local positive classical solutions by using a new approximation associated with time-evolution operators. Although the system does not equip usual comparison principle for solutions to partial differential equation, a priori bounds are derived by enclosing and renormalising arguments of solutions to the corresponding ordinary differential equations. Furthermore, time-global existence, invariant regions and asymptotic behaviours of solutions follow from such a priori bounds.

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.

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.

The problem of inverting the total divergence operator is central to finding components of a given conservation law. This might not be taxing for a low-order conservation law of a scalar partial differential equation, but integrable systems have conservation laws of arbitrarily high order that must be found with the aid of computer algebra. Even low-order conservation laws of complex systems can be hard to find and invert. This paper describes a new, efficient approach to the inversion problem. Two main tools are developed: partial Euler operators and partial scalings. These lead to a line integral formula for the inversion of a total derivative and a procedure for inverting a given total divergence concisely.

Evaporation within porous media is both a multiscale and interface-driven process, since the phase change at the evaporating interfaces within the pores generates a vapour flow and depends on the transport of vapour through the porous medium. While homogenised models of flow and chemical transport in porous media allow multiscale processes to be modelled efficiently, it is not clear how the multiscale effects impact the interface conditions required for these homogenised models. In this paper, we derive a homogenised model, including effective interface conditions, for the motion of an evaporation front through a porous medium, using a combined homogenisation and boundary layer analysis. This analysis extends previous work for a purely diffusive problem to include both gas flow and the advective–diffusive transport of material. We investigate the effect that different microscale models describing the chemistry of the evaporation have on the homogenised interface conditions. In particular, we identify a new effective parameter, $\mathcal{L}$, the average microscale interface length, which modifies the effective evaporation rate in the homogenised model. Like the effective diffusivity and permeability of a porous medium, $\mathcal{L}$ may be found by solving a periodic cell problem on the microscale. We also show that the different microscale models of the interface chemistry result in fundamentally different fine-scale behaviour at, and near, the interface.

The classic game of Battleship involves two players taking turns attempting to guess the positions of a fleet of vertically or horizontally positioned enemy ships hidden on a $10\times 10$ grid. One variant of this game, also referred to as Battleship Solitaire, Bimaru or Yubotu, considers the game with the inclusion of X-ray data, represented by knowledge of how many spots are occupied in each row and column in the enemy board. This paper considers the Battleship puzzle problem: the problem of reconstructing an enemy fleet from its X-ray data. We generate non-unique solutions to Battleship puzzles via certain reflection transformations akin to Ryser interchanges. Furthermore, we demonstrate that solutions of Battleship puzzles may be reliably obtained by searching for solutions of the associated classical binary discrete tomography problem which minimise the discrete Laplacian. We reformulate this optimisation problem as a quadratic unconstrained binary optimisation problem and approximate solutions via a simulated annealer, emphasising the future practical applicability of quantum annealers to solving discrete tomography problems with predefined structure.

The frameworks of Baikov–Gazizov–Ibragimov (BGI) and Fushchich–Shtelen (FS) approximate symmetries are used to study symmetry properties of partial differential equations with a small parameter. In general, it is shown that unlike the case of ordinary differential equations (ODEs), unstable BGI point symmetries of unperturbed partial differential equations (PDEs) do not necessarily yield local approximate symmetries for the perturbed model. While some relations between the BGI and FS approaches can be established, the two methods yield different approximate symmetry classifications. Detailed classifications are presented for two nonlinear PDE families. The second family includes a one-dimensional wave equation describing the wave motion in a hyperelastic material with a single family of fibers. For this model, approximate symmetries can be used to compute approximate closed-form solutions. Wave breaking times are found numerically and using the approximate solutions, which yield comparable results.

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.

In this contribution, we present a modelling and simulation framework for parametrised lithium-ion battery cells. We first derive a continuum model for a rather general intercalation battery cell on the basis of non-equilibrium thermodynamics. In order to efficiently evaluate the resulting parameterised non-linear system of partial differential equations, the reduced basis method is employed. The reduced basis method is a model order reduction technique on the basis of an incremental hierarchical approximate proper orthogonal decomposition approach and empirical operator interpolation. The modelling framework is particularly well suited to investigate and quantify degradation effects of battery cells. Several numerical experiments are given to demonstrate the scope and efficiency of the modelling framework.

Narrow escape and narrow capture problems which describe the average times required to stop the motion of a randomly travelling particle within a domain have applications in various areas of science. While for general domains, it is known how the escape time decreases with the increase of the trap sizes, for some specific 2D and 3D domains, higher-order asymptotic formulas have been established, providing the dependence of the escape time on the sizes and locations of the traps. Such results allow the use of global optimisation to seek trap arrangements that minimise average escape times. In a recent paper (Iyaniwura (2021) SIAM Rev. 63(3), 525–555), an explicit size- and trap location-dependent expansion of the average mean first passage time (MFPT) in a 2D elliptic domain was derived. The goal of this work is to systematically seek global minima of MFPT for $1\leq N\leq 50$ traps in elliptic domains using global optimisation techniques and compare the corresponding putative optimal trap arrangements for different values of the domain eccentricity. Further, an asymptotic formula for the average MFPT in elliptic domains with N circular traps of arbitrary sizes is derived, and sample optimal configurations involving non-equal traps are computed.

Symmetries and adjoint-symmetries are two fundamental (coordinate-free) structures of PDE systems. Recent work has developed several new algebraic aspects of adjoint-symmetries: three fundamental actions of symmetries on adjoint-symmetries; a Lie bracket on the set of adjoint-symmetries given by the range of a symmetry action; a generalised Noether (pre-symplectic) operator constructed from any non-variational adjoint-symmetry. These results are illustrated here by considering five examples of physically interesting nonlinear PDE systems – nonlinear reaction-diffusion equations, Navier-Stokes equations for compressible viscous fluid flow, surface-gravity water wave equations, coupled solitary wave equations and a nonlinear acoustic equation.