Time-averaged Wall Shear Stress Research Articles

Hemodynamic characteristics of pulsatile blood flow through bifurcated stenosed carotid artery

Purpose Carotid artery is often associated with plaque deposition because of its shape and associated flow features. The shape of stenosed bifurcation is characterised by bifurcation angle (ß), planarity angle (α) and severity of stenosis (b). In the present work, three-dimensional numerical computations have been performed to analyse the effect of these geometrical parameters of carotid bifurcation on the characteristics of flow. Design/methodology/approach Governing equations of this study were solved using ANSYS Fluent 20.1 and the blood flow was considered as laminar, pulsatile and non-Newtonian. Instantaneous flow behaviour has been illustrated using vorticity, velocity and helicity contours, whereas the time-averaged wall shear stress ( τw¯) and oscillatory shear index (OSI) quantify the time-averaged behaviour. Findings The recirculation zone and secondary flow are ascertained to be stronger for higher bifurcation angle as compared to the lower bifurcation angle. Strength of the secondary flow is found to reduce with increase in α from 0° to 10°, whereas it grows as α varies from 10° to 20°. For higher bifurcation angles, τw¯ is lower than 2 Pa and OSI is greater than 0.2 on the outer walls. Similar observations were made for τw¯ and OSI distribution on bottom wall in non-planar cases, which predicted atherogenic locations. Originality/value The values for ß were taken as 30°, 45°, 60° and 75°, whereas for α, range of 0°–20° was chosen. The stenosis was considered on the outer wall of internal carotid artery and its severity was considered within the range of 0%–60%.

Intra-Aneurysmal High-Resolution 4D MR Flow Imaging for Hemodynamic Imaging Markers in Intracranial Aneurysm Instability.

