We consider the explicit solution to the axisymmetric diffusion equation. We recast the solution in the form of a Mellin inversion formula, and outline a method to compute a formula for $u(r,t)$ as a series using the Cauchy residue theorem. As a consequence, we are able to represent the solution to the axisymmetric diffusion equation as a rapidly converging series.

We consider a model for the dynamics of growing cell populations with heterogeneous mobility and proliferation rate. The cell phenotypic state is described by a continuous structuring variable and the evolution of the local cell population density function (i.e. the cell phenotypic distribution at each spatial position) is governed by a non-local advection–reaction–diffusion equation. We report on the results of numerical simulations showing that, in the case where the cell mobility is bounded, compactly supported travelling fronts emerge. More mobile phenotypic variants occupy the front edge, whereas more proliferative phenotypic variants are selected at the back of the front. In order to explain such numerical results, we carry out formal asymptotic analysis of the model equation using a Hamilton–Jacobi approach. In summary, we show that the locally dominant phenotypic trait (i.e. the maximum point of the local cell population density function along the phenotypic dimension) satisfies a generalised Burgers’ equation with source term, we construct travelling-front solutions of such transport equation and characterise the corresponding minimal speed. Moreover, we show that, when the cell mobility is unbounded, front edge acceleration and formation of stretching fronts may occur. We briefly discuss the implications of our results in the context of glioma growth.

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.

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.

This work studies functional difference equations of the second order with a potential belonging to a special class of meromorphic functions. The equations depend on a spectral parameter. Consideration of this type of equations is motivated by applications in diffraction theory and by construction of eigenfunctions for the Laplace operator in angular domains. In particular, such eigenfunctions describe eigenoscillations of acoustic waves in angular domains with ‘semitransparent’ boundary conditions. For negative values of the spectral parameter, we study essential and discrete spectrum of the equations and describe properties of the corresponding solutions. The study is based on the reduction of the functional difference equations to integral equations with a symmetric kernel. A sufficient condition is formulated for the potential that ensures existence of the discrete spectrum. The obtained results are applied for studying the behaviour of eigenfunctions for the Laplace operator in adjacent angular domains with the Robin-type boundary conditions on their common boundary. At infinity, the eigenfunctions vanish exponentially as was expected. However, the rate of such decay depends on the observation direction. In particular, in a vicinity of some directions, the regime of decay is switched from one to another and such asymptotic behaviour is described by a Fresnel-type integral.

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.

A Doyle–Fuller–Newman (DFN) model for the charge and discharge of nano-structured lithium iron phosphate (LFP) cathodes is formulated on the basis that lithium transport within the nanoscale LFP electrode particles is much faster than cell discharge, and is therefore not rate limiting. We present some numerical solutions to the model and show that for relevant parameter values, and a variety of C-rates, it is possible for sharp discharge fronts to form and intrude into the electrode from its outer edge(s). These discharge fronts separate regions of fully utilised LFP electrode particles from those that are not. Motivated by this observation an asymptotic solution to the model is sought. The results of the asymptotic analysis of the DFN model lead to a reduced order model, which we term the reaction front model (or RFM). Favourable agreement is shown between solutions to the RFM and the full DFN model in appropriate parameter regimes. The RFM is significantly cheaper to solve than the DFN model, and therefore has the potential to be used in scenarios where computational costs are prohibitive, e.g. in optimisation and parameter estimation problems or in engineering control systems.

A hybrid asymptotic-numerical method is developed to approximate the mean first passage time (MFPT) and the splitting probability for a Brownian particle in a bounded two-dimensional (2D) domain that contains absorbing disks, referred to as “traps”, of asymptotically small radii. In contrast to previous studies that required traps to be spatially well separated, we show how to readily incorporate the effect of a cluster of closely spaced traps by adapting a recently formulated least-squares approach in order to numerically solve certain local problems for the Laplacian near the cluster. We also provide new asymptotic formulae for the MFPT in 2D spatially periodic domains where a trap cluster is centred at the lattice points of an oblique Bravais lattice. Over all such lattices with fixed area for the primitive cell, and for each specific trap set, the average MFPT is smallest for a hexagonal lattice of traps.

In this paper, we study the initial-boundary value problem of a repulsion Keller–Segel system with a logarithmic sensitivity modelling the reinforced random walk. By establishing an energy–dissipation identity, we prove the existence of classical solutions in two dimensions as well as existence of weak solutions in the three-dimensional setting. Moreover, it is shown that the weak solutions enjoy an eventual regularity property, i.e., it becomes regular after certain time T > 0. An exponential convergence rate towards the spatially homogeneous steady states is obtained as well. We adopt a new approach developed recently by the author to study the eventual regularity. The argument is based on observation of the exponential stability of constant solutions in scaling-invariant spaces together with certain dissipative property of the global solutions in the same spaces.

We propose and investigate a stage-structured SLIRM epidemic model with latent period in a spatially continuous habitat. We first show the existence of semi-travelling waves that connect the unstable disease-free equilibrium as the wave coordinate goes to − ∞, provided that the basic reproduction number $\mathcal {R}_0 > 1$ and $c > c_*$ for some positive number $c_*$. We then use a combination of asymptotic estimates, Laplace transform and Cauchy's integral theorem to show the persistence of semi-travelling waves. Based on the persistent property, we construct a Lyapunov functional to prove the convergence of the semi-travelling wave to an endemic (positive) equilibrium as the wave coordinate goes to + ∞. In addition, by the Laplace transform technique, the non-existence of bounded semi-travelling wave is also proved when $\mathcal {R}_0 > 1$ and $0 < c < c_*$. This indicates that $c_*$ is indeed the minimum wave speed. Finally simulations are given to illustrate the evolution of profiles.

