Correlation between LAA Morphological Features and Computational Fluid Dynamics Analysis for Non-Valvular Atrial Fibrillation Patients

Impact of left atrial wall motion assumptions in fluid simulations on proposed predictors of thrombus formation

Atrial fibrillation (AF) poses a significant risk of stroke due to thrombus formation, which primarily occurs in the left atrial appendage (LAA). Medical image-based computational fluid dynamics (CFD) simulations can provide valuable insight into patient-specific hemodynamics and could potentially enhance personalized assessment of thrombus risk. However, the importance of accurately representing the left atrial (LA) wall dynamics has not been fully resolved. In this study, we compared four modeling scenarios; rigid walls, a generic wall motion based on a reference motion, a semi-generic wall motion based on patient-specific motion, and patient-specific wall motion based on medical images. We considered a LA geometry acquired from 4D computed tomography during AF, systematically performed convergence tests to assess the numerical accuracy of our solution strategy, and quantified the differences between the four approaches. The results revealed that wall motion had no discernible impact on LA cavity hemodynamics, nor on the markers that indicate thrombus formation. However, the flow patterns within the LAA deviated significantly in the rigid model, indicating that the assumption of rigid walls may lead to errors in the estimated risk factors. In contrast, the generic, semi-generic, and patient-specific cases were qualitatively similar. The results highlight the crucial role of wall motion on hemodynamics and predictors of thrombus formation, and also demonstrate the potential of using a generic motion model as a surrogate for the more complex patient-specific motion. While the present study considered a single case, the employed CFD framework is entirely open-source and designed for adaptability, allowing for integration of additional models and generic motions.

Left atrial appendage occlusion: On the need of a numerical model to simulate the implant procedure

Left atrial appendage occlusion (LAAO) is a percutaneous procedure to prevent thromboembolism in patients affected by atrial fibrillation. Despite its demonstrated efficacy, the LAA morphological complexity hinders the procedure, resulting in postprocedural drawbacks (device-related thrombus and peri-device leakage). Local anatomical features may cause difficulties in the device's positioning and affect the effectiveness of the device's implant. The current work proposes a detailed FE model of the LAAO useful to investigate implant scenarios and derive clinical indications. A high-fidelity model of the Watchman FLX device and simplified parametric conduits mimicking the zone of the LAA where the device is deployed were developed. Device-conduit interactions were evaluated by looking at clinical indicators such as device-wall gap, possible cause of leakage, and device protrusion. As expected, the positioning of the crimped device before the deployment was found to significantly affect the implant outcomes: clinician's choices can be improved if FE models are used to optimize the pre-operative planning. Remarkably, also the wall mechanical stiffness plays an important role. However, this parameter value is unknown for a specific LAA, a crucial point that must be correctly defined for developing an accurate FE model. Finally, numerical simulations outlined how the device's configuration on which the clinician relies to assess the implant success (i.e., the deployed configuration with the device still attached to the catheter) may differ from the actual final device's configuration, relevant for achieving a safe intervention.

On the importance of fundamental computational fluid dynamics toward a robust and reliable model of left atrial flows

Computational fluid dynamics (CFD) studies of left atrial flows have reached a sophisticated level, for example, revealing plausible relationships between hemodynamics and stresses with atrial fibrillation. However, little focus has been on fundamental fluid modeling of LA flows. The purpose of this study was to investigate the spatiotemporal convergence, along with the differences between high- (HR) versus normal-resolution/accuracy (NR) solution strategies, respectively. Rigid wall CFD simulations were conducted on 12 patient-specific left atrial geometries obtained from computed tomography scans, utilizing a second-order accurate and space/time-centered solver. The convergence studies showed an average variability of around 30% and 55% for time averaged wall shear stress (WSS), oscillatory shear index (OSI), relative residence time (RRT), and endothelial cell activation potential (ECAP), even between intermediate spatial and temporal resolutions, in the left atrium (LA) and left atrial appendage (LAA), respectively. The comparison between HR and NR simulations showed good correlation in the LA for WSS, RRT, and ECAP (R 2 > .9 ), but not for OSI (R 2 = .63 ). However, there were poor correlations in the LAA especially for OSI, RRT, and ECAP (R 2 = .55, .63, and .61, respectively), except for WSS (R 2 = .81 ). The errors are comparable to differences previously reported with disease correlations. To robustly predict atrial hemodynamics and stresses, numerical resolutions of 10 M elements (i.e., Δ x = ∼ .5 mm) and 10 k time-steps per cycle seem necessary (i.e., one order of magnitude higher than normally used in both space and time). In conclusion, attention to fundamental numerical aspects is essential toward establishing a plausible, robust, and reliable model of LA flows.

