A mathematical model of the development of an oil and gas field is presented. The field development process is treated as sequential in nature. Completion of a well and its production are considered to be random processes. The model uses results from renewal theory where the completion of a well and failure to produce economical amount of oil or gas are analogous to the failure of a component. In principle, the theory described can give the complete probability distribution associated with a field development. Explicit expressions are given for the expected value and variance of the number of completed wells.

The following first problem is posed: is a correct ‘entropy solution’ of the Cauchy problem for the fifth-order degenerate non-linear dispersion equations (NDEs), same as for the classic Euler one ut + uux = 0, These two quasi-linear degenerate partial differential equations (PDEs) are chosen as typical representatives; so other (2m + 1)th-order NDEs of non-divergent form admit such shocks waves. As a related second problem, the opposite initial shock S+(x) = −S−(x) = sign x is shown to be a non-entropy solution creating a rarefaction wave, which becomes C∞ for any t > 0. Formation of shocks leads to non-uniqueness of any ‘entropy solutions’. Similar phenomena are studied for a fifth-order in time NDE uttttt = (uux)xxxx in normal form.

On the other hand, related NDEs, such as are shown to admit smooth compactons, as oscillatory travelling wave solutions with compact support. The well-known non-negative compactons, which appeared in various applications (first examples by Dey, 1998, Phys. Rev. E, vol. 57, pp. 4733–4738, and Rosenau and Levy, 1999, Phys. Lett. A, vol. 252, pp. 297–306), are non-existent in general and are not robust relative to small perturbations of parameters of the PDE.

We consider a family of multi-phase Stefan problems for a certain one-dimensional model of cell-to-cell adhesion and diffusion, which takes the form of a non-linear forward–backward parabolic equation. In each material phase the cell density stays either high or low, and phases are connected by jumps across an ‘unstable’ interval. We develop an existence theory for such problems which allows for the annihilation of phases and the subsequent continuation of solutions. Stability results for the long-time behaviour of solutions are also obtained, and, where necessary, the analysis is complemented by numerical simulations.

The coupled motion between shallow-water sloshing in a moving vehicle and the vehicle dynamics is considered, with the vehicle dynamics restricted to horizontal motion. The paper is motivated by Cooker's experiments and theory for water waves in a suspended container. A new derivation of the coupled problem in the Eulerian fluid representation is given. However, it is found that transformation to a Lagrangian representation leads to a formulation which has nice properties for numerical simulation. In the Lagrangian representation, a simple and fast numerical algorithm with excellent energy conservation over long times, based on the Störmer–Verlet method, is implemented. Numerical simulations of the coupled dynamics in both the linear and non-linear case are presented.

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.

We study the coupled surface and grain boundary motion in bi- and tricrystals in three-space dimensions, building on previous work by the authors on the simplified two-dimensional case. The motion of the interfaces, which in this paper are presented by two-dimensional hypersurfaces, is described by two types of normal velocities: motion by mean curvature and motion by surface diffusion. Three hypersurfaces meet at triple-junction lines, where junction conditions need to hold. Similarly, boundary conditions are prescribed where an interface meets an external boundary, and these conditions naturally give rise to contact angles. We present a variational formulation of the flows, which leads to a fully practical finite-element approximation that exhibits excellent mesh properties, with no mesh smoothing or remeshing required in practice. For the introduced parametric finite-element approximation we show well posedness and, in general, unconditional stability, i.e. there is no restriction on the chosen time-step size. Moreover, the induced discrete equations are linear and easy to solve. A generalisation to anisotropic surface energies is straightforward. Several numerical results in two- and three-space dimensions are presented, including simulations for thermal grooving and sintering. Three-dimensional simulations featuring quadruple junction points, non-standard boundary contact angles and fully anisotropic surface energies are also presented.

