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%.

Rupture point is associated with divergent hemodynamics in intracranial aneurysms.

Background: Understanding the risk factors leading to intracranial aneurysm (IA) rupture have still not been fully clarified. They are vital for proper medical guidance of patients harboring unruptured IAs. Clarifying the hemodynamics associated with the point of rupture could help could provide useful information about some of the risk factors. Thus far, few studies have studied this issue with often diverging conclusions. Methods: We identified a point of rupture in patients operated for an IAs during surgery, using a combination of preoperative computed tomography (CT) and computed tomography angiography (CTA). Hemodynamic parameters were calculated both for the aneurysm sac as a whole and the point of rupture. In two cases, the results of CFD were compared with those of the experiment using particle image velocimetry (PIV). Results: We were able to identify 6 aneurysms with a well-demarcated point of rupture. In four aneurysms, the rupture point was near the vortex with low wall shear stress (WSS) and high oscillatory shear index (OSI). In one case, the rupture point was in the flow jet with high WSS. In the last case, the rupture point was in the significant bleb and no specific hemodynamic parameters were found. The CFD results were verified in the PIV part of the study. Conclusion: Our study shows that different hemodynamic scenarios are associated with the site of IA rupture. The numerical simulations were confirmed by laboratory models. This study further supports the hypothesis that various pathological pathways may lead to aneurysm wall damage resulting in its rupture.

Association of carotid wall shear stress measured by vector flow mapping technique with ba-PWV: a pilot study.

Objective: Arterial stiffness is an important tissue biomarker of the progression of atherosclerotic diseases. Brachial-ankle pulse wave velocity (ba-PWV) is a gold standard of arterial stiffness measurement widely used in Asia. Changes in vascular wall shear stress (WSS) lead to artery wall remodeling, which could give rise to an increase in arterial stiffness. The study aimed to explore the association between ba-PWV and common carotid artery (CCA) WSS measured by a newly invented vascular vector flow mapping (VFM) technique. Methods: We included 94 subjects free of apparent cardiovascular disease (CVD) and divided them into a subclinical atherosclerosis (SA) group (N = 47) and non subclinical atherosclerosis (NSA) group (N = 47). CCA WSS was measured using the VFM technique. Bivariate correlations between CCA WSS and other factors were assessed with Pearson's, Spearman's, or Kendall's coefficient of correlation, as appropriate. Partial correlation analysis was conducted to examine the influence of age and sex. Multiple linear stepwise regression was used for the analysis of independent determinants of CCA WSS. Receiver operating characteristic (ROC) analysis was performed to find the association between CCA WSS and 10-year CVD risk. Results: The overall subjects had a mean age of 47.9 ± 11.2 years, and males accounted for 52.1%. Average systolic CCA WSS was significantly correlated with ba-PWV (r = -0.618, p < 0.001) in the SA group. Multiple linear stepwise regression analysis confirmed that ba-PWV was an independent determinant of average systolic CCA WSS (ß = -0.361, p = 0.003). The area under the curve (AUC) of average systolic CCA WSS for 10-year CVD risk ≥10% was 0.848 (p < 0.001) in the SA group. Conclusions: Average systolic CCA WSS was significantly correlated with ba-PWV and was associated with 10-year CVD risk ≥10% in the SA group. Therefore, CCA WSS measured by the VFM technique could be used for monitoring and screening subjects with potential CVD risks.

Glagov revisited: coronary artery disease phenotype on non-invasive imaging provides rationale for implementing preventive pharmacotherapy-a case report.