The role of the pulmonary veins on left atrial flow patterns and thrombus formation

Atrial fibrillation (AF) is the most common human arrhythmia, forming thrombi mostly in the left atrial appendage (LAA). However, the relation between LAA morphology, blood patterns and clot formation is not yet fully understood. Furthermore, the impact of anatomical structures like the pulmonary veins (PVs) have not been thoroughly studied due to data acquisition difficulties. In-silico studies with flow simulations provide a detailed analysis of blood flow patterns under different boundary conditions, but a limited number of cases have been reported in the literature. To address these gaps, we investigated the influence of PVs on LA blood flow patterns and thrombus formation risk through computational fluid dynamics simulations conducted on a sizeable cohort of 130 patients, establishing the largest cohort of patient-specific LA fluid simulations reported to date. The investigation encompassed an in-depth analysis of several parameters, including pulmonary vein orientation (e.g., angles) and configuration (e.g., number), LAA and LA volumes as well as their ratio, flow, and mass-less particles. Our findings highlight the total number of particles within the LAA as a key parameter for distinguishing between the thrombus and non-thrombus groups. Moreover, the angles between the different PVs play an important role to determine the flow going inside the LAA and consequently the risk of thrombus formation. The alignment between the LAA and the main direction of the left superior pulmonary vein, or the position of the right pulmonary vein when it exhibits greater inclination, had an impact to distinguish the control group vs. the thrombus group. These insights shed light on the intricate relationship between PV configuration, LAA morphology, and thrombus formation, underscoring the importance of comprehensive blood flow pattern analyses.

Personalized biomechanical insights in atrial fibrillation: opportunities & challenges

INTRODUCTION

Atrial fibrillation (AF) is an increasingly prevalent and significant worldwide health problem. Manifested as an irregular atrial electrophysiological activation, it is associated with many serious health complications. AF affects the biomechanical function of the heart as contraction follows the electrical activation, subsequently leading to reduced blood flow. The underlying mechanisms behind AF are not fully understood, but it is known that AF is highly correlated with the presence of atrial fibrosis, and with a manifold increase in risk of stroke.

AREAS COVERED

In this review, we focus on biomechanical aspects in atrial fibrillation, current and emerging use of clinical images, and personalized computational models. We also discuss how these can be used to provide patient-specific care.

EXPERT OPINION

Understanding the connection betweenatrial fibrillation and atrial remodeling might lead to valuable understanding of stroke and heart failure pathophysiology. Established and emerging imaging modalities can bring us closer to this understanding, especially with continued advancements in processing accuracy, reproducibility, and clinical relevance of the associated technologies. Computational models of cardiac electromechanics can be used to glean additional insights on the roles of AF and remodeling in heart function.

Uncertainty Quantification in the In Vivo Image-Based Estimation of Local Elastic Properties of Vascular Walls

Introduction: Patient-specific computational models are a powerful tool for planning cardiovascular interventions. However, the in vivo patient-specific mechanical properties of vessels represent a major source of uncertainty. In this study, we investigated the effect of uncertainty in the elastic module (E) on a Fluid–Structure Interaction (FSI) model of a patient-specific aorta. Methods: The image-based χ-method was used to compute the initial E value of the vascular wall. The uncertainty quantification was carried out using the generalized Polynomial Chaos (gPC) expansion technique. The stochastic analysis was based on four deterministic simulations considering four quadrature points. A deviation of about ±20% on the estimation of the E value was assumed. Results: The influence of the uncertain E parameter was evaluated along the cardiac cycle on area and flow variations extracted from five cross-sections of the aortic FSI model. Results of stochastic analysis showed the impact of E in the ascending aorta while an insignificant effect was observed in the descending tract. Conclusions: This study demonstrated the importance of the image-based methodology for inferring E, highlighting the feasibility of retrieving useful additional data and enhancing the reliability of in silico models in clinical practice.

Imaging and biophysical modelling of thrombogenic mechanisms in atrial fibrillation and stroke