We consider a dynamic frictional contact problem between an elastic-visco-plastic body and a foundation. The contact is modelled with a normal damped response condition of such a type that the normal velocity is restricted with unilateral constraint, associated with the Coulomb law in which the coefficient of friction may depend on the velocity. We derive a variational formulation of the problem which has the form of a system coupling an integro–differential equation for the stress field with an evolutionary variational inequality for the displacement field. This inequality is approximated by a variational equation using a smoothing of the friction and the penalty approximation of the unilateral condition. The existence of a weak solution to the variational equation is proved by the Galerkin method for an auxiliary problem with given viscoplastic part of the stress and a fixed point argument. The solvability of the original problem is proved by passing to the limit of the penalty parameter and the smoothing parameter. This convergence is based on a certain regularity of solutions which is verified with the use of a local rectification of the boundary and a translation method.

This paper by a criminologist explains why it makes more sense to model criminal acts than to model criminals, how many preconceptions about crime can mislead modellers and offers some simple crime modelling ideas. Many opportunities for simulation now exist, and new opportunities for real-data modelling are emerging. The author suggests mathematical models of crime, including offender foraging for crime targets, as a rich area for future research.

We study the behaviour of the steady-state voltage potentials in a material composed of a two-dimensional object surrounded by a rough thin layer and embedded in an ambient medium. The roughness of the layer is supposed to be ϵ-periodic, ϵ being the magnitude of the mean thickness of the layer. For ϵ tending to zero, we determine approximate transmission conditions in order to replace the rough thin layer by these conditions on the boundary of the interior material. This paper extends the previous works (Poignard, 2009, Math. Meth. Appl. Sci., vol. 32, pp. 435–453; Poignard et al., 2008, IEEE Trans. Magnet., vol. 44, no. 6, pp. 1154–1157) of the third author, which deal with smooth thin layers.

Proceeding from a recent construction (Fraenkel 2007, Arch. Rational Mech. Anal. 183, 187–214), this paper contains three contributions to the theory of the extreme gravity wave on water (or “wave of greatest height”) initiated by Stokes in 1880. The first is a solution, of the extreme form of the integral equation due to Nekrasov, that consists of an explicit function plus a (rigorously) estimated remainder; the remainder vanishes at the crest and at the trough of the wave and in between has a maximum value less than 0.0026 times the maximum value of the explicit part. The second contribution establishes local uniqueness for a variant of the Nekrasov equation; we prove that this equation has only one solution between upper and lower bounding functions that are moderately far apart. The third contribution consists of a table and graph of strict upper and lower bounds for the interface between air and water in the physical plane.

We are interested in the possible contributions of mathematical modelling of crime. We refer to numerous and quite recent papers that analyse and discuss empirical data in an attempt to discover stylised trends worthy of being understood through simple models. We summarise part of this literature and try to understand the reasons of important discrepancies in their conclusions. Then we present some recent modelling attempts that may help to understand the large variance in the statistical-based conclusions.

Motivated by the technology of magnetically targeted drug and gene delivery, in which a magnetic field is used to direct magnetic carrier particles from the circulation to a target site, we develop a continuum model for the motion of particles (magnetic carriers) subject to an external body force (magnetic field) in a flow of a concentrated suspension of a species of neutrally buoyant particles (blood). An advection–diffusion equation describes the evolution of the carrier particles as they advect in the flow under the action of an external body force, and diffuse as a result of random interactions with the suspension of neutrally buoyant particles (shear-induced diffusion). The model is analysed for the case in which there is steady Poiseuille flow in a cylindrical vessel, the diffusive effects are weak and there is weak carrier uptake along the walls of the vessel. The method of matched asymptotic expansions is used to show that carriers are concentrated in a boundary layer along the vessel wall and, further, that there is a carrier flux along this layer which results in a sub-layer, along one side of the vessel, in which carriers are even more highly concentrated. Three distinguished limits are identified: they correspond to cases for which (i) the force is sufficiently weak that most particles move through the vessel without entering the boundary layers along the walls of the vessel and (ii) and (iii) to a force which is sufficiently strong that a significant fraction of the particles enter the boundary layers and, depending upon the carrier absorption from the vessel walls, there is insignificant/significant axial carrier flux in these layers.

In the preceding paper by Fraenkel & Harwin, the existence theorem leaves a gap in the description of the extreme wave that is constructed there. The form of the wave in a neighbourhood of the crest is not described accurately enough. This weakness is repaired in the present paper.

