Atrial fibrillation (AF) is nowadays the most common human arrhythmia and it is considered a marker of an increased risk of embolic stroke. It is known that 99% of AF-related thrombi are generated in the left atrial appendage (LAA), an anatomical structure located within the left atrium (LA). Left atrial appendage occlusion (LAAO) has become a good alternative for nonvalvular AF patients with contraindications to anticoagulants. However, there is a non-negligible number of device-related thrombus (DRT) events, created next to the device surface. In silico fluid simulations can be a powerful tool to better understand the relation between LA anatomy, haemodynamics, and the process of thrombus formation. Despite the increasing literature in LA fluid modelling, a consensus has not been reached yet in the community on the optimal modelling choices and boundary conditions for generating realistic simulations. In this line, we have performed a sensitivity analysis of several boundary conditions scenarios, varying inlet/outlet and LA wall movement configurations, using patient-specific imaging data of six LAAO patients (three of them with DRT at follow-up). Mesh and cardiac cycle convergence were also analysed. The boundary conditions scenario that better predicted DRT cases had echocardiography-based velocities at the mitral valve outlet, a generic pressure wave from an AF patient at the pulmonary vein inlets, and a dynamic mesh approach for LA wall deformation, emphasizing the need for patient-specific data for realistic simulations. The obtained promising results need to be further validated with larger cohorts, ideally with ground truth data, but they already offer unique insights on thrombogenic risk in the left atria.

Background. Atrial fibrillation (AF) is considered the most common human arrhythmia. In nonvalvular AF, around 99 of thrombi are formed in the left atrial appendage (LAA). Nevertheless, there is not a consensus in the community about the relevant factors to stratify the AF population according to thrombogenic risk. Objective. To demonstrate the need for combining left atrial morphological and haemodynamics indices to improve the thrombogenic risk assessment in nonvalvular AF patients. Methods. A cohort of 71 nonvalvular AF patients was analysed. Statistical analysis, regression models, and random forests were used to analyse the differences between morphological and haemodynamics parameters, extracted from computational simulations built on 3D rotational angiography images, between patients with and without transient ischemic attack (TIA) or cerebrovascular accident (CVA). Results. The analysis showed that models composed of both morphological and haemodynamic factors were better predictors of TIA/CVA compared with models based on either morphological or haemodynamic factors separately. Maximum ostium diameter, length of the centreline, blood flow velocity within the LAA, oscillatory shear index, and time average wall shear stress parameters were found to be key risk factors for TIA/CVA prediction. In addition, TIA/CVA patients presented more flow stagnation within the LAA. Conclusion. Thrombus formation in the LAA is the result of multiple factors. Analyses based only on morphological or haemodynamic parameters are not precise enough to predict such a phenomenon, as demonstrated in our results; a better patient stratification can be obtained by jointly analysing morphological and haemodynamic features.

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.

Left atrial appendage occlusion devices (LAAO) are a feasible alternative for non-valvular atrial fibrillation (AF) patients at high risk of thromboembolic stroke and contraindication to antithrombotic therapies. However, optimal LAAO device configurations (i.e., size, type, location) remain unstandardized due to the large anatomical variability of the left atrial appendage (LAA) morphology, leading to a 4-6% incidence of device-related thrombus (DRT). In-silico simulations have the potential to assess DRT risk and identify the key factors, such as suboptimal device positioning. This work presents fluid simulation results computed on 20 patient-specific left atrial geometries, analysing different commercially available LAAO occluders, including plug-type and pacifier-type devices. In addition, we explored two distinct device positions: 1) the real post-LAAO intervention configuration derived from follow-up imaging; and 2) one covering the pulmonary ridge if it was not achieved during the implantation (13 out of 20). In total, 33 different configurations were analysed. In-silico indices indicating high risk of DRT (e.g., low blood flow velocities and flow complexity around the device) were combined with particle deposition analysis based on a discrete phase model. The obtained results revealed that covering the pulmonary ridge with the LAAO device may be one of the key factors to prevent DRT. Moreover, disk-based devices exhibited enhanced adaptability to various LAA morphologies and, generally, demonstrated a lower risk of abnormal events after LAAO implantation.