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.

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.

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.

A theory for the dynamics and statics of growing elastic rods is presented. First, a single growing rod is considered and the formalism of three-dimensional multiplicative decomposition of morphoelasticity is used to describe the bulk growth of Kirchhoff elastic rods. Possible constitutive laws for growth are discussed and analysed. Second, a rod constrained or glued to a rigid substrate is considered, with the mismatch between the attachment site and the growing rod inducing stress. This stress can eventually lead to instability, bifurcation, and buckling. © 2012 Elsevier Ltd. All rights reserved.

© 2015 Elsevier Inc. In this paper, three multiscale methods for coupling of mesoscopic (compartment-based) and microscopic (molecular-based) stochastic reaction-diffusion simulations are investigated. Two of the three methods that will be discussed in detail have been previously reported in the literature; the two-regime method (TRM) and the compartment-placement method (CPM). The third method that is introduced and analysed in this paper is called the ghost cell method (GCM), since it works by constructing a "ghost cell" in which molecules can disappear and jump into the compartment-based simulation. Presented is a comparison of sources of error. The convergent properties of this error are studied as the time step δ. t (for updating the molecular-based part of the model) approaches zero. It is found that the error behaviour depends on another fundamental computational parameter h, the compartment size in the mesoscopic part of the model. Two important limiting cases, which appear in applications, are considered:. (i)δt→0 and h is fixed;(ii)δt→0 and h→0 such that δt/h is fixed. The error for previously developed approaches (the TRM and CPM) converges to zero only in the limiting case (ii), but not in case (i). It is shown that the error of the GCM converges in the limiting case (i). Thus the GCM is superior to previous coupling techniques if the mesoscopic description is much coarser than the microscopic part of the model.

A cylindrical elastic tube under uniform radial external pressure will buckle circumferentially to a non-circular cross-section at a critical pressure. The buckling represents an instability of the inner or outer edge of the tube. This is a common phenomenon in biological tissues, where it is referred to as mucosal folding. Here, we investigate this buckling instability in a growing elastic tube. A change in thickness due to growth can have a dramatic impact on circumferential buckling, both in the critical pressure and the buckling pattern. We consider both single- and bi-layer tubes and multiple boundary conditions. We highlight the competition between geometric effects, i.e. the change in tube dimensions, and mechanical effects, i.e. the effect of residual stress, due to differential growth. This competition can lead to non-intuitive results, such as a tube growing to be thinner and yet buckle at a higher pressure. © 2011 Elsevier Ltd. All rights reserved.

We consider the efficient solution of the Cahn-Hilliard variational inequality using an implicit time discretization, which is formulated as an optimal control problem with pointwise constraints on the control. By applying a semi-smooth Newton method combined with a Moreau-Yosida regularization technique for handling the control constraints we show superlinear convergence in function space. At the heart of this method lies the solution of large and sparse linear systems for which we propose the use of preconditioned Krylov subspace solvers using an effective Schur complement approximation. Numerical results illustrate the competitiveness of this approach. © 2014 Elsevier Inc.

A mathematically challenging model of dynamic wetting as a process of interface formation has been, for the first time, fully incorporated into a numerical code based on the finite element method and applied, as a test case, to the problem of capillary rise. The motivation for this work comes from the fact that, as discovered experimentally more than a decade ago, the key variable in dynamic wetting flows - the dynamic contact angle - depends not just on the velocity of the three-phase contact line but on the entire flow field/geometry. Hence, to describe this effect, it becomes necessary to use the mathematical model that has this dependence as its integral part. A new physical effect, termed the 'hydrodynamic resist to dynamic wetting', is discovered where the influence of the capillary's radius on the dynamic contact angle, and hence on the global flow, is computed. The capabilities of the numerical framework are then demonstrated by comparing the results to experiments on the unsteady capillary rise, where excellent agreement is obtained. Practical recommendations on the spatial resolution required by the numerical scheme for a given set of non-dimensional similarity parameters are provided, and a comparison to asymptotic results available in limiting cases confirms that the code is converging to the correct solution. The appendix gives a user-friendly step-by-step guide specifying the entire implementation and allowing the reader to easily reproduce all presented results, including the benchmark calculations. © 2012 Elsevier Inc.