This paper is devoted to the study of the asymptotic behaviour of the high-frequency spectrum of the wave equation with periodic coefficients in a ‘thin’ elastic strip Ση=(0, 1)×(−η/2, η/2), η > 0. The main geometric assumption is that the structure period is of the order of magnitude of the strip thickness η and is chosen in such a way that η−1 is a positive large integer. On the boundary ∂Ση, we set Dirichlet (clamped) or Neumann (traction-free) boundary conditions. Aiming to describe sequences of eigenvalues of order η−2 in the above problem, which correspond to oscillations of high frequencies of order η−1, we study an appropriately rescaled limit of the spectrum. Using a suitable notion of two-scale convergence for bounded operators acting on two-scale spaces, we show that the limiting spectrum consists of two parts: the Bloch (or band) spectrum and the ‘boundary’ spectrum. The latter corresponds to sequences of eigenvectors concentrating on the vertical boundaries of Ση, and is characterised by a problem set in a semi-infinite periodic strip with either clamped or stress-free boundary conditions. Based on the observation that some of the related eigenvalues can be found by solving an appropriate periodic-cell problem, we use modal methods to investigate finite-thickness semi-infinite waveguides. We compare our results with those for finite-thickness infinite waveguides given in Adams et al. (Proc. R. Soc. Lond. A, vol. 464, 2008, pp. 2669–2692). We also study infinite-thickness semi-infinite waveguides in order to gain insight into the finite-height analogue. We develop an asymptotic algorithm making use of the unimodular property of the modal method to demonstrate that in the weak contrast limit, and when wavenumber across the guide is fixed, there is at most one surface wave per gap in the spectrum. Using the monomode property of the waveguide we can consider the gap structure for the nth mode, when doing so, for traction-free boundaries, we find exactly one surface wave in each n-band gap.

A new approach is proposed for the description of thin viscous flows near the edges of a solid surface. For a steady flow, the lubrication approximation and the no-slip condition are assumed to be valid on most of the surface, except for relatively small neighbourhoods of the edges, where a universality principle is postulated: the behaviour of the liquid in these regions is universally determined by flux, external conditions and material properties. The resulting mathematical model is formulated as an ordinary differential equation involving the height of the liquid film and the flux as unknowns, and analytical results are outlined. The form of the universal functions which describe the behaviour in the edge regions is also discussed, obtaining conditions of compatibility with lubrication theory for small fluxes. Finally, an ordinary differential equation is introduced for the description of intermediate asymptotic profiles of a liquid film which flows off a bounded solid surface.

In this paper, we argue that quantitative empirical research to explain and predict criminal and related behaviour can benefit greatly from explicit theories of action linking individual and contextual factors in the causation of crime. Such theories foster a systematic selection of causal variables for data collection and hypothesis testing instead of a more indiscriminate accumulation of ‘risk factor’ correlates. Moreover, action theory encourages statistical modelling of crime causation beyond the most common linear regression. This paper illustrates both points by estimating two empirical models – a conventional logistic model and a Rasch model – on scenario response data concerning youth violence. The findings of this study show that the extent to which young people indicate a violent response to a provocation is dependent on their (law relevant) morality and ability to exercise self-control as well as the deterrent qualities (monitoring) of the setting.

We study equilibrium liquid crystal configurations in three-dimensional geometries, within the continuum Landau-de Gennes theory. We obtain explicit bounds for the equilibrium scalar order parameters in terms of the temperature and material-dependent constants. We explicitly quantify the temperature regimes where the Landau-de Gennes predictions match and the temperature regimes where the Landau-de Gennes predictions do not match the probabilistic second-moment definition of the Q-tensor order parameter. The regime of agreement may be interpreted as the regime of validity of the Landau-de Gennes theory since the Landau-de Gennes theory predicts large values of the equilibrium scalar order parameters – larger than unity, in the low-temperature regime. We discuss a modified Landau-de Gennes energy functional which yields physically realistic values of the equilibrium scalar order parameters in all temperature regimes.