Error estimation and adaptivity for incompressible, non–linear (hyper–)elasticity

A Galerkin finite element method is developed for non–linear, incompressible (hyper) elasticity, and a posteriori error estimates are derived for both linear functionals of the solution and linear functionals of the stress on a boundary where Dirichlet boundary conditions are applied. A second, higher order method for calculating a linear functional of the stress on a Dirichlet boundary is also presented together with an a posteriori error estimator for this approach. An implementation for a 2D model problem with known solution demonstrates the accuracy of the error estimators. Finally the a posteriori error estimate is shown to provide a basis for effective mesh adaptivity

Cardiac Electromechanics: The effect of contraction model on the mathematical problem and accuracy of the numerical scheme

Models of cardiac electromechanics usually contain a contraction model determining the active tension induced at the cellular level, and the equations of nonlinear elasticity to determine tissue deformation in response to this active tension. All contraction models are dependent on cardiac electro-physiology, but can also be dependent on\ud the stretch and stretch-rate in the fibre direction. This fundamentally affects the mathematical problem being solved, through classification of the governing PDEs, which affects numerical schemes that can be used to solve the governing equations. We categorise contraction models into three types, and for each consider questions such as classification and the most appropriate choice from two numerical methods (the explicit and implicit schemes). In terms of mathematical classification, we consider the question of strong ellipticity of the total strain energy (important for precluding ‘unnatural’ material behaviour) for stretch-rate-independent contraction models; whereas for stretch-rate-dependent contraction models we introduce a corresponding third-order problem and explain how certain choices of boundary condition could lead to constraints on allowable initial condition. In terms of suitable numerical methods, we show that an explicit approach (where the contraction model is integrated in the timestep prior to the bulk deformation being computed) is: (i) appropriate for stretch-independent contraction models; (ii) only conditionally-stable, with the stability criterion independent of timestep, for contractions models which just depend on stretch (but not stretch-rate), and (iii) inappropriate for stretch-rate-dependent models

Error bounds on block Gauss Seidel solutions of coupled\ud multiphysics problems

Mathematical models in many fields often consist of coupled sub–models, each of which describe a different physical process. For many applications, the quantity of interest from these models may be written as a linear functional of the solution to the governing equations. Mature numerical solution techniques for the individual sub–models often exist. Rather than derive a numerical solution technique for the full coupled model, it is therefore natural to investigate whether these techniques may be used by coupling in a block Gauss–Seidel fashion. In this study, we derive two a posteriori bounds for such linear functionals. These bounds may be used on each Gauss–Seidel iteration to estimate the error in the linear functional computed using the single physics solvers, without actually solving the full, coupled problem. We demonstrate the use of the bound first by using a model problem from linear algebra, and then a linear ordinary differential equation example. We then investigate the effectiveness of the bound using a non–linear coupled fluid–temperature problem. One of the bounds derived is very sharp for most linear functionals considered, allowing us to predict very accurately when to terminate our block Gauss–Seidel iteration.\ud \ud Copyright c 2000 John Wiley & Sons, Ltd

Nonuniqueness in a minimal model for cell motility

Two–phase flow models have been used previously to model cell motility, however these have rapidly become very complicated, including many physical processes, and are opaque. Here we demonstrate that even the simplest one–dimensional, two–phase, poroviscous, reactive flow model displays a number of behaviours relevant to cell crawling. We present stability analyses that show that an asymmetric perturbation is required to cause a spatially uniform, stationary strip of cytoplasm to move, which is relevant to cell polarization. Our numerical simulations identify qualitatively distinct families of travelling–wave solution that co–exist at certain parameter values. Within each family, the crawling speed of the strip has a bell–shaped dependence on the adhesion strength. The model captures the experimentally observed behaviour that cells crawl quickest at intermediate adhesion strengths, when the substrate is neither too sticky nor too slippy

Knowledge and the artefact

This paper discusses ways that knowledge may be found in or through artefacts. One purpose is to suggest situations where artefacts might be central to a narrative, rather than secondary to a text. A second purpose is to suggest ways that design and production of artefacts might be instrumental in eliciting knowledge. Four general situations are proposed: (1) Simple Forms - an artefact demonstrates or describes a principle or technique. (2) Communication of Process - artefacts arising from a process make the process explicit. (3) Artefacts Within the Research - artefacts are instrumental in advancing the research by communicating ideas or information. (4) Knowledge Elicited by Artefacts - artefacts provide a stimulus or context which enables information to be uncovered. .</p

The rigidity of periodic body-bar frameworks on the three-dimensional fixed torus

We present necessary and sufficient conditions for the generic rigidity of body-bar frameworks on the three-dimensional fixed torus. These frameworks correspond to infinite periodic body-bar frameworks in $\mathbb{R}^3$ with a fixed periodic lattice.Comment: 31 pages, 12 figure

A note on heat and mass transfer from a sphere in Stokes\ud flow at low Péclet number

We consider the low Péclet number, Pe ≪ 1, asymptotic solution for steady-state heat and mass transfer from a sphere immersed in Stokes flow with a Robin boundary condition on its surface, representing Newton cooling or a first-order chemical reaction. The application of van Dyke’s rule up to terms of O(Pe3) shows that the O(Pe3 log Pe) terms in the expression for the average Nusselt/Sherwood number are double those previously derived in the literature. Inclusion of the O(Pe3) terms is shown to increase significantly the range of validity of the expansion

Comparison of bacterioneuston and bacterioplankton dynamics during a phytoplankton bloom in a fjord mesocosm

The bacterioneuston is the community of Bacteria present in surface microlayers, the thin surface film that forms the interface between aquatic environments and the atmosphere. In this study we compared bacterial cell abundance and bacterial community structure of the bacterioneuston and the bacterioplankton (from the subsurface water column) during a phytoplankton bloom mesocosm experiment. Bacterial cell abundance, determined by flow cytometry, followed a typical bacterioplankton response to a phytoplankton bloom, with Synechococcus and high nucleic acid (HNA) bacterial cell numbers initially falling, probably due to selective protist grazing. Subsequently HNA and low nucleic acid (LNA) bacterial cells increased in abundance but Synechococcus did not. There was no significant difference between bacterioneuston and bacterioplankton cell abundances during the experiment. Conversely, distinct and consistent differences between the bacterioneuston and the bacterioplankton community structure were observed. This was monitored simultaneously by Bacteria 16S rRNA gene terminal restriction fragment length polymorphism (T-RFLP) and denaturing gradient gel electrophoresis (DGGE). The conserved patterns of community structure observed in all of the mesocosms indicate that the bacterioneuston is distinctive and non-random