Numerical analysis of hemodynamic changes in the left atrium due to atrial fibrillation

Atrial fibrillation increases thrombogenicity of LVAD therapy

BACKGROUND

This study investigates the hypothesis that presence of atrial fibrillation (AF) in LVAD patients increases thrombogenicity in the left ventricle (LV) and exacerbates stroke risk.

METHODS

Using an anatomical LV model implanted with an LVAD inflow cannula, we analyze thrombogenic risk and blood flow patterns in either AF or sinus rhythm (SR) using unsteady computational fluid dynamics (CFD). To analyze platelet activation and thrombogenesis in the LV, hundreds of thousands of platelets are individually tracked to quantify platelet residence time (RT) and shear stress accumulation history (SH).

RESULTS

The irregular and chaotic mitral inflow associated with AF results in markedly different intraventricular flow patterns, with profoundly negative impact on blood flow-induced stimuli experienced by platelets as they traverse the LV. Twice as many platelets accumulated very high SH in the LVAD + AF case, resulting in a 36% increase in thrombogenic potential score, relative to the LVAD + SR case.

CONCLUSIONS

This supports the hypothesis that AF results in unfavorable blood flow patterns in the LV adding to an increased stroke risk for LVAD + AF patients. Quantification of thrombogenic risk associated with AF for LVAD patients may help guide clinical decision-making on interventions to mitigate the increased risk of thromboembolic events.

A comprehensive stroke risk assessment by combining atrial computational fluid dynamics simulations and functional patient data

Stroke, a major global health concern often rooted in cardiac dynamics, demands precise risk evaluation for targeted intervention. Current risk models, like theCHA 2 DS 2 -VASc score, often lack the granularity required for personalized predictions. In this study, we present a nuanced and thorough stroke risk assessment by integrating functional insights from cardiac magnetic resonance (CMR) with patient-specific computational fluid dynamics (CFD) simulations. Our cohort, evenly split between control and stroke groups, comprises eight patients. Utilizing CINE CMR, we compute kinematic features, revealing smaller left atrial volumes for stroke patients. The incorporation of patient-specific atrial displacement into our hemodynamic simulations unveils the influence of atrial compliance on the flow fields, emphasizing the importance of LA motion in CFD simulations and challenging the conventional rigid wall assumption in hemodynamics models. Standardizing hemodynamic features with functional metrics enhances the differentiation between stroke and control cases. While standalone assessments provide limited clarity, the synergistic fusion of CMR-derived functional data and patient-informed CFD simulations offers a personalized and mechanistic understanding, distinctly segregating stroke from control cases. Specifically, our investigation reveals a crucial clinical insight: normalizing hemodynamic features based on ejection fraction fails to differentiate between stroke and control patients. Differently, when normalized with stroke volume, a clear and clinically significant distinction emerges and this holds true for both the left atrium and its appendage, providing valuable implications for precise stroke risk assessment in clinical settings. This work introduces a novel framework for seamlessly integrating hemodynamic and functional metrics, laying the groundwork for improved predictive models, and highlighting the significance of motion-informed, personalized risk assessments.

Blood flow analysis with computational fluid dynamics in the left atrium after left atrial plication: a prospective study

OBJECTIVE

This study aimed to evaluate blood flow stagnation in an enlarged left atrium (LA) and prove that left atrial plication (LAP) could alleviate the stagnation.

METHODS

Five patients with chronic atrial fibrillation who underwent mitral valve surgery followed by LAP for an enlarged LA with a ≥ 60-mm diameter were included. We performed computational fluid dynamics (CFD) analysis using preoperative and postoperative computed tomography and four-dimensional flow magnetic resonance imaging. Additionally, computer graphics were used to create virtual left atrial appendage resection (LAAR) images. We performed CFD analysis to assess blood flow stagnation in the LA for three groups: preoperative, LAAR, and LAP.

RESULTS

When the average and constant stagnation volumes were both set to 100 preoperatively, the average stagnation volumes of the LAAR and LAP groups were 67.42 ± 18.64 and 35.88 ± 8.20, respectively. The constant stagnation volumes of these groups reduced to 45.01 ± 7.43 and 21.14 ± 7.70, respectively. The LAP group also had significantly lower average and constant stagnation volumes than those in the LAAR group (p = 0.006 and p = 0.033, respectively).

CONCLUSIONS

Blood flow stagnation was noted in the LAA and enlarged LA. CFD analysis revealed that LAP for the enlarged LA improved blood flow stagnation more than the virtual LAAR alone.

CLINICAL TRIAL REGISTRY NUMBER

UMIN000049923.

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 impact of contrast retention on thrombus formation risks in patients with atrial fibrillation: A numerical study

Background

Contrast retention (CR) is an important predictor of left atrial appendage thrombus (LAAT) and stroke in patients with non-valvular atrial fibrillation (AF). We sought to explore the underlying mechanisms of CR using computational fluid dynamic (CFD) simulations.

Methods

A total of 12 patients with AF who underwent both cardiac computed tomography angiography (CTA) and transesophageal echocardiography (TEE) before left atrial appendage occlusion (LAAO) were included in the study. The patients were allocated into the CR group or non-CR group based on left atrial appendage (LAA) angiography. Patient-specific models were reconstructed to evaluate time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), and endothelial cell activation potential (ECAP). Additionally, the incidence of thrombosis was predicted using residence time (RT) at different time-points.

Results

TAWSS was lower [median (Interquartile Range) 0.27 (0.19-0.47) vs 1.35 (0.92-1.79), p < 0.001] in LAA compared to left atrium. In contrast, RRT [1438 (409.70-13869) vs 2.23 (1.81-3.14), p < 0.001] and ECAP [122.70 (30.01-625.70) vs 0.19 (0.16-0.27), p < 0.001)] was higher in the LAA. The patients in the CR group had significantly higher RRT [(mean ± SD) 16274 ± 11797 vs 639.70 ± 595.20, p = 0.009] and ECAP [610.80 ± 365.30 vs 54.26 ± 54.38, p = 0.004] in the LAA compared to the non-CR group. Additionally, patients with CR had a wider range of thrombus-prone regions [0.44(0.27-0.66)% vs 0.05(0.03-0.27)%, p = 0.009] at the end of the 15th cardiac cycle.

Conclusions

These findings suggest that CR might be an indicator of high-risk thrombus formation in the LAA. And CT-based CFD simulation may be a feasible substitute for the evaluation of LAA thrombotic risk in patients with AF, especially in patients with CR.

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.

A comprehensive stroke risk assessment by combining atrial computational fluid dynamics simulations and functional patient data

Stroke, a major global health concern often rooted in cardiac dynamics, demands precise risk evaluation for targeted intervention. Current risk models, like the CHA2DS2-VASc score, often lack the granularity required for personalized predictions. In this study, we present a nuanced and thorough stroke risk assessment by integrating functional insights from cardiac magnetic resonance (CMR) with patient-specific computational fluid dynamics (CFD) simulations. Our cohort, evenly split between control and stroke groups, comprises eight patients. Utilizing CINE CMR, we compute kinematic features, revealing smaller left atrial volumes for stroke patients. The incorporation of patient-specific atrial displacement into our hemodynamic simulations unveils the influence of atrial compliance on the flow fields, emphasizing the importance of LA motion in CFD simulations and challenging the conventional rigid wall assumption in hemodynamics models. Standardizing hemodynamic features with functional metrics enhances the differentiation between stroke and control cases. While standalone assessments provide limited clarity, the synergistic fusion of CMR-derived functional data and patient-informed CFD simulations offers a personalized and mechanistic understanding, distinctly segregating stroke from control cases. Specifically, our investigation reveals a crucial clinical insight: normalizing hemodynamic features based on ejection fraction fails to differentiate between stroke and control patients. Differently, when normalized with stroke volume, a clear and clinically significant distinction emerges and this holds true for both the left atrium and its appendage, providing valuable implications for precise stroke risk assessment in clinical settings. This work introduces a novel framework for seamlessly integrating hemodynamic and functional metrics, laying the groundwork for improved predictive models, and highlighting the significance of motion-informed, personalized risk assessments.

