Investigation of Atrial Vortices Using a Novel Right Heart Model and Possible Implications for Atrial Thrombus Formation

Investigating the importance of left atrial compliance on fluid dynamics in a novel mock circulatory loop

The left atrium (LA) hemodynamic indices hold prognostic value in various cardiac diseases and disorders. To understand the mechanisms of these conditions and to assess the performance of cardiac devices and interventions, in vitro models can be used to replicate the complex physiological interplay between the pulmonary veins, LA, and left ventricle. In this study, a comprehensive and adaptable in vitro model was created. The model includes a flexible LA made from silicone and allows distinct control over the systolic and diastolic functions of both the LA and left ventricle. The LA was mechanically matched with porcine LAs through expansion tests. Fluid dynamic measures were validated against the literature and pulmonary venous flows recorded on five healthy individuals using magnetic resonance flow imaging. Furthermore, the fluid dynamic measures were also used to construct LA pressure-volume loops. The in vitro pressure and flow recordings expressed a high resemblance to physiological waveforms. By decreasing the compliance of the LA, the model behaved realistically, elevating the a- and v-wave peaks of the LA pressure from 12 to 19 mmHg and 22 to 26 mmHg, respectively, while reducing the S/D ratio of the pulmonary venous flowrate from 1.5 to 0.3. This model provides a realistic platform and framework for developing and evaluating left heart procedures and interventions.

Reduced Left Atrial Rotational Flow Is Independently Associated With Embolic Brain Infarcts

BACKGROUND

Up to 25% of embolic strokes occur in individuals without atrial fibrillation (AF) or other identifiable mechanisms.

OBJECTIVES

This study aims to assess whether left atrial (LA) blood flow characteristics are associated with embolic brain infarcts, independently of AF.

METHODS

The authors recruited 134 patients: 44 with a history of ischemic stroke and 90 with no history of stroke but CHA2DS2VASc score ≥1. Cardiac magnetic resonance (CMR) evaluated cardiac function and LA 4-dimensional flow parameters, including velocity and vorticity (a measure of rotational flow), and brain magnetic resonance imaging (MRI) was performed to detect large noncortical or cortical infarcts (LNCCIs) (likely embolic), or nonembolic lacunar infarcts.

RESULTS

Patients (41% female; age 70 ± 9 years) had moderate stroke risk (median CHA2DS2VASc = 3, Q1-Q3: 2-4). Sixty-eight (51%) had diagnosed AF, of whom 58 (43%) were in AF during CMR. Thirty-nine (29%) had ≥1 LNCCI, 20 (15%) had ≥1 lacunar infarct without LNCCI, and 75 (56%) had no infarct. Lower LA vorticity was significantly associated with prevalent LNCCIs after adjustment for AF during CMR, history of AF, CHA2DS2VASc score, LA emptying fraction, LA indexed maximum volume, left ventricular ejection fraction, and indexed left ventricular mass (OR: 2.06 [95% CI: 1.08-3.92 per SD]; P = 0.027). By contrast, LA flow peak velocity was not significantly associated with LNCCIs (P = 0.21). No LA parameter was associated with lacunar infarcts (all P > 0.05).

CONCLUSIONS

Reduced LA flow vorticity is significantly and independently associated with embolic brain infarcts. Imaging LA flow characteristics may aid identification of individuals who would benefit from anticoagulation for embolic stroke prevention, regardless of heart rhythm.

On non-Kolmogorov turbulence in blood flow and its possible role in mechanobiological stimulation