Background: Glagov et al. showed that no reduction in vessel lumen occurred until the atherosclerotic plaque burden exceeded 40% of the vessel area. Most major adverse cardiac events occurring in the first 4 years after a myocardial infarction arise from untreated angiographically mild, non-flow-limiting lesions at the time of the index event. We report how computed tomography (CT) coronary angiography (CCTA) can be used to non-invasively risk stratify a patient with non-obstructive coronary artery disease (CAD) and guide further management. Case summary: A 69-year-old non-smoking female with hypertension, dyslipidaemia, and hypothyroidism presented with atypical chest pain. Electrocardiogram and left ventricular ejection fraction were normal. Her lipidic profile was normal. CCTA showed a lipid-rich plaque with very low attenuation (<30 HU) in the left main stem (LMS) extending into the proximal left anterior descending (LAD) and in the mid LAD artery. The maximum plaque burden in the LMS was 67% with a remodelling index of 1.375, and an area stenosis of 22%. Tissue characterization showed a lipid-rich plaque with a thin fibrous cap. The perivascular fat attenuation index (FAI) in the proximal LAD was suggestive of (-69 HU) inflamed perivascular fat. Shear stress analysis of the LMS plaque showed normal wall shear stress (WSS); however, the axial plaque stress was high. Her medications were intensified to rosuvastatin 20 mg once daily (OD) and ezetimibe 10 mg OD. The patient remained asymptomatic at 6 months follow-up. Discussion: Our case exemplifies the value of CCTA as a diagnostic 'one-stop shop' (CCTA, finite element analysis, computed tomographic density [CTD], tissue characterization analysis, FAI analysis, WSS and wall strain, and etc.) when stratifying a patient with non-obstructive CAD. With further development of novel potent anti-lipidaemic and anti-inflammatory drugs, non-obstructive lesions with adverse plaque and haemodynamic parameters will have the opportunity to be treated with additional preventive pharmacological therapy.

Pathophysiology of intracranial aneurysms in monozygotic twins: A rare case study from hemodynamic perspectives.

Hemodynamic mechanisms of the formation and growth of intracranial aneurysms (IA) in monozygotic twins (MTs) are still under-reported. To partially fill such knowledge gap, this study employed an experimentally validated numerical model to compare hemodynamics in 3 anatomical and 5 ablation study neurovascular models from a rare pair of MTs in terms of 7 critical hemodynamic parameters. Numerical results showed significant differences in hemodynamics between the MTs, although they share the same genes, indicating that genetic mutation and environmental factors might affect neurovascular morphologies and cause hemodynamic changes. After virtual removals of IAs in the ablation study, the locations where the aneurysmal sac/bleb generated in bifurcated anterior cerebral arteries (ACAs) register a locally high instantaneous wall shear stress (IWSS) of 52.9 and 70.1 Pa at the systolic peak in twin A and twin B, respectively. Same scenario can be observed in the distribution of instantaneous wall shear stress gradient (IWSSG), with 571.1 Pa/mm for twin A and 301.3 Pa/mm for twin B due to aggressive blood impingements, leading to IA generation. The fenestrated complex approaching ACA bifurcations in twin A may assist IA growth and rupture, via. Causing abnormal IWSS of 116.3 Pa, IWSSG of 832.5 Pa/mm, and oscillatory shear index (OSI) of 0.49. The bleb in twin B has high risks of progression and possible rupture as the IA suffers relatively low IWSS and high OSI. Additionally, IA generation can change blood flow rates in each connected artery, then affecting blood supplies to associated tissues and organs.

Validation of hemodynamic stress calculation in coronary computed tomography angiography versus intravascular ultrasound.

Background: Development in computational fluid dynamics and 3D construction could facilitate the calculation of hemodynamic stresses in coronary computed tomography angiography (CCTA). However, the agreement between CCTA derived stresses and intravascular ultrasound/intravascular coronary angiography (IVUS/ICA)-derived stresses remains undetermined. Thus, the purpose of this study is to investigate if CCTA can serve as alternative to IVUS/ICA for hemodynamic evaluation. Methods: In this retrospective study, 13 patients (14 arteries) with unstable angina who underwent both CCTA and IVUS/ICA at an interval of less than 7 days were consecutively included at the Chinese PLA General Hospital within the year of 2021. Slice-level minimal lumen area (MLA), percent area stenosis, velocity, pressure, Reynolds number, wall shear stress (WSS) and axial plaque stress (APS) were determined by both modalities. The agreement between CCTA and IVUS/ICA was assessed using the intraclass correlation coefficient (ICC), Pearson's correlation coefficient and Bland-Altman analysis. Results: CCTA overestimated the degree of area stenosis (50.22%±16.15% vs. 36.41%±19.37%, P=0.004) with the MLA showing no significant difference (5.81±2.24 vs. 6.72±2.04 mm2, P=0.126). No statistical difference was observed in WSS (6.57±6.26 vs. 5.98±5.55 Pa, P=0.420) and APS (16.03±1,159.45 vs. -1.27±890.39 Pa, P=0.691) between CCTA and IVUS. Good correlation was found in velocity (ICC: 0.796, 95% CI: 0.752-0.833), Reynolds number (ICC: 0.810, 95% CI: 0.768-0.844) and WSS (ICC: 0.769, 95% CI: 0.718-0.810), while the ICC of APS was (ICC: 0.341, 95% CI: 0.197-0.458), indicating a relatively poor correlation. Conclusions: CCTA can serve as a satisfactory alternative to the reference standard, IVUS/ICA in morphology simulation and hemodynamic stress calculation, especially in the calculation of WSS.