Pulmonary vein flow split effects in patient-specific simulations of left atrial flow

Disruptions to left atrial (LA) blood flow, such as those caused by atrial fibrillation (AF), can lead to thrombosis in the left atrial appendage (LAA) and an increased risk of systemic embolism. LA hemodynamics are influenced by various factors, including LA anatomy and function, and pulmonary vein (PV) inflow conditions. In particular, the PV flow split can vary significantly among and within patients depending on multiple factors. In this study, we investigated how changes in PV flow split affect LA flow transport, focusing for the first time on blood stasis in the LAA, using a high-fidelity patient-specific computational fluid dynamics (CFD) model. We use an Immersed Boundary Method, simulating the flow in a fixed, uniform Cartesian mesh and imposing the movement of the LA walls with a moving Lagrangian mesh generated from 4D Computerized Tomography images. We analyzed LA anatomies from eight patients with varying atrial function, including three with AF and either a LAA thrombus or a history of Transient Ischemic Attacks (TIAs). Using four different flow splits (60/40% and 55/45% through right and left PVs, even flow rate, and same velocity through each PV), we found that flow patterns are sensitive to PV flow split variations, particularly in planes parallel to the mitral valve. Changes in PV flow split also had a significant impact on blood stasis and could contribute to increased risk for thrombosis inside the LAA, particularly in patients with AF and previous LAA thrombus or a history of TIAs. Our study highlights the importance of considering patient-specific PV flow split variations when assessing LA hemodynamics and identifying patients at increased risk for thrombosis and stroke. This knowledge is relevant to planning clinical procedures such as AF ablation or the implementation of LAA occluders.

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.

Machine Learning and the Conundrum of Stroke Risk Prediction

Stroke is a leading cause of death worldwide. With escalating healthcare costs, early non-invasive stroke risk stratification is vital. The current paradigm of stroke risk assessment and mitigation is focused on clinical risk factors and comorbidities. Standard algorithms predict risk using regression-based statistical associations, which, while useful and easy to use, have moderate predictive accuracy. This review summarises recent efforts to deploy machine learning (ML) to predict stroke risk and enrich the understanding of the mechanisms underlying stroke. The surveyed body of literature includes studies comparing ML algorithms with conventional statistical models for predicting cardiovascular disease and, in particular, different stroke subtypes. Another avenue of research explored is ML as a means of enriching multiscale computational modelling, which holds great promise for revealing thrombogenesis mechanisms. Overall, ML offers a new approach to stroke risk stratification that accounts for subtle physiologic variants between patients, potentially leading to more reliable and personalised predictions than standard regression-based statistical associations.

Subject-specific factors affecting particle residence time distribution of left atrial appendage in atrial fibrillation: A computational model-based study

Background

Atrial fibrillation (AF) is a prevalent arrhythmia, that causes thrombus formation, ordinarily in the left atrial appendage (LAA). The conventional metric of stroke risk stratification, CHA2DS2-VASc score, does not account for LAA morphology or hemodynamics. We showed in our previous study that residence time distribution (RTD) of blood-borne particles in the LAA and its associated calculated variables (i.e., mean residence time, tm, and asymptotic concentration, C∞) have the potential to improve CHA2DS2-VASc score. The purpose of this research was to investigate the effects of the following potential confounding factors on LAA tm and C∞: (1) pulmonary vein flow waveform pulsatility, (2) non-Newtonian blood rheology and hematocrit level, and (3) length of the simulation.

Methods

Subject-Specific data including left atrial (LA) and LAA cardiac computed tomography, cardiac output (CO), heart rate, and hematocrit level were gathered from 25 AF subjects. We calculated LAA tm and C∞ based on series of computational fluid dynamics (CFD) analyses.

Results

Both LAA tm and C∞ are significantly affected by the CO, but not by temporal pattern of the inlet flow. Both LAA tm and C∞ increase with increasing hematocrit level and both calculated indices are higher for non-Newtonian blood rheology for a given hematocrit level. Further, at least 20,000 s of CFD simulation is needed to calculate LAA tm and C∞ values reliably.

Conclusions

Subject-specific LA and LAA geometries, CO, and hematocrit level are essential to quantify the subject-specific proclivity of blood cell tarrying inside LAA in terms of the RTD function.

Medical record data-enabled machine learning can enhance prediction of left atrial appendage thrombosis in nonvalvular atrial fibrillation

BACKGROUND

As a major complication of non-valvular atrial fibrillation (NVAF), left atrial appendage (LAA) thrombosis is associated with cerebral ischemic strokes, as well as high morbidity. Due to insufficient incorporation of risk factors, most current scoring methods are limited to the analysis of relationships between clinical characteristics and LAA thrombosis rather than detecting potential risk. Therefore, this study proposes a clinical data-driven machine learning method to predict LAA thrombosis of NVAF.

METHODS

Patients with NVAF from January 2014 to June 2022 were enrolled from Southwest Hospital. We selected 40 variables for analysis, including demographic data, medical history records, laboratory results, and the structure of LAA. Three machine learning algorithms were adopted to construct classifiers for the prediction of LAA thrombosis risk. The most important variables related to LAA thrombosis and their influences were recognized by SHapley Addictive exPlanations method. In addition, we compared our model with CHADS2 and CHADS2-VASc scoring methods.

RESULTS

A total of 713 participants were recruited, including 127 patients with LAA thrombosis and 586 patients with no obvious thrombosis. The consensus models based on Random Forest and eXtreme Gradient Boosting LAA thrombosis prediction (RXTP) achieved the best accuracy of 0.865, significantly outperforming CHADS2 score and CHA2DS2-VASc score (0.757 and 0.754, respectively). The SHAP results showed that B-type natriuretic peptide, left atrial appendage width, C-reactive protein, Fibrinogen and estimated glomerular filtration rate are closely related to the risk of LAA thrombosis in nonvalvular atrial fibrillation.

CONCLUSIONS

The RXTP-NVAF model is the most effective model with the greatest ROC value and recall rate. The summarized risk factors obtained from SHAP enable the optimization of the treatment strategy, thereby preventing thromboembolism events and the occurrence of cardiogenic ischemic stroke.

Prognostic impact of plasma xanthine oxidoreductase activity in patients with heart failure with atrial fibrillation

BACKGROUND

Xanthine oxidoreductase (XOR) is a rate-limiting enzyme for uric acid (UA) production and plays an important role in generating reactive oxygen species (ROS). Overproduction of ROS is reported to contribute to the pathophysiology of atrial fibrillation (AF), however, the prognostic impact of plasma XOR activity in patients with heart failure (HF) with AF is undetermined.

METHODS

We measured plasma XOR activity in 475 HF patients, including those with sinus rhythm (HF-SR, n = 211), and those with AF (HF-AF, n = 264). The type of AF included paroxysmal (n = 128) and persistent (n = 136) AF. All patients were prospectively followed up for a median period of 804 days.

RESULTS

HF-AF patients had significantly higher plasma XOR activity and serum UA levels compared with HF-SR patients. Both plasma XOR activity and serum UA levels were higher in patients with persistent AF than in those with SR and with paroxysmal AF. Multivariate linear regression analysis showed that persistent AF was independently associated with increased XOR activity. During the follow-up period, there were 79 major adverse cardiovascular events (MACEs). HF-AF patients with MACEs had higher plasma XOR activity compared with those without MACEs, while there were no significant differences in serum UA levels. Multivariate Cox proportional analysis showed that high XOR activity was an independent risk factor for MACEs after adjustment for confounding factors. Kaplan-Meier analysis revealed that the high XOR activity group had a higher risk of MACEs than the low XOR activity group. The prediction model was significantly improved by the addition of XOR activity to the basic predictors.

