Computation of forces from deformed visco-elastic biological tissues

We present a least-squares based inverse analysis of visco-elastic biological tissues. The proposed method computes the set of contractile forces (dipoles) at the cell boundaries that induce the observed and quantified deformations. We show that the computation of these forces requires the regularisation of the problem functional for some load configurations that we study here. The functional measures the error of the dynamic problem being discretised in time with a second-order implicit time-stepping and in space with standard finite elements. We analyse the uniqueness of the inverse problem and estimate the regularisation parameter by means of an L-curved criterion. We apply the methodology to a simple toy problem and to an in vivo set of morphogenetic deformations of the Drosophila embryo.Peer ReviewedPostprint (author's final draft

Modelling Rod-like Flexible Biological Tissues for Medical Training

This paper outlines a framework for the modelling of slender rod-like biological tissue structures in both global and local scales. Volumetric discretization of a rod-like structure is expensive in computation and therefore is not ideal for applications where real-time performance is essential. In our approach, the Cosserat rod model is introduced to capture the global shape changes, which models the structure as a one-dimensional entity, while the local deformation is handled separately. In this way a good balance in accuracy and efficiency is achieved. These advantages make our method appropriate for the modelling of soft tissues for medical training applications

Time-gated transillumination of biological tissues and tissuelike phantoms

The applicability and limits of time-resolved transillumination to determine the internal details of biological tissues are investigated by phantom experiments. By means of line scans across a sharp edge, the spatial resolution (Ax) and its dependence on the time-gate width (At) can be determined. Additionally, measurements of completely absorbing bead pairs embedded in a turbid medium demonstrate the physical resolution in a more realistic case. The benefit of time resolution is especially high for a turbid medium with a comparatively small reduced scattering coefficient of approximately pL,' = 0.12 mm-1. Investigations with partially absorbing beads and filled plastic tubes demonstrate the high sensitivity of time-resolving techniques with respect to spatial variations in scattering or absorption coefficients that are due to the embedded disturber. In particular, it is shown that time gating is sensitive to variations in scattering coefficients. Key words: Time-resolved transillumination, turbid media, light scattering, streak camera

Classical and all-floating FETI methods for the simulation of arterial tissues

High-resolution and anatomically realistic computer models of biological soft tissues play a significant role in the understanding of the function of cardiovascular components in health and disease. However, the computational effort to handle fine grids to resolve the geometries as well as sophisticated tissue models is very challenging. One possibility to derive a strongly scalable parallel solution algorithm is to consider finite element tearing and interconnecting (FETI) methods. In this study we propose and investigate the application of FETI methods to simulate the elastic behavior of biological soft tissues. As one particular example we choose the artery which is - as most other biological tissues - characterized by anisotropic and nonlinear material properties. We compare two specific approaches of FETI methods, classical and all-floating, and investigate the numerical behavior of different preconditioning techniques. In comparison to classical FETI, the all-floating approach has not only advantages concerning the implementation but in many cases also concerning the convergence of the global iterative solution method. This behavior is illustrated with numerical examples. We present results of linear elastic simulations to show convergence rates, as expected from the theory, and results from the more sophisticated nonlinear case where we apply a well-known anisotropic model to the realistic geometry of an artery. Although the FETI methods have a great applicability on artery simulations we will also discuss some limitations concerning the dependence on material parameters.Comment: 29 page

Low-Dimensional Stochastic Modeling of the Electrical Properties of Biological Tissues

Uncertainty quantification plays an important role in biomedical engineering as measurement data is often unavailable and literature data shows a wide variability. Using state-of-the-art methods one encounters difficulties when the number of random inputs is large. This is the case, e.g., when using composite Cole-Cole equations to model random electrical properties. It is shown how the number of parameters can be significantly reduced by the Karhunen-Loeve expansion. The low-dimensional random model is used to quantify uncertainties in the axon activation during deep brain stimulation. Numerical results for a Medtronic 3387 electrode design are given.Comment: 4 pages, 5 figure

Role of cell deformability in the two-dimensional melting of biological tissues

The size and shape of a large variety of polymeric particles, including biological cells, star polymers, dendrimes, and microgels, depend on the applied stresses as the particles are extremely soft. In high-density suspensions these particles deform as stressed by their neighbors, which implies that the interparticle interaction becomes of many-body type. Investigating a two-dimensional model of cell tissue, where the single particle shear modulus is related to the cell adhesion strength, here we show that the particle deformability affects the melting scenario. On increasing the temperature, stiff particles undergo a first-order solid/liquid transition, while soft ones undergo a continuous solid/hexatic transition followed by a discontinuous hexatic/liquid transition. At zero temperature the melting transition driven by the decrease of the adhesion strength occurs through two continuous transitions as in the Kosterlitz, Thouless, Halperin, Nelson, and Young scenario. Thus, there is a range of adhesion strength values where the hexatic phase is stable at zero temperature, which suggests that the intermediate phase of the epithelial-to-mesenchymal transition could be hexatic type