Fetal XCMR: a numerical phantom for fetal cardiovascular magnetic resonance imaging.

Validating new techniques for fetal cardiovascular magnetic resonance (CMR) is challenging due to random fetal movement that precludes repeat measurements. Consequently, fetal CMR development has been largely performed using physical phantoms or postnatal volunteers. In this work, we present an open-source simulation designed to aid in the development and validation of new approaches for fetal CMR. Our approach, fetal extended Cardiac-Torso cardiovascular magnetic resonance imaging (Fetal XCMR), builds on established methods for simulating CMR acquisitions but is tailored toward the dynamic physiology of the fetal heart and body. We present comparisons between the Fetal XCMR phantom and data acquired in utero, resulting in image quality, anatomy, tissue signals and contrast. Existing extended Cardiac-Torso models are modified to create maternal and fetal anatomy, combined according to simulated motion, mapped to CMR contrast, and converted to CMR data. To provide a comparison between the proposed simulation and experimental fetal CMR images acquired in utero, images from a typical scan of a pregnant woman are included and simulated acquisitions were generated using matching CMR parameters, motion and noise levels. Three reconstruction (static, real-time, and CINE), and two motion estimation methods (translational motion, fetal heart rate) from data acquired in transverse, sagittal, coronal, and short-axis planes of the fetal heart were performed to compare to in utero acquisitions and demonstrate feasibility of the proposed simulation framework. Overall, CMR contrast, morphologies, and relative proportions of the maternal and fetal anatomy are well represented by the Fetal XCMR images when comparing the simulation to static images acquired in utero. Additionally, visualization of maternal respiratory and fetal cardiac motion is comparable between Fetal XCMR and in utero real-time images. Finally, high quality CINE image reconstructions provide excellent delineation of fetal cardiac anatomy and temporal dynamics for both data types. The fetal CMR phantom provides a new method for evaluating fetal CMR acquisition and reconstruction methods by simulating the underlying anatomy and physiology. As the field of fetal CMR continues to grow, new methods will become available and require careful validation. The fetal CMR phantom is therefore a powerful and convenient tool in the continued development of fetal cardiac imaging

Assessing the performance of ultrafast vector flow imaging in the neonatal heart via multiphysics modeling and In vitro experiments

Ultrafast vector flow imaging would benefit newborn patients with congenital heart disorders, but still requires thorough validation before translation to clinical practice. This paper investigates 2-D speckle tracking (ST) of intraventricular blood flow in neonates when transmitting diverging waves at ultrafast frame rate. Computational and in vitro studies enabled us to quantify the performance and identify artifacts related to the flow and the imaging sequence. First, synthetic ultrasound images of a neonate's left ventricular flow pattern were obtained with the ultrasound simulator Field II by propagating point scatterers according to 3-D intraventricular flow fields obtained with computational fluid dynamics (CFD). Noncompounded diverging waves (opening angle of 60 degrees) were transmitted at a pulse repetition frequency of 9 kHz. ST of the B-mode data provided 2-D flow estimates at 180 Hz, which were compared with the CFD flow field. We demonstrated that the diastolic inflow jet showed a strong bias in the lateral velocity estimates at the edges of the jet, as confirmed by additional in vitro tests on a jet flow phantom. Furthermore, ST performance was highly dependent on the cardiac phase with low flows (< 5 cm/s), high spatial flow gradients, and out-of-plane flow as deteriorating factors. Despite the observed artifacts, a good overall performance of 2-D ST was obtained with a median magnitude underestimation and angular deviation of, respectively, 28% and 13.5 degrees during systole and 16% and 10.5 degrees during diastole

The LifeV library: engineering mathematics beyond the proof of concept

LifeV is a library for the finite element (FE) solution of partial differential equations in one, two, and three dimensions. It is written in C++ and designed to run on diverse parallel architectures, including cloud and high performance computing facilities. In spite of its academic research nature, meaning a library for the development and testing of new methods, one distinguishing feature of LifeV is its use on real world problems and it is intended to provide a tool for many engineering applications. It has been actually used in computational hemodynamics, including cardiac mechanics and fluid-structure interaction problems, in porous media, ice sheets dynamics for both forward and inverse problems. In this paper we give a short overview of the features of LifeV and its coding paradigms on simple problems. The main focus is on the parallel environment which is mainly driven by domain decomposition methods and based on external libraries such as MPI, the Trilinos project, HDF5 and ParMetis. Dedicated to the memory of Fausto Saleri.Comment: Review of the LifeV Finite Element librar

Evolution of spiral and scroll waves of excitation in a mathematical model of ischaemic border zone

