We propose and study a class of parabolic-ordinary differential equation models involving chemotaxis and haptotaxis of a species following signals indirectly produced by another, non-motile one. The setting is motivated by cancer invasion mediated by interactions with the tumour microenvironment, but has much wider applicability, being able to comprise descriptions of biologically quite different problems. As a main mathematical feature constituting a core difference to both classical Keller–Segel chemotaxis systems and Chaplain–Lolas type chemotaxis–haptotaxis systems, the considered model accounts for certain types of indirect signal production mechanisms. The main results assert unique global classical solvability under suitably mild assumptions on the system parameter functions in associated spatially two-dimensional initial-boundary value problems. In particular, this rigorously confirms that at least in two-dimensional settings, the considered indirectness in signal production induces a significant blow-up suppressing tendency also in taxis systems substantially more general than some particular examples for which corresponding effects have recently been observed.

The Neumann–Poincaré (NP) operator, a singular integral operator on the boundary of a domain, naturally appears when one solves a conductivity transmission problem via the boundary integral formulation. Recently, a series expression of the NP operator was developed in two dimensions based on geometric function theory [34]. In this paper, we investigate geometric properties of composite materials using this series expansion. In particular, we obtain explicit formulas for the polarisation tensor and the effective conductivity for an inclusion or a periodic array of inclusions of arbitrary shape with extremal conductivity, in terms of the associated exterior conformal mapping. Also, we observe by numerical computations that the spectrum of the NP operator has a monotonic behaviour with respect to the shape deformation of the inclusion. Additionally, we derive inequality relations of the coefficients of the Riemann mapping of an arbitrary Lipschitz domain using the properties of the polarisation tensor corresponding to the domain.

When ceramics are heated inside a microwave cavity, a well-known phenomenon is the occurrence of hot spots – localised regions of high temperature. This phenomenon was modelled by Kriegsmann ((1997), IMA J. Appl. Math. 59(2), pp. 123–146; (2001), IMA J. Appl. Math. 66(1), pp. 1–32) using a non-local evolution PDE. We investigate profile and the stability of hot spots in one and two dimensions by using Kriegsmann's model with exponential non-linearity. The linearised problem associated with hot-spot-type solutions possesses two classes of eigenvalues. The first type is the large eigenvalues associated with the stability of the hot-spot profile and in this particular model there cannot be instability associated with these eigenvalues. The second type is the small eigenvalues associated with translation invariance. We show that the hot spots can become unstable due to the presence of small eigenvalues, and we characterise the instability thresholds. In particular, we show that for the material with low heat conductivity (such as ceramics), and in the presence of a variable electric field, the hot spots are typically stable inside a plate (in two dimensions) but can become unstable for a slab (in one dimension) provided that the microwave power is sufficiently large. On the other hand, for materials with high heat conductivity, the interior hot spots are unstable and move to the boundary of the domain in either one or two dimensions. For materials with moderate heat conductivity, the stability of hot spots is determined by both the geometry and the electric field inside the microwave cavity.

Recently we have proposed a monostable reaction-diffusion system to explain the Neolithic transition from hunter-gatherer life to farmer life in Europe. The system is described by a three-component system for the populations of hunter-gatherer (H), sedentary farmer (F1) and migratory one (F2). The conversion between F1 and F2 is specified by such a way that if the total farmers F1 + F2 are overcrowded, F1 actively changes to F2, while if it is less crowded, the situation is vice versa. In order to include this property in the system, the system incorporates a critical parameter (say F0) depending on the development of farming technology in a monotonically increasing way. It determines whether the total farmers are either over crowded (F1 + F2 >F0) or less crowded (F1 + F2 <F0) ( [9, 20]). Previous numerical studies indicate that the structure of travelling wave solutions of the system is qualitatively similar to the one of the Fisher-KPP equation, that the asymptotically expanding velocity of farmers is equal to the minimal velocity (say cm(F0)) of travelling wave solutions, and that cm(F0) is monotonically decreasing as F0 increases. The latter result suggests that the development of farming technology suppresses the expanding velocity of farmers. As a partial analytical result to this property, the purpose of this paper is to consider the two limiting cases where F0 = 0 and F0 → ∞, and to prove cm(0)>cm(∞).

In this paper we prove the existence of a solution for a system of nonlinear parabolic partial differential equations arising from thermoelectric modelling of metallurgical electrodes undergoing a phase change. The model consists of an electromagnetic problem for eddy current computation coupled with a Stefan problem for temperature. The proof uses a regularized problem obtained by truncating the source term in temperature equation. Passing to the limit requires fine a priori estimates leading to compactness.

The use of both higher-order amplitude corrections and complex rays to increase the accuracy of the approximation of the eigenvalues of Helmholtz equation on an annulus is demonstrated.

We study a system of two coupled parabolic equations that models the spreading of a drop of an insoluble surfactant on a thin liquid film. Unlike the previously known results, the surface diffusion coefficient is not assumed constant and depends on the surfactant concentration. We obtain sufficient conditions for finite speed of support propagation and for waiting-time phenomenon by application of an extension of Stampacchia's lemma for a system of functional equations.