CONCLUSIONS

HF-AF patients had significantly higher plasma XOR activity compared with HF-SR patients. Plasma XOR activity proved to be a reliable indicator for MACEs in HF-AF patients.

Comparative Analysis of Temperature Rise between Convective Heat Transfer Method and Computational Fluid Dynamics Method in an Anatomy-Based Left Atrium Model during Pulsed Field Ablation: A Computational Study

The non-thermal effects are considered one of the prominent advantages of pulsed field ablation (PFA). However, at higher PFA doses, the temperature rise in the tissue during PFA may exceed the thermal damage threshold, at which time intracardiac pulsatile blood flow plays a crucial role in suppressing this temperature rise. This study aims to compare the effect of heat dissipation of the different methods in simulating the pulsatile blood flow during PFA. This study first constructed an anatomy-based left atrium (LA) model and then applied the convective heat transfer (CHT) method and the computational fluid dynamics (CFD) method to the model, respectively, and the thermal convective coefficients used in the CHT method are 984 (W/m²*K) (blood-myocardium interface) and 4372 (W/m²*K) (blood-catheter interface), respectively. Then, it compared the effect of the above two methods on the maximum temperature of myocardium and blood, as well as the myocardial ablation volumes caused by irreversible electroporation (IRE) and hyperthermia under different PFA parameters. Compared with the CFD method, the CHT method underestimates the maximum temperature of myocardium and blood; the differences in the maximum temperature of myocardium and blood between the two methods at the end of the last pulse are significant (>1 °C), and the differences in the maximum temperature of blood at the end of the last pulse interval are significant (>1 °C) only at a pulse amplitude greater than 1000 V or pulse number greater than 10. Under the same pulse amplitude and different heat dissipation methods, the IRE ablation volumes are the same. Compared with the CFD method, the CHT method underestimates the hyperthermia ablation volume; the differences in the hyperthermia ablation volume are significant (>1 mm³) only at a pulse amplitude greater than 1000 V, a pulse interval of 250 ms, or a pulse number greater than 10. Additionally, the hyperthermia ablation isosurfaces are completely wrapped by the IRE ablation isosurfaces in the myocardium. Thus, during PFA, compared with the CFD method, the CHT method cannot accurately simulate the maximum myocardial temperature; however, except at the above PFA parameters, the CHT method can accurately simulate the maximum blood temperature and the myocardial ablation volume caused by IRE and hyperthermia. Additionally, within the range of the PFA parameters used in this study, the temperature rise during PFA may not lead to the appearance of additional hyperthermia ablation areas beyond the IRE ablation area in the myocardium.

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.

Preliminary findings on left atrial appendage occlusion simulations applying different endocardial devices

Atrial fibrillation (AF) is one of the most investigated arrhythmias since it is associated with a five-fold increase in the risk of strokes. Left atrium dilation and unbalanced and irregular contraction caused by AF favour blood stasis and, consequently, stroke risk. The left atrial appendage (LAA) is the site of the highest clots formation, increasing the incidence of stroke in AF population. For many years oral anticoagulation therapy has been the most used AF treatment option available to decrease stroke risk. Unfortunately, several contraindications including bleeding risk increase, interference with other drugs and with multiorgan functioning, might outweigh its remarkable benefits on thromboembolic events. For these reasons, in recent years, other approaches have been designed, including LAA percutaneous closure. Unfortunately, nowadays, LAA occlusion (LAAO) is restricted to small subgroups of patients and require a certain level of expertise and training to successfully complete the procedure without complications. The most critical clinical problems associated with LAAO are represented by peri-device leaks and device related thrombus (DRT). The anatomical variability of the LAA plays a key role in the choice of the correct LAA occlusion device and in its correct positioning with respect to the LAA ostium during the implant. In this scenario, computational fluid dynamics (CFD) simulations could have a crucial role in improving LAAO intervention. The aim of this study was to simulate the fluid dynamics effects of LAAO in AF patients to predict hemodynamic changes due to the occlusion. LAAO was simulated by applying two different types of closure devices based on the plug and the pacifier principles on 3D LA anatomical models derived from real clinical data in five AF patients. CFD simulations were performed on the left atrium model before and after the LAAO intervention with each device. Blood velocity, particle washout and endothelial damage were computed to quantify flow pattern changes after the occlusion in relation to the thrombogenic risk. Our preliminary results confirmed an improved blood washout after the simulated implants and the capability of foreseeing thrombogenic risk based on endothelial damage and maximum blood velocities in different scenarios. This tool may help to identify effective device configurations in limiting stroke risk for patient-specific LA morphologies.

Cerebral hemodynamics during atrial fibrillation: Computational fluid dynamics analysis of lenticulostriate arteries using 7 T high-resolution magnetic resonance imaging

Atrial fibrillation (AF) is the most common cardiac arrhythmia, inducing irregular and faster heart beating. Aside from disabling symptoms-such as palpitations, chest discomfort, and reduced exercise capacity-there is growing evidence that AF increases the risk of dementia and cognitive decline, even in the absence of clinical strokes. Among the possible mechanisms, the alteration of deep cerebral hemodynamics during AF is one of the most fascinating and least investigated hypotheses. Lenticulostriate arteries (LSAs)-small perforating arteries perpendicularly departing from the anterior and middle cerebral arteries and supplying blood flow to basal ganglia-are especially involved in silent strokes and cerebral small vessel diseases, which are considered among the main vascular drivers of dementia. We propose for the first time a computational fluid dynamics analysis to investigate the AF effects on the LSAs hemodynamics by using 7 T high-resolution magnetic resonance imaging (MRI). We explored different heart rates (HRs)-from 50 to 130 bpm-in sinus rhythm and AF, exploiting MRI data from a healthy young male and internal carotid artery data from validated 0D cardiovascular-cerebral modeling as inflow condition. Our results reveal that AF induces a marked reduction of wall shear stress and flow velocity fields. This study suggests that AF at higher HR leads to a more hazardous hemodynamic scenario by increasing the atheromatosis and thrombogenesis risks in the LSAs region.

Impact of atrial fibrillation on left atrium haemodynamics: A computational fluid dynamics study

We analyse the haemodynamics of the left atrium, highlighting differences between healthy individuals and patients affected by atrial fibrillation. The computational study is based on patient-specific geometries of the left atria to simulate blood flow dynamics. We design a novel procedure to compute the boundary data for the 3D haemodynamic simulations, which are particularly useful in absence of data from clinical measurements. With this aim, we introduce a parametric definition of atrial displacement, and we use a closed-loop lumped parameter model of the whole cardiovascular circulation conveniently tuned on the basis of the patient's characteristics. We evaluate several fluid dynamics indicators for atrial haemodynamics, validating our numerical results in terms of clinical measurements; we investigate the impact of geometric and clinical characteristics on the risk of thrombosis. To highlight the correlation of thrombus formation with atrial fibrillation, according to medical evidence, we propose a novel indicator: age stasis. It arises from the combination of Eulerian and Lagrangian quantities. This indicator identifies regions where slow flow cannot properly rinse the chamber, accumulating stale blood particles, and creating optimal conditions for clots formation.

Effect of Anisotropic Electrical Conductivity Induced by Fiber Orientation on Ablation Characteristics of Pulsed Field Ablation in Atrial Fibrillation Treatment: A Computational Study

