Search

Purpose: Personalized computational simulations of the heart could open up new improved approaches to diagnosis and surgery assistance systems. While it is fully recognized that myocardial fiber orientation is central for the construction of realistic computational models of cardiac electromechanics, the role of its overall architecture and connectivity remains unclear. Morphological studies show that the distribution of cardiac muscular fibers at the basal ring connects epicardium and endocardium. However, computational models simplify their distribution and disregard the basal loop. This work explores the influence in computational simulations of fiber distribution at different short-axis cuts. Methods: We have used a highly parallelized computational solver to test different fiber models of ventricular muscular connectivity. We have considered two rule-based mathematical models and an own-designed method preserving basal connectivity as observed in experimental data. Simulated cardiac functional scores (rotation, torsion and longitudinal shortening) were compared to experimental healthy ranges using generalized models (rotation) and Mahalanobis distances (shortening, torsion). Results: The probability of rotation was significantly lower for ruled-based models [95% CI (0.13, 0.20)] in comparison with experimental data [95% CI (0.23, 0.31)]. The Mahalanobis distance for experimental data was in the edge of the region enclosing 99% of the healthy population. Conclusions: Cardiac electromechanical simulations of the heart with fibers extracted from experimental data produce functional scores closer to healthy ranges than rule-based models disregarding architecture connectivity.

In this work, we present a fully coupled fluid-electro-mechanical model of a 50th percentile human heart. The model is implemented on Alya, the BSC multi-physics parallel code, capable of running efficiently in supercomputers. Blood in the cardiac cavities is modeled by the incompressible Navier-Stokes equations and an arbitrary Lagrangian-Eulerian (ALE) scheme. Electrophysiology is modeled with a monodomain scheme and the O'Hara-Rudy cell model. Solid mechanics is modeled with a total Lagrangian formulation for discrete strains using the Holzapfel-Ogden cardiac tissue material model. The three problems are simultaneously and bidirectionally coupled through an electromechanical feedback and a fluid-structure interaction scheme. In this paper, we present the scheme in detail and propose it as a computational cardiac workbench.

Alya is a multi-physics simulation code developed at Barcelona Supercomputing Center (BSC). From its inception Alya code is designed using advanced High Performance Computing programming techniques to solve coupled problems on supercomputers efficiently. The target domain is engineering, with all its particular features: complex geometries and unstructured meshes, coupled multi-physics with exotic coupling schemes and physical models, ill-posed problems, flexibility needs for rapidly including new models, etc. Since its beginnings in 2004, Alya has scaled well in an increasing number of processors when solving single-physics problems such as fluid mechanics, solid mechanics, acoustics, etc. Over time, we have made a concerted effort to maintain and even improve scalability for multi-physics problems. This poses challenges on multiple fronts, including: numerical models, parallel implementation, physical coupling models, algorithms and solution schemes, meshing process, etc. In this paper, we introduce Alya's main features and focus particularly on its solvers. We present Alya's performance up to 100.000 processors in Blue Waters, the NCSA supercomputer with selected multi-physics tests that are representative of the engineering world. The tests are incompressible flow in a human respiratory system, low Mach combustion problem in a kiln furnace, and coupled electro-mechanical contraction of the heart. We show scalability plots for all cases and discuss all aspects of such simulations, including solver convergence., The authors would like to thank the following fellow researchers and institutions: • The Private Sector Program at NCSA and the BlueWaters sustained-petascale computing project-supported by the National Science Foundation (award number OCI 07-25070) and the state of Illinois. • Denis Doorly and Alister Bates (Imperial College London, UK), collaborators of the airways study. Part of this work was financed by European PRACE Type B/C projects. • The heart geometry was provided by Dr. A. Berruezo (Hospital Clinic de Barcelona) in collaboration with R. Sebastian (UVEG) and O. Camara (UPF), partially financed through project TIN2011-28067 from MINECO, Spain. • Part of the cardiac model development was financed by the grant SEV-2011-00067 of Severo Ochoa Program, awarded by the Spanish Government. • Part of the kiln model development was financed by the European Commission in the framework of the FP7 Collaborative project “Advanced Technologies for the Production of Cement and Clean Aggregates from Construction and Demolition Waste (C2CA)”, Grant Agreement No 265189., Peer Reviewed, Postprint (author's final draft)

Alya is a multi-physics simulation code developed at Barcelona Supercomputing Center (BSC). From its inception Alya code is designed using advanced High Performance Computing programming techniques to solve coupled problems on supercomputers efficiently. The target domain is engineering, with all its particular features: complex geometries and unstructured meshes, coupled multi-physics with exotic coupling schemes and physical models, ill-posed problems, flexibility needs for rapidly including new models, etc. Since its beginnings in 2004, Alya has scaled well in an increasing number of processors when solving single-physics problems such as fluid mechanics, solid mechanics, acoustics, etc. Over time, we have made a concerted effort to maintain and even improve scalability for multi-physics problems. This poses challenges on multiple fronts, including: numerical models, parallel implementation, physical coupling models, algorithms and solution schemes, meshing process, etc. In this paper, we introduce Alya's main features and focus particularly on its solvers. We present Alya's performance up to 100.000 processors in Blue Waters, the NCSA supercomputer with selected multi-physics tests that are representative of the engineering world. The tests are incompressible flow in a human respiratory system, low Mach combustion problem in a kiln furnace, and coupled electro-mechanical contraction of the heart. We show scalability plots for all cases and discuss all aspects of such simulations, including solver convergence. © 2016 Elsevier B.V. All rights reserved.

Peer Reviewed