This paper presents the current state of mathematical modelling of the electrochemical behaviour of lithium-ion batteries (LIBs) as they are charged and discharged. It reviews the models developed by Newman and co-workers, both in the cases of dilute and moderately concentrated electrolytes and indicates the modelling assumptions required for their development. Particular attention is paid to the interface conditions imposed between the electrolyte and the active electrode material; necessary conditions are derived for one of these, the Butler–Volmer relation, in order to ensure physically realistic solutions. Insight into the origin of the differences between various models found in the literature is revealed by considering formulations obtained by using different measures of the electric potential. Materials commonly used for electrodes in LIBs are considered and the various mathematical models used to describe lithium transport in them discussed. The problem of upscaling from models of behaviour at the single electrode particle scale to the cell scale is addressed using homogenisation techniques resulting in the pseudo-2D model commonly used to describe charge transport and discharge behaviour in lithium-ion cells. Numerical solution to this model is discussed and illustrative results for a common device are computed.

This paper investigates global dynamics of an infection age-space structured cholera model. The model describes the vibrio cholerae transmission in human population, where infection-age structure of vibrio cholerae and infectious individuals are incorporated to measure the infectivity during the different stage of disease transmission. The model is described by reaction–diffusion models involving the spatial dispersal of vibrios and the mobility of human populations in the same domain Ω ⊂ ℝn. We first give the well-posedness of the model by converting the model to a reaction–diffusion model and two Volterra integral equations and obtain two constant equilibria. Our result suggest that the basic reproduction number determines the dichotomy of disease persistence and extinction, which is achieved by studying the local stability of equilibria, disease persistence and global attractivity of equilibria.

The study of spherically symmetric motion is important for the theory of explosion waves. In this paper, we consider a ‘spherical piston’ problem for the relativistic Euler equations, which describes the wave motion produced by a sphere expanding into an infinite surrounding medium. We use the reflected characteristics method to construct a global piecewise smooth solution with a single shock of this spherical piston problem, provided that the speed of the sphere is a small perturbation of a constant speed.

A non-linear reaction–diffusion system with cross-diffusion describing the COVID-19 outbreak is studied using the Lie symmetry method. A complete Lie symmetry classification is derived and it is shown that the system with correctly specified parameters admits highly non-trivial Lie symmetry operators, which do not occur for all known reaction–diffusion systems. The symmetries obtained are also applied for finding exact solutions of the system in the most interesting case from applicability point of view. It is shown that the exact solutions derived possess typical properties for describing the pandemic spread under 1D approximation in space and lead to the distributions, which qualitatively correspond to the measured data of the COVID-19 spread in Ukraine.

The existence of a symmetric mode in an elastic solid wedge for all admissible values of the Poisson ratio and arbitrary interior angles close to π has been proven.

This paper studies the periodic trajectories of a novel age-structured prey–predator system with Michaelis–Menten functional response including delays and asymmetric diffusion. To begin with, the system is turned into an abstract non-densely defined Cauchy problem, and a time-lag effect in their interaction is investigated. Next, we acquire that this system appears a periodic orbit near the positive steady state by employing the method of integrated semigroup and the Hopf bifurcation theory for semilinear equations with non-dense domain, which is also the main result of this article. Finally, in order to illustrate our theoretical analysis more vividly, we make some numerical simulations and give some discussions.

When an explosive burns, gaseous products are formed as a result. The interaction of the burning solid and gas is not well understood. More specifically, the process of the gaseous product heating the explosive is yet to be explored in detail. The present work sets out to fill some of that gap using mathematical modelling: this aims to track the temperature profile in the explosive. The work begins by modelling single-step reactions using a simple Arrhenius model. The model is then extended to include three-step reaction. An alternative asymptotic approach is also employed. There is close agreement between results from the full reaction-diffusion problem and the asymptotic problem.

We study the Muskat problem describing the vertical motion of two immiscible fluids in a two-dimensional homogeneous porous medium in an Lp-setting with p ∈ (1, ∞). The Sobolev space $W_p^s(\mathbb R)$ with s = 1+1/p is a critical space for this problem. We prove, for each s ∈ (1+1/p, 2) that the Rayleigh–Taylor condition identifies an open subset of $W_p^s(\mathbb R)$ within which the Muskat problem is of parabolic type. This enables us to establish the local well-posedness of the problem in all these subcritical spaces together with a parabolic smoothing property.

In this paper, we study a mathematical model for an infectious disease caused by a virus such as Cholera without lifetime immunity. Due to the different mobility for susceptible, infected human and recovered human hosts, the diffusion coefficients are assumed to be different. The resulting system is governed by a strongly coupled reaction–diffusion system with different diffusion coefficients. Global existence and uniqueness are established under certain assumptions on known data. Moreover, global asymptotic behaviour of the solution is obtained when some parameters satisfy certain conditions. These results extend the existing results in the literature. The main tool used in this paper comes from the delicate theory of elliptic and parabolic equations. Moreover, the energy method and Sobolev embedding are used in deriving a priori estimates. The analysis developed in this paper can be employed to study other epidemic models in biological, ecological and health sciences.

