Fast multiscale modeling of cardiac electrophysiology including Purkinje system

In this paper, we present a modeling methodology to couple the cardiac conduction system to cardiac myocytes through a model of Purkinje-ventricular junctions to yield fast and realistic electrical activation of the ventricles. A patient-specific biventricular geometry is obtained from processing computed tomography scan data. A one-manifold implementation of the fast marching method based on Eikonal-type equations is used for modeling heart electrophysiology, which facilitates the multiscale 1-D-3-D coupling at very low computational costs. The method is illustrated in in-silico experiments where we analyze and compare alternative pacing strategies on the same patient-specific anatomy. We also show very good agreement between the results from the proposed approach and more detailed and comprehensive biophysical models for modeling cardiac electrophysiology. The effect of atrioventricular delay on the distribution of activation time in myocardium is studied with two experiments. Given the reasonable computational times and realistic activation sequences provided by our method, it can have an important clinical impact on the selection of optimal implantation sites of pacing leads or placement of ablation catheter's tip in the context of cardiac rhythm management therapies.     

Reconstructing the 3D shape and bone mineral density distribution of the proximal femur from dual-energy X-ray absorptiometry

The accurate diagnosis of osteoporosis has gained increasing importance due to the aging of our society. Areal bone mineral density (BMD) measured by dual-energy X-ray absorptiometry (DXA) is an established criterion in the diagnosis of osteoporosis. This measure, however, is limited by its two-dimensionality. This work presents a method to reconstruct both the 3D bone shape and 3D BMD distribution of the proximal femur from a single DXA image used in clinical routine. A statistical model of the combined shape and BMD distribution is presented, together with a method for its construction from a set of quantitative computed tomography (QCT) scans. A reconstruction is acquired in an intensity based 3D-2D registration process whereby an instance of the model is found that maximizes the similarity between its projection and the DXA image. Reconstruction experiments were performed on the DXA images of 30 subjects, with a model constructed from a database of QCT scans of 85 subjects. The accuracy was evaluated by comparing the reconstructions with the same subject QCT scans. The method presented here can potentially improve the diagnosis of osteoporosis and fracture risk assessment from the low radiation dose and low cost DXA devices currently used in clinical routine.     

Solid damping in micro electro mechanical systems

This paper focuses on the problem of the numerical evaluation of dissipation induced by thermoelastic coupling in microelectromechanical systems. An ad hoc conceived, FE based, numerical procedure for the evaluation of the thermoelastic dissipation is proposed and the numerical results are compared with analytical solutions. In order to introduce in the numerical response a dependence on the size of the resonating devices, which is experimentally observed at very small dimensions, a new enhanced non-local coupled thermoelastic model is proposed and the first results are discussed. An erratum to this article can be found at     

Propagation of measurement noise through backprojection reconstruction in electrical impedance tomography

A framework to analyze the propagation of measurement noise through backprojection reconstruction algorithms in electrical impedance tomography (EIT) is presented. Two measurement noise sources were considered: noise in the current drivers and in the voltage detectors. The influence of the acquisition system architecture (serial/semi-parallel) is also discussed. Three variants of backprojection reconstruction are studied: basic (unweighted), weighted and exponential backprojection. The results of error propagation theory have been compared with those obtained from simulated and experimental data. This comparison shows that the approach provides a good estimate of the reconstruction error variance. It is argued that the reconstruction error in EIT images obtained via backprojection can be approximately modeled as a spatially nonstationary Gaussian distribution. This methodology allows us to develop a spatial characterization of the reconstruction error in EIT images.     

Finite element modelling of a rotating piezoelectric ultrasonic motor

The evaluation of the performance of ultrasonic motors as a function of input parameters such as the driving frequency, voltage input and pre-load on the rotor is of key importance to their development and is here addressed by means of a finite element three-dimensional model. First the stator is simulated as a fully deformable elastic body and the travelling wave dynamics is accurately reproduced; secondly the interaction through contact between the stator and the rotor is accounted for by assuming that the rotor behaves as a rigid surface. Numerical results for the whole motor are finally compared to available experimental data.