The asymptotic behavior of solutions to a family of Dirichlet boundary value problems, involving differential operators in divergence form, on a domain equipped with a Finsler metric is investigated. Solutions are shown to converge uniformly to the distance function to the boundary of the domain, which takes into account the Finsler norm involved in the equation. This implies that a well-known result in the analysis of problems modeling torsional creep continues to hold in this more general setting.

We study the existence of entropy solutions by assuming the right-hand side function f to be an integrable function for some elliptic nonlocal p-Laplacian type problems. Moreover, the existence of weak solutions for the corresponding parabolic cases is also established. The main aim of this paper is to provide some positive answers for the two questions proposed by Chipot and de Oliveira (Math. Ann., 2019, 375, 283-306).

We prove that any simple planar travelling wave solution to the membrane equation in spatial dimension $d\geqslant 3$ with bounded spatial extent is globally nonlinearly stable under sufficiently small compactly supported perturbations, where the smallness depends on the size of the support of the perturbation as well as on the initial travelling wave profile. The main novelty of the argument is the lack of higher order peeling in our vector-field-based method. In particular, the higher order energies (in fact, all energies at order $2$ or higher) are allowed to grow polynomially (but in a controlled way) in time. This is in contrast with classical global stability arguments, where only the ‘top’ order energies used in the bootstrap argument exhibit growth, and reflects the fact that the background travelling wave solution has ‘infinite energy’ and the coefficients of the perturbation equation are not asymptotically Lorentz invariant. Nonetheless, we can prove that the perturbation converges to zero in $C^{2}$ by carefully analysing the nonlinear interactions and exposing a certain ‘vestigial’ null structure in the equations.

The work in this paper concerns the study of different approximations for one-dimensional one-phase Stefan-like problems with a space-dependent latent heat. It is considered two different problems, which differ from each other in their boundary condition imposed at the fixed face: Dirichlet and Robin conditions. The approximate solutions are obtained by applying the heat balance integral method (HBIM), the modified HBIM and the refined integral method (RIM). Taking advantage of the exact analytical solutions, we compare and test the accuracy of the approximate solutions. The analysis is carried out using the dimensionless generalised Stefan number (Ste) and Biot number (Bi). It is also studied the case when Bi goes to infinity in the problem with a convective condition, recovering the approximate solutions when a temperature condition is imposed at the fixed face. Some numerical simulations are provided in order to assert which of the approximate integral methods turns out to be optimal. Moreover, we pose an approximate technique based on minimising the least-squares error, obtaining also approximate solutions for the classical Stefan problem.

This paper is concerned with the existence results for generalized transition waves of space periodic and time heterogeneous lattice Fisher-KPP equations. By constructing appropriate subsolutions and supersolutions, we show that there is a critical wave speed such that a transition wave solution exists as soon as the least mean of wave speed is above this critical speed. Moreover, the critical speed we construct is proved to be minimal in some particular cases, such as space-time periodic or space independent.

In this article, we establish a Lyapunov-type inequality for the following extremal Pucci’s equation:

$$\begin{eqnarray}\left\{\begin{array}{@{}ll@{}}{\mathcal{M}}_{\unicode[STIX]{x1D706},\unicode[STIX]{x1D6EC}}^{+}(D^{2}u)+b(x)|Du|+a(x)u=0 & \text{in}~\unicode[STIX]{x1D6FA},\\ u=0 & \text{on}~\unicode[STIX]{x2202}\unicode[STIX]{x1D6FA},\end{array}\right.\end{eqnarray}$$

$\unicode[STIX]{x1D6FA}$

$\mathbb{R}^{N}$

$N\geq 2$

In this article, we consider diffusive transport of a reactive substance in a saturated porous medium including variable porosity. Thereby, the evolution of the microstructure is caused by precipitation of the transported substance. We are particularly interested in analysing the model when the equations degenerate due to clogging. Introducing an appropriate weighted function space, we are able to handle the degeneracy and obtain analytical results for the transport equation. Also the decay behaviour of this solution with respect to the porosity is investigated. There a restriction on the decay order is assumed, that is, besides low initial concentration also dense precipitation leads to possible high decay. We obtain nonnegativity and boundedness for the weak solution to the transport equation. Moreover, we study an ordinary differential equation (ODE) describing the change of porosity. Thereby, the control of an appropriate weighted norm of the gradient of the porosity is crucial for the analysis of the transport equation. In order to obtain global in time solutions to the overall coupled system, we apply a fixed point argument. The problem is solved for substantially degenerating hydrodynamic parameters.

In this paper we study a general eigenvalue problem for the so called (p, 2)-Laplace operator on a smooth bounded domain Ω ⊂ ℝN under a nonlinear Steklov type boundary condition, namely

\[\left\{ \begin{aligned} -\Delta_pu-\Delta u & =\lambda a(x)u \quad {\rm in}\ \Omega,\\ (|\nabla u|^{p-2}+1)\dfrac{\partial u}{\partial\nu} & =\lambda b(x)u \quad {\rm on}\ \partial\Omega . \end{aligned} \right.\]

In this paper, we study the nonlinear diffusion equation associated with a particle system where the common drift depends on the rate of absorption of particles at a boundary. We provide an interpretation of this equation, which is also related to the supercooled Stefan problem, as a structural credit risk model with default contagion in a large interconnected banking system. Using the method of heat potentials, we derive a coupled system of Volterra integral equations for the transition density and for the loss through absorption. An approximation by expansion is given for a small interaction parameter. We also present a numerical solution algorithm and conduct computational tests.