Pulsed field ablation (PFA) is a promising new ablation modality for the treatment of atrial fibrillation (AF); however, the effect of fiber orientation on the ablation characteristics of PFA in AF treatment is still unclear, which is likely an essential factor in influencing the ablation characteristics. This study constructed an anatomy-based left atrium (LA) model incorporating fiber orientation and selected various electrical conductivity and ablation targets to investigate the effect of anisotropic electrical conductivity (AC), compared with isotropic electrical conductivity (IC), on the ablation characteristics of PFA in AF treatment. The results show that the percentage differences in the size of the surface ablation area between AC and IC are greater than 73.71%; the maximum difference in the size of the ablation isosurface between AC and IC at different locations in the atrial wall is 3.65 mm (X-axis), 3.65 mm (Z-axis), and 4.03 mm (X-axis), respectively; and the percentage differences in the size of the ablation volume are greater than 6.9%. Under the condition of the pulse, the amplitude is 1000 V, the total PFA duration is 1 s, and the pulse train interval is 198.4 ms; the differences in the temperature increase between AC and IC in LA are less than 2.46 °C. Hence, this study suggests that in further exploration of the computational study of PFA in AF treatment using the same or similar conditions as those used here (myocardial electrical conductivity, pulse parameters, and electric field intensity damage threshold), to obtain more accurate computational results, it is necessary to adopt AC rather than IC to investigate the size of the surface ablation area, the size of the ablation isosurface, or the size of the ablation volume generated by PFA in LA. Moreover, if only investigating the temperature increase generated by PFA in LA, adopting IC instead of AC for simplifying the model construction process is reasonable.

Computational Fluid Dynamic Model of Left Atrium to Analyze Hemodynamic Manifestation during Atrial Fibrillation

In this paper, we present a computational fluid dynamic (CFD) model of left atrium (LA) to analyze the manifestation and progression of atrial fibrillation (AF) in terms of hemodynamic metrics. We propose a coupled lumped-CFD (0d-3d) pipeline to model and predict the pulsatile flow and pressure fields of three-dimensional cardiac chamber under the influence of sinus rhythm, high frequency AF (HF-AF) and LA remodeled AF, considering the interactions between the heart and the arterial system through a separately modeled 0d lumped hemodynamic cardiac model. A novel rhythm generator is modeled to generate modulated cardiac chamber compliance and decoupled auricular and ventricular contraction rate to synthesize variation in sinus rhythm and subsequent AF generation. CFD simulation were solved using subject specific CT scan. Systemic and pulmonary flow and pressure along with metrics related to wall shear stress in LA were derived. Left ventricular (LV) hemodynamic parameters associated with global cardio vascular evaluation like ejection fraction, stroke volume, cardiac output, etc. were also generated for all the rhythmic disturbance under consideration. The proposed 0d-3d coupled hemodynamic model of the LA can provide useful insights on the dynamics of AF manifestation and predict vulnerable regions in the cardiac chambers as well as arterial vasculature for probable thrombogenic plaque formation that leads to stroke and infraction, leading to heart failure.

Inversion of Left Atrial Appendage Will Cause Compressive Stresses in the Tissue: Simulation Study of Potential Therapy

In atrial fibrillation (AF), thromboembolic events can result from the particular conformation of the left atrial appendage (LAA) bearing increased clot formation and accumulation. Current therapies to reduce the risk of adverse events rely on surgical exclusion or percutaneous occlusion, each of which has drawbacks limiting application and efficacy. We sought to quantify the hemodynamic and structural loads of a novel potential procedure to partially invert the “dead” LAA space to eliminate the auricle apex where clots develop. A realistic left atrial geometry was first achieved from the heart anatomy of the Living Heart Human Model (LHHM) and then the left atrial appendage inversion (LAAI) was simulated by finite-element analysis. The LAAI procedure was simulated by pulling the elements at the LAA tip and prescribing a displacement motion along a predefined path. The deformed configuration was then used to develop a computational flow analysis of LAAI. Results demonstrated that the inverted LAA wall undergoes a change in the stress distribution from tensile to compressive in the inverted appendage, and this can lead to resorption of the LAA tissue as per a reduced stress/resorption relationship. Computational flow analyses highlighted a slightly nested low-flow velocity pattern for the inverted LAA with minimal differences from that of a model without inversion of the LAA apex. Our study revealed important insights into the biomechanics of LAAI and demonstrated the inversion of the stress field (from tensile to compressive), which &can ultimately lead the long-term resorption of the LAA.

Fractal Geometry Illustrated Left Atrial Appendage Morphology That Predicted Thrombosis and Stroke in Patients With Atrial Fibrillation

Background

This study aims to correlate the morphological complexity of left atrial appendage (LAA) with thrombosis and stroke in patients with atrial fibrillation (AF).

Methods

The training cohort consisted of 46 patients with AF (age 55.8 ± 7.2 years, 73.9% men) who were referred for radiofrequency catheter ablation. An independent validation cohort consisting of 443 patients with AF was enrolled for further verification. All patients in the training cohort underwent both transesophageal echocardiography (TEE) and enhanced computed tomography (CT). Fractal dimension (FD) analysis was performed to evaluate the morphological complexity of LAAs quantitatively. Clinical and imaging manifestations, FD of LAAs, and diagnostic accuracy were investigated and compared between patients with AF in both training and validation cohorts.

Results

In the training cohort, LAAs (n = 22) with thrombi had significantly higher FD than those without thrombi (n = 24) h 0.44 ± 0.07 vs. 2.35 ± 0.11, p = 0.003). Receiver-operating characteristic (ROC) analysis suggested that the diagnostic accuracy of FD combined with a CHA₂DS₂-VaSc score was significantly higher than that of the CHA₂DS₂-VaSc score alone in low- to moderate-risk patients with AF (area under the curve 0.8479 vs. 0.6958, p = 0.009). The results were also validated in an independent external validation cohort and demonstrated that increased FD was associated with stroke. Hemodynamic analysis revealed that LAAs with thrombi and high FD were prone to blood stasis and lower blood flow rate.

Conclusion

LAA morphological complexity is closely associated with thrombosis and stroke in patients with paroxysmal AF. A new risk assessment system combining CHA₂DS₂-VaSc score and FD has a higher diagnostic accuracy in predicting LAA thrombosis.

Evaluation of Stroke Risk in Patients With Atrial Fibrillation Using Morphological and Hemodynamic Characteristics

Background

It is well known that the thrombus triggering stroke in patients with atrial fibrillation (AF) mainly comes from the left atrial appendage (LAA). This study aims to characterize the morphological and hemodynamic parameters and evaluate their differences between AF patients with and without a stroke history.

Methods

Cardiac CT images were obtained from AF patients with (n = 10) and without a history of stroke (n = 10). 3D models of the left atrium (LA) were reconstructed by processing the CT image, and the LA/LAA morphological parameters were measured. Computational fluid dynamics (CFD) simulations were performed to calculate the hemodynamic parameters in LA. The species-transport model and discrete phase model (DPM) were applied to analyze blood residual ratio and particle residual ratio, two qualitative parameters for thrombus formation and flow-out potential, respectively.

Results

There were significant differences in LAA actual depth (p = 0.002), and direct length (p = 0.049) between the non-stroke and stroke groups. Significant differences were also found in certain hemodynamic parameters. The blood residual ratio in LAA was significantly smaller in the stroke group than in the non-stroke group (p < 0.05). Moreover, the particle residual ratio within LAA was significantly smaller in the stroke groups than in the non-stroke group (p < 0.05).

Conclusion

There are significant differences in both morphological and hemodynamic parameters between AF patients with and without a stroke history. A high blood residual ratio in LAA confirms that thrombus is more likely to form in AF patients. A significantly smaller particle residual ratio in the stroke group may suggest the thrombus formed with LAA is more likely to flow out of LAA, leading to a higher risk of stroke. The proposed morphological and hemodynamic parameters may be additional risk factors that can be used to better risk stratify AF patients.