Prediction of aneurysm instability is crucial to guide treatment decisions and to select appropriate patients with unruptured intracranial aneurysms (IAs) for preventive treatment. High-resolution 4D MR flow imaging and 3D quantification of aneurysm morphology could offer insights and new imaging markers for aneurysm instability. In this cross-sectional study, we aim to identify 4D MR flow imaging markers for aneurysm instability by relating hemodynamics in the aneurysm sac to 3D morphologic proxy parameters for aneurysm instability. In 35 patients with 37 unruptured IAs, a 3T MRA and a 7T 4D MRI flow scan were performed. Five hemodynamic parameters-peak-systolic wall shear stress (WSSMAX) and time-averaged wall shear stress (WSSMEAN), oscillatory shear index (OSI), mean velocity, and velocity pulsatility index-were correlated to 6 3D morphology proxy parameters of aneurysm instability-major axis length, volume, surface area (all 3 size parameters), flatness, shape index, and curvedness-by Pearson correlation with 95% CI. Scatterplots of hemodynamic parameters that correlated with IA size (major axis length) were created. WSSMAX and WSSMEAN correlated negatively with all 3 size parameters (strongest for WSSMEAN with volume (r = -0.70, 95% CI -0.83 to -0.49) and OSI positively (strongest with major axis length [r = 0.87, 95% CI 0.76-0.93]). WSSMAX and WSSMEAN correlated positively with shape index (r = 0.61, 95% CI 0.36-0.78 and r = 0.49, 95% CI 0.20-0.70, respectively) and OSI negatively (r = -0.82, 95% CI -0.9 to -0.68). WSSMEAN and mean velocity correlated negatively with flatness (r = -0.35, 95% CI -0.61 to -0.029 and r = -0.33, 95% CI -0.59 to 0.007, respectively) and OSI positively (r = 0.54, 95% CI 0.26-0.74). Velocity pulsatility index did not show any statistically relevant correlation. Out of the 5 included hemodynamic parameters, WSSMAX, WSSMEAN, and OSI showed the strongest correlation with morphologic 3D proxy parameters of aneurysm instability. Future studies should assess these promising new imaging marker parameters for predicting aneurysm instability in longitudinal cohorts of patients with IA.

Assessing the impact of fetal-type posterior cerebral artery variations on cerebral hemodynamics

The circle of Willis (CoW) is a critical, arterial structure that ensures balanced, cerebral-blood supply. The fetal-type posterior cerebral artery (f-PCA) is a CoW variant that can significantly affect hemodynamics and elevate the risk of cerebrovascular diseases. This study used computational fluid dynamics simulations and a patient-specific, three-dimensional model to evaluate the hemodynamic effects of the f-PCA variants on cerebral-blood flow and key hemodynamic indices—such as time-averaged wall-shear stress (TAWSS), oscillatory shear index (OSI), pulsatility index, and resistive index. The fetal ratio (FR) is defined as the ratio of the diameter of the posterior communicating artery (PCoA) to that of the first segment (P1) of the PCA. Our findings indicate that as the FR increases, the contribution of the basilar artery to the second segment (P2) of PCA decreases significantly. Specifically, the flow rate through ipsilateral P1 decreased by 40.0% for FR = 1 and 70.9% for FR = 2, with the internal carotid artery (ICA) compensating for this reduction. Moreover, variations in f-PCA led to significant increases in TAWSS and OSI in key arterial segments (including the ipsilateral P1, PCoA, and the anterior communicating artery), which are associated with a higher risk of aneurysm initiation and growth. Under conditions of unilateral stenosis in the ipsilateral ICA, f-PCA models exhibit a more complex and pronounced impact on blood flow than models without f-PCA, emphasizing the need for detailed hemodynamic assessments in clinical evaluations and preoperative planning to mitigate the risks associated with CoW anatomical variations.

Simulation of plaque formation in a realistic geometry of a human aorta: effects of endothelial layer properties, heart rate, and hypertension.

Nowadays, cardiovascular diseases are the most common cause of death worldwide. Besides, atherosclerosis is a cardiovascular disease that occurs with persistent narrowing of arteries, especially medium and large-sized arteries. Atherosclerosis begins with a local elevation in the permeability of the arterial wall as a result of endothelial inflammation. Subsequently, excess LDL permeates into the arterial wall. Then, through several chemical responses and reactions, foam cells are produced. These foam cells serve as a crucial indicator for assessing the development of atherosclerosis within the arteries. In this study, the effect of endothelial layer modeling, heart rate (HR) and hypertension on the foam cell accumulation is numerically investigated in a patient-specific geometry of the human thoracic aorta. Navier-Stokes, Darcy, and mass transfer equations are used to obtain the velocity and concentration field within the domain. Regarding the dependence of endothelial cell properties on time-averaged wall shear stress, it is observed that foam cells are mainly concentrated in the outer curvature of the aortic arch, downstream of the left subclavian artery. However, considering oscillatory-shear-rate as the determinant of endothelial cell properties leads to the accumulation of foam cells in the inner curvature of the descending aorta. Regarding the HR, with the increase of HR, the volume average concentration of the foam cell decreases. However, there is no substantial difference between the cases of different HRs. Moreover, foam cell concentration significantly increases in the hypertension case. This result implies that a slight increase in the blood pressure may induce irreparable problems in the circulatory system.

Computational Model for Early-Stage Aortic Valve Calcification Shows Hemodynamic Biomarkers.

Heart disease is a leading cause of mortality, with calcific aortic valve disease (CAVD) being the most prevalent subset. Being able to predict this disease in its early stages is important for monitoring patients before they need aortic valve replacement surgery. Thus, this study explored hydrodynamic, mechanical, and hemodynamic differences in healthy and very mildly calcified porcine small intestinal submucosa (PSIS) bioscaffold valves to determine any notable parameters between groups that could, possibly, be used for disease tracking purposes. Three valve groups were tested: raw PSIS as a control and two calcified groups that were seeded with human valvular interstitial and endothelial cells (VICs/VECs) and cultivated in calcifying media. These two calcified groups were cultured in either static or bioreactor-induced oscillatory flow conditions. Hydrodynamic assessments showed metrics were below thresholds associated for even mild calcification. Young's modulus, however, was significantly higher in calcified valves when compared to raw PSIS, indicating the morphological changes to the tissue structure. Fluid-structure interaction (FSI) simulations agreed well with hydrodynamic results and, most notably, showed a significant increase in time-averaged wall shear stress (TAWSS) between raw and calcified groups. We conclude that tracking hemodynamics may be a viable biomarker for early-stage CAVD tracking.

Relative Residence Time Can Account for Half of the Anatomical Variation in Fatty Streak Prevalence Within the Right Coronary Artery.

The patchy anatomical distribution of atherosclerosis has been attributed to variation in haemodynamic wall shear stress (WSS). The consensus is that low WSS and a high Oscillatory Shear Index (OSI) trigger the disease. We found that atherosclerosis at aortic branch sites correlates threefold better with transverse WSS (transWSS), a metric which quantifies multidirectional near-wall flow. Coronary artery disease has greater clinical significance than aortic disease but computation of WSS metrics is complicated by the substantial vessel motion occurring during each cardiac cycle. Here we present the first comparison of the distribution of atherosclerosis with WSS metrics computed for moving coronary arteries. Maps of WSS metrics were computed using dynamic geometries reconstructed from angiograms of ten non-stenosed human right coronary arteries (RCAs). They were compared with maps of fatty streak prevalence derived from a previous study of 1852 RCAs. Time average WSS (TAWSS), OSI, transWSS and the cross-flow index (CFI), a non-dimensional form of the transWSS, gave non-significant or significant but low spatial correlations with lesion prevalence. The highest correlation coefficient (0.71) was for the relative residence time (RRT), a metric that decreases with TAWSS and increases with OSI. The coefficient was not changed if RRT was calculated using CFI, which captures multidirectional WSS only, rather than OSI, which encompasses both multidirectional and oscillatory WSS. Contrary to our earlier findings in the aorta, low WSS in combination with highly multidirectional flow correlates best with lesion location in the RCA, explaining approximately half of its anatomical variation.

Intracranial bypass for giant aneurysms treatment assessed by computational fluid dynamics (CFD) analysis

Unruptured giant intracranial aneurysms (GIA) are those with diameters of 25 mm or greater. As aneurysm size is correlated with rupture risk, GIA natural history is poor. Parent artery occlusion or trapping plus bypass revascularization should be considered to encourage intra-aneurysmal thrombosis when other treatment options are contraindicated. The mechanistic background of these methods is poorly studied. Thus, we assessed the potential of computational fluid dynamics (CFD) and fluid–structure interaction (FSI) analyses for clinical use in the preoperative stage. A CFD investigation in three patient-specific flexible models of whole arterial brain circulation was performed. A C6 ICA segment GIA model was created based on CT angiography. Two models were then constructed that simulated a virtual bypass in combination with proximal GIA occlusion, but with differing middle cerebral artery (MCA) recipient vessels for the anastomosis. FSI and CFD investigations were performed in three models to assess changes in flow pattern and haemodynamic parameters alternations (wall shear stress (WSS), oscillatory shear index (OSI), maximal time averaged WSS (TAWSS), and pressure). General flow splitting across the entire domain was affected by virtual bypass procedures, and any deficiency was partially compensated by a specific configuration of the circle of Willis. Following the implementation of bypass procedures, a reduction in haemodynamic parameters was observed within the aneurysm in both cases under analysis. In the case of the temporal MCA branch bypass, the decreases in the studied parameters were slightly greater than in the frontal MCA branch bypass. The reduction in the magnitude of the chosen area-averaged parameters (averaged over the aneurysm wall surface) was as follows: WSS 35.7%, OSI 19.0%, TAWSS 94.7%, and pressure 24.2%. FSI CFD investigation based on patient-specific anatomy models with subsequent stimulation of virtual proximal aneurysm occlusion in conjunction with bypass showed that this method creates a pro-thrombotic favourable environment whilst reducing intra-aneurysmal pressure leading to shrinking. MCA branch recipient selection for optimum haemodynamic conditions should be evaluated individually in the preoperative stage.

Wall Shear Stress (WSS) Analysis in Atherosclerosis in Partial Ligated Apolipoprotein E Knockout Mouse Model through Computational Fluid Dynamics (CFD).

Atherosclerosis involves an inflammatory response due to plaque formation within the arteries, which can lead to ischemic stroke and heart disease. It is one of the leading causes of death worldwide, with various contributing factors such as hyperlipidemia, hypertension, obesity, diabetes, and smoking. Wall shear stress (WSS) is also known as a contributing factor of the formation of atherosclerotic plaques. Since the causes of atherosclerosis cannot be attributed to a single factor, clearly understanding the mechanisms and causes of its occurrence is crucial for preventing the disease and developing effective treatment strategies. To better understand atherosclerosis and define the correlation between various contributing factors, computational fluid dynamics (CFD) analysis is primarily used. CFD simulates WSS, the frictional force caused by blood flow on the vessel wall with various hemodynamic changes. Using apolipoprotein E knockout (ApoE-KO) mice subjected to partial ligation and a high-fat diet at 1-week, 2-week, and 4-week intervals as an atherosclerosis model, CFD analysis was conducted along with the reconstruction of carotid artery blood flow via magnetic resonance imaging (MRI) and compared to the inflammatory factors and pathological staining. In this experiment, a comparative analysis of the effects of high WSS and low WSS was conducted by comparing the standard deviation of time-averaged wall shear stress (TAWSS) at each point within the vessel wall. As a novel approach, the standard deviation of TAWSS within the vessel was analyzed with the staining results and pathological features. Since the onset of atherosclerosis cannot be explained by a single factor, the aim was to find the correlation between the thickness of atherosclerotic plaques and inflammatory factors through standard deviation analysis. As a result, the gap between low WSS and high WSS widened as the interval between weeks in the atherosclerosis mouse model increased. This finding not only linked the occurrence of atherosclerosis to WSS differences but also provided a connection to the causes of vulnerable plaques.

Impact of blood viscosity on hemodynamics of large intracranial aneurysms

BackgroundHemodynamic factors play an important role in the formation and rupture of intracranial aneurysms. Blood viscosity has been recognized as a potential factor influencing the hemodynamics of aneurysms. Computational fluid dynamics (CFD) is one of the main methods to study aneurysm hemodynamics. However, current CFD studies often set the viscosity to a standard value, neglecting the effect of individualized viscosity on hemodynamics. We investigate the impact of blood viscosity on hemodynamics in large intracranial aneurysm (IA) and assess the potential implications for aneurysm growth and rupture risk. MethodsCFD simulations of 8 unruptured large internal carotid artery aneurysms were conducted using pulsatile inlet conditions. For each aneurysm, CFD simulations were performed at 5 different viscosity levels (0.004, 0.006, 0.008, 0.010, and 0.012 Pa·s). Differences in hemodynamic parameters across viscosity levels were compared using paired t-tests, and the correlation between viscosity and hemodynamic parameters was analyzed. ResultsIncreasing blood viscosity leads to significant decrease in blood flow velocity within aneurysms. Time-averaged wall shear stress (WSS) showed significant positive correlation with viscosity, particularly at the aneurysm neck. Oscillatory shear index (OSI) showed general decreasing trend with increased viscosity, while it displayed an irregular pattern in a few cases. ConclusionsVariations in viscosity markedly influence velocity, WSS, and OSI in aneurysms, suggesting a role in modulating aneurysm growth and rupture risk. Incorporating patient-specific viscosity values in CFD simulations is vital for accurate and reliable outcomes.

A new method for scaling inlet flow waveform in hemodynamic analysis of aortic dissection.

Computational fluid dynamics (CFD) simulations have shown great potentials in cardiovascular disease diagnosis and postoperative assessment. Patient-specific and well-tuned boundary conditions are key to obtaining accurate and reliable hemodynamic results. However, CFD simulations are usually performed under non-patient-specific flow conditions due to the absence of in vivo flow and pressure measurements. This study proposes a new method to overcome this challenge by tuning inlet boundary conditions using data extracted from electrocardiogram (ECG). Five patient-specific geometric models of type B aortic dissection were reconstructed from computed tomography (CT) images. Other available data included stoke volume (SV), ECG, and 4D-flow magnetic resonance imaging (MRI). ECG waveforms were processed to extract patient-specific systole to diastole ratio (SDR). Inlet boundary conditions were defined based on a generic aortic flow waveform tuned using (1) SV only, and (2) with ECG and SV (ECG + SV). 4D-flow MRI derived inlet boundary conditions were also used in patient-specific simulations to provide the gold standard for comparison and validation. Simulations using inlet flow waveform tuned with ECG + SV not only successfully reproduced flow distributions in the descending aorta but also provided accurate prediction of time-averaged wall shear stress (TAWSS) in the primary entry tear (PET) and abdominal regions, as well as maximum pressure difference, ∆Pmax, from the aortic root to the distal false lumen. Compared with simulations with inlet waveform tuned with SV alone, using ECG + SV in the tuning method significantly reduced the error in false lumen ejection fraction at the PET (from 149.1% to 6.2%), reduced errors in TAWSS at the PET (from 54.1% to 5.7%) and in the abdominal region (from 61.3% to 11.1%), and improved ∆Pmax prediction (from 283.1% to 18.8%) However, neither of these inlet waveforms could be used for accurate prediction of TAWSS in the ascending aorta. This study demonstrates the importance of SDR in tailoring inlet flow waveforms for patient-specific hemodynamic simulations. A well-tuned flow waveform is essential for ensuring that the simulation results are patient-specific, thereby enhancing the confidence and fidelity of computational tools in future clinical applications.

Hemodynamics in the treatment of pseudoaneurysm caused by extreme constriction of aortic arch with coated stent.

To evaluate the changes in distal vascular morphology and hemodynamics in patients with extremely severe aortic coarctation (CoA) after covered palliative (CP) stent dilation with different surgical strategies. Perioperative computed tomography angiography and digital subtraction angiography were utilized to construct three aortic models with varying stenosis rates and one follow-up model in a patient with extremely severe CoA. The models included: an idealized non-stenosed model (A: 0%), a model post initial stent deployment (B: 28%), a model post balloon expansion (C: 39%), and a model 18 months after post-balloon expansion (D: 39%). Consistent boundary conditions were applied to all models, and hemodynamic simulation was conducted using the pure fluid method. The narrowest and distal diameter of the stent increased by 34.71% and 59.29%, respectively, from model B to C. Additionally, the distal diameter of the stent increased by -13.80% and +43.68% compared to the descending aorta diameter, respectively. Furthermore, the ellipticity of the maximum cross-section of the aneurysm region in model A to D continued to increase. The oscillatory shear index at the stenosis to the region of the aneurysm were found to be higher in Models A and B, and lower in Models C and D. At the moment of maximum flow velocity, the blood flow distribution in models A and B was more uniform in the widest section of the blood vessels at the distal end of the stenosis, whereas models C and D exhibited disturbed blood flow with more than 2 eddy currents. The time-averaged wall shear stress (TAWSS) decreased in the distal and basal aneurysms, while it significantly increased at the step position. The aneurysmal region exhibited an endothelial cell activation potential value lower than 0.4 Pa-1. In patients with extremely severe CoA, it is crucial to ensure that the expanded diameter at both ends of the CP stent does not exceed the native vascular diameter during deployment. Our simulation results demonstrate that overdilation leads to a decrease in the TAWSS above the injured vessel, creating an abnormal hemodynamic environment that may contribute to the development and enlargement of false aneurysms in the early postoperative period. ClinicalTrials.gov, (NCT02917980).

Numerical assessment of using various outlet boundary conditions on the hemodynamics of an idealized left coronary artery model.

Vascular diseases are greatly influenced by the hemodynamic parameters and the accuracy of determining these parameters depends on the use of correct boundary conditions. The present work carries out a two-way fluid-structure interaction (FSI) simulation to investigate the effects of outlet pressure boundary conditions on the hemodynamics through the left coronary artery bifurcation with moderate stenosis (50%) in the left anterior descending (LAD) branch. The Carreau viscosity model is employed to characterise the shear-thinning behaviour of blood. The results of the study reveal that the employment of zero pressure at the outlet boundaries significantly overestimates the values of hemodynamic variables like wall shear stress (WSS), and time-averaged wall shear stress (TAWSS) compared with human healthy and pulsatile pressure outlet conditions. However, the difference between these variables is marginally low for human healthy and pulsatile pressure outlets. The oscillatory shear index (OSI) remains the same across all scenarios, indicating independence from the outlet boundary condition. Furthermore, the magnitude of negative axial velocity and pressure drop across the plaque are found to be higher at the zero pressure outlet boundary condition.

Haemodynamic study of left nonthrombotic iliac vein lesions: a preliminary report

Nonthrombotic iliac vein lesions (NIVLs) are significant causes of chronic venous insufficiency (CVI) in the left lower limb and symptom recurrence following left lower limb varicose vein treatment. The goal of this study was to explore the haemodynamic and morphological characteristics of iliac veins in patients with NIVLs. Pressure at the caudal end of the stenotic left common iliac vein (LCIV) segment, local blood flow velocity, and time-averaged wall shear stress in the stenotic segment exhibited positive correlations with the clinical CVI classification (R = 0.92, p < 0.001; R = 0.94, p < 0.001; R = 0.87, p < 0.001), while the relative retention time showed a negative correlation (R = -0.94, p < 0.001). The pressure difference (∆P) between the two ends of the stenotic segment and the velocity difference (∆V) between the stenotic segment and the caudal end were positively correlated with the clinical classification (R = 0.92, p < 0.001; R = 0.9, p < 0.001). The cross-sectional area stenosis rate and length of the stenotic LCIV segment were positively correlated with the clinical classification (R = 0.93, p < 0.001; R = 0.63, p < 0.001). The results suggest that haemodynamic assessment of the iliac vein could effectively portray blood flow disturbances in stenotic segments of the LCIV, potentially reflecting the degree of iliac vein stenosis. Haemodynamic indicators are correlated with the severity of clinical CVI symptoms.

Validation of classification system for isolated superior mesenteric artery dissections using image-based computational flow analysis

Background and ObjectiveThe isolated superior mesenteric artery dissection (ISMAD) is a rare but potentially fatal vascular disorder. Classifications for ISMAD were previously proposed based on morphometric features. However, the classification systems were not standardized and verified yet. This study conducted computational flow analysis to validate the latest classification system of ISMAD and aid clinical decision-making based on hemodynamic parameters. Methods62 patients with ISMAD were included and classified into different types according to false lumen structures (five types, Type I-V) and true lumen patency (two types, Type P and Type S) according to Qiu classification system. Computational fluid dynamics and three-dimensional structural analyses were conducted on the basis of computed tomography angiography datasets. Quantitative and qualitative functional analyses were performed via parameters of interest including volume flow of each minute, pressure drop, pressure gradient, the derivative parameters of wall shear stress such as time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and the relative residence time (RRT). Statistical analyses were conducted among different ISMAD types. ResultsTAWSS, OSI and RRT showed significant difference among different types when classified using false lumen structures. In detail, Type IV showed significantly higher TAWSS than other types (p = 0.007). OSI was obviously higher in Type II (p = 0.015). Type IV also presented the lowest RRT (p = 0.005). The pressure drop, pressure gradient, OSI and RRT showed higher value in Type S than that in Type P, demonstrating a statistical significance with p values of 0.017, 0.041, 0.001 and 0.012, respectively. While Type P had larger volume flow than Type S (p = 0.041). ConclusionsThe notable differences in hemodynamic features among different types demonstrated the feasibility of Qiu classification system. The evaluation based on hemodynamic simulation might also provide insights into risk identification and guide therapeutic decisions for ISMAD.

Role of aneurysmal hemodynamic changes in pathogenesis of headaches associated with unruptured cerebral aneurysms.

One symptom commonly associated with the presence of unruptured intracranial aneurysms is headache. In this study, the authors aimed to analyze factors associated with headaches among patients with intracranial aneurysms, with special consideration of hemodynamic parameters. The authors prospectively included 96 patients with 122 unruptured intracranial aneurysms. The authors obtained detailed medical history including current diseases and medications, as well as blood pressure values taken during hospitalization from the patients' medical records. The short-form McGill Pain Questionnaire was administered to each patient at admission and 3-6 months after the procedure to assess type and severity of headache. Based on imaging data, the authors obtained 3D reconstruction of each patients' aneurysm dome with feeding artery. The authors performed computational fluid dynamics analysis of blood flow through prepared models using OpenFOAM. Blood was modeled as Newtonian fluid, using the incompressible transient solver. Patient-specific internal carotid artery (ICA) blood velocity waves obtained with Doppler ultrasound were set as inlet boundary conditions. After performing simulation, the authors calculated the hemodynamic parameters of the aneurysm dome. A total of 30 patients (31.25%) reported having headaches. In multivariate logistic regression analysis, female sex (OR 2.81, 95% CI 2.51-4.86; p < 0.01), ICA aneurysm location (OR 7.93, 95% CI 5.51-8.52; p < 0.01), multiple aneurysms (OR 6.05, 95% CI 1.83-11.83; p = 0.02), mean dome blood velocity (OR 3.10, 95% CI 2.01-3.30; p < 0.01) and time-averaged wall shear stress (OR 1.18, 95% CI 1.47-2.72; p = 0.04) were independently associated with the presence of headache. Additionally, 17 patients (56.67%) reported complete relief of symptoms after the procedure. In multivariate logistic regression analysis, the mean blood flow in the ICA was independently associated with complete resolution of headaches after aneurysm treatment (OR 2.32, 95% CI 1.57-3.28; p < 0.01). Hemodynamic parameters of intracranial aneurysms might be associated with headaches and their relief after aneurysm treatment.

In-silico analysis of hemodynamic indicators in idealized stented coronary arteries for varying stent indentation

In this work, we investigate the effects of stent indentation on hemodynamic indicators in stented coronary arteries. Our aim is to assess in-silico risk factors for in-stent restenosis (ISR) and thrombosis after stent implantation. The proposed model is applied to an idealized artery with Xience V stent for four indentation percentages and three mesh refinements. We analyze the patterns of hemodynamic indicators arising from different stent indentations and propose an analysis of time-averaged WSS (TAWSS), topological shear variation index (TSVI), oscillatory shear index (OSI), and relative residence time (RRT). We observe that higher indentations display higher frequency of critically low TAWSS, high TSVI, and non-physiological OSI and RRT. Furthermore, an appropriate mesh refinement is needed for accurate representation of hemodynamics in the stent vicinity. The results suggest that disturbed hemodynamics could play a role in the correlation between high indentation and ISR.

Blood vessel wall shear stress determines regions of liposome accumulation in angiogenic vasculature.

Nanoparticles used for drug delivery often require intravenous administration exposing them to fluid forces within the vasculature, yet the impact of blood flow on nanoparticle delivery remains incompletely understood. Here, we utilized transgenic zebrafish embryos to investigate the relationship between the accumulation of fluorescently labeled PEGylated liposomes and various hemodynamic factors (such as flow velocity, wall shear stress (WSS), and flow pattern) across a wide range of angiogenic blood vessels. We reconstructed 3D models of vascular structures from confocal images and used computational fluid dynamics to calculate local WSS, velocities, and define flow patterns. The spatial distribution of fluorescently labeled liposomes was subsequently mapped within the same 3D space and correlated with local hemodynamic parameters. Through the integration of computational fluid dynamics and in vivo experimentation, we show that liposomes accumulated in vessel regions with WSS between 0.1-0.8 Pa, displaying an inverse linear correlation (R2 > 0.85) between time-averaged wall shear stress and liposome localization in vivo. Interestingly, flow pattern did not appear to impact liposome accumulation. Collectively, our findings suggest the potential of stealth liposomes for passive targeting of low-flow vasculature, including capillaries and intricate angiogenic vasculature resembling that of tumor vessel networks.

A comprehensive MRI-based computational model of blood flow in compliant aorta using radial basis function interpolation

BackgroundProperly understanding the origin and progression of the thoracic aortic aneurysm (TAA) can help prevent its growth and rupture. For a better understanding of this pathogenesis, the aortic blood flow has to be studied and interpreted in great detail. We can obtain detailed aortic blood flow information using magnetic resonance imaging (MRI) based computational fluid dynamics (CFD) with a prescribed motion of the aortic wall.MethodsWe performed two different types of simulations—static (rigid wall) and dynamic (moving wall) for healthy control and a patient with a TAA. For the latter, we have developed a novel morphing approach based on the radial basis function (RBF) interpolation of the segmented 4D-flow MRI geometries at different time instants. Additionally, we have applied reconstructed 4D-flow MRI velocity profiles at the inlet with an automatic registration protocol.ResultsThe simulated RBF-based movement of the aorta matched well with the original 4D-flow MRI geometries. The wall movement was most dominant in the ascending aorta, accompanied by the highest variation of the blood flow patterns. The resulting data indicated significant differences between the dynamic and static simulations, with a relative difference for the patient of 7.47±14.18% in time-averaged wall shear stress and 15.97±43.32% in the oscillatory shear index (for the whole domain).ConclusionsIn conclusion, the RBF-based morphing approach proved to be numerically accurate and computationally efficient in capturing complex kinematics of the aorta, as validated by 4D-flow MRI. We recommend this approach for future use in MRI-based CFD simulations in broad population studies. Performing these would bring a better understanding of the onset and growth of TAA.

Understanding the role of the left atrial appendage on the flow in the atrium.

The exclusion/occlusion of the left atrial appendage (LAA) is a treatment option for atrial fibrillation (AF) patients who are at high risk of stroke and high risk of bleeding. As the role of the LAA is not well understood or explored, this study aims to assess its role on flow dynamics in the left atrium. Computational fluid dynamics (CFD) simulations were carried out for nine AF patients before and after LAA exclusion. The flow parameters investigated included the LA velocities, Time Averaged Wall Shear Stress (TAWSS), Oscillatory Shear Index (OSI), Relative Residence Time (RRT), and Pressure in the LA. This study shows that, on average, a decrease in TAWSS (1.82 ± 1.85 Pa to 1.27 ± 0.96 Pa, p < 0.05) and a slight increase in OSI (0.16 ± 0.10 to 0.17 ± 0.10, p < 0.05), RRT (1.87 ± 1.84 Pa-1 to 2.11 ± 1.78 Pa-1, p < 0.05), and pressure (-19.2 ± 6.8 mmHg to -15.3 ± 8.3 mmHg, p < 0.05) were observed in the LA after the exclusion of the LAA, with a decrease in low-magnitude velocities. The exclusion of the LAA seems to be associated with changes in LA flow dynamics. Further studies are needed to elucidate the clinical implications of these changes.

Non-invasive fractional flow reserve estimation in coronary arteries using angiographic images

Coronary artery disease is the leading global cause of mortality and Fractional Flow Reserve (FFR) is widely regarded as the gold standard for assessing coronary artery stenosis severity. However, due to the limitations of invasive FFR measurements, there is a pressing need for a highly accurate virtual FFR calculation framework. Additionally, it’s essential to consider local haemodynamic factors such as time-averaged wall shear stress (TAWSS), which play a critical role in advancement of atherosclerosis. This study introduces an innovative FFR computation method that involves creating five patient-specific geometries from two-dimensional coronary angiography images and conducting numerical simulations using computational fluid dynamics with a three-element Windkessel model boundary condition at the outlet to predict haemodynamic distribution. Furthermore, four distinct boundary condition methodologies are applied to each geometry for comprehensive analysis. Several haemodynamic features, including velocity, pressure, TAWSS, and oscillatory shear index are investigated and compared for each case. Results show that models with average boundary conditions can predict FFR values accurately and observed errors between invasive FFR and virtual FFR are found to be less than 5%.