Fluid–structure interaction analysis of eccentricity and leaflet rigidity on thrombosis biomarkers in bioprosthetic aortic valve replacements

This work intends to study the effect of aortic annulus eccentricity and leaflet rigidity on the performance, thrombogenic risk and calcification risk in bioprosthetic aortic valve replacements (BAVRs). To address these questions, a two-way immersed fluid–structure interaction (FSI) computational model was implemented in a high-performance computing (HPC) multi-physics simulation software, and validated against a well-known FSI benchmark. The aortic valve bioprosthesis model is qualitatively contrasted against experimental data, showing good agreement in closed and open states. Regarding the performance of BAVRs, the model predicts that increasing eccentricities yield lower geometric orifice areas (GOAs) and higher normalized transvalvular pressure gradients (TPGs) for healthy cardiac outputs during systole, agreeing with in vitro experiments. Regions with peak values of residence time are observed to grow with eccentricity in the sinus of Valsalva, indicating an elevated risk of thrombus formation for eccentric configurations. In addition, the computational model is used to analyze the effect of varying leaflet rigidity on both performance, thrombogenic and calcification risks with applications to tissue-engineered prostheses. For more rigid leaflets it predicts an increase in systolic and diastolic TPGs, and decrease in systolic GOA, which translates to decreased valve performance. The peak shear rate and residence time regions increase with leaflet rigidity, but their volume-averaged values were not significantly affected. Peak solid stresses are also analyzed, and observed to increase with rigidity, elevating risk of valve calcification and structural failure. To the authors' knowledge this is the first computational FSI model to study the effect of eccentricity or leaflet rigidity on thrombogenic biomarkers, providing a novel tool to aid device manufacturers and clinical practitioners.This project has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 713673. The research leading to these results has also received funding from “la Caixa” Foundation, with fellowship ID: LCF/BQ/DI18/11660044, and has been co-funded by the project CompBioMed2 (H2020-EU.1.4.1.3. Grant No. 823712)Peer ReviewedPostprint (published version

Effects of detailed ventricular anatomy on the blood flow

The presented study is a preliminary test and analysis of the role of trabeculae and papillary muscles in the hemodynamics of the left ventricle (LV)

Mind the gap: quantification of incomplete ablation patterns after pulmonary vein isolation using minimum path search

Pulmonary vein isolation (PVI) is a common procedure for the treatment of atrial fibrillation (AF). A successful isolation produces a continuous lesion (scar) completely encircling the veins that stops activation waves from propagating to the atrial body. Unfortunately, the encircling lesion is often incomplete, becoming a combination of scar and gaps of healthy tissue. These gaps are potential causes of AF recurrence, which requires a redo of the isolation procedure. Late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) is a non-invasive method that may also be used to detect gaps, but it is currently a time-consuming process, prone to high inter-observer variability. In this paper, we present a method to semi-automatically identify and quantify ablation gaps. Gap quantification is performed through minimum path search in a graph where every node is a scar patch and the edges are the geodesic distances between patches. We propose the Relative Gap Measure (RGM) to estimate the percentage of gap around a vein, which is defined as the ratio of the overall gap length and the total length of the path that encircles the vein. Additionally, an advanced version of the RGM has been developed to integrate gap quantification estimates from different scar segmentation techniques into a single figure-of-merit. Population-based statistical and regional analysis of gap distribution was performed using a standardised parcellation of the left atrium. We have evaluated our method on synthetic and clinical data from 50 AF patients who underwent PVI with radiofrequency ablation. The population-based analysis concluded that the left superior PV is more prone to lesion gaps while the left inferior PV tends to have less gaps (p<0.05 in both cases), in the processed data. This type of information can be very useful for the optimization and objective assessment of PVI interventions

Left Ventricular Trabeculations Decrease the Wall Shear Stress and Increase the Intra-Ventricular Pressure Drop in CFD Simulations