A Detailed Study to Discover the Trade between Left Atrial Blood Flow, Expression of Calcium-Activated Potassium Channels and Valvular Atrial Fibrillation

Background: The present study aimed to explore the correlation between calcium-activated potassium channels, left atrial flow field mechanics, valvular atrial fibrillation (VAF), and thrombosis. The process of transforming mechanical signals into biological signals has been revealed, which offers new insights into the study of VAF. Methods: Computational fluid dynamics simulations use numeric analysis and algorithms to compute flow parameters, including turbulent shear stress (TSS) and wall pressure in the left atrium (LA). Real-time PCR and western blotting were used to detect the mRNA and protein expression of IKCa2.3/3.1, ATK1, and P300 in the left atrial tissue of 90 patients. Results: In the valvular disease group, the TSS and wall ressure in the LA increased, the wall pressure increased in turn in all disease groups, mainly near the mitral valve and the posterior portion of the LA, the increase in TSS was the most significant in each group near the mitral valve, and the middle and lower part of the back of the LA and the mRNA expression and protein expression levels of IKCa2.3/3.1, AKT1, and P300 increased (p < 0.05) (n = 15). The present study was preliminarily conducted to elucidate whether there might be a certain correlation between IKCa2.3 and LA hemodynamic changes. Conclusions: The TSS and wall pressure changes in the LA are correlated with the upregulation of mRNA and protein expression of IKCa2.3/3.1, AKT1, and P300.

Blood Flow Simulation to Determine the Risk of Thrombosis in the Fontan Circulation: Comparison between Atriopulmonary and Total Cavopulmonary Connections

Three-dimensional computational fluid dynamics (CFD) simulations were performed in the anastomotic region of the Fontan route between the venae cava and pulmonary arteries to investigate the risk of thrombosis due to blood stasis in the Fontan circulation. The finite volume method based on the time-averaged continuity and Navier–Stokes equations combined with the k-ω SST turbulent model was used in the CFD simulations. Low shear rate (SR) and SR on the wall (WSR) of <10 s−1 were used as markers to assess blood stasis as a cause of blood coagulation. Simulated blood flow velocity and both SR and WSR were reduced in the right atrium (RA) as the cavity of a flow channel in the atriopulmonary connection (APC) Fontan model, whereas the values increased in the total cavopulmonary connection (TCPC) Fontan model, which has no cavity. The volume of SR <10 s−1 and wall surface area of WSR <10 s−1 were, respectively, 4.6–261.8 cm3 and 1.2–38.3 cm2 in the APC Fontan model, and 0.1–0.3 cm3 and 0.1–0.6 cm2 in the TCPC Fontan model. The SR and WSR increased in the APC model with a normal-sized RA and the TCPC model as the flow rate of blood from the inferior vena cava increased with exercise; however, the SR and WSR in the RA decreased in the APC model with a dilated RA owing to the development of a recirculating flow. These findings suggest that the APC Fontan has a higher risk of thrombosis due to blood stasis than the TCPC Fontan and a higher RA dilation is associated with a higher risk of thrombosis from a fluid mechanics perspective.

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.

Computational fluid dynamics of blood flow in an idealized left human heart

We construct an idealized computational model of the left human heart for the study of the blood flow dynamics in the left atrium and ventricle. We solve the Navier-Stokes equations in the ALE formulation and we prescribe the left heart wall displacement based on physiological data; moreover, we consider the presence of both the mitral and aortic valves through the resistive method. We simulate the left heart hemodynamics by means of the finite element method and we consider the variational multiscale large eddy simulation (LES) formulation to account for the transitional and nearly turbulent regimes of the blood flow in physiological conditions. The main contribution of this paper is the characterization of the blood flow in an idealized configuration of the left heart aiming at reproducing function in normal conditions. Our assessment is based on the analysis of instantaneous and phase averaged velocity fields, blood pressure, and other clinically meaningful fluid dynamics indicators. Finally, we show that our idealized computational model can be suitably used to study and critically discuss pathological scenarios like that of a regurgitant mitral valve.

Sensitivity Analysis of In Silico Fluid Simulations to Predict Thrombus Formation after Left Atrial Appendage Occlusion

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.

Presence of Left Atrial Fibrosis May Contribute to Aberrant Hemodynamics and Increased Risk of Stroke in Atrial Fibrillation Patients

Atrial fibrillation (AF) patients are at high risk of stroke, with the left atrial appendage (LAA) found to be the most common site of clot formation. Presence of left atrial (LA) fibrosis has also been associated with higher stroke risk. However, the mechanisms for increased stroke risk in patients with atrial fibrotic remodeling are poorly understood. We sought to explore these mechanisms using fluid dynamic analysis and to test the hypothesis that the presence of LA fibrosis leads to aberrant hemodynamics in the LA, contributing to increased stroke risk in AF patients. We retrospectively collected late-gadolinium-enhanced MRI (LGE-MRI) images of eight AF patients (four persistent and four paroxysmal) and reconstructed their 3D LA surfaces. Personalized computational fluid dynamic simulations were performed, and hemodynamics at the LA wall were quantified by wall shear stress (WSS, friction of blood), oscillatory shear index (OSI, temporal directional change of WSS), endothelial cell activation potential (ECAP, ratio of OSI and WSS), and relative residence time (RRT, residence time of blood near the LA wall). For each case, these hemodynamic metrics were compared between fibrotic and non-fibrotic portions of the wall. Our results showed that WSS was lower, and OSI, ECAP, and RRT was higher in the fibrotic region as compared to the non-fibrotic region, with ECAP (p = 0.001) and RRT (p = 0.002) having significant differences. Case-wise analysis showed that these differences in hemodynamics were statistically significant for seven cases. Furthermore, patients with higher fibrotic burden were exposed to larger regions of high ECAP, which represents regions of low WSS and high OSI. Consistently, high ECAP in the vicinity of the fibrotic wall suggest that local blood flow was slow and oscillating that represents aberrant hemodynamic conditions, thus enabling prothrombotic conditions for circulating blood. AF patients with high LA fibrotic burden had more prothrombotic regions, providing more sites for potential clot formation, thus increasing their risk of stroke.

Subject-Specific Calculation of Left Atrial Appendage Blood-Borne Particle Residence Time Distribution in Atrial Fibrillation

Atrial fibrillation (AF) is the most common arrhythmia that leads to thrombus formation, mostly in the left atrial appendage (LAA). The current standard of stratifying stroke risk, based on the CHA2DS2-VASc score, does not consider LAA morphology, and the clinically accepted LAA morphology-based classification is highly subjective. The aim of this study was to determine whether LAA blood-borne particle residence time distribution and the proposed quantitative index of LAA 3D geometry can add independent information to the CHA2DS2-VASc score. Data were collected from 16 AF subjects. Subject-specific measurements included left atrial (LA) and LAA 3D geometry obtained by cardiac computed tomography, cardiac output, and heart rate. We quantified 3D LAA appearance in terms of a novel LAA appearance complexity index (LAA-ACI). We employed computational fluid dynamics analysis and a systems-based approach to quantify residence time distribution and associated calculated variable (LAA mean residence time, t m) in each subject. The LAA-ACI captured the subject-specific LAA 3D geometry in terms of a single number. LAA t m varied significantly within a given LAA morphology as defined by the current subjective method and it was not simply a reflection of LAA geometry/appearance. In addition, LAA-ACI and LAA t m varied significantly for a given CHA2DS2-VASc score, indicating that these two indices of stasis are not simply a reflection of the subjects' clinical status. We conclude that LAA-ACI and LAA t m add independent information to the CHA2DS2-VASc score about stasis risk and thereby can potentially enhance its ability to stratify stroke risk in AF patients.

A comprehensive comparison of various patient-specific CFD models of the left atrium for atrial fibrillation patients

BACKGROUND

