Gliomas are very aggressive brain tumours, in which tumour cells gain the ability to penetrate the surrounding normal tissue. The invasion mechanisms of this type of tumour remain to be elucidated. Our work is motivated by the migration/proliferation dichotomy (go-or-grow) hypothesis, i.e. the antagonistic migratory and proliferating cellular behaviours in a cell population, which may play a central role in these tumours. In this paper, we formulate a simple go-or-grow model to investigate the dynamics of a population of glioma cells for which the switch from a migratory to a proliferating phenotype (and vice versa) depends on the local cell density. The model consists of two reaction-diffusion equations describing cell migration, proliferation and a phenotypic switch. We use a combination of numerical and analytical techniques to characterize the development of spatio-temporal instabilities and travelling wave solutions generated by our model. We demonstrate that the density-dependent go-or-grow mechanism can produce complex dynamics similar to those associated with tumour heterogeneity and invasion.

In-plane vibrations of an elastic rod clamped at both extremities are studied. The rod is modeled as an extensible planar Kirchhoff elastic rod under large displacements and rotations. Equilibrium configurations and vibrations around these configurations are computed analytically in the incipient post-buckling regime. Of particular interest is the variation of the first mode frequency as the load is increased through the buckling threshold. The loading type is found to have a crucial importance as the first mode frequency is shown to behave singularly in the zero thickness limit in the case of prescribed axial displacement, whereas a regular behavior is found in the case of prescribed axial load. © 2013 Elsevier Ltd.

We investigate the deformation of a thin elastic sheet floating on a liquid surface and subject to a uniaxial compression. We show that at a critical compression the sheet delaminates from the liquid over a finite region forming a delamination "blister." This blistering regime adds to the wrinkling and localized folding regimes that have been studied previously. The transition from wrinkled to blistered states occurs when delamination becomes energetically favorable compared with wrinkling. We determine the initial blister size and the evolution of blister size with continuing compression before verifying our theoretical results with experiments at a macroscopic scale. © 2011 American Physical Society.

Prior research into the effect of convection on steady-state mass transfer from a spherical particle embedded in a porous medium has used the Darcy model to describe the flow. However, a limitation of the Darcy model is that it does not account for viscous effects near boundaries. Brinkman modified the Darcy model to include these effects by introducing an extra viscous term. Here we investigate the impact of this extra viscous term on the steady-state mass transfer from a sphere at low Péclet number, Pe 1. We use singular perturbation techniques to find the approximate asymptotic solution for the concentration profile. Mass-release from the surface of the sphere is described by a Robin boundary condition, which represents a first-order chemical reaction. We find that a larger Brinkman viscous boundary layer renders mass transport by convection less effective, and reduces the asymmetry in the peri-sphere concentration profiles. We provide simple analytical expressions that can be used to calculate the concentration profiles, as well as the local and average Sherwood numbers; and comparison to numerical simulations verifies the order of magnitude of the error in the asymptotic expansions. In the appropriate limits, the asymptotic results agree with solutions previously obtained for Stokes and Darcy flow. The solution for Darcy flow with a Robin boundary condition has not been considered previously in the literature and is a new result. Whilst the article has been formulated in terms of mass transfer, the analysis is also applicable to heat transfer, with concentration replaced by temperature and the Sherwood number by the Nusselt number. © 2013 Elsevier Ltd. All rights reserved.

© 2012 Society for Industrial and Applied Mathematics. Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration of stochastic differential equations. This approach is illustrated with the constructions of new methods of weak order two, in particular, semi-implicit integrators well suited for stiff (meansquare stable) stochastic problems, and implicit integrators that exactly conserve all quadratic first integrals of a stochastic dynamical system. Numerical examples confirm the theoretical results and show the versatility of our methodology.

The centrifugally driven viscous fingering problem arises when two immiscible fluids of different densities flow in a rotating Hele-Shaw cell. In this conventional setting an interplay between capillary and centrifugal forces makes the fluid-fluid interface unstable, leading to the formation of fingered structures that compete dynamically and reach different lengths. In this context, it is known that finger competition is very sensitive to changes in the viscosity contrast between the fluids. We study a variant of such a rotating flow problem where the fluids react and produce a gellike phase at their separating boundary. This interface is assumed to be elastic, presenting a curvature-dependent bending rigidity. A perturbative weakly nonlinear approach is used to investigate how the elastic nature of the interface affects finger competition events. Our results unveil a very different dynamic scenario, in which finger length variability is not regulated by the viscosity contrast, but rather determined by two controlling quantities: a characteristic radius and a rigidity fraction parameter. By properly tuning these quantities one can describe a whole range of finger competition behaviors even if the viscosity contrast is kept unchanged. © 2014 American Physical Society.