Aortic wall shear stress in bicuspid aortic valve disease-10-year follow-up.

Background: Bicuspid aortic valve (BAV) disease leads to deviant helical flow patterns especially in the mid-ascending aorta (AAo), potentially causing wall alterations such as aortic dilation and dissection. Among others, wall shear stress (WSS) could contribute to the prediction of long-term outcome of patients with BAV. 4D flow in cardiovascular magnetic resonance (CMR) has been established as a valid method for flow visualization and WSS estimation. The aim of this study is to reevaluate flow patterns and WSS in patients with BAV 10 years after the initial evaluation. Methods: Fifteen patients (median age 34.0 years) with BAV were re-evaluated 10 years after the initial study from 2008/2009 using 4D flow by CMR. Our particular patient cohort met the same inclusion criteria as in 2008/2009, all without enlargement of the aorta or valvular impairment at that time. Flow patterns, aortic diameters, WSS and distensibility were calculated in different aortic regions of interest (ROI) with dedicated software tools. Results: Indexed aortic diameters in the descending aorta (DAo), but especially in the AAo did not change in the 10-year period. Median difference 0.05 cm/m2 (95% CI: 0.01 to 0.22; P=0.06) for AAo and median difference -0.08 cm/m2 (95% CI: -0.12 to 0.01; P=0.07) for DAo. WSS values were lower in 2018/2019 at all measured levels. Aortic distensibility decreased by median 25.6% in the AAo, while stiffness increased concordantly (median +23.6%). Conclusions: After a ten years' follow-up of patients with isolated BAV disease, indexed aortic diameters did not change in this patient cohort. WSS was lower compared to values generated 10 years earlier. Possibly a drop of WSS in BAV could serve as a marker for a benign long-term course and implementation of more conservative treatment strategies.

First insights into the role of wall shear stress in the development of a distal stent graft induced new entry through computational fluid dynamics simulations.

Background: Distal stent graft induced new entry (dSINE) is an emerging complication after frozen elephant trunk (FET) procedure. The aim of this computational fluid dynamics (CFD) study was to investigate the role of wall shear stress (WSS) after the development of dSINE based on hemodynamic changes. Methods: Aortic diameter and WSS of five patients who developed a dSINE after FET procedure were retrospectively analyzed before and after the occurrence of dSINE. Patient-specific 3-dimentional surface models of the aortic lumen were reconstructed from computed tomography angiographic datasets (pre dSINE: n=5, dSINE: n=5) to perform steady-state CFD simulations with laminar blood flow and zero pressure outlet conditions. WSS was calculated at the level of the stent graft (SG), the landing zone of the SG and at a location further distal to the SG, as well as on the outer and inner curvature of the aorta from SG center to its distal end. Results: Post dSINE occurrence, median WSS increased significantly from 0.87 [interquartile range (IQR): 0.83-1.03] to 1.55 (IQR: 1.09-2.70) Pa, (P=0.043) within the SG and from 1.22 (IQR: 0.81-1.44) to 1.76 (IQR: 1.55-3.60) Pa, (P=0.043) at the landing zone of the SG. A non-significant increase from 1.22 (IQR: 0.59-3.50) to 2.58 (IQR: 1.16-3.78) Pa, (P=0.686) further downstream was observed. WSS at the outer curvature of the SG was significantly higher compared to WSS at the inner curvature for dSINE. Conclusions: Adverse hemodynamic conditions in the form of elevated WSS consist inside and at the distal end of the SG as well as at the outer curvature of the aorta, which may contribute to weakening of the aortic wall. These new findings emphasize the relevance and potential of WSS in dSINE for additional adverse events, such as aortic rupture. Further prospective studies are warranted to explore if the combination of clinical parameters with WSS might be useful to decide which patients require an urgent reintervention in terms of a SG extension.

