In this special issue of the European Journal of Applied Mathematics we celebrate the 75th birthday of its founding editor, Professor John Ockendon.

John started his research career with his D.Phil., awarded in 1965, working with Alan Tayler at the University of Oxford. John's background was in fluid mechanics, covering a wide range of topics, including, lubrication flows, flows in porous media, hypersonic aerodynamics, sloshing phenomena, and flow separation. He soon became involved in industrial research which led, in particular, to his interest in the field of free boundary problems, an area in which he became a driving force, through papers, conferences, and particularly through the fostering of regular international meetings and personal interactions. His pioneering work in the field was particularly concentrated in diffusion-controlled moving boundary problems that arise in various applications, such as in phase change problems, elastic-contact problems and the fundamental Hele-Shaw free boundary problem. The industrial collaborations also led to mathematical studies of lens design, fibre extrusion, fluidised beds, glass manufacture and semiconductor devices, to name but a few, and to problems in scattering and wave theory, non-linear wave propagation and non-linear oscillations, as well as mechanics and heat and mass transfer. He continues to work on PDEs and asymptotics applied to new well-founded models from the real world, with a particular focus in the physical and earth sciences.

We consider Hele-Shaw flows driven by injection of a highly shear-thinning power-law fluid (of exponent n) in the absence of surface tension. We formulate the problem in terms of the streamfunction ψ, which satisfies the p-Laplacian equation ∇·(|∇ψ|p−2∇ψ) = 0 (with p = (n+1)/n) and use the method of matched asymptotic expansions in the large n (extreme-shear-thinning) limit to find an approximate solution. The results show that significant flow occurs only in (I) segments of a (single) circle centred on the injection point, whose perimeters comprise the portion of free boundary closest to the injection point and (II) an exponentially small region around the injection point and (III) a transition region to the rest of the fluid: while the flow in the latter is exponentially slow it can be characterised in detail.

Asymptotic homogenisation via the method of multiple scales is considered for problems in which the microstructure comprises inclusions of one material embedded in a matrix formed from another. In particular, problems are considered in which the interface conditions include a global balance law in the form of an integral constraint; this may be zero net charge on the inclusion, for example. It is shown that for such problems care must be taken in determining the precise location of the interface; a naive approach leads to an incorrect homogenised model. The method is applied to the problems of perfectly dielectric inclusions in an insulator, and acoustic wave propagation through a bubbly fluid in which the gas density is taken to be negligible.

The superheating that usually occurs when a solid is melted by volumetric heating can produce irregular solid–liquid interfaces. Such interfaces can be visualised in ice, where they are sometimes known as Tyndall stars. This paper describes some of the experimental observations of Tyndall stars and a mathematical model for the early stages of their evolution. The modelling is complicated by the strong crystalline anisotropy, which results in an anisotropic kinetic undercooling at the interface; it leads to an interesting class of free boundary problems that treat the melt region as infinitesimally thin.

We consider a coating flow of nematic liquid crystal (NLC) fluid film on an inclined substrate. Exploiting the small aspect ratio in the geometry of interest, a fourth-order nonlinear partial differential equation is used to model the free surface evolution. Particular attention is paid to the interplay between the bulk elasticity and the anchoring conditions at the substrate and free surface. Previous results have shown that there exist two-dimensional travelling wave solutions that translate down the substrate. In contrast to the analogous Newtonian flow, such solutions may be unstable to streamwise perturbations. Extending well-known results for Newtonian flow, we analyse the stability of the front with respect to transverse perturbations. Using full numerical simulations, we validate the linear stability theory and present examples of downslope flow of nematic liquid crystal in the presence of both transverse and streamwise instabilities.