This review concentrates on the fluctuations of the velocities of sedimenting spheres, and on the structural instability of a suspension of settling fibers. For many years, theoretical estimates and numerical simulations predicted the fluctuations of the velocities of spheres to increase with the size of the container, whereas experiments found no such variation. Two ideas have increased our understanding. First, the correlation length of the velocity fluctuations was found experimentally to be 20 interparticle separations. Second, in dilute suspensions, a vertical variation in the concentration due to the spreading of the front with the clear fluid can inhibit the velocity fluctuations. In a very dilute regime, a homogeneous suspension of fibers suffers a spontaneous instability in which fast descending fiber-rich columns are separated by rising fiber-sparse columns. In a semidilute regime, the settling is hindered, more so than for spheres. © 2011 by Annual Reviews. All rights reserved.

Mathematical models in biology and physiology are often represented by large systems of non-linear ordinary differential equations. In many cases, an observed behaviour may be written as a linear functional of the solution of this system of equations. A technique is presented in this study for automatically identifying key terms in the system of equations that are responsible for a given linear functional of the solution. This technique is underpinned by ideas drawn from a posteriori error analysis. This concept has been used in finite element analysis to identify regions of the computational domain and components of the solution where a fine computational mesh should be used to ensure accuracy of the numerical solution. We use this concept to identify regions of the computational domain and components of the solution where accurate representation of the mathematical model is required for accuracy of the functional of interest. The technique presented is demonstrated by application to a model problem, and then to automatically deduce known results from a cell-level cardiac electrophysiology model. © 2010 Elsevier Inc.

The beat patterns of mammalian sperm flagella can be categorised into two different types. The first involves symmetric waves propagating down the flagellum with a net linear propulsion of the sperm cell. The second, hyperactive, waveform is classified by vigorous asymmetric waves of higher amplitude, lower wavenumber and frequency propagating down the flagellum resulting in highly curved trajectories. The latter beat pattern is part of the capacitation process whereby sperm prepare for the prospective penetration of the zona pellucida and fusion with the egg. Hyperactivation is often observed to initiate as sperm escape from epithelial and ciliary bindings formed within the isthmic regions of the female oviducts, leading to a conjecture in the literature that this waveform is mechanically important for sperm escape. Hence, we explore the mechanical effects of hyperactivation on a tethered sperm, focussing on a Newtonian fluid. Using a resistive force theory model we demonstrate that hyperactivation can indeed generate forces that pull the sperm away from a tethering point and consequently a hyperactivated sperm cell bound to an epithelial surface need not always be pushed by its flagellum. More generally, directions of the forces generated by tethered flagella are insensitive to reductions in beat frequency and the detailed flagellar responses depend on the nature of the binding at the tethering point. Furthermore, waveform asymmetry and amplitude increases enhance the tendency for a tethered flagellum to start tugging on its binding. The same is generally predicted to be true for reductions in the wavenumber of the flagellum beat, but not universally so, emphasising the dynamical complexity of flagellar force generation. Finally, qualitative observations drawn from experimental data of human sperm bound to excised female reproductive tract are also presented and are found to be consistent with the theoretical predictions. © 2012 Elsevier Ltd.

The invasion of pest insects often changes or destroys a native ecosystem, and can result in food shortages and disease endemics. Issues such as the environmental effects of chemical control methods, the economic burden of maintaining control strategies and the risk of pest resistance still remain, and mosquito-borne diseases such as malaria and dengue fever prevail in many countries, infecting over 100 million worldwide in 2010. One environmentally friendly method for mosquito control is the Sterile Insect Technique (SIT). This species-specific method of insect control relies on the mass rearing, sterilization and release of large numbers of sterile insects. An alternative transgenic method is the Release of Insects carrying a Dominant Lethal (RIDL). Our objective is to consider contrasting control strategies for two invasive scenarios via SIT and RIDL: an endemic case and an emerging outbreak. We investigate how the release rate and size of release region influence both the potential for control success and the resources needed to achieve it, under a range of conditions and control strategies, and we discuss advantageous strategies with respect to reducing the release resources and strategy costs (in terms of control mosquito numbers) required to achieve complete eradication of wild-type mosquitoes. © 2013 Elsevier Ltd.