Clinical Application of High-Frame-Rate Vector Flow Imaging in Evaluation of Carotid Atherosclerotic Stenosis.

OBJECTIVE: This study seeks to evaluate the value of the high-frame-rate vector flow imaging technique in assessing the hemodynamic changes of carotid atherosclerotic stenosis in aging people (>60 years old). METHODS: Aging patients diagnosed with carotid atherosclerotic stenosis who underwent carotid high-frame-rate vector flow imaging examination were prospectively enrolled. A Mindray Resona7s ultrasound machine equipped with high-frame-rate vector flow function was used for ultrasound evaluation. First, B mode ultrasound and color Doppler flow imaging were used to evaluate carotid stenosis. Then, the vector arrows and flow streamline detected by V Flow were analyzed and the wall shear stress values (Pa) at the carotid stenosis site were measured. All patients were divided into symptomatic and asymptomatic groups according to whether they had acute/subacute stroke or other clinical symptoms within 2 weeks before ultrasound examination. The results of digital subtraction angiography or computed tomography angiography were used as the gold standard. The stenosis rate was calcified, according to North American Symptomatic Carotid Endarterectomy Trial criteria. The diagnostic values of wall shear stress, conventional ultrasound, and the combined diagnosis in carotid atherosclerotic stenosis were compared. RESULTS: Finally, 88 patients with carotid atherosclerotic plaque were enrolled (71 males (80.7%), mean age 67.6 ± 5.4 years). The success rate of high-frame-rate vector flow imaging was 96.7% (88/91). The WSS value of symptomatic carotid stenosis (1.4 ± 0.15 Pa) was significantly higher than that of asymptomatic carotid stenosis (0.80 ± 0.08 Pa) (p < 0.05). Taking the wall shear stress value > 0.78 Pa as the diagnostic criteria for symptomatic carotid atherosclerotic plaque, the area under receiver operating characteristic curves was 0.79 with 87.1% sensitivity and 69.6% specificity. The area under receiver operating characteristic curves of the combined diagnosis (0.966) for differentiating severe carotid atherosclerotic stenosis was significantly higher than that of conventional ultrasound and WSS value, with 89.7% sensitivity and 93.2% specificity (p < 0.05). CONCLUSION: As a non-invasive imaging method, the high-frame-rate vector flow imaging technique showed potential value in the preoperative assessment of the symptomatic carotid atherosclerotic stenosis and diagnosing carotid atherosclerotic stenosis in aging patients.

Numerical investigation of different viscosity models on pulsatile blood flow of thoracic aortic aneurysm (TAA) in a patient-specific model.

Aortic aneurysm is one of the most common aortic diseases that can lead to unfortunate consequences. Numerical simulations have an important role in the prediction of the aftereffects of vascular diseases including aneurysm. In this research, numerical simulation of pulsatile blood flow is performed for a 3-dimensional patient-specific model of a thoracic aortic aneurysm (TAA). Since the choice of blood viscosity model may have a significant impact on the simulation results, the effects of four non-Newtonian models of blood viscosity namely Carreau, Casson, Herschel-Bulkley, power low, and the Newtonian model on the wall shear stress (WSS) distribution, shear rate, and oscillatory shear index (OSI) have been analyzed. Simulation results showed that all the non-Newtonian and Newtonian models generally, predict similar patterns for blood flow and shear rate. At high flow rates in the cardiac cycle, the WSS value for all the models are similar to each other except for the power-law model due to the shear thinning behavior. All models predict high values of OSI on the inner wall of the ascending aorta and broad areas of the inner wall of the aneurysm sac. However, the Newtonian model predicts the OSI less than the non-Newtonian models in some areas of the aneurysm sac. Results indicated that the Newtonian model generally can predict the hemodynamic parameters of the blood flow similar to the non-Newtonian but for more precise analysis and to predict the regions prone to rupture and atherosclerosis, choosing a proper non-Newtonian model is recommended.

