In this paper we prove some inequalities for finite sums and infinite series with positive terms. As an application of these results we obtain some inequalities for entropies of discrete probability distributions.

To assess rotational deformity in a broken forearm, an orthopaedic surgeon needs to determine the amount of rotation of the radius from one or more two-dimensional x-rays of the fracture. This requires only simple first-year university mathematics — rotational transformations of ellipses plus a little differential calculus — which yields a general formula giving the rotation angle from information obtained from an x-ray. Preliminary comparisons with experimental results are excellent. This is a practical problem that may be useful to motivate the teaching of conic sections.

An advanced pantograph-type partial differential equation, supplemented with initial and boundary conditions, arises in a model of asymmetric cell division. Methods for solving such problems are limited owing to functional (nonlocal) terms. The separation of variables entails an eigenvalue problem that involves a nonlocal ordinary differential equation. We discuss plausible eigenvalues that may yield nontrivial solutions to the problem for certain choices of growth and division rates of cells. We also consider the asymmetric division of cells with linear growth rate which corresponds to “exponential growth” and exponential rate of cell division, and show that the solution to the problem is a certain Dirichlet series. The distribution of the first moment of the biomass is shown to be unimodal.

A perturbation model is used to predict the distance jumped by a long-jumper for a range of tailwinds and headwinds. The zeroth-order approximation is based on gravity being the only force present, the effects of drag and lift only being included in the first-order corrections. The difference in predicted distances produced by the zeroth and first-order approximations is less than 2% for headwinds or tailwinds upto 4 ms−1. Most increases or decreases due to wind are caused by changes in the run-up speed, and consequently the take-off angle and speed.

Wave reflection by a wedge in deep water is examined, where the wedge can represent a breakwater of finite length or the bow of a ship heading directly into the waves. In addition, the form of the solution allows the results to apply to ships heading at an angle into the waves. We consider a deep-water wavetrain approaching the wedge head on from infinity and being reflected. Far from the wedge there is a field of progressive waves (the incident wavetrain) while close to the wedge there is a short-crested wavefield (the incident and reflected wavetrains). A weakly-nonlinear slowly-varying averaged Lagrangian theory is used to describe the problem (see Whitham [16]) as the theory includes the nonlinear interaction between the incident and reflected wavetrains. This modelling of a short-crested wavefield allows the nonlinear wavefield to be found for broad wedges, as opposed to previous theories which are applicable to thin wedges only.

It is shown that the governing partial differential equations are hyperbolic and that the solution comprises two regions, within which the wave properties are constant separated by a wave jump. Given the wedge angle and the incident wavefield, the jump angle and the wave steepness and wavenumber of the short-crested wave-field behind the wave jump can be determined. Two solution branches are found to exist: one corresponds to regular reflection, while for small amplitudes the other is similar to Mach-reflection and so it is called near Mach-reflection. Results are presented describing both solution branches and the transition between them.

Slow catalyst poisoning can result in the sudden failure of a chemical reactor operating isothermally with substrate-inhibited kinetics. At failure, a satisfactory steady state is exchanged for one of low conversion. The method of matched asymptotic expansions is used to give a detailed description of the exchange process in the phase plane. The structure of the jump is ascertained by separate asymptotic expansions across two adjoining transition regions in which the independent variables contain unknown shifts.

In this article, we generalise Newton's diagram method for finding small solutions ξ(λ) of equations f (ξ,λ) = 0 (0,0) = 0 with f analytic (see [1, 2, 4, 6]) to the case of a multi-dimensional function f, unknown variable ζ and small parameter λ. This method was briefly described in [1]. The method has many different applications and allows one to solve some inflexible problems. In particular, the method can be used in very difficult bifurcation problems, for example, for systems with small imperfections.

Iterative methods for solving systems of linear equations may be accelerated by coarse mesh rebalance techniques. The iterative technique, the Method of Implicit Non-stationary Iteration (MINI), is examined through a local-mode Fourier analysis and compared to relaxation techniques as a potential candidate for such acceleration. Results of a global-mode Fourier analysis for MINI, relaxation methods, and the conjugate gradient method are reported for two test problems.

An approximate solution is determined for the problem of scattering of water waves by a nearly vertical plate, by reducing it to two mixed boundary-value problems (BVP) for Laplace's equation, using perturbation techniques. While the solution of one of these BVP is well-known, the other BVPs is reduced to the problem of solving two uncoupled problems, and the complete solution of the problem under consideration up to first-order accuracy is derived with a special assumption on the shape of the plate and a related approximation. Known results involving the reflection and transmission coefficients are reproduced.

The main object of present paper is to obtain a finite summation of Srivastava's general triple hypergeometric series in terms of Kampé de Fériet's double hypergeometric series. A number of finite sums of Kampé de Fériet's double hypergeometric polynomials in terms of different kinds of single hypergeometric polynomials of higher order, are obtained. Some known results of Manocha and Sharma [9], [10], Munot [11], Pathan [12], Qureshi [15], Qureshi and Pathan [16] and Srivastava [26] are deduced as special cases. A result of Pathan [13, page 316 (1.2)] is also corrected here.