The study of turbulence in physiologic blood flow is important due to its strong relevance to endothelial mechanobiology and vascular disease. Recently, Saqr et al. (Sci Rep 10, 15,492, 2020) discovered non-Kolmogorov turbulence in physiologic blood flow in vivo, traced its origins to the Navier–Stokes equation and demonstrated some of its properties using chaos and hydrodynamic-stability theories. The present work extends these findings and investigates some inherent characteristics of non-Kolmogorov turbulence in monoharmonic and multiharmonic pulsatile flow under ideal physiologic conditions. The purpose of this work is to propose a conjecture for the origins for picoNewton forces that are known to regulate endothelial cells’ functions. The new conjecture relates these forces to physiologic momentum-viscous interactions in the near-wall region of the flow. Here, we used high-resolution large eddy simulation (HRLES) to study pulsatile incompressible flow in a straight pipe of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L/D=20$$\end{document}L/D=20. The simulations presented Newtonian and Carreau–Yasuda fluid flows, at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R{e}_{m}\approx 250$$\end{document}Rem≈250, each represented by one, two and three boundary harmonics. Comparison was established based on maintaining constant time-averaged mass flow rate in all simulations. First, we report the effect of primary harmonics on the global power budget using primitive variables in phase space. Second, we describe the non-Kolmogorov turbulence in frequency domain. Third, we investigate the near-wall coherent structures in time and space domains. Finally, we propose a new conjecture for the role of turbulence in endothelial cells’ mechanobiology. The proposed conjecture correlates near-wall turbulence to a force field of picoNewton scale, suggesting possible relevance to endothelial cells mechanobiology.

Impact of Altered Vena Cava Flow Rates on Right Atrium Flow Characteristics

The right atrium (RA) combines the superior (SVC) and inferior vena cava (IVC) flows. Treatments like extracorporeal membrane oxygenation (ECMO) and hemodialysis by catheter alter IVC/SVC flows. Here we assess how altered IVC/SVC flow contributions impact RA flow. Four healthy volunteers were imaged with CT, reconstructed and combined into a patient-averaged model. Large Eddy Simulations (LES) were performed for a range of IVC/SVC flow contributions (30-70% each, increments of 5%) and common flow metrics were recorded. Model sensitivity to reconstruction domain extent, constant/pulsatile inlets and hematocrit was also assessed. Consistent with literature, a single vortex occupied the central RA across all flowrates with a smaller counter-rotating vortex, not previously reported, in the auricle. Vena cava flow was highly helical. RA turbulent kinetic energy (TKE) (P=0.027) and time-averaged wall shear stress (WSS) (P<0.001) increased with SVC flow. WSS was lower in the auricle (2 Pa, P<0.001). WSS in the vena cava were equal at IVC/SVC =65/35%. The model was highly sensitive to the reconstruction domain with cropped geometries lacking helicity in the vena cavae, altering RA flow. RA flow was not significantly affected by constant inlets or hematocrit. The rotational flow conventionally described in the RA is confirmed however a new, smaller vortex was also recorded in the auricle. When IVC flow dominates, as is normal, TKE in the RA is reduced and WSS in the vena cavae equalize. Significant helicity exists in the vena cava, a result of distal geometry and this geometry appears crucial to accurately simulating RA flow.

The impact of atrial fibrillation and stroke risk factors on left atrial blood flow characteristics

AIMS

Altered left atrial (LA) blood flow characteristics account for an increase in cardioembolic stroke risk in atrial fibrillation (AF). Here, we aimed to assess whether exposure to stroke risk factors is sufficient to alter LA blood flow even in the presence of sinus rhythm (SR).

METHODS AND RESULTS

We investigated 95 individuals: 37 patients with persistent AF, who were studied before and after cardioversion [Group 1; median CHA2DS2-VASc = 2.0 (1.5-3.5)]; 35 individuals with no history of AF but similar stroke risk to Group 1 [Group 2; median CHA2DS2-VASc = 3.0 (2.0-4.0)]; and 23 low-risk individuals in SR [Group 3; median CHA2DS2-VASc = 0.0 (0.0-0.0)]. Cardiac function and LA flow characteristics were evaluated using cardiac magnetic resonance. Before cardioversion, Group 1 displayed impaired left ventricular (LV) and LA function, reduced LA flow velocities and vorticity, and a higher normalized vortex volume (all P < 0.001 vs. Groups 2 and 3). After restoration of SR at ≥4-week post-cardioversion, LV systolic function and LA flow parameters improved significantly (all P < 0.001 vs. pre-cardioversion) and were no longer different from those in Group 2. However, in the presence of SR, LA flow peak and mean velocity, and vorticity were lower in Groups 1 and 2 vs. Group 3 (all P < 0.01), and were associated with impaired LA emptying fraction (LAEF) and LV diastolic dysfunction.

CONCLUSION

Patients at moderate-to-high stroke risk display altered LA flow characteristics in SR in association with an LA myopathic phenotype and LV diastolic dysfunction, regardless of a history of AF.