Application of Biomechanical Stress from Coronary Computed Tomography Angiography in Coronary Plaque Assessment and Prediction of Adverse Cardiovascular Events / 医用生物力学

Long-term exposure to risk factors will lead to coronary atherosclerosis, which will lead to the formation and progression of coronary plaque. Early identification of high-risk plaque characteristics will help prevent plaque rupture or erosion, thus avoiding the occurrence of acute cardiovascular events. Biomechanical stress plays an important role in progression and rupture of atherosclerotic plaques. In recent years, non-invasive coronary computed tomography angiography (CCTA) computational fluid dynamics (CFD) modeling has made it possible to acquire the corresponding biomechanical stress parameters. These coronary biomechanical stress parameters, especially wall shear stress (WSS), will aid in the development of a more accurate clinical model for predicting plaque progression and major adverse cardiovascular events ( MACE ). In this review, the biomechanical stress and the role of WSS from CCTA in atherosclerosis were introduced, and the researches on the relationship between biomechanical stress from CCTA and coronary artery diseases were discussed.

Fluid-Structure Interaction Simulation of Carotid Plaque / 医用生物力学

Objective To propose a two-way fluid-structure interaction (FSI) method based on real patients with carotid artery stenosis, and analyze the hemodynamic parameters of carotid plaques with different types at the lesion as well as deformation and stress changes of the plaque itself. Methods Three-dimensional ( 3D) modeling was performed based on computed tomography angiography ( CTA) data of patients with moderate carotid artery stenosis. The carotid artery wall model and plaque model were separated, and transient fluid structure coupling calculation was performed. The situation from early stage of carotid atherosclerosis to formation of the plaque was simulated. The plaque types were divided into thickened plaques, lipid plaques, mixed plaques and calcified plaques, among which thickened plaques were regarded as non-plaque conditions for representing the thickening of vascular intima-media. The stenotic carotid arteries with different plaque types were compared and analyzed. Results The plaques with different types had little effect on the overall blood flow, but the wall shear stress of lipid plaques at the lesion was lower than that of other plaques. With thickened plaques as a control, concurrence of the plaque would inhibit artery expansion, and lipid plaques were the most obvious. Calcified plaques had the highest average plaque structure stress, while lipid plaques had the lowest average plaque structure stress. Conclusions The method proposed in this study can analyze fluid area and solid area at the same time. The results can contribute to better understanding the influence of different plaque types on carotid artery diseases.

Investigation on the differences of hemodynamics in normal common carotid, subclavian, and common femoral arteries using the vector flow technique.

Objectives: To investigate the feasibility of the vector flow imaging (V Flow) technique to measure peripheral arterial hemodynamic parameters, including wall shear stress (WSS) and turbulence index (Tur) in healthy adults, and compare the results in different arteries. Materials and methods: Fifty-two healthy adult volunteers were recruited in this study. The maximum and mean values of WSS, and the Tur values at early-systole, mid-systole, late-systole, and early diastole for total 156 normal peripheral arteries [common carotid arteries (CCA), subclavian arteries (SCA), and common femoral arteries (CFA)] were assessed using the V Flow technique. Results: The mean WSS values for CCA, SCA, and CFA were (1.66 ± 0.68) Pa, (0.62 ± 0.30) Pa, and (0.56 ± 0.27) Pa, respectively. The mean Tur values for CCA, SCA, and CFA were (0.46 ± 1.09%), (20.7 ± 9.06%), and (24.63 ± 17.66%), respectively. The CCA and SCA, as well as the CCA and CFA, showed statistically significant differences in the mean WSS and the mean Tur (P < 0.01). The mean Tur values had a negative correlation with the mean WSS; the correlation coefficient between log(Tur) and WSS is -0.69 (P < 0.05). Conclusion: V Flow technique is a simple, practical, and feasible quantitative imaging approach for assessing WSS and Tur in peripheral arteries. It has the potential to be a useful tool for evaluating atherosclerotic plaques in peripheral arteries. The results provide a new quantitative foundation for future investigations into diverse arterial hemodynamic parameters.

