We consider heat or mass transport from a circular cylinder under a uniform crossflow at small Reynolds numbers, $\mathrm{Re}\ll 1$. This problem has been thwarted in the past by limitations inherent in the classical analyses of the singular flow problem, which have used asymptotic expansions in inverse powers of $\log \mathrm{Re}$. We here make use of the hybrid approximation of Kropinski, Ward & Keller [(1995) SIAM J. Appl. Math. 55, 1484], based upon a robust asymptotic expansion in powers of $\mathrm{Re}$. In that approximation, the “inner” streamfunction is provided by the product of a pre-factor $S$, a slowly varying function of $\mathrm{Re}$, with a $\mathrm{Re}$-independent “canonical” solution of a simple mathematical form. The pre-factor, in turn, is determined as an implicit function of $\log \mathrm{Re}$ via asymptotic matching with a numerical solution of the nonlinear single-scaled “outer” problem, where the cylinder appears as a point singularity. We exploit the hybrid approximation to analyse the transport problem in the limit of large Péclet number, $\mathrm{Pe}\gg 1$. In that limit, transport is restricted to a narrow boundary layer about the cylinder surface – a province contained within the inner region of the flow problem. With $S$ appearing as a parameter, a similarity solution is readily constructed for the boundary-layer problem. It provides the Nusselt number as $0.5799(S\,\mathrm{Pe})^{1/3}$. This asymptotic prediction is in remarkably close agreement with that of the numerical solution of the exact problem [Dennis, Hudson & Smith (1968) Phys. Fluids 11, 933] even for moderate $\mathrm{Re}$-values.

Mathematical models of polyelectrolyte gels are often simplified by assuming the gel is electrically neutral. The rationale behind this assumption is that the thickness of the electric double layer (EDL) at the free surface of the gel is small compared to the size of the gel. Hence, the thin-EDL limit is taken, in which the thickness of the EDL is set to zero. Despite the widespread use of the thin-EDL limit, the solutions in the EDL are rarely computed and shown to match to the solutions for the electrically neutral bulk. The aims of this paper are to study the structure of the EDL and establish the validity of the thin-EDL limit. The model for the gel accounts for phase separation, which gives rise to diffuse interfaces with a thickness described by the Kuhn length. We show that the solutions in the EDL can only be asymptotically matched to the solutions for an electrically neutral bulk, in general, when the Debye length is much smaller than the Kuhn length. If the Debye length is similar to or larger than the Kuhn length, then phase separation can be initiated in the EDL. This phase separation spreads into the bulk of the gel and gives rise to electrically charged layers with different degrees of swelling. Thus, the thin-EDL limit and the assumption of electroneutrality only generally apply when the Debye length is much smaller than the Kuhn length.

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.

We consider two-dimensional mass transport to a finite absorbing strip in a uniform shear flow as a model of surface-based biosensors. The quantity of interest is the Sherwood number Sh, namely the dimensionless net flux onto the strip. Considering early-time absorption, it is a function of the Péclet number Pe and the Damköhler number Da, which, respectively, represent the characteristic magnitude of advection and reaction relative to diffusion. With a view towards modelling nanoscale biosensors, we consider the limit Pe«1. This singular limit is handled using matched asymptotic expansions, with an inner region on the scale of the strip, where mass transport is diffusively dominated, and an outer region at distances that scale as Pe-1/2, where advection enters the dominant balance. At the inner region, the mass concentration possesses a point-sink behaviour at large distances, proportional to Sh. A rescaled concentration, normalised using that number, thus possesses a universal logarithmic divergence; its leading-order correction represents a uniform background concentration. At the outer region, where advection by the shear flow enters the leading-order balance, the strip appears as a point singularity. Asymptotic matching with the concentration field in that region provides the Sherwood number as

$${\rm{Sh}} = {\pi \over {\ln (2/{\rm{P}}{{\rm{e}}^{1/2}}) + 1.0559 + \beta }},$$

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 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.

A phase field approach for structural topology optimization which allows for topologychanges and multiple materials is analyzed. First order optimality conditions arerigorously derived and it is shown via formally matched asymptoticexpansions that these conditions converge to classical first order conditions obtained inthe context of shape calculus. We also discuss how to deal with triple junctions wheree.g. two materials and the void meet. Finally, we present severalnumerical results for mean compliance problems and a cost involving the least square errorto a target displacement.