Left atrial 4D flow cardiovascular magnetic resonance: a reproducibility study in sinus rhythm and atrial fibrillation

BACKGROUND

Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) allows sophisticated quantification of left atrial (LA) blood flow, and could yield novel biomarkers of propensity for intra-cardiac thrombus formation and embolic stroke. As reproducibility is critically important to diagnostic performance, we systematically investigated technical and temporal variation of LA 4D flow in atrial fibrillation (AF) and sinus rhythm (SR).

METHODS

Eighty-six subjects (SR, n = 64; AF, n = 22) with wide-ranging stroke risk (CHA2DS2VASc 0-6) underwent LA 4D flow assessment of peak and mean velocity, vorticity, vortex volume, and stasis. Eighty-five (99%) underwent a second acquisition within the same session, and 74 (86%) also returned at 30 (27-35) days for an interval scan. We assessed variability attributable to manual contouring (intra- and inter-observer), and subject repositioning and reacquisition of data, both within the same session (same-day scan-rescan), and over time (interval scan). Within-subject coefficients of variation (CV) and bootstrapped 95% CIs were calculated and compared.

RESULTS

Same-day scan-rescan CVs were 6% for peak velocity, 5% for mean velocity, 7% for vorticity, 9% for vortex volume, and 10% for stasis, and were similar between SR and AF subjects (all p > 0.05). Interval-scan variability was similar to same-day scan-rescan variability for peak velocity, vorticity, and vortex volume (all p > 0.05), and higher for stasis and mean velocity (interval scan CVs of 14% and 8%, respectively, both p < 0.05). Longitudinal changes in heart rate and blood pressure at the interval scan in the same subjects were associated with significantly higher variability for LA stasis (p = 0.024), but not for the remaining flow parameters (all p > 0.05). SR subjects showed significantly greater interval-scan variability than AF patients for mean velocity, vortex volume, and stasis (all p < 0.05), but not peak velocity or vorticity (both p > 0.05).

CONCLUSIONS

LA peak velocity and vorticity are the most reproducible and temporally stable novel LA 4D flow biomarkers, and are robust to changes in heart rate, blood pressure, and differences in heart rhythm.

Computational fluid dynamics of the right atrium: Assessment of modelling criteria for the evaluation of dialysis catheters

Central venous catheters are widely used in haemodialysis therapy, having to respect design requirements for appropriate performance. These are placed within the right atrium (RA); however, there is no prior computational study assessing different catheter designs while mimicking their native environment. Here, a computational fluid dynamics model of the RA, based on realistic geometry and transient physiological boundary conditions, was developed and validated. Symmetric, split and step catheter designs were virtually placed in the RA and their performance was evaluated by: assessing their interaction with the RA haemodynamic environment through prediction of flow vorticity and wall shear stress (WSS) magnitudes (1); and quantifying recirculation and tip shear stress (2). Haemodynamic predictions from our RA model showed good agreement with the literature. Catheter placement in the RA increased average vorticity, which could indicate alterations of normal blood flow, and altered WSS magnitudes and distribution, which could indicate changes in tissue mechanical properties. All designs had recirculation and elevated shear stress values, which can induce platelet activation and subsequently thrombosis. The symmetric design, however, had the lowest associated values (best performance), while step design catheters working in reverse mode were associated with worsened performance. Different tip placements also impacted on catheter performance. Our findings suggest that using a realistically anatomical RA model to study catheter performance and interaction with the haemodynamic environment is crucial, and that care needs to be given to correct tip placement within the RA for improved recirculation percentages and diminished shear stress values.

Histological properties of oscillating intracardiac masses associated with cardiac implantable electric devices

Background

There have been a few cases of echogenic cardiac implantable electric device (CIED) lead‐associated oscillating intracardiac masses (ICMs) in leads imaged by echocardiography. The histological properties of ICMs could help clarify the etiological diagnosis. Although there is extensive literature on mass size, the histological properties of such masses have not been characterized. The aim of this research was to clarify the histological features of oscillating ICMs in CIED patients.

Methods