We study a mathematical model proposed in the literature with the aim of describing the interactions between tumor cells and the immune system, when a periodic treatment of immunotherapy is applied. Combining some techniques from non-linear analysis (degree theory, lower and upper solutions, and theory of free-homeomorphisms in the plane), we give a detailed global analysis of the model. We also observe that for certain therapies, the maximum level of aggressiveness of a cancer, for which the treatment works (or does not work), can be computed explicitly. We discuss some strategies for designing therapies. The mathematical analysis is completed with numerical results and conclusions.

The self-interaction force of dislocation curves in metals depends on the local arrangement of the atoms and on the non-local interaction between dislocation curve segments. While these non-local segment–segment interactions can be accurately described by linear elasticity when the segments are further apart than the atomic scale of size $\varepsilon$ , this model breaks down and blows up when the segments are $O(\varepsilon)$ apart. To separate the non-local interactions from the local contribution, various models depending on $\varepsilon$ have been constructed to account for the non-local term. However, there are no quantitative comparisons available between these models. This paper makes such comparisons possible by expanding the self-interaction force in these models in $\varepsilon$ beyond the O(1)-term. Our derivation of these expansions relies on asymptotic analysis. The practical use of these expansions is demonstrated by developing numerical schemes for them, and by – for the first time – bounding the corresponding discretisation error.

An emerging technique in image segmentation, semi-supervised learning and general classification problems concerns the use of phase-separating flows defined on finite graphs. This technique was pioneered in Bertozzi and Flenner (2012, Multiscale Modeling and Simulation10(3), 1090–1118), which used the Allen–Cahn flow on a graph, and was then extended in Merkurjev et al. (2013, SIAM J. Imaging Sci.6(4), 1903–1930) using instead the Merriman–Bence–Osher (MBO) scheme on a graph. In previous work by the authors, Budd and Van Gennip (2020, SIAM J. Math. Anal.52(5), 4101–4139), we gave a theoretical justification for this use of the MBO scheme in place of Allen–Cahn flow, showing that the MBO scheme is a special case of a ‘semi-discrete’ numerical scheme for Allen–Cahn flow. In this paper, we extend this earlier work, showing that this link via the semi-discrete scheme is robust to passing to the mass-conserving case. Inspired by Rubinstein and Sternberg (1992, IMA J. Appl. Math.48, 249–264), we define a mass-conserving Allen–Cahn equation on a graph. Then, with the help of the tools of convex optimisation, we show that our earlier machinery can be applied to derive the mass-conserving MBO scheme on a graph as a special case of a semi-discrete scheme for mass-conserving Allen–Cahn. We give a theoretical analysis of this flow and scheme, proving various desired properties like existence and uniqueness of the flow and convergence of the scheme, and also show that the semi-discrete scheme yields a choice function for solutions to the mass-conserving MBO scheme.

Silicon is produced in submerged arc furnaces which are heated by electric currents passing through the furnace. It is important to understand the distribution of heating within the furnace in order to accurately model the silicon production process, yet many existing studies neglect aspects of this current flow. In the present paper, we formulate a model that couples the electrical current to thermal, material flow and chemical processes in the furnace. We then exploit disparate timescales to homogenise the model over the timescale of the alternating current, deriving averaged equations for the slow evolution of the system. Our numerical simulations predict a minimum applied current that is required in order to obtain steady-state solutions of the homogenised model and show that for high enough applied currents, two spatially heterogeneous steady-state solutions exist, with distinct crater sizes. We show that the system evolves to the steady state with a larger crater radius and explain this behaviour in terms of the overall power balance typically found within a furnace. We find that the industrial practice of stoking furnaces increases the overall rate of material consumption in the furnace, thereby improving the efficiency of silicon production.

Quasiperiodic media is a class of almost periodic media which is generated from periodic media through a ‘cut and project’ procedure. Quasiperiodic media displays some extraordinary optical, electronic and conductivity properties which call for the development of methods to analyse their microstructures and effective behaviour. In this paper, we develop the method of Bloch wave homogenisation for quasiperiodic media. Bloch waves are typically defined through a direct integral decomposition of periodic operators. A suitable direct integral decomposition is not available for almost periodic operators. To remedy this, we lift a quasiperiodic operator to a degenerate periodic operator in higher dimensions. Approximate Bloch waves are obtained for a regularised version of the degenerate operator. Homogenised coefficients for quasiperiodic media are obtained from the first Bloch eigenvalue of the regularised operator in the limit of regularisation parameter going to zero. A notion of quasiperiodic Bloch transform is defined and employed to obtain homogenisation limit for an equation with highly oscillating quasiperiodic coefficients.