Abnormal electrical activity from the boundaries of ischemic cardiac tissue is recognized as one of the major causes in generation of ischemia-reperfusion arrhythmias. Here we present theoretical analysis of the waves of electrical activity that can rise on the boundary of cardiac cell network upon its recovery from ischaemia-like conditions. The main factors included in our analysis are macroscopic gradients of the cell-to-cell coupling and cell excitability and microscopic heterogeneity of individual cells. The interplay between these factors allows one to explain how spirals form, drift together with the moving boundary, get transiently pinned to local inhomogeneities, and finally penetrate into the bulk of the well-coupled tissue where they reach macroscopic scale. The asymptotic theory of the drift of spiral and scroll waves based on response functions provides explanation of the drifts involved in this mechanism, with the exception of effects due to the discreteness of cardiac tissue. In particular, this asymptotic theory allows an extrapolation of 2D events into 3D, which has shown that cells within the border zone can give rise to 3D analogues of spirals, the scroll waves. When and if such scroll waves escape into a better coupled tissue, they are likely to collapse due to the positive filament tension. However, our simulations have shown that such collapse of newly generated scrolls is not inevitable and that under certain conditions filament tension becomes negative, leading to scroll filaments to expand and multiply leading to a fibrillation-like state within small areas of cardiac tissue.Comment: 26 pages, 13 figures, appendix and 2 movies, as accepted to PLoS ONE 2011/08/0

Negative tension of scroll wave filaments and turbulence in three-dimensional excitable media and application in cardiac dynamics

Scroll waves are vortices that occur in three-dimensional excitable media. Scroll waves have been observed in a variety of systems including cardiac tissue, where they are associated with cardiac arrhythmias. The disorganization of scroll waves into chaotic behavior is thought to be the mechanism of ventricular fibrillation, whose lethality is widely known. One possible mechanism for this process of scroll wave instability is negative filament tension. It was discovered in 1987 in a simple two variables model of an excitable medium. Since that time, negative filament tension of scroll waves and the resulting complex, often turbulent dynamics was studied in many generic models of excitable media as well as in physiologically realistic models of cardiac tissue. In this article, we review the work in this area from the first simulations in FitzHugh-Nagumo type models to recent studies involving detailed ionic models of cardiac tissue. We discuss the relation of negative filament tension and tissue excitability and the effects of discreteness in the tissue on the filament tension. Finally, we consider the application of the negative tension mechanism to computational cardiology, where it may be regarded as a fundamental mechanism that explains differences in the onset of arrhythmias in thin and thick tissue

A comparative study fourth order runge kutta-tvd Scheme and fluent software case of inlet flow problems

Inlet as part of aircraft engine plays important role in controlling the rate of airflow entering to the engine. The shape of inlet has to be designed in such away to make the rate of airflow does not change too much with angle of attack and also not much pressure losses at the time, the airflow entering to the compressor section. It is therefore understanding on the flow pattern inside the inlet is important. The present work presents on the use of the Fourth Order Runge Kutta – Harten Yee TVD scheme for the flow analysis inside inlet. The flow is assumed as an inviscid quasi two dimensional compressible flow. As an initial stage of computer code development, here uses three generic inlet models. The first inlet model to allow the problem in hand solved as the case of inlet with expansion wave case. The second inlet model will relate to the case of expansion compression wave. The last inlet model concerns with the inlet which produce series of weak shock wave and end up with a normal shock wave. The comparison result for the same test case with Fluent Software [1, 2] indicates that the developed computer code based on the Fourth Order Runge Kutta – Harten – Yee TVD scheme are very close to each other. However for complex inlet geometry, the problem is in the way how to provide an appropriate mesh model

The Evolution of Reaction-diffusion Controllers for Minimally Cognitive Agents

No description supplie

Modeling the Heart as a Communication System

Electrical communication between cardiomyocytes can be perturbed during arrhythmia, but these perturbations are not captured by conventional electrocardiographic metrics. We developed a theoretical framework to quantify electrical communication using information theory metrics in 2-dimensional cell lattice models of cardiac excitation propagation. The time series generated by each cell was coarse-grained to 1 when excited or 0 when resting. The Shannon entropy for each cell was calculated from the time series during four clinically important heart rhythms: normal heartbeat, anatomical reentry, spiral reentry, and multiple reentry. We also used mutual information to perform spatial profiling of communication during these cardiac arrhythmias. We found that information sharing between cells was spatially heterogeneous. In addition, cardiac arrhythmia significantly impacted information sharing within the heart. Entropy localized the path of the drifting core of spiral reentry, which could be an optimal target of therapeutic ablation. We conclude that information theory metrics can quantitatively assess electrical communication among cardiomyocytes. The traditional concept of the heart as a functional syncytium sharing electrical information cannot predict altered entropy and information sharing during complex arrhythmia. Information theory metrics may find clinical application in the identification of rhythm-specific treatments which are currently unmet by traditional electrocardiographic techniques.Comment: 26 pages (including Appendix), 6 figures, 8 videos (not uploaded due to size limitation