Determining possible thrombus sites in an extracorporeal device, using computational fluid dynamics-derived relative residence time

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

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

Longitudinal wall shear stress evaluation using centerline projection approach in the numerical simulations of the patient-based carotid artery

In this numerical study, areas of the carotid bifurcation and of a distal stenosis in the internal carotid artery are closely observed to evaluate the patient's current risks of ischemic stroke. An indicator for the vessel wall defects is the stress exerted by blood on the vessel tissue, typically expressed by the amplitude of the wall shear stress vector (WSS) and its oscillatory shear index. To detect negative shear stresses corresponding with reversal flow, we perform orientation-based shear evaluation. We investigate the longitudinal component of the wall shear vector, where tangential vectors aligned longitudinally with the vessel are necessary. However, resulting from imaging segmentation resolution of patients' computed tomography angiography scans and stenotic regions, the geometry model's mesh is non-smooth on its surface areas and the automatically generated tangential vector field is discontinuous and multi-directional, making an interpretation of our orientation-based risk indicators unreliable. We improve the evaluation of longitudinal shear stress by applying the projection of the vessel's centerline to the surface to construct smooth tangential field aligned longitudinally with the vessel. We validate our approach for the longitudinal WSS component and the corresponding oscillatory index by comparing them to results obtained using automatically generated tangents in both rigid and elastic vessel modeling and to amplitude-based indicators. We present the major benefit of our longitudinal WSS evaluation based on its directionality for the cardiovascular risk assessment, which is the detection of negative WSS indicating persistent reversal or transverse flow. This is impossible in the case of the amplitude-based WSS.

A computational model-based study on the feasibility of predicting post-splenectomy thrombosis using hemodynamic metrics

For portal hypertensive patients with splenomegaly and hypersplenism, splenectomy is an effective surgery to relieve the complications. However, patients who have undergone splenectomy often suffer from portal venous system thrombosis, a sequela that requires prophylaxis and timely treatment to avoid deterioration and death. The aim of this study is to investigate the feasibility of predicting post-splenectomy thrombosis using hemodynamic metrics based on computational models. First, 15 portal hypertensive patients who had undergone splenectomy were enrolled, and their preoperative clinical data and postoperative follow-up results were collected. Next, computational models of the portal venous system were constructed based on the preoperative computed tomography angiography images and ultrasound-measured flow velocities. On this basis, splenectomy was mimicked and the postoperative area of low wall shear stress (ALWSS) was simulated for each patient-specific model. Finally, model-simulated ALWSS was statistically compared with the patient follow-up results to investigate the feasibility of predicting post-splenectomy thrombosis using hemodynamic metrics. Results showed that ALWSS could predict the occurrence of post-splenectomy thrombosis with the area under the receiver operating characteristic curve (AUC) equal to 0.75. Moreover, statistical analysis implied that the diameter of the splenic vein is positively correlated with ALWSS (r = 0.883, p < 0.0001), and the anatomical structures of the portal venous system also influence the ALWSS. These findings demonstrated that the computational model-based hemodynamic metric ALWSS, which is associated with the anatomorphological features of the portal venous system, is capable of predicting the occurrence of post-splenectomy thrombosis, promoting better prophylaxis and postoperative management for portal hypertensive patients receiving splenectomy.

In-silico enhanced animal study of pulmonary artery pressure sensors: assessing hemodynamics using computational fluid dynamics

