Atrial fibrillation (AF) is a common arrhythmia mainly affecting the elderly population, which can lead to serious complications such as stroke, ischaemic attack and vascular dementia. These problems are caused by thrombi which mostly originate in the left atrial appendage (LAA), a small muscular sac protruding from left atrium. The abnormal heart rhythm associated with AF results in alterations in the heart muscle contractions and in some reshaping of the cardiac chambers. This study aims to verify if and how these physiological changes can establish hemodynamic conditions in the LAA promoting thrombus formation, by means of computational fluid dynamic (CFD) analyses. In particular, sinus and fibrillation contractility was replicated by applying wall velocity/motion to models based on healthy and dilated idealized shapes of the left atrium with a common LAA morphology. The models were analyzed and compared in terms of shear strain rate (SSR) and vorticity, which are hemodynamic parameters directly associated with thrombogenicity. The study clearly indicates that the alterations in contractility and morphology associated with AF pathologies play a primary role in establishing hemodynamic conditions which promote higher incidence of ischaemic events, consistently with the clinical evidence. In particular, in the analyzed models, the impairment in contractility determined a decrease in SSR of about 50%, whilst the chamber pathological dilatation contributed to a 30% reduction, indicating increased risk of clot formation. The equivalent rigid wall model was characterized by SSR values about one order of magnitude smaller than in the contractile models, and substantially different vortical behavior, suggesting that analyses based on rigid chambers, although common in the literature, are inadequate to provide realistic results on the LAA hemodynamics.

No abstract

No abstract

BackgroundA large majority of thrombi causing ischemic complications under atrial fibrillation (AF) originate in the left atrial appendage (LAA), an anatomical structure departing from the left atrium, characterized by a large morphological variability between individuals. This work analyses the hemodynamics simulated for different patient-specific models of LAA by means of computational fluid–structure interaction studies, modeling the effect of the changes in contractility and shape resulting from AF.MethodsThree operating conditions were analyzed: sinus rhythm, acute atrial fibrillation, and chronic atrial fibrillation. These were simulated on four patient-specific LAA morphologies, each associated with one of the main morphological variants identified from the common classification: chicken wing, cactus, windsock, and cauliflower. Active contractility of the wall muscle was calibrated on the basis of clinical evaluations of the filling and emptying volumes, and boundary conditions were imposed on the fluid to replicate physiological and pathological atrial pressures, typical of the various operating conditions.ResultsThe LAA volume and shear strain rates were analyzed over time and space for the different models. Globally, under AF conditions, all models were well aligned in terms of shear strain rate values and predicted levels of risk. Regions of low shear rate, typically associated with a higher risk of a clot, appeared to be promoted by sudden bends and focused at the trabecule and the lobes. These become substantially more pronounced and extended with AF, especially under acute conditions.ConclusionThis work clarifies the role of active and passive contraction on the healthy hemodynamics in the LAA, analyzing the hemodynamic effect of AF that promotes clot formation. The study indicates that local LAA topological features are more directly associated with a thromboembolic risk than the global shape of the appendage, suggesting that more effective classification criteria should be identified.

An efficient and accurate method is proposed to solve the incompressible flow momentum and continuity equations in computational domains partitioned into subdomains in the framework of the smoothed particle hydrodynamics method. The procedure does not require any overlap of the subdomains, which would result in the increase of the computational effort. Perfectly matching solutions are obtained at the surfaces separating neighboring blocks. The block interfaces can be both planar and curved surfaces allowing to easily decompose even geometrically complex domains.The smoothing length of the kernel function is maintained constant in each subdomain, while changing between blocks where a different resolution is required. Particles leaving each block through the interfaces are deactivated and correspondingly new particles are generated at the neighboring block using a dynamically adaptive procedure to control their frequency of release. No splitting and coalescing method is thus employed to take into account the different size and mass of the particles going through the interfaces. Mass conservation is guaranteed during the procedure, which is a challenging task in a Lagrangian method based on the domain decomposition.The test cases in both 2D and 3D approximation show the accuracy of the method and its ability to strongly reduce the computational efforts through a multi-resolution approach.