Recently, advances in medical imaging, segmentation techniques, and high-performance computing have supported the use of patient-specific computational fluid dynamics (CFD) simulations. At present, CFD-compatible atrium geometries can be easily reconstructed from atrium images, providing important insight into the atrial fibrillation (AF) phenomenon, and assistance during therapy selection and surgical procedures. However, the hypothesis assumed for such CFD models should be adequately validated.

AIM

This work aims to perform an extensive study of the different hypotheses that are commonly assumed when performing atrial simulations for AF patients, as well as to evaluate and compare the range of indices that are usually applied to assess thrombus formation within the left atrium appendage (LAA).

METHODS

The atrial geometries of two AF patients have been segmented. The resulting geometries have been registered and interpolated to construct a dynamic mesh, which has been employed to compare the rigid and flexible models. Two families of hemodynamic indices have been calculated and compared: wall shear-based and blood age distribution-based.

RESULTS

The findings of this study illustrate the importance of validating the rigid atrium hypothesis when utilizing an AF CFD model. In particular, the absence of the A-wave contraction does not avoid a certain degree of passive atrial contraction, making the rigid model a poor approximation in some cases. Moreover, a new thrombosis predicting index has been proposed, i.e., M4, which has been shown to predict stasis more effectively than other indicators.

Fluid-structure interaction in a fully coupled three-dimensional mitral-atrium-pulmonary model

This paper aims to investigate detailed mechanical interactions between the pulmonary haemodynamics and left heart function in pathophysiological situations (e.g. atrial fibrillation and acute mitral regurgitation). This is achieved by developing a complex computational framework for a coupled pulmonary circulation, left atrium and mitral valve model. The left atrium and mitral valve are modelled with physiologically realistic three-dimensional geometries, fibre-reinforced hyperelastic materials and fluid–structure interaction, and the pulmonary vessels are modelled as one-dimensional network ended with structured trees, with specified vessel geometries and wall material properties. This new coupled model reveals some interesting results which could be of diagnostic values. For example, the wave propagation through the pulmonary vasculature can lead to different arrival times for the second systolic flow wave (S2 wave) among the pulmonary veins, forming vortex rings inside the left atrium. In the case of acute mitral regurgitation, the left atrium experiences an increased energy dissipation and pressure elevation. The pulmonary veins can experience increased wave intensities, reversal flow during systole and increased early-diastolic flow wave (D wave), which in turn causes an additional flow wave across the mitral valve (L wave), as well as a reversal flow at the left atrial appendage orifice. In the case of atrial fibrillation, we show that the loss of active contraction is associated with a slower flow inside the left atrial appendage and disappearances of the late-diastole atrial reversal wave (AR wave) and the first systolic wave (S1 wave) in pulmonary veins. The haemodynamic changes along the pulmonary vessel trees on different scales from microscopic vessels to the main pulmonary artery can all be captured in this model. The work promises a potential in quantifying disease progression and medical treatments of various pulmonary diseases such as the pulmonary hypertension due to a left heart dysfunction.

Effects of Ageing on Aortic Circulation During Atrial Fibrillation; a Numerical Study on Different Aortic Morphologies

Atrial fibrillation (AF) can alter intra-cardiac flow and cardiac output that subsequently affects aortic flow circulation. These changes may become more significant where they occur concomitantly with ageing. Aortic ageing is accompanied with morphological changes such as dilation, lengthening, and arch unfolding. While the recognition of AF mechanism has been the subject of numerous studies, less focus has been devoted to the aortic circulation during the AF and there is a lack of such investigation at different ages. The current work aims to address the present gap. First, we analyse aortic flow distribution in three configurations, which attribute to young, middle and old people, using geometries constructed via clinical data. We then introduce two transient inlet flow conditions representative of key AF-associated defects. Results demonstrate that both AF and ageing negatively affect flow circulation. The main consequence of concomitant occurrence is enhancement of endothelial cell activation potential (ECAP) throughout the vascular domain, mainly at aortic arch and descending thoracic aorta, which is consistent with some clinical observations. The outcome of the current study suggests that AF exacerbates the vascular defects occurred due to the ageing, which increases the possibility of cardiovascular diseases per se.

Demonstration of Patient-Specific Simulations to Assess Left Atrial Appendage Thrombogenesis Risk

Atrial fibrillation (AF) alters left atrial (LA) hemodynamics, which can lead to thrombosis in the left atrial appendage (LAA), systemic embolism and stroke. A personalized risk-stratification of AF patients for stroke would permit improved balancing of preventive anticoagulation therapies against bleeding risk. We investigated how LA anatomy and function impact LA and LAA hemodynamics, and explored whether patient-specific analysis by computational fluid dynamics (CFD) can predict the risk of LAA thrombosis. We analyzed 4D-CT acquisitions of LA wall motion with an in-house immersed-boundary CFD solver. We considered six patients with diverse atrial function, three with either a LAA thrombus (removed digitally before running the simulations) or a history of transient ischemic attacks (LAAT/TIA-pos), and three without a LAA thrombus or TIA (LAAT/TIA-neg). We found that blood inside the left atrial appendage of LAAT/TIA-pos patients had marked alterations in residence time and kinetic energy when compared with LAAT/TIA-neg patients. In addition, we showed how the LA conduit, reservoir and booster functions distinctly affect LA and LAA hemodynamics. Finally, fixed-wall and moving-wall simulations produced different LA hemodynamics and residence time predictions for each patient. Consequently, fixed-wall simulations risk-stratified our small cohort for LAA thrombosis worse than moving-wall simulations, particularly patients with intermediate LAA residence time. Overall, these results suggest that both wall kinetics and LAA morphology contribute to LAA blood stasis and thrombosis.

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

The left atrial appendage (LAA) is a complex cardiovascular structure which can yield to thrombi formation in patients with non-valvular atrial fibrillation (AF). The study of LAA fluid dynamics together with morphological features should be investigated in order to evaluate the possible connection of geometrical and hemodynamics indices with the stroke risk. To reach this goal, we conducted a morphological analysis of four different LAA shapes considering their variation during the cardiac cycle and computational fluid dynamics (CFD) simulations in AF conditions were carried out. The analysis of main geometrical LAA parameters showed a huger ostium and a reduced motility for the cauliflower and cactus shapes, as well as a lower velocity values from the CFD analysis. Such findings are in line with literature and highlight the importance of coupling dynamics imaging data with CFD calculations for providing information not available at clinical level.

Numerical Study of Atrial Fibrillation Effects on Flow Distribution in Aortic Circulation

Atrial fibrillation (AF) is the most common type of arrhythmia, which undermines cardiac function. Atrial fibrillation is a multi-facet malady and it may occur as a result of other diseases or it may trigger other problems. One of the main complications of AF is stroke due to the possibility of clot formation inside the atrium. However, the possibility of stroke occurrence due to the AF and the location from which an embolus dispatches are subject of debate. Another hypothesis about the embolus formation during AF is thrombus formation in aorta and carotid arteries, embolus detachment and its movement. To investigate the possibility of the latter postulation, the current work suggests a parametric study to quantify the sensitivity of aortic flow to four common AF traits including lack of atrial kick, atrial remodelling, left ventricle systolic dysfunction, and high frequency fibrillation. The simulation was carried out by coupling several in-house codes and ANSYS-CFX module. The results reveal that AF traits lower flow rate at left ventricular outflow tract, which in general lowers blood perfusion to systemic, cerebral and coronary circulations. Consequently, it leads to endothelial cell activation potential (ECAP) increase and variation of flow structure that both suggest predisposed areas to atherogenesis and thrombus formation in different regions in ascending aorta, aortic arch and descending thoracic aorta.

Analysis of a coupled fluid-structure interaction model of the left atrium and mitral valve