Distribution and regional variation of wall shear stress in the curved middle cerebral artery using four-dimensional flow magnetic resonance imaging.

Background: To investigate the distribution and regional variation of wall shear stress (WSS) in the curved middle cerebral artery (MCA) in healthy individuals using four-dimensional (4D) flow magnetic resonance imaging (MRI). Methods: A total of 44 healthy participants (18 males; mean ages: 27.16±5.69 years) were included in this cross-sectional study. The WSS parameters of mean, minimum, and maximum values, the coefficient of variation of time-averaged WSS (TAWSSCV), and the maximum values of the oscillatory shear index (OSI) were calculated and compared in the curved proximal (M1) segments. Three cross-sectional planes were selected: the location perpendicular to the beginning of the long axis of the curved M1 segment of the MCA (proximal section), the most curved M1 location (curved M1 section), and the location before the insular (M2) segment bifurcation (distal section). The WSS and OSI parameters of the proximal, curved, and distal sections of the curved M1 segment were compared, including the inner and outer curvatures of the curved M1 section. Results: Of the curved M1 segments, the curved M1 section had significantly lower minimum TAWSS values than the proximal (P=0.031) and distal sections (P=0.002), and the curved M1 section had significantly higher maximum OSI values than the distal section (P=0.001). The TAWSSCV values at the curved M1 section were significantly higher than the proximal (P=0.001) and distal sections (P<0.001). At the curved M1 section, the inner curvature showed a significantly lower minimum TAWSS (P=0.013) and higher maximum OSI values (P=0.002) than the outer curvature. Conclusions: There are distribution variation of WSS and OSI parameters at the curved M1 section of the curved MCA, and the inner curvature of the curved M1 section has the lowest WSS and highest OSI distribution. The local hemodynamic features of the curved MCA may be related to the predilection for atherosclerotic plaque development.

Hemodynamic Analysis Shows High Wall Shear Stress Is Associated with Intraoperatively Observed Thin Wall Regions of Intracranial Aneurysms.

BACKGROUND: Studying the relationship between hemodynamics and local intracranial aneurysm (IA) pathobiology can help us understand the natural history of IA. We characterized the relationship between the IA wall appearance, using intraoperative imaging, and the hemodynamics from CFD simulations. METHODS: Three-dimensional geometries of 15 IAs were constructed and used for CFD. Two-dimensional intraoperative images were subjected to wall classification using a machine learning approach, after which the wall type was mapped onto the 3D surface. IA wall regions included thick (white), normal (purple-crimson), and thin/translucent (red) regions. IA-wide and local statistical analyses were performed to assess the relationship between hemodynamics and wall type. RESULTS: Thin regions of the IA sac had significantly higher WSS, Normalized WSS, WSS Divergence and Transverse WSS, compared to both normal and thick regions. Thicker regions tended to co-locate with significantly higher RRT than thin regions. These trends were observed on a local scale as well. Regression analysis showed a significant positive correlation between WSS and thin regions and a significant negative correlation between WSSD and thick regions. CONCLUSION: Hemodynamic simulation results were associated with the intraoperatively observed IA wall type. We consistently found that elevated WSS and WSSNorm were associated with thin regions of the IA wall rather than thick and normal regions.

Effects of stent shape on focal hemodynamics in intracranial atherosclerotic stenosis: A simulation study with computational fluid dynamics modeling.