In this paper, a bilayer model is derived to simulate the evolution of a thin film flow over water. This model is derived from the incompressible Navier–Stokes equations together with suitable boundary conditions including friction and capillary effects. The derivation is based on the different properties of the fluids; thus, we perform a multiscale analysis in space and time, and a different asymptotic analysis to derive a system coupling two different models: the Reynolds lubrication equation for the upper layer and the shallow water model for the lower one. We prove that the model verifies a dissipative entropy inequality up to a second-order term. Moreover, we propose a correction of the model – by taking into account the second-order extension for the pressure – that admits an exact dissipative entropy inequality. Two numerical tests are presented. In the first test, we compare the numerical results with the viscous bilayer shallow water model proposed in Narbona-Reina et al. (Comput. Model. Eng. Sci., 2009, Vol. 43, pp. 27–71). In the second test, the objective is to show some of the characteristic situations that can be studied with the proposed model. We simulate a problem of pollutant dispersion near the coast. For this test, the influence of the friction coefficient on the coastal area affected by the pollutant is studied.

This paper is concerned with some mathematical aspects of magnetic resonance imaging (MRI) of the beating human heart. In particular, we investigate the so-called retrospective gating technique which is a non-triggered technique for data acquisition and reconstruction of (approximately) periodically changing organs like the heart. We formulate the reconstruction problem as a moment problem in a Hilbert space and give the solution method. The stability of the solution is investigated and various error estimates are given. The reconstruction method consists of temporal interpolation followed by spatial Fourier inversion. Different choices for the Hilbert space ℋ of interpolating functions are possible. In particular, we study the case where ℋ is (i) the space of bandlimited functions, or (ii) the space of spline functions of odd degree. The theory is applied to reconstructions from synthetic data as well as real MRI data.

The dynamics of moving solids with unilateral contacts are often modelled by assuming rigidity, point contacts, and Coulomb friction. The canonical example of a rigid rod with one endpoint slipping in two dimensions along a fixed surface (sometimes referred to as Painlevé rod) has been investigated thoroughly by many authors. The generic transitions of that system include three classical transitions (slip-stick, slip reversal, and liftoff) as well as a singularity called dynamic jamming, i.e., convergence to a codimension 2 manifold in state space, where rigid body theory breaks down. The goal of this paper is to identify similar singularities arising in systems with multiple point contacts, and in a broader setting to make initial steps towards a comprehensive list of generic transitions from slip motion to other types of dynamics. We show that – in addition to the classical transitions – dynamic jamming remains a generic phenomenon. We also find new forms of singularity and solution indeterminacy, as well as generic routes from sliding to self-excited microscopic or macroscopic oscillations.

The high Reynolds-number structure of the laminar, chaotic and turbulent attractors is investigated in a two-dimensional Kolmogorov flow. The laminar attractors include the families of multi-phased travelling waves and quasi-periodic standing waves both of which form the backbone of the transition to a turbulent flow. At leading order, each laminar attractor under study is obtained by solving the Euler equations on a manifold subject to the appropriate periodicity and symmetry conditions. The manifold is determined by a finite number of vorticity equations, these being required to suppress the secular terms at the next order. Our results show that, for the multi-phased travelling waves, the first phase velocity can be determined by an integral conservation law for kinetic energy and the subsequent phase velocities can be evaluated by a non-linear eigenvalue problem. The results also reveal that whereas viscosity determines the smallest scales and controls the amplitude of the flow, the inertial terms govern the shape and form of the flow. The comparison of our analytical predictions for evaluating the stable single-phased travelling wave with the direct numerical simulation of the Navier–Stokes equations has been undertaken, the agreement being excellent. For sufficiently high Reynolds number, after the bifurcation to chaotic flow, all of the multi-phased travelling waves and quasi-periodic standing waves become unstable non-wandering sets. Based on the above new findings for these unstable non-wandering sets and other travelling and standing waves of this kind in phase space, necessary conditions for the invariant manifolds of the chaotic and turbulent attractors are obtained, these necessary conditions being conjectured to be also sufficient.

We study the homogenisation of the Stokes system with a non-homogeneous Fourier boundary condition on the boundary of the holes, depending on a parameter γ. Such systems arise in the modelling of the flow of an incompressible viscous fluid through a porous medium under the influence of body forces. At the limit, by using the periodic unfolding method in perforated domains, we obtain, following the values of γ, different Darcy's laws of type Mu = −N∇p + F with suitable matrices M and N with F depending on the right-hand side in the bulk term and in the boundary condition.

