Hemodynamics of the heart's left atrium based on a Variational Multiscale-LES numerical model

In this paper, we investigate the hemodynamics of a left atrium (LA) by proposing a computational model suitable to provide physically meaningful fluid dynamics indications and detailed blood flow characterization. In particular, we consider the incompressible NavierStokes equations in Arbitrary Lagrangian Eulerian (ALE) formulation to deal with the LA domain under prescribed motion. A Variational Multiscale (VMS) method is adopted to obtain a stable formulation of the Navier-Stokes equations discretized by means of the Finite Element method and to account for turbulence modeling based on Large Eddy Simulation (LES). The aim of this paper is twofold: on one hand to improve the general understanding of blood flow in the human LA in normal conditions; on the other, to analyse the effects of the turbulence VMS-LES method on a situation of blood flow which is neither laminar, nor fully turbulent, but rather transitional as in LA. Our results suggest that if relatively coarse meshes are adopted, the additional stabilization terms introduced by the VMS-LES method allow to better predict transitional effects and cycle-to-cycle blood flow variations than the standard SUPG stabilization method.

[1]  Tomoyuki Yambe,et al.  Numerical analysis of hemodynamic changes in the left atrium due to atrial fibrillation. , 2015, Journal of biomechanics.

[2]  J. Halleux,et al.  An arbitrary lagrangian-eulerian finite element method for transient dynamic fluid-structure interactions , 1982 .

[3]  Alfio Quarteroni,et al.  A Patient-Specific Computational Fluid Dynamics Model of the Left Atrium in Atrial Fibrillation: Development and Initial Evaluation , 2017, FIMH.

[4]  M. O'Neill,et al.  Modeling Left Atrial Flow, Energy, Blood Heating Distribution in Response to Catheter Ablation Therapy , 2018, Front. Physiol..

[5]  I. Marshall,et al.  Comparative Study of Magnetic Resonance Imaging and Image-Based Computational Fluid Dynamics for Quantification of Pulsatile Flow in a Carotid Bifurcation Phantom , 2003, Annals of Biomedical Engineering.

[6]  Tayfun E. Tezduyar,et al.  Mesh update strategies in parallel finite element computations of flow problems with moving boundaries and interfaces , 1994 .

[7]  M. T. Ngo,et al.  Four-Dimensional Flow Magnetic Resonance Imaging for Assessment of Velocity Magnitudes and Flow Patterns in The Human Carotid Artery Bifurcation: Comparison with Computational Fluid Dynamics , 2019, Diagnostics.

[8]  A. Quarteroni.  Numerical Models for Differential Problems , 2009 .

[9]  Jiafeng Zhang,et al.  Computational characterization of flow and blood damage potential of the new maglev CH-VAD pump versus the HVAD and HeartMate II pumps , 2020, The International journal of artificial organs.

[10]  T. Hughes,et al.  Variational and Multiscale Methods in Turbulence , 2005 .

[11]  C. Peskin,et al.  Modelling cardiac fluid dynamics and diastolic function , 2001, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[12]  Luca Dedè,et al.  Semi-implicit BDF time discretization of the Navier–Stokes equations with VMS-LES modeling in a High Performance Computing framework , 2015 .

[13]  Christian Vergara,et al.  Computational Fluid-Dynamic Analysis after Carotid Endarterectomy: Patch Graft versus Direct Suture Closure. , 2017, Annals of vascular surgery.

[14]  A. Hazel,et al.  Spatial comparison between wall shear stress measures and porcine arterial endothelial permeability. , 2004, American journal of physiology. Heart and circulatory physiology.

[15]  Alfio Quarteroni,et al.  Fluid dynamics of an idealized left ventricle: the extended Nitsche's method for the treatment of heart valves as mixed time varying boundary conditions , 2017 .