The stability of a one-spike solution to a general class of reaction-diffusion (RD) system with both regular and anomalous diffusion is analyzed. The method of matched asymptotic expansions is used to construct a one-spike equilibrium solution and to derive a nonlocal eigenvalue problem (NLEP) that determines the stability of this solution on an O(1) time-scale. For a particular sub-class of the reaction kinetics, it is shown that the discrete spectrum of this NLEP is determined in terms of the roots of certain simple transcendental equations that involve two key parameters related to the choice of the nonlinear kinetics. From a rigorous analysis of these transcendental equations by using a winding number approach and explicit calculations, sufficient conditions are given to predict the occurrence of Hopf bifurcations of the one-spike solution. Our analysis determines explicitly the number of possible Hopf bifurcation points as well as providing analytical formulae for them. The analysis is implemented for the shadow limit of the RD system defined on a finite domain and for a one-spike solution of the RD system on the infinite line. The theory is illustrated for two specific RD systems. Finally, in parameter ranges for which the Hopf bifurcation is unique, it is shown that the effect of sub-diffusion is to delay the onset of the Hopf bifurcation.

The slow dynamics and linearized stability of a two-spike quasi-equilibrium solution to ageneral class of reaction-diffusion (RD) system with and without sub-diffusion isanalyzed. For both the case of regular and sub-diffusion, the method of matched asymptoticexpansions is used to derive an ODE characterizing the spike locations in the absence ofany 𝒪(1) time-scale instabilities of the two-spike quasi-equilibrium profile. These fastinstabilities result from unstable eigenvalues of a certain nonlocal eigenvalue problem(NLEP) that is derived by linearizing the RD system around the two-spike quasi-equilibriumsolution. For a particular sub-class of the reaction kinetics, it is shown that thediscrete spectrum of this NLEP is determined by the roots of some simple transcendentalequations. From a rigorous analysis of these transcendental equations, explicit sufficientconditions are given to predict the occurrence of either Hopf bifurcations or competitioninstabilities of the two-spike quasi-equilibrium solution. The theory is illustrated forseveral specific choices of the reaction kinetics.

In this article, we consider a domain consisting of two cavities linked by a hole of small size. We derive a numerical method to compute an approximation of the eigenvalues of an elliptic operator without refining in the neighborhood of the hole. Several convergence rates are obtained and illustrated by numerical simulations.

In this paper, we will show how to obtain asymptotic solutions for the problem of pricing Asian options. Under the assumption that the underlying follows geometric Brownian motion, we will derive Taylor expansion series for the fixed and floating strike Asian options. While there will be no analytical formulae for calculating expansion coefficients, we will provide relatively simple algorithms for calculating them. The methodology is particularly effective for the case of continuously sampled fixed-strike Asian calls where it takes only seconds to obtain constants for the Taylor expansion series that can converge beyond 10 significant digits. It is needless to say that we need to calculate Taylor expansion constants only once and the option price would be an analytical expression constructed from a cumulative normal distribution function, an exponential function and finite sums.

Recent attention has focused on deriving localised pulse solutions to various systems of reaction–diffusion equations. In this paper, we consider the evolution of localised pulses in the Brusselator activator–inhibitor model, long considered a paradigm for the study of non-linear equations, in a finite one-dimensional domain with feed of the inhibitor through the boundary and global feed of the activator. We employ the method of matched asymptotic expansions in the limit of small activator diffusivity and small activator and inhibitor feeds. The disparity of diffusion lengths between the activator and inhibitor leads to pulse-type solutions in which the activator is localised while the inhibitor varies on an O(1) length scale. In the asymptotic limit considered, the pulses become spikes described by Dirac delta functions and evolve slowly in time until equilibrium is reached. Such quasi-equilibrium solutions with N activator pulses are constructed and a differential-algebraic system of equations (DAE) is derived, characterising the slow evolution of the locations and the amplitudes of the pulses. We find excellent agreement for the pulse evolution between the asymptotic theory and the results of numerical computations. An algebraic system for the equilibrium pulse amplitudes and locations is derived from the equilibrium points of the DAE system. Both symmetric equilibria, corresponding to a common pulse amplitude, and asymmetric pulse equilibria, for which the pulse amplitudes are different, are constructed. We find that for a positive boundary feed rate, pulse spacing of symmetric equilibria is no longer uniform, and that for sufficiently large boundary flux, pulses at the edges of the pattern may collide with and remain fixed at the boundary. Lastly, stability of the equilibrium solutions is analysed through linearisation of the DAE, which, in contrast to previous approaches, provides a quick way to calculate the small eigenvalues governing weak translation-type instabilities of equilibrium pulse patterns.