Popis: |
In this paper, we propose a parallel fluid-structure interaction solution method for the simulation of the diastolic process of the human left ventricle. The three dimensional geometry of a patient-specific left ventricle is first segmented from a medical image, and then based on which an unstructured mesh is generated. The coupled fluid-structure interaction equations are spatially discretized by a finite element method on the unstructured mesh and a fully implicit scheme is introduced for the temporal discretization. At each timestep of the fully implicit method, a scalable parallel method is applied for solving the large scale nonlinear system. Numerical experiments show that the proposed method scales up to 1920 processor cores on the Tianhe-2A supercomputer for solving a problem with over 7.9 million unstructured tetrahedrons, which shows a potential to do high-fidelity simulation of human heart for clinical applications. |