We consider a spherical particle levitating above a liquid bath owing to the Leidenfrost effect, where the vapour of either the bath or sphere forms an insulating film whose pressure supports the sphere’s weight. Starting from a reduced formulation based on a lubrication-type approximation, we use matched asymptotics to describe the morphology of the vapour film assuming that the sphere is small relative to the capillary length (small Bond number) and that the densities of the bath and sphere are comparable. We find that this regime is comprised of two formally infinite sequences of distinguished limits which meet at an accumulation point, the limits being defined by the smallness of an intrinsic evaporation number relative to the Bond number. These sequences of limits reveal a surprisingly intricate evolution of the film morphology with increasing sphere size. Initially, the vapour film transitions from a featureless morphology, where the thickness profile is parabolic, to a neck–bubble morphology, which consists of a uniform pressure bubble bounded by a narrow and much thinner annular neck. Gravity effects then become important in the bubble leading to sequential formation of increasingly smaller neck–bubble pairs near the symmetry axis. This process terminates when the pairs closest to the symmetry axis become indistinguishable and merge. Subsequently, the inner section of that merger transitions into a uniform-thickness film that expands radially, gradually squishing increasingly larger neck–bubble pairs into a region of localised oscillations sandwiched between the uniform film and what remains of the bubble whose radial extent is presently comparable to the uniform film; the neck–bubble pairs farther from the axis remain essentially intact. Ultimately, the uniform film gobbles up the largest outermost bubble, whereby the morphology simplifies to a uniform film bounded by localised oscillations. Overall, the asymptotic analysis describes the continuous evolution of the vapour film from a neck–bubble morphology typical of a Leidenfrost drop levitating above a flat solid substrate to a uniform-film morphology which resembles that in the case of a large liquid drop levitating above a liquid bath.

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.

We study of a new type of multi-bump blowup solutions of the Ginzburg–Landau equation. Multi-bump blowup solutions have previously been found in numeric simulations, asymptotic analysis and were proved to exist via geometric construction. In the geometric construction of the solutions, the existence of two types of multi-bump solutions was shown. One type is exponentially small at ξ=0, and the other type of solutions is algebraically small at ξ=0. So far, the first type of solutions were studied asymptotically. Here, we analyse the solutions which are algebraically small at ξ=0 by using asymptotic methods. This construction is essentially different from the existing one, and ideas are obtained from the geometric construction. Hence, this is a good example of where both asymptotic analysis and geometric methods are needed for the overall picture.

We consider heat transfer processes where radiation in a large number of frequency bands plays a dominant role. In the simulation of such processes, the radiative transfer equation has to be solved repeatedly. To obtain an efficient and accurate solution method, we propose a new hybrid algorithm which combines fast solvers for the radiative transfer equation in the low and high absorption regime, respectively. A key role is played by an exact, residual based error formula. The algorithm is applied to a cooling problem of high quality optical glass.

Basic methods from bifurcation theory are applied to derive a criterion that predicts when a symmetric tricritical point may occur in a transition between condensed ordered phases described by any finite number of scalar order parameters. At such a point, a change of order takes place in the phase transition, which passes from first to second order, or vice versa.

We consider the following system of equations: where the spatial average ⟨ B ⟩ = 0 and μ > σ > 0. This system plays an important role as a Ginzburg-Landau equation with a mean field in several areas of the applied sciences and the steady-states of this system extend to periodic steady-states with period L on the real line which are observed in experiments. Our approach is by combining methods of nonlinear functional analysis such as nonlocal eigenvalue problems and the variational characterization of eigenvalues with Jacobi elliptic integrals. This enables us to give a complete classification of all stable steady-states for any positive L.

The gradient flow structure of the model introduced in Cermelli & Gurtin (1999, The motion of screw dislocations in crystalline materials undergoing antiplane shear: glide, cross-slip, fine cross-slip. Arch. Rational Mech. Anal.148(1), 3–52) for the dynamics of screw dislocations is investigated by means of a generalised minimising-movements scheme approach. The assumption of a finite number of available glide directions, together with the “maximal dissipation criterion” that governs the equations of motion, results into solving a differential inclusion rather than an ODE. This paper addresses how the model in Cermelli & Gurtin is connected to a time-discrete evolution scheme which explicitly confines dislocations to move at each time step along a single glide direction. It is proved that the time-continuous model in Cermelli & Gurtin is the limit of these time-discrete minimising-movement schemes when the time step converges to 0. The study presented here is a first step towards a generalisation of standard gradient flow theory that allows for dissipations which cannot be described by a metric.

This special issue is one of the very first dedicated to crime modelling in a journal of applied mathematics. It emphasizes one of the new areas at the Social Science frontier, where modelling and mathematical tools are put to use with a view to shed light on phenomena previously thought to be outside of their reach. Pioneering research is increasingly being carried out in many different areas in the life sciences or social sciences, often under the heading of the study of complex systems. When addressing issues regarding society, individuals or the collective behaviours of humans, several questions naturally arise about the modelling enterprise. What is the nature and role of modelling in social sciences? What is one to expect from these new approaches? The case of economics, which has relied on mathematics for a very long time now, can serve as a paradigm for what is happening in other social sciences.