We present a coupled left atrium-mitral valve model based on computed tomography scans with fibre-reinforced hyperelastic materials. Fluid-structure interaction is realised by using an immersed boundary-finite element framework. Effects of pathological conditions, eg, mitral valve regurgitation and atrial fibrillation, and geometric and structural variations, namely, uniform vs non-uniform atrial wall thickness and rule-based vs atlas-based fibre architectures, on the system are investigated. We show that in the case of atrial fibrillation, pulmonary venous flow reversal at late diastole disappears, and the filling waves at the left atrial appendage orifice during systole have reduced magnitude. In the case of mitral regurgitation, a higher atrial pressure and disturbed flows are seen, especially during systole, when a large regurgitant jet can be found with the suppressed pulmonary venous flow. We also show that both the rule-based and atlas-based fibre defining methods lead to similar flow fields and atrial wall deformations. However, the changes in wall thickness from non-uniform to uniform tend to underestimate the atrial deformation. Using a uniform but thickened wall also lowers the overall strain level. The flow velocity within the left atrial appendage, which is important in terms of appendage thrombosis, increases with the thickness of the left atrial wall. Energy analysis shows that the kinetic and dissipation energies of the flow within the left atrium are altered differently by atrial fibrillation and mitral valve regurgitation, providing a useful indication of the atrial performance in pathological situations.

A Proof of Concept for Computational Fluid Dynamic Analysis of the Left Atrium in Atrial Fibrillation on a Patient-Specific Basis

Atrial fibrillation (AF) is associated with a fivefold increase in the risk of cerebrovascular events, being responsible of 15–18% of all strokes. The morphological and functional remodeling of the left atrium (LA) caused by AF favors blood stasis and, consequently, stroke risk. In this context, several clinical studies suggest that the stroke risk stratification could be improved by using hemodynamic information on the LA and the left atrial appendage (LAA). The goal of this study was to develop a personalized computational fluid dynamics (CFD) model of the LA which could clarify the hemodynamic implications of AF on a patient-specific basis. In this paper, we present the developed model and its application to two AF patients as a preliminary advancement toward an optimized stroke risk stratification pipeline.

Cloud Database Construction for the Expressway Design by the use of the Medical Information

Your heart rate and blood pressure are respond to the curve, slope, lane width, and road surface friction coefficient of the expressway design. However, no report was reported concerning about the Expressway design from the viewpoint of medical information of the driver until now. To prevent the traffic accident, human factor is of course one of the most important factors. In this study, the Cloud Database Construction for the Expressway Design by the use of the Medical Information had been tried to carry out. HR response and PWV responses had been tried to be analyzed by the sensors in the car during driving. LF, HF and LF/HF of Heart rate variability had been calculated and tagged with expressway information including left and right curve, slope, lane width, and road surface friction coefficient. Furthermore, pulse of the descending aorta had been tried to be recorded from the sensor in a driver seat, so, the pulse wave velocity and blood pressure could be evaluated. Recording system of an Eye movement, pupil diameter, cerebral blood flow, and EEG are now under construction. So, all human driver's data will be combined in the Cloud of the Central office. this method will be useful for the development of the designing method the Expressway in near future.

In silico Optimization of Left Atrial Appendage Occluder Implantation Using Interactive and Modeling Tools

According to clinical studies, around one third of patients with atrial fibrillation (AF) will suffer a stroke during their lifetime. Between 70 and 90% of these strokes are caused by thrombus formed in the left atrial appendage. In patients with contraindications to oral anticoagulants, a left atrial appendage occluder (LAAO) is often implanted to prevent blood flow entering in the LAA. A limited range of LAAO devices is available, with different designs and sizes. Together with the heterogeneity of LAA morphology, these factors make LAAO success dependent on clinician's experience. A sub-optimal LAAO implantation can generate thrombi outside the device, eventually leading to stroke if not treated. The aim of this study was to develop clinician-friendly tools based on biophysical models to optimize LAAO device therapies. A web-based 3D interactive virtual implantation platform, so-called VIDAA, was created to select the most appropriate LAAO configurations (type of device, size, landing zone) for a given patient-specific LAA morphology. An initial LAAO configuration is proposed in VIDAA, automatically computed from LAA shape features (centreline, diameters). The most promising LAAO settings and LAA geometries were exported from VIDAA to build volumetric meshes and run Computational Fluid Dynamics (CFD) simulations to assess blood flow patterns after implantation. Risk of thrombus formation was estimated from the simulated hemodynamics with an index combining information from blood flow velocity and complexity. The combination of the VIDAA platform with in silico indices allowed to identify the LAAO configurations associated to a lower risk of thrombus formation; device positioning was key to the creation of regions with turbulent flows after implantation. Our results demonstrate the potential for optimizing LAAO therapy settings during pre-implant planning based on modeling tools and contribute to reduce the risk of thrombus formation after treatment.

The Impact of Left Atrium Appendage Morphology on Stroke Risk Assessment in Atrial Fibrillation: A Computational Fluid Dynamics Study

Atrial fibrillation (AF) carries out a 5-fold increase in stroke risk, related to embolization of thrombi clotting in left atrium (LA). Left atrial appendage (LAA) is the site with the highest blood stasis which causes thrombus formation. About 90 % of the intracardiac thrombi in patients with cardioembolic events originally develop in the LAA. Recent studies have been focused on the association between LAA anatomical features and stroke risk and provided conflicting results. Haemodynamic and fluid dynamic information on the LA and mostly on the LAA may improve stroke risk stratification. Therefore, the aim of this study was the design and development of a workflow to quantitatively define the influence of the LAA morphology on LA hemodynamics. Five 3D LA anatomical models, obtained from real clinical data, which were clearly different as regard to LAA morphology were used. For each LAA we identified and computed several parameters describing its geometry. Then, one LA chamber model was chosen and a framework was developed to connect the different LAAs belonging to the other four patients to this model. These new anatomical models represented the computational domain for the computational fluid dynamics (CFD) study; simulations of the hemodynamics within the LA and LAA were performed in order to evaluate the interplay of the LAA shape on the blood flow characteristics in AF condition. CFD simulations were carried out for five cardiac cycles. Blood velocity, vorticity, LAA orifice velocity, residence time computed in the five models were compared and correlated with LAA morphologies. Results showed that not only complex morphologies were characterized by low velocities, low vorticity and consequently could carry a higher thrombogenic risk; even qualitatively simple morphologies showed a thrombogenic risk equal, or even higher, than more complex auricles. CFD results supported the hypothesis that LAA geometric characteristics plays a key-role in defining thromboembolic risk. LAA geometric parameters could be considered, coupled with the morphological characteristics, for a comprehensive evaluation of the regional blood stasis. The proposed procedure might address the development of a tool for patient-specific stroke risk assessment and preventive treatment in AF patients, relying on morpho-functional defintion of each LAA type.

Modeling Left Atrial Flow, Energy, Blood Heating Distribution in Response to Catheter Ablation Therapy

Introduction: Atrial fibrillation (AF) is a widespread cardiac arrhythmia that commonly affects the left atrium (LA), causing it to quiver instead of contracting effectively. This behavior is triggered by abnormal electrical impulses at a specific site in the atrial wall. Catheter ablation (CA) treatment consists of isolating this driver site by burning the surrounding tissue to restore sinus rhythm (SR). However, evidence suggests that CA can concur to the formation of blood clots by promoting coagulation near the heat source and in regions with low flow velocity and blood stagnation.

Methods: A patient-specific modeling workflow was created and applied to simulate thermal-fluid dynamics in two patients pre- and post-CA. Each model was personalized based on pre- and post-CA imaging datasets. The wall motion and anatomy were derived from SSFP Cine MRI data, while the trans-valvular flow was based on Doppler ultrasound data. The temperature distribution in the blood was modeled using a modified Pennes bioheat equation implemented in a finite-element based Navier-Stokes solver. Blood particles were also classified based on their residence time in the LA using a particle-tracking algorithm.