To assess whether in-silico models can be used to predict the risk of thrombus formation in pulmonary artery pressure sensors (PAPS), a chronic animal study using pigs was conducted. Computed tomography (CT) data was acquired before and immediately after implantation, as well as one and three months after the implantation. Devices were implanted into 10 pigs, each one in the left and right pulmonary artery (PA), to reduce the required number of animal experiments. The implantation procedure aimed at facilitating optimal and non-optimal positioning of the devices to increase chances of thrombus formation. Eight devices were positioned non-optimally. Three devices were positioned in the main PA instead of the left and right PA. Pre-interventional PA geometries were reconstructed from the respective CT images, and the devices were virtually implanted at the exact sites and orientations indicated by the follow-up CT after one month. Transient intra-arterial hemodynamics were calculated using computational fluid dynamics. Volume flow rates were modelled specifically matching the animals body weights. Wall shear stresses (WSS) and oscillatory shear indices (OSI) before and after device implantation were compared. Simulations revealed no relevant changes in any investigated hemodynamic parameters due to device implantation. Even in cases, where devices were implanted in a non-optimal manner, no marked differences in hemodynamic parameters compared to devices implanted in an optimal position were found. Before implantation time and surface-averaged WSS was 2.35 ± 0.47  Pa, whereas OSI was 0.08 ± 0.17 , respectively. Areas affected by low WSS magnitudes were 2.5 ± 2.7  cm2 , whereas the areas affected by high OSI were 18.1 ± 6.3  cm2 . After device implantation, WSS and OSI were 2.45 ± 0.49  Pa and 0.08 ± 0.16 , respectively. Surface areas affected by low WSS and high OSI were 2.9 ± 2.7  cm2 , and 18.4 ± 6.1  cm2 , respectively. This in-silico study indicates that no clinically relevant differences in intra-arterial hemodynamics are occurring after device implantation, even at non-optimal positioning of the sensor. Simultaneously, no embolic events were observed, suggesting that the risk for thrombus formation after device implantation is low and independent of the sensor position.

Hemodynamics of thrombus formation in intracranial aneurysms: An in silico observational study

How prevalent is spontaneous thrombosis in a population containing all sizes of intracranial aneurysms? How can we calibrate computational models of thrombosis based on published data? How does spontaneous thrombosis differ in normo- and hypertensive subjects? We address the first question through a thorough analysis of published datasets that provide spontaneous thrombosis rates across different aneurysm characteristics. This analysis provides data for a subgroup of the general population of aneurysms, namely, those of large and giant size (>10 mm). Based on these observed spontaneous thrombosis rates, our computational modeling platform enables the first in silico observational study of spontaneous thrombosis prevalence across a broader set of aneurysm phenotypes. We generate 109 virtual patients and use a novel approach to calibrate two trigger thresholds: residence time and shear rate, thus addressing the second question. We then address the third question by utilizing this calibrated model to provide new insight into the effects of hypertension on spontaneous thrombosis. We demonstrate how a mechanistic thrombosis model calibrated on an intracranial aneurysm cohort can help estimate spontaneous thrombosis prevalence in a broader aneurysm population. This study is enabled through a fully automatic multi-scale modeling pipeline. We use the clinical spontaneous thrombosis data as an indirect population-level validation of a complex computational modeling framework. Furthermore, our framework allows exploration of the influence of hypertension in spontaneous thrombosis. This lays the foundation for in silico clinical trials of cerebrovascular devices in high-risk populations, e.g., assessing the performance of flow diverters in aneurysms for hypertensive patients.

Predicting the risk of postsplenectomy thrombosis in patients with portal hypertension using computational hemodynamics models: A proof-of-concept study

BACKGROUND

The high incidence of thrombosis in the portal venous system following splenectomy (a frequently adopted surgery for treating portal hypertension in patients with splenomegaly and hypersplenism) is a critical clinical issue. The aim of this study was to address whether quantification of postsplenectomy hemodynamics has potential value for assessing the risk of postsplenectomy thrombosis.

METHODS

Computational models were constructed for three portal hypertensive patients treated with splenectomy based on their preoperative clinical data to quantify hemodynamics in the portal venous system before and after splenectomy, respectively. Each patient was followed up for three or five months after surgery and examined with CT to screen potential thrombosis.

FINDINGS

The area ratio of wall regions exposed to low wall shear stress was small before splenectomy in all patients, which increased markedly after splenectomy and exhibited enlarged inter-patient differences. The largest area ratio of low wall shear stress and most severe flow stagnation after splenectomy were predicted for the patient suffering from postsplenectomy thrombosis, with the wall regions exposed to low wall shear stress corresponding well with the CT-detected distribution of thrombus. Further analyses revealed that postoperative hemodynamic characteristics were considerably influenced by the anatomorphological features of the portal venous system.

