Application of hpDGFEM to mechanisms at channel microband electrodes

We extend our earlier work (Harriman et al., Oxford University Computing Laboratory Technical Report NA04/19) on hp-DGFEM for disc electrodes to the case of reaction mechanisms to the increasingly popular channel microband electrode configuration. We present results for the simple E reaction mechanism (convection-diffusion equation), for the ECE and EC2E reaction mechanisms (linear and nonlinear systems of reaction-convection- diffusion equations, respectively) and for the DISP1 and DISP2 reaction mechanisms (linear and nonlinear coupled systems of reaction-convection-diffusion equations, respectively). In all cases we demonstrate excellent agreement with previous results using relatively coarse meshes and without the need for streamline-diffusion stabilisation, even at high flow rates

Prioritizing Healthy Child Development Could Prevent Child Prostitution

The CHILDREN AT RISK documentaries “Domestic Minor Sex Trafficking” and “International Human Trafficking” draw attention to the dire consequences of our failure as a society to ensure that all children are raised with healthy experiences in safe and loving environments. It is our collective responsibility to put policies and services into place to prevent child prostitution from happening in the first place, while also providing treatment and care for the victims of prostitution. We must embed the prevention of child prostitution into a broader vision for healthy child development and encourage our national, state, and local policymakers to prioritize the development and implementation of a comprehensive and coordinated strategy for children

Approximation of linear functionals using an hp-adaptive discontinuous Galerkin finite element method

We consider the problem of computing a linear functional of the solution of an elliptic partial differential equation to within a given tolerance. We drive an a posteriori error bound for the linear functional and use this as the basis of an hp-adaptive discontinuous Galerkin finite element algorithm to deliver the functional to within a prescribed error tolerance

Models for pattern formation in somitogenesis: a marriage of cellular and molecular biology

Somitogenesis, the process by which a bilaterally symmetric pattern of cell aggregations is laid down in a cranio-caudal sequence in early vertebrate development, provides an excellent model study for the coupling of interactions at the molecular and cellular level. Here, we review some of the key experimental results and theoretical models related to this process. We extend a recent chemical pre-pattern model based on the cell cycle Journal of Theoretical Biology 207 (2000) 305-316, by including cell movement and show that the resultant model exhibits the correct spatio-temporal dynamics of cell aggregation. We also postulate a model to account for the recently observed spatio-temporal dynamics at the molecular level

Validity of the Cauchy-Born rule applied to discrete cellular-scale models of biological tissues.

The development of new models of biological tissues that consider cells in a discrete manner is becoming increasingly popular as an alternative to continuum methods based on partial differential equations, although formal relationships between the discrete and continuum frameworks remain to be established. For crystal mechanics, the discrete-to-continuum bridge is often made by assuming that local atom displacements can be mapped homogeneously from the mesoscale deformation gradient, an assumption known as the Cauchy-Born rule (CBR). Although the CBR does not hold exactly for noncrystalline materials, it may still be used as a first-order approximation for analytic calculations of effective stresses or strain energies. In this work, our goal is to investigate numerically the applicability of the CBR to two-dimensional cellular-scale models by assessing the mechanical behavior of model biological tissues, including crystalline (honeycomb) and noncrystalline reference states. The numerical procedure involves applying an affine deformation to the boundary cells and computing the quasistatic position of internal cells. The position of internal cells is then compared with the prediction of the CBR and an average deviation is calculated in the strain domain. For center-based cell models, we show that the CBR holds exactly when the deformation gradient is relatively small and the reference stress-free configuration is defined by a honeycomb lattice. We show further that the CBR may be used approximately when the reference state is perturbed from the honeycomb configuration. By contrast, for vertex-based cell models, a similar analysis reveals that the CBR does not provide a good representation of the tissue mechanics, even when the reference configuration is defined by a honeycomb lattice. The paper concludes with a discussion of the implications of these results for concurrent discrete and continuous modeling, adaptation of atom-to-continuum techniques to biological tissues, and model classification

A Comparison of Numerical Methods used for\ud Finite Element Modelling of Soft Tissue\ud Deformation

Soft tissue deformation is often modelled using incompressible nonlinear elasticity, with solutions computed using the finite element method. There are a range of options available when using the finite element method, in particular, the polynomial degree of the basis functions used for interpolating position and pressure, and the type of element making up the mesh. We investigate the effect of these choices on the accuracy of the computed solution, using a selection of model problems motivated by typical deformations seen in soft tissue modelling. We set up model problems with discontinuous material properties (as is the case for the breast), steeply changing gradients in the body force (as found in contracting cardiac tissue), and discontinuous first derivatives in the solution at the boundary, caused by a discontinuous applied force (as in the breast during mammography). We find that the choice of pressure basis functions are vital in the presence of a material interface, higher-order schemes do not perform as well as may be expected when there are sharp gradients, and in general that it is important to take the expected regularity of the solution into account when choosing a numerical scheme

Validity of the Cauchy-Born rule applied to discrete cellular-scale models of biological tissues

The development of new models of biological tissues that consider cells in a discrete manner is becoming increasingly popular as an alternative to PDE-based continuum methods, although formal relationships between the discrete and continuum frameworks remain to be established. For crystal mechanics, the discrete-to-continuum bridge is often made by assuming that local atom displacements can be mapped homogeneously from the mesoscale deformation gradient, an assumption known as the Cauchy-Born rule (CBR). Although the CBR does not hold exactly for non-crystalline materials, it may still be used as a first order approximation for analytic calculations of effective stresses or strain energies. In this work, our goal is to investigate numerically the applicability of the CBR to 2-D cellular-scale models by assessing the mechanical behaviour of model biological tissues, including crystalline (honeycomb) and non-crystalline reference states. The numerical procedure consists in precribing an affine deformation on the boundary cells and computing the position of internal cells. The position of internal cells is then compared with the prediction of the CBR and an average deviation is calculated in the strain domain. For centre-based models, we show that the CBR holds exactly when the deformation gradient is relatively small and the reference stress-free configuration is defined by a honeycomb lattice. We show further that the CBR may be used approximately when the reference state is perturbed from the honeycomb configuration. By contrast, for vertex-based models, a similar analysis reveals that the CBR does not provide a good representation of the tissue mechanics, even when the reference configuration is defined by a honeycomb lattice. The paper concludes with a discussion of the implications of these results for concurrent discrete/continuous modelling, adaptation of atom-to-continuum (AtC) techniques to biological tissues and model classification

Challenges of ultra large scale integration of biomedical computing systems

The NCRI Informatics Initiative is overseeing the implementation of an informatics framework for the UK cancer research community. The framework advocates an integrated multidisciplinary method of working between scientific and medical communities. Key to this process is community adoption of high quality acquisition, storage, sharing and integration of diverse data elements to improve knowledge of the causes, prevention and treatment of cancer. The integration of the complex data and meta-data used by these multiple communities is a significant challenge and there are technical, resource-based and sociological issues to be addressed. In this paper we review progress aimed at establishing the framework and outline key challenges in ultra large scale integration of biomedical computing systems

The importance of adjoint consistency in the approximation of linear functionals using the discontinuous Galerkin finite element method

We describe how a discontinuous Galerkin finite element method with interior penalty can be used to compute the solution to an elliptic partial differential equation and a linear functional of this solution can be evaluated. We show that, in order to have an adjoint consistent method and thus obtain optimal rates of convergence of the functional, a symmetric interior penalty Galerkin method must be used and, when the functional depends on the derivative of the solution of the PDE, an equivalent formulation of the functional must be used