Background and aims: The shape of a stent could influence focal hemodynamics and subsequently plaque growth or in-stent restenosis in intracranial atherosclerotic stenosis (ICAS). In this preliminary study, we aim to investigate the associations between stent shapes and focal hemodynamics in ICAS, using computational fluid dynamics (CFD) simulations with manually manipulated stents of different shapes. Methods: We built an idealized artery model, and reconstructed four patient-specific models of ICAS. In each model, three variations of stent geometry (i.e., enlarged, inner-narrowed, and outer-narrowed) were developed. We performed static CFD simulation on the idealized model and three patient-specific models, and transient CFD simulation of three cardiac cycles on one patient-specific model. Pressure, wall shear stress (WSS), and low-density lipoprotein (LDL) filtration rate were quantified in the CFD models, and compared between models with an inner- or outer-narrowed stent vs. an enlarged stent. The absolute difference in each hemodynamic parameter was obtained by subtracting values from two models; a normalized difference (ND) was calculated as the ratio of the absolute difference and the value in the enlarged stent model, both area-averaged throughout the arterial wall. Results: The differences in focal pressure in models with different stent geometry were negligible (ND<1% for all cases). However, there were significant differences in the WSS and LDL filtration rate with different stent geometry, with ND >20% in a static model. Observable differences in WSS and LDL filtration rate mainly appeared in area adjacent to and immediately distal to the stent. In the transient simulation, the LDL filtration rate had milder temporal fluctuations than WSS. Conclusions: The stent geometry might influence the focal WSS and LDL filtration rate in ICAS, with negligible effect on pressure. Future studies are warranted to verify the relevance of the changes in these hemodynamic parameters in governing plaque growth and possibly in-stent restenosis in ICAS.

Derivation of vascular wall shear stress from 1000 fps high-speed angiography (HSA) velocity distributions.

Pathological changes in blood flow lead to altered hemodynamic forces, which are responsible for a number of conditions related to the remodeling and regeneration of the vasculature. More specifically, wall shear stress (WSS) has been shown to be a significant hemodynamic parameter with respect to aneurysm growth and rupture, as well as plaque activation leading to increased risk of stroke. In-vivo measurement of shear stress is difficult due to the stringent requirements on spatial resolution near the wall boundaries, as well as the deviation from the commonly assumed parabolic flow behavior at the wall. In this work, we propose an experimental method of in-vitro WSS calculations from high-temporal resolution velocity distributions, which are derived from 1000 fps high-speed angiography (HSA). The high-spatial and temporal resolution of our HSA detector makes such high-resolution velocity gradient measurements feasible. Presented here is the methodology for calculation of WSS in the imaging plane, as well as initial results for a variety of vascular geometries at physiologically realistic flow rates. Further, the effect of spatial resolution on the gradient calculation is explored using CFD-derived velocity data. Such angiographic-based analysis with HSA has the potential to provide critical hemodynamic feedback in an interventional setting, with the overarching objective of supporting clinical decision-making and improving patient outcomes.

Influence of a single hemodialysis on left ventricular energy loss and wall shear stress in patients with uremic cardiomyopathy assessed with vector flow mapping.