Results: SR simulations showed multiple short-lived vortices with an average blood velocity of 0.2-0.22 m/s. In contrast, AF patients presented a slower vortex and stagnant flow in the LA appendage, with the average blood velocity reduced to 0.08–0.14 m/s. Restoration of SR also increased the blood kinetic energy and the viscous dissipation due to the presence of multiple vortices. Particle tracking showed a dramatic decrease in the percentage of blood remaining in the LA for longer than one cycle after CA (65.9 vs. 43.3% in patient A and 62.2 vs. 54.8% in patient B). Maximum temperatures of 76° and 58°C were observed when CA was performed near the appendage and in a pulmonary vein, respectively.

Conclusion: This computational study presents novel models to elucidate relations between catheter temperature, patient-specific atrial anatomy and blood velocity, and predict how they change from SR to AF. The models can quantify blood flow in critical regions, including residence times and temperature distribution for different catheter positions, providing a basis for quantifying stroke risks.

Sensitivity analysis of geometrical parameters to study haemodynamics and thrombus formation in the left atrial appendage

The left atrial appendage (LAA) is a complex and heterogeneous protruding structure of the left atrium (LA). In atrial fibrillation patients, it is the location where 90% of the thrombi are formed. However, the role of the LAA in thrombus formation is not fully known yet. The main goal of this work is to perform a sensitivity analysis to identify the most relevant LA and LAA morphological parameters in atrial blood flow dynamics. Simulations were run on synthetic ellipsoidal left atria models where different parameters were individually studied: pulmonary veins and mitral valve dimensions; LAA shape; and LA volume. Our computational analysis confirmed the relation between large LAA ostia, low blood flow velocities and thrombus formation. Additionally, we found that pulmonary vein configuration exerted a critical influence on LAA blood flow patterns. These findings contribute to a better understanding of the LAA and to support clinical decisions for atrial fibrillation patients.

Computational Fluid Dynamic Analysis of the Left Atrial Appendage to Predict Thrombosis Risk

During Atrial Fibrillation (AF) more than 90% of the left atrial thrombi responsible for thromboembolic events originate in the left atrial appendage (LAA), a complex small sac protruding from the left atrium (LA). Current available treatments to prevent thromboembolic events are oral anticoagulation, surgical LAA exclusion, or percutaneous LAA occlusion. However, the mechanism behind thrombus formation in the LAA is poorly understood. The aim of this work is to analyse the hemodynamic behaviour in four typical LAA morphologies - "Chicken wing", "Cactus", "Windsock" and "Cauliflower" - to identify potential relationships between the different shapes and the risk of thrombotic events. Computerised tomography (CT) images from four patients with no LA pathology were segmented to derive the 3D anatomical shape of LAA and LA. Computational Fluid Dynamic (CFD) analyses based on the patient-specific anatomies were carried out imposing both healthy and AF flow conditions. Velocity and shear strain rate (SSR) were analysed for all cases. Residence time in the different LAA regions was estimated with a virtual contrast agent washing out. CFD results indicate that both velocity and SSR decrease along the LAA, from the ostium to the tip, at each instant in the cardiac cycle, thus making the LAA tip more prone to fluid stagnation, and therefore to thrombus formation. Velocity and SSR also decrease from normal to AF conditions. After four cardiac cycles, the lowest washout of contrast agent was observed for the Cauliflower morphology (3.27% of residual contrast in AF), and the highest for the Windsock (0.56% of residual contrast in AF). This suggests that the former is expected to be associated with a higher risk of thrombosis, in agreement with clinical reports in the literature. The presented computational models highlight the major role played by the LAA morphology on the hemodynamics, both in normal and AF conditions, revealing the potential support that numerical analyses can provide in the stratification of patients under risk of thrombus formation, towards personalised patient care.

Computational Flow Dynamic Analysis of Right and Left Atria in Patent Foramen Ovale: Potential Links with Atrial Fibrillation

Background

An impairment of the left atrial function similar to that usually observed in atrial fibrillation (AF) has been observed also in patients with patent foramen ovale (PFO) and permanent right-to-left shunting (RLS).

Methods

We reconstructed the geometrical model of right atrium (RA), PFO, left atrium (LA) and left atrial appendage (LAA) of 65 patients with mild (36 patients mean age 45.5±6.8 years, 24 females) or permanent (29 patients, mean age 45.1±5.3 years, 21 females) RLS using anatomical data obtained by transoesophageal echocardiography (TEE) and cardiac MRI, performed as a part of our institutional screening protocol for paradoxical embolism. Using computational fluid dynamic analysis (CFD) we assessed the vorticity magnitude in both the LA and LAA to analyse a possible rheological relationship between PFO and AF.

Results

The anatomical models, in terms of dimensions, were comparable among the patients with mild and permanent RLS. A higher vorticity magnitude was observed in the mild shunt both in the LA (101.12±21.3 vs 88.3±22.6, p=0.02) and LAA (62±14.4 vs 32.4±12.3, p<0.01) when compared to the permanent R-L shunting.

Conclusion

The lower vorticity magnitude across the LA and LAA in patients with permanent RLS suggests a possible higher stagnation of the blood in these anatomical sites, similarly as previously observed in patients with AF.

Experimental and CFD flow studies in an intracranial aneurysm model with Newtonian and non-Newtonian fluids

BACKGROUND

According to the clinical data, flow conditions play a major role in the genesis of intracranial aneurysms. The disorder of the flow structure is the cause of damage of the inner layer of the vessel wall, which leads to the development of cerebral aneurysms. Knowledge of the alteration of the flow field in the aneurysm region is important for treatment.

OBJECTIVE

The aim is to study quantitatively the flow structure in an patient-specific aneurysm model of the internal carotid artery using both experimental and computational fluid dynamics (CFD) methods with Newtonian and non-Newtonian fluids.

METHODS

A patient-specific geometry of aneurysm of the internal carotid artery was used. Patient data was segmented and smoothed to obtain geometrical model. An elastic true-to-scale silicone model was created with stereolithography. For initial investigation of the blood flow, the flow was visualized by adding particles into the silicone model. The precise flow velocity measurements were done using 1D Laser Doppler Anemometer with a spatial resolution of 50 μ m and a temporal resolution of 1 ms. The local velocity measurements were done at a distance of 4 mm to each other. A fluid with non-Newtonian properties was used in the experiment. The CFD simulations for unsteady-state problem were done using constructed hexahedral mesh for Newtonian and non-Newtonian fluids.

RESULTS

Using 1D laser Doppler Anemometer the minimum velocity magnitude at the end of systole -0.01 m/s was obtained in the aneurysm dome while the maximum velocity 1 m/s was at the center of the outlet segment. On central cross section of the aneurysm the maximum velocity value is only 20% of the average inlet velocity. The average velocity on the cross-section is only 11% of the inlet axial velocity. Using the CFD simulation the wall shear stresses for Newtonian and non-Newtonian fluid at the end of systolic phase (t= 0.25 s) were computed. The wall shear stress varies from 3.52 mPa (minimum value) to 10.21 Pa (maximum value) for the Newtonian fluid. For the non-Newtonian fluid the wall shear stress minimum is 2.94 mPa; the maximum is 9.14 Pa. The lowest value of the wall shear stress for both fluids was obtained at the dome of the aneurysm while the highest wall shear stress was at the beginning of the outlet segment. The vortex in the aneurysm region is unstable during the cardiac cycle. The clockwise rotation of the streamlines at the inlet segment for Newtonian and non-Newtonian fluid is shown.

CONCLUSION

The results of the present study are in agreement with the hemodynamics theory of aneurysm genesis. Low value of wall shear stress is observed at the aneurysm dome which can cause a rupture of an aneurysm.