[16]  Tsuyoshi Murata,et al.  {m , 1934, ACML.

[17]  A. Quarteroni,et al.  A Proof of Concept for Computational Fluid Dynamic Analysis of the Left Atrium in Atrial Fibrillation on a Patient-Specific Basis , 2019, Journal of Biomechanical Engineering.

[18]  Thomas J. R. Hughes,et al.  Large eddy simulation of turbulent channel flows by the variational multiscale method , 2001 .

[19]  T. Hughes.  Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , 1995 .

[20]  Alfio Quarteroni,et al.  Complex blood flow patterns in an idealized left ventricle: A numerical study. , 2017, Chaos.

[21]  Mark D. Huffman,et al.  AHA Statistical Update Heart Disease and Stroke Statistics — 2012 Update A Report From the American Heart Association WRITING GROUP MEMBERS , 2010 .

[22]  Numerical simulation of plane Poiseuille turbulence , 1994 .

[23]  Katharine H Fraser,et al.  The use of computational fluid dynamics in the development of ventricular assist devices. , 2011, Medical engineering & physics.

[24]  Alfio Quarteroni,et al.  Isogeometric approximation of cardiac electrophysiology models on surfaces: An accuracy study with application to the human left atrium , 2017 .

[25]  高等学校計算数学学報編輯委員会編.  高等学校計算数学学報 = Numerical mathematics , 1979 .

[26]  Yuri Bazilevs,et al.  Aerodynamic and FSI Analysis of Wind Turbines with the ALE-VMS and ST-VMS Methods , 2014 .

[27]  Long Chen.  FINITE ELEMENT METHOD , 2013 .

[28]  Eduardo J. Godoy,et al.  Detailed Anatomical and Electrophysiological Models of Human Atria and Torso for the Simulation of Atrial Activation , 2015, PloS one.

[29]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

[30]  T. Hughes,et al.  Large Eddy Simulation and the variational multiscale method , 2000 .

[31]  Gianni Pedrizzetti,et al.  Intraventricular vortex flow changes in the infarcted left ventricle: numerical results in an idealised 3D shape , 2011, Computer methods in biomechanics and biomedical engineering.

[32]  Alessandro Veneziani,et al.  The LifeV library: engineering mathematics beyond the proof of concept , 2017, 1710.06596.

[33]  Alfio Quarteroni,et al.  A numerical study of isotropic and anisotropic constitutive models with relevance to healthy and unhealthy cerebral arterial tissues , 2016 .

[34]  Jung Hee Seo,et al.  Computational modeling and analysis of intracardiac flows in simple models of the left ventricle , 2012 .

[35]  I. Mokhov,et al.  Assessment of the impact of the planetary scale on the decay of blocking and the use of phase diagrams and enstrophy as a diagnostic , 2007 .

[36]  A P Yoganathan,et al.  Left ventricular blood flow patterns in normal subjects: a quantitative analysis by three-dimensional magnetic resonance velocity mapping. , 1995, Journal of the American College of Cardiology.

[37]  Laurent Younes,et al.  Hemodynamics in the Left Atrium and Its Effect on Ventricular Flow Patterns. , 2015, Journal of biomechanical engineering.

[38]  P. Moin,et al.  Eddies, streams, and convergence zones in turbulent flows , 1988 .

[39]  D. Wilcox.  Turbulence modeling for CFD , 1993 .

[40]  M. Fortin,et al.  Mixed Finite Element Methods and Applications , 2013 .

[41]  Thomas J. R. Hughes,et al.  Patient-specific isogeometric fluid–structure interaction analysis of thoracic aortic blood flow due to implantation of the Jarvik 2000 left ventricular assist device , 2009 .

[42]  F. Menter,et al.  Ten Years of Industrial Experience with the SST Turbulence Model , 2003 .

[43]  D. Ku,et al.  Pulsatile Flow and Atherosclerosis in the Human Carotid Bifurcation: Positive Correlation between Plaque Location and Low and Oscillating Shear Stress , 1985, Arteriosclerosis.

[44]  Danna Zhou,et al.  d. , 1840, Microbial pathogenesis.

[45]  T. Korakianitis,et al.  Numerical simulation of cardiovascular dynamics with healthy and diseased heart valves. , 2006, Journal of biomechanics.

[46]  Alfio Quarteroni,et al.  Comparisons between reduced order models and full 3D models for fluid-structure interaction problems in haemodynamics , 2014, J. Comput. Appl. Math..

[47]  Samuele Rubino,et al.  Numerical comparisons of finite element stabilized methods for a 2D vortex dynamics simulation at high Reynolds number , 2019, Computer Methods in Applied Mechanics and Engineering.

[48]  Simone Scacchi,et al.  Parallel multilevel solvers for the cardiac electro-mechanical coupling , 2015 .

[49]  Giorgio Galanti,et al.  Comparative numerical study on left ventricular fluid dynamics after dilated cardiomyopathy. , 2013, Journal of biomechanics.

[50]  T. Hughes,et al.  Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows , 2007 .

[51]  László Szilágyi,et al.  A patient specific electro-mechanical model of the heart , 2011, Comput. Methods Programs Biomed..

[52]  Bradley D. Allen,et al.  4D flow imaging with MRI. , 2014, Cardiovascular diagnosis and therapy.

[53]  Alfio Quarteroni,et al.  Integrated Heart—Coupling multiscale and multiphysics models for the simulation of the cardiac function , 2017 .

[54]  David A. Steinman,et al.  An image-based modeling framework for patient-specific computational hemodynamics , 2008, Medical & Biological Engineering & Computing.

[55]  Jiun-Jr Wang,et al.  Expanding application of the Wiggers diagram to teach cardiovascular physiology. , 2014, Advances in physiology education.

[56]  Alfio Quarteroni,et al.  Mathematical Modelling of the Human Cardiovascular System , 2019 .

[57]  Franck Nicoud,et al.  Image-based large-eddy simulation in a realistic left heart , 2014 .

[58]  Einar Heiberg,et al.  Left ventricular fluid kinetic energy time curves in heart failure from cardiovascular magnetic resonance 4D flow data , 2015, Journal of Cardiovascular Magnetic Resonance.

[59]  Toshiaki Hisada,et al.  Multiphysics simulation of left ventricular filling dynamics using fluid-structure interaction finite element method. , 2004, Biophysical journal.

[60]  Gianni Pedrizzetti,et al.  Three-dimensional filling flow into a model left ventricle , 2005, Journal of Fluid Mechanics.

[61]  Deepak K Gupta,et al.  Left atrial structure and function in atrial fibrillation: ENGAGE AF-TIMI 48. , 2014, European heart journal.

[62]  Thomas J. R. Hughes,et al.  Multiscale and Stabilized Methods , 2007 .

[63]  M. Behr,et al.  The variational multiscale formulation for the fully‐implicit log‐morphology equation as a tensor‐based blood damage model , 2019, International journal for numerical methods in biomedical engineering.

[64]  Alfio Quarteroni,et al.  The Impact of Left Atrium Appendage Morphology on Stroke Risk Assessment in Atrial Fibrillation: A Computational Fluid Dynamics Study , 2019, Front. Physiol..