The aim of the present study is to characterize the hemodynamics of left ventricular (LV) geometries to examine the impact of trabeculae and papillary muscles (PMs) on blood flow using high performance computing (HPC). Five pairs of detailed and smoothed LV endocardium models were reconstructed from high-resolution magnetic resonance images (MRI) of ex-vivo human hearts. The detailed model of one LV pair is characterized only by the PMs and few big trabeculae, to represent state of art level of endocardial detail. The other four detailed models obtained include instead endocardial structures measuring ≥1 mm2 in cross-sectional area. The geometrical characterizations were done using computational fluid dynamics (CFD) simulations with rigid walls and both constant and transient flow inputs on the detailed and smoothed models for comparison. These simulations do not represent a clinical or physiological scenario, but a characterization of the interaction of endocardial structures with blood flow. Steady flow simulations were employed to quantify the pressure drop between the inlet and the outlet of the LVs and the wall shear stress (WSS). Coherent structures were analyzed using the Q-criterion for both constant and transient flow inputs. Our results show that trabeculae and PMs increase the intra-ventricular pressure drop, reduce the WSS and disrupt the dominant single vortex, usually present in the smoothed-endocardium models, generating secondary small vortices. Given that obtaining high resolution anatomical detail is challenging in-vivo, we propose that the effect of trabeculations can be incorporated into smoothed ventricular geometries by adding a porous layer along the LV endocardial wall. Results show that a porous layer of a thickness of 1.2·10−2 m with a porosity of 20 kg/m2 on the smoothed-endocardium ventricle models approximates the pressure drops, vorticities and WSS observed in the detailed models.This paper has been partially funded by CompBioMed project, under H2020-EU.1.4.1.3 European Union’s Horizon 2020 research and innovation programme, grant agreement n◦ 675451. FS is supported by a grant from Severo Ochoa (n◦ SEV-2015-0493-16-4), Spain. CB is supported by a grant from the Fundació LaMarató de TV3 (n◦ 20154031), Spain. TI and PI are supported by the Institute of Engineering in Medicine, USA, and the Lillehei Heart Institute, USA.Peer ReviewedPostprint (published version

Ventricular anatomical complexity and sex differences impact predictions from electrophysiological computational models

The aim of this work was to analyze the influence of sex hormones and anatomical details (trabeculations and false tendons) on the electrophysiology of healthy human hearts. Additionally, sex- and anatomy-dependent effects of ventricular tachycardia (VT) inducibility are presented. To this end, four anatomically normal, human, biventricular geometries (two male, two female), with identifiable trabeculations, were obtained from high-resolution, ex-vivo MRI and represented by detailed and smoothed geometrical models (with and without the trabeculations). Additionally one model was augmented by a scar. The electrophysiology finite element model (FEM) simulations were carried out, using O’Hara-Rudy human myocyte model with sex phenotypes of Yang and Clancy. A systematic comparison between detailed vs smooth anatomies, male vs female normal hearts was carried out. The heart with a myocardial infarction was subjected to a programmed stimulus protocol to identify the effects of sex and anatomical detail on ventricular tachycardia inducibility. All female hearts presented QT-interval prolongation however the prolongation interval in comparison to the male phenotypes was anatomy-dependent and was not correlated to the size of the heart. Detailed geometries showed QRS fractionation and increased T-wave magnitude in comparison to the corresponding smoothed geometries. A variety of sustained VTs were obtained in the detailed and smoothed male geometries at different pacing locations, which provide evidence of the geometry-dependent differences regarding the prediction of the locations of reentry channels. In the female phenotype, sustained VTs were induced in both detailed and smooth geometries with RV apex pacing, however no consistent reentry channels were identified. Anatomical and physiological cardiac features play an important role defining risk in cardiac disease. These are often excluded from cardiac electrophysiology simulations. The assumption that the cardiac endocardium is smooth may produce inaccurate predictions towards the location of reentry channels in in-silico tachycardia inducibility studiesJA-S, FS, GH and MV are supported by the European Union’s Horizon 2020 research and innovation programme under grant agreements No 675451 (Compbiomed project phase 1) and No 823712 (CompBioMed project, phase 2) and project No 777204 (SilicoFCM project). Part of the simulation computing hours were provided by the CompBioMed project phase 1. JA-S was awarded computation time from Red Espanola de Supercomputacion (RES). (Activity IDs: FI-2018-2-0049 and BCV-2019-2-0014) JA-S is funded by a Ramon y Cajal fellowship (RYC-2017-22532), Ministerio de Ciencia e Innovacion, Spain; and by Plan Estatal de Investigacion Cientifica y Tecnica y de Innovacion 2017-2020 from the Ministerio de Ciencia e Innovacion y Universidades (PID2019-104356RBC41/AEI/10.13039/501100011033): meHeart ME PID2019-104356RB-C44. CB is funded by the Torres Quevedo Program (PTQ2018-010290), Ministerio de Ciencia e Innovacion, Spain. MV, GH and CB are funded by the Spanish Neotec project EXP - 00123159/SNEO-20191113 Generador de corazones virtuales. LKGM was funded by Fundacion Carolina-BBVA. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. There was no additional external funding received for this study.Peer Reviewed"Article signat per 11 autors/es: Pablo Gonzalez-Martin,Federica Sacco,Constantine Butakoff,Ruben Doste,Carlos Bederian,Lilian K. Gutierrez Espinosa de los Monteros,Guillaume Houzeaux,Paul A. Iaizzo,Tinen L. Iles,Mariano Vazquez,Jazmin Aguado-Sierra"Postprint (published version

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

BACKGROUND: Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. METHODS: The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow , maximum aortic root flow , average LVAD flow , and maximum LVAD flow ); (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. RESULTS: Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. Conclusions This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.This project was funded in part by the FDA Critical Path Initiative and by an appointment to the Research Participation Program at the Division of Biomedical Physics, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, U.S. Food and Drug Administration, administered by the Oak Ridge Institute for Science, and Education through an interagency agreement between the U.S. Department of Energy and FDA to RAG. MV and AS acknowledge the funding from the project CompBioMed2 (H2020-EU.1.4.1.3. Grant number: 823712), SilicoFCM (H2020-EU.3.1.5. Grant number: 777204), and NEOTEC 2019 - "Generador de Corazones Virtuales" (“Ministerio de Economía y competititvidad”, EXP - 00123159 / SNEO-20191113). AS salary is partially funded by the “Ministerio de Economía y competititvidad” under the Torres Quevedo Program (grant number: PTQ2019-010528). CB salary is partially funded by the Torres Quevedo Program (grant number: PTQ2018-010290). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Peer ReviewedPostprint (published version

Evaluating the roles of detailed endocardial structures on right ventricular haemodynamics by means of CFD simulations

Computational modelling plays an important role in right ventricular (RV) haemodynamic analysis. However, current approaches use smoothed ventricular anatomies. The aim of this study is to characterise RV haemodynamics including detailed endocardial structures like trabeculae, moderator band, and papillary muscles. Four paired detailed and smoothed RV endocardium models (2 male and 2 female) were reconstructed from ex vivo human hearts high‐resolution magnetic resonance images. Detailed models include structures with ≥1 mm2 cross‐sectional area. Haemodynamic characterisation was done by computational fluid dynamics simulations with steady and transient inflows, using high‐performance computing. The differences between the flows in smoothed and detailed models were assessed using Q‐criterion for vorticity quantification, the pressure drop between inlet and outlet, and the wall shear stress. Results demonstrated that detailed endocardial structures increase the degree of intra‐ventricular pressure drop, decrease the wall shear stress, and disrupt the dominant vortex creating secondary small vortices. Increasingly turbulent blood flow was observed in the detailed RVs. Female RVs were less trabeculated and presented lower pressure drops than the males. In conclusion, neglecting endocardial structures in RV haemodynamic models may lead to inaccurate conclusions about the pressures, stresses, and blood flow behaviour in the cavity.The DICOMdatasetswere provided by the Visible Heart R Laboratory, obtained byMRI scanning of perfusion fixed hearts that were graciously donated by the organ donors and their families through LifeSource. Part of the simulation hours were provided by the CompBioMed project in the Archer supercomputer, EPCC, UK.Peer ReviewedPostprint (author's final draft

Quantitative analysis of lead position vs. correction of electrical dyssynchrony in an experimental model of LBBB/CRT

International audienceCardiac resynchronization therapy (CRT) is a recommended treatment in patients with electrical dyssynchrony such as left bundle branch block (LBBB). The determination of the optimal leads position, and the quantification of the changes in electrical activation are two current major challenges. In this paper, we investigate these aspects through electroanatomical data from a controlled experimental protocol, which studied pigs with no structural disease under LBBB and CRT conditions. We propose to use a quasi-conformal mapping technique to standardize electroanatomical maps of endo- and epi-cardial walls of both ventricles to a common reference geometry, in which simple quantitative indices can be computed. Then, we investigate the relation between leads and simple surrogates of the recovery of the electrical activation based either on total activation times or on the spatial distribution of the patterns. Our methodology allows a better understanding of the complex electrical activation patterns in LBBB and CRT, and confirms hypotheses about the optimal leads position from previous studies