Preoperative echocardiography was performed in all candidates for CIED removal. In the patients with ICMs, specimens were obtained by 3 methods: direct tissue collection during open‐heart surgery; tissue collection together with the CIED lead during transvenous extraction; and tissue collection by catheter vacuum during transvenous CIED removal. A standard histopathological examination of ICM tissue was performed.

Results

A total of 106 patients underwent lead removal in our institute (April 2009‐March 2018); 14 patients had an ICM (13.2%), and 7 specimens were obtained in patients with CIED lead‐related ICM. Following histological examination, 2 types of ICM were identified: one mainly composed of thickened endocardium (EN type; 3 patients), and the other mainly an aggregate of inflammatory cells as a neutrophil cell (NC type; 4 patients).

Conclusions

Two histological types of intracardiac masses, including a thickened endocardium type and a neutrophil cell type, were identified. These classifications might help make an accurate histological diagnosis of lead‐associated intracardiac masses.

A Novel Approach for 3D-Structural Identification through Video Recording: Magnified Tracking

Advancements in optical imaging devices and computer vision algorithms allow the exploration of novel diagnostic techniques for use within engineering systems. A recent field of application lies in the adoption of such devices for non-contact vibrational response recordings of structures, allowing high spatial density measurements without the burden of heavy cabling associated with conventional technologies. This, however, is not a straightforward task due to the typically low-amplitude displacement response of structures under ambient operational conditions. A novel framework, namely Magnified Tracking (MT), is proposed herein to overcome this limitation through the synergistic use of two computer vision techniques. The recently proposed phase-based motion magnification (PBMM) framework, for amplifying motion in a video within a defined frequency band, is coupled with motion tracking by means of particle tracking velocimetry (PTV). An experimental campaign was conducted to validate a proof-of-concept, where the dynamic response of a shear frame was measured both by conventional sensors as well as a video camera setup, and cross-compared to prove the feasibility of the proposed non-contact approach. The methodology was explored both in 2D and 3D configurations, with PTV revealing a powerful tool for the measurement of perceptible motion. When MT is utilized for tracking “imperceptible” structural responses (i.e., below PTV sensitivity), via the use of PBMM around the resonant frequencies of the structure, the amplified motion reveals the operational deflection shapes, which are otherwise intractable. The modal results extracted from the magnified videos, using PTV, demonstrate MT to be a viable non-contact alternative for 3D modal identification with the benefit of a spatially dense measurement grid.

Hemodynamic Changes in the Right Ventricle Induced by Variations of Cardiac Output: A Possible Mechanism for Arrhythmia Occurrence in the Outflow Tract

The rationale of this paper is to investigate right ventricular (RV) hemodynamics in relation to changes in cardiac output, and in particular to study exercise-induced stresses at the RV outflow tract (RVOT), which is a common site of ventricular arrhythmias in the athlete’s heart. We hypothesize that the thin-walled RVOT is exposed to high wall shear stresses (WSS) during physiological states associated with high cardiac output such as exercise, and therefore, may be particularly prone to substrate formation leading to ventricular tachyarrhythmias. 3D Particle Tracking Velocimetry (3D-PTV), an optical imaging method, has been performed in a novel anatomically accurate compliant silicone right heart model derived from a high resolution MRI heart scan of a healthy male proband. RV and RVOT flow patterns at resting conditions were obtained from two healthy athletic male proband’s hearts and two patients with arrhythmogenic right ventricular cardiomyopathy/dysplasia (ARVC/D) via phase contrast magnetic resonance imaging (PC-MRI). The healthy case was used as a reference for validating the in vitro flow patterns of the silicone model, while the diseased cases were used to generalize our findings and investigate possible changes in hemodynamic stresses with RV morphological remodelling. Our results showed that both healthy and diseased geometries consistently displayed an increased WSS in the RVOT relative to the rest of the RV. We found that increases in cardiac output may lead to increases of mean kinetic energy (MKE), laminar viscous dissipation and WSS at the RVOT. Furthermore, higher peak WSS magnitudes were found for the diseased cases. The identified high WSS regions may correlate with the common site of RVOT ventricular tachycardia in athletes and patients with ARVC/D. Our results imply that exercise, as well as anatomical and functional remodeling might alter RV wall shear stress both in magnitude and spatial distribution, leading to increased hemodynamic stresses in the RVOT.

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.