A synthesis is presented of two recent studies on modelling the nonlinear neuro-mechanical hearing processes in mosquitoes and in mammals. In each case, a hierarchy of models is considered in attempts to understand data that shows nonlinear amplification and compression of incoming sound signals. The insect’s hearing is tuned to the vicinity of a single input frequency. Nonlinear response occurs via an arrangement of many dual capacity neuro-mechanical units called scolopidia within the Johnston’s organ. It is shown how the observed data can be captured by a simple nonlinear oscillator model that is derived from homogenization of a more complex model involving a radial array of scolopidia. The physiology of the mammalian cochlea is much more complex, with hearing occurring via a travelling wave along a tapered, compartmentalized tube. Waves travel a frequency-dependent distance along the tube, at which point they are amplified and “heard”. Local models are reviewed for the pickup mechanism, within the outer hair cells of the organ of Corti. The current debate in the literature is elucidated, on the relative importance of two possible nonlinear mechanisms: active hair bundles and somatic motility. It is argued that the best experimental agreement can be found when the nonlinear terms include longitudinal coupling, the physiological basis of which is described. A discussion section summarizes the lessons learnt from both studies and attempts to shed light on the more general question of what constitutes a good mathematical model of a complex physiological process.

The effects on adsorption of the geometry of the solid may be studied through calculations based on a (distance)−ε (ε> 3) intermolecular potential. This paper establishes the result that the potential due to an infinitely long polygonal homogeneous solid prism, at position r in the plane of its right section, is – . Here ρi = ∣ r − ri ∣, where the ri are the position vectors of the n vertices of the polygon, and θij are the angles r − ri makes with the two sides of the polygon which meet at vertex ri. The g's are exact functions of θij. They are, in general, integrals of associated Legendre functions, but they are elementary for ε an even integer. A similar result holds for the potential within an infinitely long polygonal prismatic cavity. The analysis involves a systematic superposition schema and the concept of a supplementary potential with datum within the solid at infinity. The cases ε = 6 and ε = 4 are treated in detail and illustrative solutions given for the following configurations: semi-infinite laminae, deep rectangular cracks, square prisms, square prismatic cavities and regular n-gonal prismatic cavities.

We consider an optimal control problem with, possibly time-dependent, constraints on state and control variables, jointly. Using only elementary methods, we derive a sufficient condition for optimality. Although phrased in terms reminiscent of the necessary condition of Pontryagin, the sufficient condition is logically independent, as can be shown by a simple example.

This is a review of thin-body and slender-body theories, with indications of some new applications. Topics discussed include bodies with near-constant surface pressure, subsonic and supersonic aerodynamics, ship hydrodynamics, slender bodies in Stokes flow, slender footings in elastic media, and slender moonpools. Mathematical features of the thin- and slender-body approximations are also discussed, especially nonlocal convolution terms modelling three-dimensionality in the otherwise two-dimensional near field, end effects, and the role of the logarithm of the slenderness ratio. This review was presented by the first author as the IMA Lighthill Memorial Lecture at the British Applied Mathematics Colloquium (BAMC) 2004.

We present new exact solutions for the flow of liquid during constant-rate expression from a finite thickness of liquid-saturated porous material with nonlinear properties. By varying a single nonlinearity parameter and a dimensionless expression rate, we systematically investigate the effect of nonlinearity and of an impermeable barrier (e.g. a piston). We illustrate the water profile shape and the water ratio deficit at the expression surface (e.g. a filter membrane) as a function of time.

In this paper, we study the existence of positive solutions for the one-dimensional p-Laplacian differential equation, subject to the multipoint boundary condition by applying a monotone iterative method.

Rayleigh–Taylor instability occurs when a heavier fluid overlies a lighter fluid, and the two seek to exchange positions under the effect of gravity. We present a linearized theory for arbitrary three-dimensional (3D) initial disturbances that grow in time, and calculate the evolution of the interface for early times. A new spectral method is introduced for the fully 3D nonlinear problem in a Boussinesq fluid, where the interface between the light and heavy fluids is approximated with a smooth but rapid density change in the fluid. The results of large-scale numerical calculation are presented in fully 3D geometry, and compared and contrasted with the early-time linearized theory.

Spatial data sets can be analysed by counting the number of objects in equally sized bins. The bin counts are related to the Pólya urn process, where coloured balls (for example, white or black) are removed from the urn at random. If there are insufficient white or black balls for the prescribed number of trials, the Pólya urn process becomes untenable. In this case, we modify the Pólya urn process so that it continues to describe the removal of volume within a spatial distribution of objects. We determine when the standard formula for the variance of the standard Pólya distribution gives a good approximation to the true variance. The variance quantifies an index for assessing whether a spatial point data set is at its most randomly distributed state, called the complete spatial randomness (CSR) state. If the bin size is an order of magnitude larger than the size of the objects, then the standard formula for the CSR limit is indicative of when the CSR state has been attained. For the special case when the object size divides the bin size, the standard formula is in fact exact.

New polynomial solutions of the Navier-Stokes equations for steady uni-directional flow of viscous incompressible fluid, with a free surface, down inclined channels of specialized cross-section are considered. An inverse method is uded to obtain the geometrical shape of the channel by equating the polynomials solution to zero (i.e. the no-slip condition) and thence determining the boundary shape.

A small gas bubble in a liquid, when driven by intense ultrasound, collapses and emits light in a process called Single-Bubble Sonoluminescence (SBSL). While the dynamics of driven bubbles are well studied, less is known of the physical conditions in the gas or whether it is necessary to include ionisation in simpler studies of bubble dynamics. In this study, a model was derived from Rayleigh-Plesset dynamics, a van der Waals equation of state and the first law of thermodynamics (including interfacial heat transfer and ionisation). Stronger model ionisation reduced the maximum collapse temperature, and altered other collapse characteristics. Chaotic parameter regions are proximal to, but not coincident with, known stable SL regions. Resonant behaviour was only markedly affected by ionisation close to these chaotic regions.