INTERPRETATION

Postoperative hemodynamic conditions in the portal venous system are highly patient-specific and have a potential link to postsplenectomy thrombosis, which indicates that patient-specific hemodynamic studies may serve as a complement to routine clinical assessments for refining risk stratification and postoperative patient management.

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.

Influences of Anatomorphological Features of the Portal Venous System on Postsplenectomy Hemodynamic Characteristics in Patients With Portal Hypertension: A Computational Model-Based Study

Splenectomy, as an effective surgery for relieving complications caused by portal hypertension, is often accompanied by a significantly increased incidence of postoperative thrombosis in the portal venous system (PVS). While the underlying mechanisms remain insufficiently understood, the marked changes in hemodynamic conditions in the PVS following splenectomy have been suggested to be a potential contributing factor. The aim of this study was to investigate the influences of the anatomorphological features of the PVS on hemodynamic characteristics before and after splenectomy, with emphasis on identifying the specific anatomorphological features that make postoperative hemodynamic conditions more clot-promoting. For this purpose, idealized computational hemodynamics models of the PVS were constructed based on general anatomical structures and population-averaged geometrical parameters of the PVS. In the models, we incorporated various anatomorphological variations to represent inter-patient variability. The analyses of hemodynamic data were focused on the spatial distribution of wall shear stress (WSS) and the area ratio of wall regions exposed to low WSS (ALS). Obtained results showed that preoperative hemodynamic conditions were comparable among different models in terms of space-averaged WSS and ALS (all were small) irrespective of the considerable differences in spatial distribution of WSS, whereas, the inter-model differences in ALS were significantly augmented after splenectomy, with the value of ALS reaching up to over 30% in some models, while being smaller than 15% in some other models. Postoperative ALS was mainly determined by the anatomical structure of the PVS, followed by some morphogeometrical parameters, such as the diameter and curvature of the splenic vein, and the distance between the inferior mesenteric vein and splenoportal junction. Relatively, the angles between tributary veins and trunk veins only had mild influences on ALS. In addition, a marked increase in blood viscosity was predicted after splenectomy, especially in regions with low WSS, which may play an additive role to low WSS in initiating thrombosis. These findings suggest that the anatomical structure and some morphogeometrical features of the PVS are important determinants of hemodynamic conditions following splenectomy, which may provide useful clues to assessing the risk of postsplenectomy thrombosis based on medical imaging data.

Arteriovenous fistula maturation and the influence of fluid dynamics

Arteriovenous fistula creation is the preferred vascular access for haemodialysis therapy, but has a large failure rate in the maturation period. This period generally lasts 6 to 8 weeks after surgical creation, in which the vein and artery undergo extensive vascular remodelling. In this review, we outline proposed mechanisms for both arteriovenous fistula maturation and arteriovenous fistula failure. Clinical, animal and computational studies have not yet shown a definitive link between any metric and disease development, although a number of theories based on wall shear stress metrics have been suggested. Recent work allowing patient-based longitudinal studies may hold the key to understanding arteriovenous fistula maturation processes.

Haemodynamic impacts of myocardial bridge length: A congenital heart disease

OBJECTIVES

There is an association between long and thick myocardial bridging (MB), haemodynamic perturbations and increased risk of myocardial infarction. This study aims to investigate the alteration in coronary haemodynamics with increasing the length of MB.

METHODS

Angiography and intravascular ultrasound were performed in 10 patients with varying length of MB in the left anterior descending (LAD) artery. In silico models of MB were developed based on the reconstructed three-dimensional model of the LAD. The entire LAD was divided into 3 segments, proximal (pre-bridge), bridge and distal (post-bridge). Transient computational fluid dynamics simulations were performed to derive distribution of blood residence time and wall shear stress (WSS) over entire vessel including proximal, bridge and distal segments.

RESULTS