Atrial fibrillation (AF) underlies almost one third of all ischaemic strokes, with the left atrial appendage (LAA) identified as the primary thromboembolic source. Current stroke risk stratification approaches, such as the CHA₂DS₂-VASc score, rely mostly on clinical comorbidities, rather than thrombogenic mechanisms such as blood stasis, hypercoagulability and endothelial dysfunction-known as Virchow's triad. While detection of AF-related thrombi is possible using established cardiac imaging techniques, such as transoesophageal echocardiography, there is a growing need to reliably assess AF-patient thrombogenicity prior to thrombus formation. Over the past decade, cardiac imaging and image-based biophysical modelling have emerged as powerful tools for reproducing the mechanisms of thrombogenesis. Clinical imaging modalities such as cardiac computed tomography, magnetic resonance and echocardiographic techniques can measure blood flow velocities and identify LA fibrosis (an indicator of endothelial dysfunction), but imaging remains limited in its ability to assess blood coagulation dynamics. In-silico cardiac modelling tools-such as computational fluid dynamics for blood flow, reaction-diffusion-convection equations to mimic the coagulation cascade, and surrogate flow metrics associated with endothelial damage-have grown in prevalence and advanced mechanistic understanding of thrombogenesis. However, neither technique alone can fully elucidate thrombogenicity in AF. In future, combining cardiac imaging with in-silico modelling and integrating machine learning approaches for rapid results directly from imaging data will require development under a rigorous framework of verification and clinical validation, but may pave the way towards enhanced personalised stroke risk stratification in the growing population of AF patients. This Review will focus on the significant progress in these fields.

Left atrial appendage shape impacts on the left atrial flow hemodynamics: A numerical hypothesis generating study on two cases

BACKGROUND AND OBJECTIVES

The left atrial appendage (LAA) is the most common region for thrombus formation in atrial fibrillation (AF). Morphological parameters such as shape, size, and LAA volume can cause insufficient effectiveness of available therapeutic options. This study aimed to examine blood flow inside LAA and its removal effects. Computational fluid dynamic (CFD) simulations were carried out on two patients with different morphologies.

METHODS

Two patients' CT was used to reconstruct the 3D geometries of the left atrium (LA) and left atrial appendage (LAA). Then, the geometries were refined in the mentioned software, and the LAA in some models was removed. Next, in generated 3D volume mesh, sinus rhythm (SR) and atrial fibrillation (AF) outflow velocity were imposed at the mitral valve as boundary conditions. Finally, CFD simulation was conducted to analyzing blood flow within LA with/without LA.

RESULTS

The results confirmed that velocity and vorticity decreased under AF conditions inside the LA domain for both patients. However, removing LAA may cause unpredictable consequences, due to different shape and volume of LAA. LAA removal had insignificant effects on velocity and vorticity within LA in SR-mitral outflow. However, removing LAA increased the blood flow rate by 9.15% and vorticity by 7.27% for patient one under AF rhythm (SR)-outflow. In contrast, for patient two, LAA removal in both AF and SR decreased velocity and vorticity within the LA domain. In SR-mitral outflow, velocity dropped by 18.8 %, and vorticity by 13.2%. Also, under AF velocity and vorticity decreased by 23.33% and 18.6% respectively. Meanwhile, the results indicated that the vorticity magnitude increased inside the LAA under AF associated with the risk of thrombus formation, particularly for patient one under AF. The distal part of LAA in both patients was the most common region for blood stasis because of the lowest velocity magnitude.

CONCLUSION

Overall, the morphology of LAA could be the critical parameter to determine the possibility of thrombosis formation, particularly under AF conditions. High volume, low blood flow velocity and two-lobe-appendage are more likely to have blood stasis. Furthermore, the morphology difference can affect the LAA removal result and make it more complicated. So, it could be challenging to generalize LAA removal as a therapeutic option for different patients. The implication of this CFD observation needs more investigation.

3D Printing in Modern Cardiology

BACKGROUND

3D printing represents an emerging technology in the field of cardiovascular medicine. 3D printing can help to perform a better analysis of complex anatomies to optimize intervention planning.

METHODS

A systematic review was performed to illustrate the 3D printing technology and to describe the workflow to obtain 3D printed models from patient-specific images. Examples from our laboratory of the benefit of 3D printing in planning interventions were also reported.

RESULTS

3D printing technique is reliable when applied to high-quality 3D image data (CTA, CMR, 3D echography), but it still needs the involvement of expert operators for image segmentation and mesh refinement. 3D printed models could be useful in interventional planning, although prospective studies with comprehensive and clinically meaningful endpoints are required to demonstrate the clinical utility.

CONCLUSION

3D printing can be used to improve anatomy understanding and surgical planning.