Background: The influence of hemodialysis (HD) on hydromechanics of the left ventricle has not been reported. This study evaluated the left ventricular summation of energy loss (EL-SUM), average energy loss (EL-AVE), and wall shear stress (WSS) before and after HD using vector flow mapping (VFM) in patients with end-stage renal disease (ESRD). Methods: We prospectively recruited 40 patients receiving long-term HD and excluded those with structural cardiac disease. Echocardiography was performed before and within 24 hours after HD. Conventional echocardiographic parameters, summation, and average energy loss (EL-SUM, EL-AVE, EL-base, EL-mid and EL-apex), and WSS in each segment were compared. Results: A total of 40 patients with uremia were recruited. After HD, left ventricular EL-AVE-total, and EL-SUM-total decreased significantly in the early diastolic [29.43 (18.76 to 46.28) vs. 17.70 (10.76 to 95.60) N/(m2·s) and 12 (6 to 17) vs. 5 (3 to 11) e-2 J; P<0.001, respectively], mid-diastolic [17.07 (10.38 to 24.35) vs. 10.29 (5.86 to 16.30) N/(m2·s) and 7 (3 to 10) vs. 4 (2 to 6) e-2 J; P<0.001, respectively], and early systolic [17.82 (12.79 to 24.77) vs.14.90 (10.23 to 19.05) N/(m2·s) P=0.011 and 8 (5 to 11) vs. 5 (4 to 8) e-2 J, P=0.002, respectively] phases. It was revealed that HD did not change EL-AVE-total and EL-SUM-total in the late diastolic and late systolic phases. The EL-AVE decreased after HD in the left ventricular (LV) basal [50.70 (24.19 to 77.92) vs. 26.00 (11.50 to 47.68) N/(m2·s); P<0.001] and mid [15.52 (8.88 to 20.90) vs. 9.47 (6.41 to 14.21) N/(m2·s); P=0.001] segments during the early diastolic phase; in the LV basal [18.64 (10.33 to 29.80) vs. 10.25 (6.98 to 19.43) N/(m2·s); P<0.001), mid (15.70 (9.93 to 23.08) vs. 9.99 (6.03 to 16.25) N/(m2·s); P<0.001), and apical [9.78 (4.06 to 15.77) vs. 4.52 (3.14 to 10.36) N/(m2·s); P=0.001) segments during the mid-diastolic phase; in the LV mid [14.34 (8.34 to 23.88) vs. 9.36 (6.48 to 17.05) N/(m2·s); P=0.013] and apex [11.25 (6.37 to 21.88) vs. 6.60 (5.33 to 12.17) N/(m2·s); P=0.016] segments during the late diastolic phase; and in the apical [10.28 (6.05 to 17.01) vs. 7.59 (3.73 to 13.20) N/(m2·s) P=0.025] segment during the early systolic phase. After HD, WSS significantly reduced in the mid-diastolic [0.51 (0.32 to 0.69) vs. 0.38 (0.30 to 0.46) Pa, P=0.001] and early systolic [0.60 (0.45 to 0.81) vs. 0.57 (0.42 to 0.68) Pa, P=0.029] phases. There was no change in WSS during the early diastolic, late diastolic, and late systolic phases. Conclusions: After HD, EL and WSS of LV decrease during the systolic and diastolic phases. The VFM can reflect the LV hemodynamics in patients undergoing HD under different fluid loads.

Effects of Pulsatile Flow Rate and Shunt Ratio in Bifurcated Distal Arteries on Hemodynamic Characteristics Involved in Two Patient-Specific Internal Carotid Artery Sidewall Aneurysms: A Numerical Study.

The pulsatile flow rate (PFR) in the cerebral artery system and shunt ratios in bifurcated arteries are two patient-specific parameters that may affect the hemodynamic characteristics in the pathobiology of cerebral aneurysms, which needs to be identified comprehensively. Accordingly, a systematic study was employed to study the effects of pulsatile flow rate (i.e., PFR-I, PFR-II, and PFR-III) and shunt ratio (i.e., 75:25 and 64:36) in bifurcated distal arteries, and transient cardiac pulsatile waveform on hemodynamic patterns in two internal carotid artery sidewall aneurysm models using computational fluid dynamics (CFD) modeling. Numerical results indicate that larger PFRs can cause higher wall shear stress (WSS) in some local regions of the aneurysmal dome that may increase the probability of small/secondary aneurysm generation than under smaller PFRs. The low WSS and relatively high oscillatory shear index (OSI) could appear under a smaller PFR, increasing the potential risk of aneurysmal sac growth and rupture. However, the variances in PFRs and bifurcated shunt ratios have rare impacts on the time-average pressure (TAP) distributions on the aneurysmal sac, although a higher PFR can contribute more to the pressure increase in the ICASA-1 dome due to the relatively stronger impingement by the redirected bloodstream than in ICASA-2. CFD simulations also show that the variances of shunt ratios in bifurcated distal arteries have rare impacts on the hemodynamic characteristics in the sacs, mainly because the bifurcated location is not close enough to the sac in present models. Furthermore, it has been found that the vortex location plays a major role in the temporal and spatial distribution of the WSS on the luminal wall, varying significantly with the cardiac period.