Discrete cilia modelling with singularity distributions

We discuss in detail techniques for modelling flows due to finite and infinite arrays of beating cilia. An efficient technique, based on concepts from previous ‘singularity models’ is described, that is accurate in both near and far-fields. Cilia are modelled as curved slender ellipsoidal bodies by distributing Stokeslet and potential source dipole singularities along their centrelines, leading to an integral equation that can be solved using a simple and efficient discretisation. The computed velocity on the cilium surface is found to compare favourably with the boundary condition. We then present results for two topics of current interest in biology. 1) We present the first theoretical results showing the mechanism by which rotating embryonic nodal cilia produce a leftward flow by a ‘posterior tilt,’ and track particle motion in an array of three simulated nodal cilia. We find that, contrary to recent suggestions, there is no continuous layer of negative fluid transport close to the ciliated boundary. The mean leftward particle transport is found to be just over 1 μm/s, within experimentally measured ranges. We also discuss the accuracy of models that represent the action of cilia by steady rotlet arrays, in particular, confirming the importance of image systems in the boundary in establishing the far-field fluid transport. Future modelling may lead to understanding of the mechanisms by which morphogen gradients or mechanosensing cilia convert a directional flow to asymmetric gene expression. 2) We develop a more complex and detailed model of flow patterns in the periciliary layer of the airway surface liquid. Our results confirm that shear flow of the mucous layer drives a significant volume of periciliary liquid in the direction of mucus transport even during the recovery stroke of the cilia. Finally, we discuss the advantages and disadvantages of the singularity technique and outline future theoretical and experimental developments required to apply this technique to various other biological problems, particularly in the reproductive system

The Metapopulation Dynamics of an Infectious Disease: Tuberculosis in Possums

An SEI metapopulation model is developed for the spread of an infectious agent by migration. The model portrays two age classes on a number of patches connected by migration routes which are used as host animals mature. A feature of this model is that the basic reproduction ratio may be computed directly, using a scheme that separates topography, demography, and epidemiology. We also provide formulas for individual patch basic reproduction numbers and discuss their connection with the basic reproduction ratio for the system. The model is applied to the problem of spatial spread of bovine tuberculosis in a possum population. The temporal dynamics of infection are investigated for some generic networks of migration links, and the basic reproduction ratio is computed—its value is not greatly different from that for a homogeneous model. Three scenarios are considered for the control of bovine tuberculosis in possums where the spatial aspect is shown to be crucial for the design of disease management operation

The metapopulation dynamics of an infectious disease: a case study of tuberculosis in possums

An SEI metapopulation model is developed for the spread of an infectious agent by migration. The model portrays two age classes on a number of patches connected by migration routes which are used as host animals mature. A feature of this model is that the basic reproduction ratio may be computed directly, using a scheme that separates topography, demography, and epidemiology. We also provide formulas for individual patch basic reproduction numbers and discuss their connection with the basic reproduction ratio for the system. The model is applied to the problem of spatial spread of bovine tuberculosis in a possum population. The temporal dynamics of infection are investigated for some generic networks of migration links, and the basic reproduction ratio is computed—its value is not greatly different from that for a homogeneous model. Three scenarios are considered for the control of bovine tuberculosis in possums where the spatial aspect is shown to be crucial for the design of disease management operations

A stochastic model for MRSA transmission within a hospital ward incorporating environmental contamination

Methicillin-resistant Staphylococcus aureus (MRSA) transmission in hospital wards is associated with adverse outcomes for patients and increased costs for hospitals. The transmission process is inherently stochastic and the randomness emphasized by the small population sizes involved. As such, a stochastic model was proposed to describe the MRSA transmission process, taking into account the related contribution and modelling of the associated microbiological environmental contamination. The model was used to evaluate the performance of five common interventions and their combinations on six potential outcome measures of interest under two hypothetical disease burden settings. The model showed that the optimal intervention combination varied depending on the outcome measure and burden setting. In particular, it was found that certain outcomes only required a small subset of targeted interventions to control the outcome measure, while other outcomes still reported reduction in the outcome distribution with up to all five interventions included. This study describes a new stochastic model for MRSA transmission within a ward and highlights the use of the generalized Mann–Whitney statistic to compare the distribution of the outcome measures under different intervention combinations to assist in planning future interventions in hospital wards under different potential outcome measures and disease burden

Dynamics of tear film deposition and draining

This paper investigates the deposition of the tear film on the cornea of the human eye. The tear film is laid down by the motion of the upper eyelid and then subsequently flows and thins. Of particular interest is the stability of the tear layer and the development of dry patches on the cornea. While there has been significant research on the behaviour of tear films between blinks, this paper focuses on understanding the mechanisms which control the shape and thickness of the deposited film and how this affects the subsequent film behaviour. Numerical and analytical methods are applied to a lubrication model which includes the effects of surface tension, viscosity, gravity and evaporation. The model reveals the importance of the eyelid velocity, motion of the surface lipid layer and the storage of tear film between blinks. Keywords: eye; tear film; film deposition; film rupture