Predicting TAVI paravalvular regurgitation outcomes based on numerical simulation of the aortic annulus eccentricity and perivalvular areas

Trans-catheter aortic valve implantation (TAVI) is an increasingly adopted technique which provides a minimal invasive solution for patients who suffer from severe aortic stenosis. Some complications of the procedure could be annular rupture or paravalvular leakage, both related with adverse outcome. In TAVI with balloon expandable devices, a mismatch between those two factors leads to a conflict situation, where improving one worsens the other. The presented research proposes a methodology that uses numerical simulation to obtain certain TAVI outcomes related with aortic regurgitation due to paravalvular leakage, such as perivalvular area, aortic eccentricity or annular pressure. The application of the methodology for two patients shows the possibility of predicting those quantities. The highest stress values are distributed along the contact area. Results also show that a great deformation on the aortic annulus does not necessarily imply a higher stress; pressure can either be converted into root reshape or into root stretching. Validation of the results was done using scientific publications, clinical guidelines and clinical reports. Numerical simulation provides a suitable tool that could possibly contribute to optimize the planification procedure adjusting the mismatch between size and pressure.

Impact of left atrial appendage location on risk of thrombus formation in patients with atrial fibrillation

Most strokes in patients with atrial fibrillation (AF) are thought to arise from thrombus formation in the left atrial appendage (LAA). Assessing the hemodynamics in LAA and left atrium (LA) may provide some insights in the evaluation of the risk of thrombus formation. This study aims to find out the impact of different LAA locations with respect of LA on the risk of thrombus formation within LAA in patients with AF. Three different LAA locations at LA were modeled and a fully coupled fluid-structure interaction analysis was performed. A discrete phase method was used for particle residence analysis to evaluate risk of the thrombus formation. The results showed that LAA positions on the LA affected the LAA flow velocity distribution, passive contraction ability, and particle residence. In particular, the left pulmonary veins (PVs) had a greater influence on the LAA hemodynamics than the right PVs. The LAA had the lowest contractibility when it was located between left superior and left inferior PVs, and in this case, a larger number of particles were resided, which indicated a higher risk of thrombus formation. The present work provides a quantitative way to evaluate the risk of thrombus formation within LAA in patients with AF.

Correlation between micro and macrostructural biaxial behavior of ascending thoracic aneurysm: a novel experimental technique

Mechanical properties and microstructural modifications of vessel tissues are strongly linked, as established in the state of the art of cardiovascular diseases. Techniques to obtain both mechanical and structural information are reported, but the possibility to obtain real-time microstructural and macrostructural data correlated is still lacking. An experimental approach to characterize the aortic tissue is presented. A setup integrating biaxial traction and Small Angle Light Scattering (SALS) analysis is described. The system was adopted to test ex-vivo aorta specimens from healthy and aneusymatic (aTAA) cases. A significant variation of the fiber dispersion with respect to the unloaded state was encountered during the material traction. The corresponding microstructural and mechanical data were successfully used to fit a given anisotropic constitutive model, with satisfactory R² values (0.97±0.11 and 0.96±0.17, for aTAA and healthy population, respectively) and fiber dispersion parameters variations between the aTAA and healthy populations (0.39±0.23 and 0.15±0.10). The method integrating the biaxial/SALS technique was validated, allowing for real-time synchronization between mechanical and microstructural analysis of anisotropic biological tissues.

A novel formulation for the study of the ascending aortic fluid dynamics with in vivo data

Numerical simulations to evaluate thoracic aortic hemodynamics include a computational fluid dynamic (CFD) approach or fluid-structure interaction (FSI) approach. While CFD neglects the arterial deformation along the cardiac cycle by applying a rigid wall simplification, on the other side the FSI simulation requires a lot of assumptions for the material properties definition and high computational costs. The aim of this study is to investigate the feasibility of a new strategy, based on Radial Basis Functions (RBF) mesh morphing technique and transient simulations, able to introduce the patient-specific changes in aortic geometry during the cardiac cycle. Starting from medical images, aorta models at different phases of cardiac cycle were reconstructed and a transient shape deformation was obtained by proper activating incremental RBF solutions during the simulation process. The results, in terms of main hemodynamic parameters, were compared with two performed CFD simulations for the aortic model at minimum and maximum volume. Our implemented strategy copes the actual arterial variation during cardiac cycle with high accuracy, capturing the impact of geometrical variations on fluid dynamics, overcoming the complexity of a standard FSI approach.