With increasing the length of MB, a decreasing trend was observed in the WSS over proximal segment whereas an increasing trend was found in the WSS over bridge segment. When patients were divided into 2 groups based on the average length of MB in the whole cohort (L_ave = 23.92 mm), patients with bridges longer than L_ave had smaller WSS and higher residence time in the proximal segment compared to those with bridges shorter than L_ave (0.59 ± 0.31 vs 0.21 ± 0.14 Pa and 0.0021 ± 0.0015 vs 0.0045 ± 0.0021 s). In contrast, patients with bridges longer than L_ave had greater WSS in the bridge segment compared to those with bridges shorter than L_ave (1.37 ± 1.66 vs 2.53 ± 3.14 Pa). No significant difference was found in the distal WSS of patients with short and long bridges.

CONCLUSION

Our findings revealed a direct relationship between the length of MB and haemodynamic perturbations in the proximal segment such that the increased length of MB is associated with decreased WSS and increased residence time.

Successful use of retrograde branched extension limb assembling technique in endovascular repair of pararenal abdominal aortic aneurysm

Surgeon-modified retrograde branched extension limb assembling technique and bridged endografts were successfully used to exclude an asymptomatic pararenal abdominal aortic aneurysm and to reconstruct the superior mesenteric artery and bilateral renal arteries in a case with high-grade celiac artery stenosis, nondilated aorta above the superior mesenteric artery, and large lumen below the renal arteries. In patient-specific models for hemodynamics analysis, enhanced flow diversion to visceral arteries up to 6-month follow-up confirmed treatment feasibility; however, endograft configurations could be improved to avoid sharp corners at bifurcations, thereby ensuring smooth flow transport and possibly reducing risk for endograft narrowing or the development of thrombosis.

Hemodynamic parameters that may predict false-lumen growth in type-B aortic dissection after endovascular repair: A preliminary study on long-term multiple follow-ups

Thoracic endovascular aortic repair (TEVAR) is commonly applied in type-B aortic dissection. For patients with dissection affects descending aorta and extends downward to involve abdominal aorta and possibly iliac arteries, false lumen (FL) expansion might occur post-TEVAR. Predictions of dissection development may assist in medical decision on re-intervention or surgery. In this study, two patients are selected with similar morphological features at initial presentation but with different long-term FL development post-TEVAR (stable and enlarged FL). Patient-specific models are established for each of the follow-ups. Flow boundaries and computational validations are obtained from Doppler ultrasound velocimetry. By analyzing the hemodynamic parameters, the false-to-true luminal pressure difference (PDiff) and particle relative residence time (RRT) are found related to FL remodeling. It is found that (i) the position of the first FL flow entry is the watershed of negative-and-positive PDiff and, in long-term follow-ups, and the position of largest PDiff is consistent with that of the greatest increase of FL width; (ii) high RRT occurs at the FL proximal tip and similar magnitude of RRT is found in both stable and enlarged cases; (iii) comparing to the RRT at 7days post-TEVAR, an increase of RRT afterwards in short-term is found in the stable case while a slight decrease of this parameter is found in the enlarged case, indicating that the variation of RRT in short-term post-TEVAR might be potential to predict long-term FL remodeling.

Spiral blood flow in aorta-renal bifurcation models

The presence of a spiral arterial blood flow pattern in humans has been widely accepted. It is believed that this spiral component of the blood flow alters arterial haemodynamics in both positive and negative ways. The purpose of this study was to determine the effect of spiral flow on haemodynamic changes in aorta-renal bifurcations. In this regard, a computational fluid dynamics analysis of pulsatile blood flow was performed in two idealised models of aorta-renal bifurcations with and without flow diverter. The results show that the spirality effect causes a substantial variation in blood velocity distribution, while causing only slight changes in fluid shear stress patterns. The dominant observed effect of spiral flow is on turbulent kinetic energy and flow recirculation zones. As spiral flow intensity increases, the rate of turbulent kinetic energy production decreases, reducing the region of potential damage to red blood cells and endothelial cells. Furthermore, the recirculation zones which form on the cranial sides of the aorta and renal artery shrink in size in the presence of spirality effect; this may lower the rate of atherosclerosis development and progression in the aorta-renal bifurcation. These results indicate that the spiral nature of blood flow has atheroprotective effects in renal arteries and should be taken into consideration in analyses of the aorta and renal arteries.