The aim of the paper is to explore non-local reverse-space matrix non-linear Schrödinger equations and their inverse scattering transforms. Riemann–Hilbert problems are formulated to analyse the inverse scattering problems, and the Sokhotski–Plemelj formula is used to determine Gelfand–Levitan–Marchenko-type integral equations for generalised matrix Jost solutions. Soliton solutions are constructed through the reflectionless transforms associated with poles of the Riemann–Hilbert problems.

Semi-supervised and unsupervised machine learning methods often rely on graphs to model data, prompting research on how theoretical properties of operators on graphs are leveraged in learning problems. While most of the existing literature focuses on undirected graphs, directed graphs are very important in practice, giving models for physical, biological or transportation networks, among many other applications. In this paper, we propose a new framework for rigorously studying continuum limits of learning algorithms on directed graphs. We use the new framework to study the PageRank algorithm and show how it can be interpreted as a numerical scheme on a directed graph involving a type of normalised graph Laplacian. We show that the corresponding continuum limit problem, which is taken as the number of webpages grows to infinity, is a second-order, possibly degenerate, elliptic equation that contains reaction, diffusion and advection terms. We prove that the numerical scheme is consistent and stable and compute explicit rates of convergence of the discrete solution to the solution of the continuum limit partial differential equation. We give applications to proving stability and asymptotic regularity of the PageRank vector. Finally, we illustrate our results with numerical experiments and explore an application to data depth.

We extend previous weak well-posedness results obtained in Frigeri et al. (2017, Solvability, Regularity, and Optimal Control of Boundary Value Problems for PDEs, Vol. 22, Springer, Cham, pp. 217–254) concerning a non-local variant of a diffuse interface tumour model proposed by Hawkins-Daarud et al. (2012, Int. J. Numer. Method Biomed. Engng.28, 3–24). The model consists of a non-local Cahn–Hilliard equation with degenerate mobility and singular potential for the phase field variable, coupled to a reaction–diffusion equation for the concentration of a nutrient. We prove the existence of strong solutions to the model and establish some high-order continuous dependence estimates, even in the presence of concentration-dependent mobilities for the nutrient variable in two spatial dimensions. Then, we apply the new regularity results to study an inverse problem identifying the initial tumour distribution from measurements at the terminal time. Formulating the Tikhonov regularised inverse problem as a constrained minimisation problem, we establish the existence of minimisers and derive first-order necessary optimality conditions.

We derive and numerically implement various asymptotic approximations for the lowest or principal eigenvalue of the Laplacian with a periodic arrangement of localised traps of small \[\mathcal{O}(\varepsilon )\] spatial extent that are centred at the lattice points of an arbitrary Bravais lattice in \[{\mathbb{R}^2}\]. The expansion of this principal eigenvalue proceeds in powers of \[\nu \equiv - 1/\log (\varepsilon {d_c})\], where dc is the logarithmic capacitance of the trap set. An explicit three-term approximation for this principal eigenvalue is derived using strong localised perturbation theory, with the coefficients in this series evaluated numerically by using an explicit formula for the source-neutral periodic Green’s function and its regular part. Moreover, a transcendental equation for an improved approximation to the principal eigenvalue, which effectively sums all the logarithmic terms in powers of v, is derived in terms of the regular part of the periodic Helmholtz Green’s function. By using an Ewald summation technique to first obtain a rapidly converging infinite series representation for this regular part, a simple Newton iteration scheme on the transcendental equation is implemented to numerically evaluate the improved ‘log-summed’ approximation to the principal eigenvalue. From a numerical computation of the PDE eigenvalue problem defined on the fundamental Wigner–Seitz (WS) cell for the lattice, it is shown that the three-term asymptotic approximation for the principal eigenvalue agrees well with the numerical result only for a rather small trap radius. In contrast, the log-summed asymptotic result provides a very close approximation to the principal eigenvalue even when the trap radius is only moderately small. For a circular trap, the first few transcendental correction terms that further improves the log-summed approximation for the principal eigenvalue are derived. Finally, it is shown numerically that, amongst all Bravais lattices with a fixed area of the primitive cell, the principal eigenvalue is maximised for a regular hexagonal arrangement of traps.

Existence of non-negative weak solutions is shown for a full curvature thin-film model of a liquid thin film flowing down a vertical fibre. The proof is based on the application of a priori estimates derived for energy-entropy functionals. Long-time behaviour of these weak solutions is analysed and, under some additional constraints for the model parameters and initial values, convergence towards a travelling wave solution is obtained. Numerical studies of energy minimisers and travelling waves are presented to illustrate analytical results.

Based on the fact that the incubation periods of epidemic disease in asymptomatically infected and infected individuals are inevitable and different, we propose a diffusive susceptible, asymptomatically infected, symptomatically infected and vaccinated (SAIV) epidemic model with delays in this paper. To see whether epidemic disease can propagate spatially with a constant speed, we focus on the travelling wave solution for this model. When the basic reproduction number of the corresponding spatial-homogenous delayed differential system is greater than one and the wave speed is greater than or equal to the critical speed, we prove that this model admits nontrivial positive travelling wave solutions. Our theoretical results are of benefit to the prevention and control of epidemic.