We investigate the effect of mass transfer on the evolution of a thin, two-dimensional, partially wetting drop. While the effects of viscous dissipation, capillarity, slip and uniform mass transfer are taken into account, other effects, such as gravity, surface tension gradients, vapour transport and heat transport, are neglected in favour of mathematical tractability. Our focus is on a matched-asymptotic analysis in the small-slip limit, which reveals that the leading-order outer formulation and contact-line law depend delicately on both the sign and the size of the mass transfer flux. This leads, in particular, to novel generalisations of Tanner's law. We analyse the resulting evolution of the drop on the timescale of mass transfer and validate the leading-order predictions by comparison with preliminary numerical simulations. Finally, we outline the generalisation of the leading-order formulations to prescribed non-uniform rates of mass transfer and to three dimensions.

We consider the steady two-dimensional thin-film version of a problem concerning a weightless non-isothermal free fluid film subject to thermocapillarity, proposed and analysed by Pukhnachev and co-workers. Specifically, we extend and correct the paper by Karabut and Pukhnachev (J. Appl. Mech. Tech. Phys. 49, 568–579, 2008), in which the problem is solved numerically, and in which it is claimed that there exists a unique solution for any value of a prescribed heat-flux parameter in the model. We present a closed-form (parametric) solution of the problem, and from this show that, on the contrary, solutions exist only when the heat-flux parameter is less than a critical value found numerically by Karabut and Pukhnachev, and that when this condition is satisfied there are in fact two solutions, one of which recovers that obtained numerically by Karabut and Pukhnachev, the other being new.

In this paper, we consider the straining flow of a weakly interacting polymer–surfactant solution below a free surface, with the bulk surfactant concentration above the critical micelle concentration. We formulate a set of coupled differential equations describing the concentration of monomers, micelles, polymer, and polymer–micelle aggregates in the flow. We analyse the model in several asymptotic limits, and make predictions about the distribution of each of the species. In particular, in the large-reaction-rate limit we find that the model predicts a region near the free surface where no micelles or aggregates are present, and beneath this a region where the concentration of surfactant is constant, across which the concentration of aggregates increases until all the free polymer is consumed. For certain parameter regimes, a maximum in the concentration of the polymer–micelle complex occurs within the bulk fluid. In the finite-reaction-rate limit, micelles, and aggregates are present right up to the free surface, and the plateau in the concentration of surfactant in the bulk is no longer present. Results from the asymptotic theory compare favorably with full numerical solutions.

The hybrid numerical-asymptotic (HNA) approach aims to reduce the computational cost of conventional numerical methods for high-frequency wave scattering problems by enriching the numerical approximation space with oscillatory basis functions, chosen based on partial knowledge of the high-frequency solution asymptotics. In this paper, we propose a new methodology for the treatment of shadow boundary effects in HNA boundary element methods, using the classical geometrical theory of diffraction phase functions combined with mesh refinement. We develop our methodology in the context of scattering by a class of sound-soft non-convex polygons, presenting a rigorous numerical analysis (supported by numerical results) which proves the effectiveness of our HNA approximation space at high frequencies. Our analysis is based on a study of certain approximation properties of the Fresnel integral and related functions, which govern the shadow boundary behaviour.

The theoretical investigation here of a three-dimensional array of jets of fluid (air guns) and their interference is motivated by applications to the food sorting industry especially. Three-dimensional motion without symmetry is addressed for arbitrary jet cross-sections and incident velocity profiles. Asymptotic analysis based on the comparatively long axial length scale of the configuration leads to a reduced longitudinal vortex system providing a slender flow model for the complete array response. Analytical and numerical studies, along with comparisons and asymptotic limits or checks, are presented for various cross-sectional shapes of nozzle and velocity inputs. The influences of swirl and of unsteady jets are examined. Substantial cross-flows are found to occur due to the interference. The flow solution is non-periodic in the cross-plane even if the nozzle array itself is periodic. The analysis shows that in general the bulk of the three-dimensional motion can be described simply in a cross-plane problem but the induced flow in the cross-plane is sensitively controlled by edge effects and incident conditions, a feature which applies to any of the array configurations examined. Interference readily alters the cross-flow direction and misdirects the jets. Design considerations centre on target positioning and jet swirling.