Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Neurol., 08 October 2025

Sec. Applied Neuroimaging

Volume 16 - 2025 | https://doi.org/10.3389/fneur.2025.1686189

CTA and DSA based computational fluid dynamics models for morphological and hemodynamic assessment of intracranial atherosclerotic stenosis

  • 1Department of Neurology, The First Affiliated Hospital of Soochow University, Suzhou, Jiangsu, China
  • 2Institute of Stroke Research, Soochow University, Suzhou, China
  • 3ArteryFlow Technology Co., Ltd., Hangzhou, China
  • 4Department of Interventional, The First Affiliated Hospital of Soochow University, Suzhou, Jiangsu, China

Background: Intracranial atherosclerotic stenosis (ICAS) is a primary cause of ischemic stroke. Accurate assessment of anatomical and hemodynamic characteristics is crucial for treatment planning, yet current clinical evaluation primarily relies on luminal stenosis.

Objective: This study aims to compare computational fluid dynamics (CFD) models based on digital subtraction angiography (DSA), computed tomography angiography (CTA) and CTA model incorporating DSA hemodynamic information (CMD) integrating DSA flow data with CTA morphological structure, evaluating their differences and consistency in ICAS assessment.

Methods: 40 ICAS patients who underwent CTA and DSA were retrospectively included. Patient-specific CFD simulations were performed using standardized boundary conditions to assess morphological data and hemodynamic parameters, including pressure ratio, wall shear stress ratio, and high shear stress areas. Statistical analyses included paired comparisons, intraclass correlation coefficients (ICC), and Bland–Altman analysis.

Results: CTA-based models demonstrated excellent consistency with DSA in anatomical measurements (ICC > 0.90). The CMD approach enhanced consistency in functional metrics, with CMD-derived PR and WSSR highly concordant with DSA results. When using CTA alone, WSSR was slightly underestimated, particularly in middle artery lesions. Subgroup analysis indicated that lesion location significantly influences flow and shear stress patterns.

Conclusion: CTA-based CFD modeling serves as a reliable non-invasive alternative to DSA for morphological ICAS assessment. The CMD method further improves the accuracy of functional evaluation by integrating flow data. These findings support the integration of anatomical imaging with hemodynamic modeling to enhance the clinical potential for stroke risk stratification.

1 Introduction

Ischemic stroke remains one of the leading causes of death and long-term disability among adults worldwide, imposing a substantial medical and economic burden on healthcare systems and society as a whole (1, 2). Intracranial atherosclerotic stenosis (ICAS) is one of the most common causes of ischemic stroke or transient ischemic attack (TIA) worldwide (3). Despite advances in medical management, patients with severe arterial stenosis (≥70%) continue to face a significantly elevated risk of stroke recurrence, particularly in the territory supplied by the stenotic vessel (4, 5). As a result, the degree of vascular stenosis has become a key determinant in clinical decision-making, particularly in guiding the need for endovascular interventions.

However, mounting evidence suggests that luminal stenosis alone may not be sufficient to accurately predict the risk of stroke recurrence (5, 6). Stroke pathogenesis is increasingly recognized as a multifactorial process that involves not only the severity of luminal narrowing but also other lesion-specific characteristics, including plaque composition and vulnerability, collateral circulation status, hemodynamic alterations, and downstream perfusion deficits (711). These additional factors offer the potential for more comprehensive and precise risk stratification.

In this context, the concept of functional assessment, originally developed for coronary artery disease using fractional flow reserve (FFR) (1214) has been adapted to cerebrovascular disease through the development of fractional flow (FF) assessments for ICAS (15, 16). Concurrently, the application of computational fluid dynamics (CFD) in ICAS has emerged as a promising tool for noninvasive (17), patient-specific hemodynamic evaluation based on anatomical imaging modalities such as digital subtraction angiography (DSA) (18, 19), computed tomography angiography (CTA) (9, 2022) and magnetic resonance angiography (MRA) (2325).

Although many CFD studies have evaluated hemodynamic features in atherosclerotic lesions, most are limited to a single imaging modality, restricting cross-modal comparisons. As DSA and CTA are both widely used in clinical practice, a systematic comparison of CFD-derived morphological and hemodynamic parameters is necessary to understand their respective strengths. This study aimed to compare patient-specific CFD models based on DSA and CTA in ICAS, focusing on key parameters such as pressure gradient, wall shear stress, and flow patterns, to evaluate the consistency and reliability of functional assessments across modalities.

2 Materials and methods

2.1 Subjects and clinical data

A systematic review was conducted on the clinical data of 40 patients with ischemic stroke due to intracranial atherosclerotic stenosis, admitted to the First Affiliated Hospital of Soochow University between June 2021 and September 2022. The inclusion criteria were: (1) Computed tomography angiography (CTA) was performed within 1 week prior to Digital subtraction angiography (DSA); (2) Cerebral infarction caused by luminal stenosis of intracranial arteries, including the intracranial carotid arteries (ICA), the middle cerebral arteries (MCA), and the intracranial segments of the vertebrobasilar arteries (VA). Exclusion criteria included: (1) Ischemic stroke caused by factors other than atherosclerosis, such as Moyamoya disease, vasculitis, dissection, tumors, or cardioembolism; (2) Patients who had previously undergone interventional and/or surgical treatment of intracranial and extracranial arteries; (3) Cases of complete arterial occlusion (Figure 1). This study has received approval from the local Institutional Review Board (IRB) of the participating center. The data have been anonymized, and the requirement for informed consent has been waived.

Figure 1
Flowchart detailing a study involving patient screening and enrollment (n=40) for CTA within one week before DSA. It includes imaging acquisition, CFD analysis, boundary conditions, solution strategy, mesh generation, and morphological and hemodynamic parameter extraction. Statistical analysis follows for morphological and hemodynamic differences, consistency, and stenosis subgroup analysis.

Figure 1. Flowchart of the present study. Flowchart of the present study. Involving the screening and enrollment of 40 patients, covering aspects such as image acquisition, CFD analysis, boundary condition setting, solution strategies, mesh generation, and extraction of morphological and hemodynamic parameters. Statistical analysis is subsequently conducted to explore morphological and hemodynamic differences, consistency, and subgroup analysis of stenosis.

2.2 Computed tomography angiography

The CT scanner selected is the Revolution 256 spiral CT machine from GE Company in the United States. CTA examination was conducted by injecting 40 mL of iodinated contrast agent and 50 mL of saline, both injected at a rate of 5 mL/s. Arterial and venous phase images were obtained with scanning parameters of 0.625 mm layer thickness, 512 × 512 matrix, and a scan range from the aortic arch to the top of the skull. The resolution of CTA is approximately 0.488 mm for each direction.

2.3 Digital subtraction angiography operation

We used Siemens Artis Zeego or Toshiba INFX-8000 V DSA machine as interventional imaging equipment. The right femoral artery was punctured using Seldinger modified technique, and a 4F H1 arterial sheath was inserted. Then, 2,000 U heparin was injected intravenously. A 5F pigtail catheter and a single curved catheter were used to perform angiography of the aortic arch, bilateral carotid arteries, and bilateral vertebral arteries, and contrast agent (Ioversol 320 mg/mL) was injected. The angiography parameters were as follows: injection rate 4-5 mL/s, injection dose 6-8 mL, shooting speed 6 frames/s. The average resolution of CTA is about 0.37 mm (0.35 mm ~ 0.39 mm).

2.4 Model reconstruction

The DICOM files of DSA and CTA of head vessels analyzed in this paper are all from the Neusoft Imaging Workstation of the First Affiliated Hospital of Soochow University. First, the image data in DICOM format was imported into MIMICS, Version 20.0, for image segmentation and reconstruction (26). Removing side branches (radius < 50% of parent vessel) has negligible impact on PR measurement accuracy (27), but unknown effect on other parameters, so retained stenotic peripheral branches where possible.

2.5 Hemodynamic simulation of CFD

2.5.1 Mesh generation and Independence verification

Polyhedral volume meshs with prism layers near the wall were generated for each model. To ensure the accuracy of the simulation results was independent of mesh density, a mesh independence study was conducted on a representative model with 62% MCA diameter stenosis. The translesional pressure ratio (PR) was monitored as the key metric. The solution was deemed mesh-independent when the PR value varied by less than 2% upon further refinement. Based on this study, maximum cell size of 0.15 mm was adopted for all simulations, resulting in approximately 1 to 2 million cells per model.

2.5.2 Boundary conditions

For models derived from DSA images, the patient-specific time-averaged flow velocity at the inlet boundary was determined using the TIMI frame count method (27) and converted to a flow rate. A constant velocity-inlet boundary condition was prescribed at the model inlet. For models derived from CTA images, the patient-specific time-averaged flow velocity at the inlet boundary was obtained from Transcranial Doppler (TCD) ultrasound measurements of the relevant artery (28). TCD examinations were conducted within 48 h of patient admission by two certified neurosonographers, utilizing a transcranial Doppler device (EMS-9 PB; Shenzhen Delica Medical Equipment Co., Ltd., Shenzhen, China). The mean flow velocities of the middle cerebral artery, intracranial segments of the internal carotid artery, and the proximal portion of the vertebrobasilar artery were measured, with three independent measurements averaged to minimize operational error. The resulting patient-specific mean flow velocity values were employed as inlet boundary conditions for the model. All outlets were modeled as split outlet boundary conditions with flow rates proportional to the cube of their diameters. The CTA model incorporating DSA hemodynamic information (CMD) used the same geometry and outlet conditions as the CTA group, but with the inlet flow rate replaced by the DSA-derived TIMI frame count value.

2.5.3 Solution strategy

Patient-specific models derived from DSA and CTA were simulated using the finite volume method (FVM) in STAR-CCM + v.20.02 (Siemens Digital Industries Software, Plano, TX, USA). Blood was modeled as an incompressible Newtonian fluid (density = 1,060 kg/m3, viscosity = 0.0035 Pa·s) with rigid, no-slip vessel walls (9, 20, 21). Steady-state simulations under mean-flow conditions were performed to assess translesional pressure drop, neglecting pulsatility effects. The SIMPLE algorithm coupled pressure and velocity, with second-order spatial discretization. Convergence was achieved when normalized residuals for continuity and momentum fell below 10−5. The blood flow simulation method similar to DSA has recently been validated by comparing the simulated values with pressure wire measurements in the literature (19).

2.5.4 Definition of relevant parameters

The comparative evaluation of DSA and CTA for stenotic vessels focuses on two aspects. Firstly, the morphological evaluation mainly includes six anatomical parameters, percentage of diameter stenosis (DS%) (29), percentage of area stenosis (AS%), narrow proximal diameter, narrow distal diameter, minimum lumen diameter (30) and stenosis length. Secondly, the hemodynamic evaluation includes four parameters, Narrow flow rate, translesional pressure ratio (PR): defined as the ratio of mean pressure distal to the stenosis to that proximal to the stenosis (PR = P_distal / P_proximal) (9, 22), wall shear stress ratio (WSSR): defined as the ratio of mean wall shear stress within the stenotic throat to the mean wall shear stress in the proximal reference segment (WSSR = WSS_stenosis / WSS_proximal) (9, 23) and the area of high wall shear stress (High-WSS area): defined as the luminal surface area exposed to wall shear stress values greater than 3.0 Pa, a threshold widely used to indicate abnormally elevated shear stress in intracranial and carotid arteries (20).

2.6 Statistical analysis

DSA is recognized as the gold standard for detecting vascular stenosis (30). Consequently, in this study, the DSA-based model was established as the “standard model,” and the CTA-based model along with the CMD were compared against the “standard model DM.” A comprehensive statistical analysis was conducted to evaluate the morphological and hemodynamic consistency among DSA, CTA, and CMD. For continuous variables, the Shapiro–Wilk test was employed to assess normality. Depending on the data distribution, paired t-tests or Wilcoxon signed-rank tests were used for comparisons between two groups, while the Kruskal-Wallis test was applied for comparisons among three groups with non-normally distributed data. To assess the agreement between imaging modalities, Bland–Altman plots were utilized to evaluate consistency, and the absolute agreement ICC were calculated using a two-way random effects model. All statistical tests were two-sided, with a p-value < 0.05 considered statistically significant.

3 Results

3.1 Baseline characteristics

A total of 40 patients with intracranial arterial stenosis were included. The most frequently affected vascular territory was the MCA, accounting for 77.5% (31/40) of cases, followed by the ICA (15%) and the vertebrobasilar artery VA (7.5%). The mean fasting blood glucose level was 6.45 ± 1.93 mmol/L, with a mean triglyceride level of 1.42 ± 0.76 mmol/L and a mean HbA1c of 6.38 ± 1.07%. The levels of high-density lipoprotein (HDL) and low-density lipoprotein cholesterol (LDL-C) were 1.12 ± 0.31 mmol/L and 2.40 ± 0.77 mmol/L, respectively. Regarding neurological status, the median modified Rankin Scale (mRS) score prior to admission was 1 (interquartile range [IQR]: 0–2), and the median National Institutes of Health Stroke Scale (NIHSS) score at admission was 0 (IQR: 0–2). In terms of medical history, 29 patients (72.5%) had a history of hypertension, 16 (40.0%) had dyslipidemia, and 9 (22.5%) had diabetes. Six patients (15.0%) were current or former smokers, 5 (12.5%) had a history of ischemic heart disease, and 23 patients (57.5%) had a prior history of stroke or TIA.

3.2 Morphology data analysis of DSA and CTA

3.2.1 Analysis of morphological assessment differences

To assess the concordance between DSA and CTA in morphological measurements (Table 1), we conducted a paired statistical comparison of key vascular indices, including diameter stenosis, area stenosis, proximal and distal reference diameters of narrow, minimal lumen diameter, and stenosis length. Across all parameters, the measurements from DSA and CTA exhibited no statistically significant differences. Specifically, the median diameter stenosis rate in the DSA group was 59.95% [IQR: 49.35–76.50], whereas in the CTA group, it was 58.95% ± 19.84 (p = 0.618). For area stenosis, the DSA yielded 82.98% [71.52–94.00], compared to 83.56% [71.97–92.94] from CTA (p = 0.510). Similarly, no significant discrepancies were observed in proximal diameter (DSA: 2.29 [1.90–2.74] mm; CTA: 2.25 [1.86–2.85] mm; p = 0.154) or distal diameter (DSA: 2.23 ± 0.72 mm; CTA: 2.26 ± 0.73 mm; p = 0.068). Likewise, minimal lumen diameter (DSA: 1.02 ± 0.52 mm; CTA: 1.05 ± 0.55 mm; p = 0.509) and stenosis length (DSA: 4.28 ± 2.18 mm; CTA: 4.43 ± 2.27 mm; p = 0.164) did not demonstrate significant differences. These findings indicate a high level of agreement between DSA and CTA in evaluating the morphological aspects of intracranial stenosis, suggesting that CTA can serve as a reliable non-invasive alternative for morphological assessment in clinical practice.

Table 1
www.frontiersin.org

Table 1. Morphological differences analysis between DSA and CTA.

3.2.2 Morphological assessment consistency analysis

To evaluate the consistency between DSA and CTA in assessing vascular morphology, we performed Bland–Altman and ICC analyses across six paired anatomical parameters. Bland–Altman plots demonstrated that for all parameters (Figure 2A), the mean differences between DSA and CTA measurements were small, with most values falling within the 95% limits of agreement (LoA), indicating no significant systematic bias. Specifically, the mean differences (DSA – CTA) were: Percentage of Diameter Stenosis: 0.17, with LoA ranging from −10.71 to 11.05; Percentage of Area Stenosis: 0.25, with LoA from −9.57 to 10.06; Proximal Diameter: 0.04 mm, LoA from −0.62 to 0.71 mm; Distal Diameter: −0.03 mm, LoA from −0.60 to 0.54 mm; Minimum Lumen Diameter: −0.03 mm, LoA from −0.52 to 0.46 mm; Stenosis Length: −0.15 mm, LoA from −2.67 to 2.36 mm. These results suggest that individual-level variation between DSA and CTA was acceptable and showed no systematic measurement shift. Furthermore, intraclass correlation analysis using a two-way random-effects model revealed excellent agreement between DSA and CTA (Figure 2B): Percentage of Diameter Stenosis: ICC = 0.996; Percentage of Area Stenosis: ICC = 0.995; Proximal Diameter: ICC = 0.976; Distal Diameter: ICC = 0.991; Minimum Lumen Diameter: ICC = 0.902; Stenosis Length: ICC = 0.985; All ICC values exceeded 0.90, indicating excellent consistency between DSA and CTA across all morphological measurements. These findings confirm that CTA-based vascular modeling can serve as a reliable alternative to DSA for quantifying anatomical features relevant to intracranial arterial stenosis.

Figure 2
Panel A displays Bland-Altman plots for six parameters: Diameter Stenosis, Area Stenosis, Narrow Proximal Diameter, Narrow Distal Diameter, Minimal Lumen Diameter, and Stenosis Length. Each plot compares differences (DSA-CTA) against the mean, with mean difference and limits of agreement. Panel B is a bar chart showing the Intraclass Correlation Coefficient (ICC) for each parameter, with a reference line for good reliability at 0.75.

Figure 2. Consistency Analysis of Morphological Parameters between DSA and CTA. (A) Bland-Altman plots were employed to compare six paired anatomical measurements obtained from digital subtraction angiography (DSA) and computed tomography angiography (CTA), encompassing percent diameter stenosis, percent area stenosis, proximal diameter, distal diameter, minimal lumen diameter, and stenosis length. The blue dashed lines represent the mean differences, while the red and green dashed lines denote the Upper and Lower 95% limits of agreement (LoA), respectively. (B) Intraclass correlation coefficient (ICC) analysis results for the same six parameters.

3.2.3 Morphological subgroup analysis

Subgroup analysis was conducted to compare morphological characteristics of intracranial atherosclerotic lesions located in the MCA (n = 31), ICA (n = 6), and VA (n = 3), using both DSA and CTA (Table 2). For the percentage of diameter stenosis, both DSA and CTA showed no statistically significant differences across vascular territories (DSA: MCA 60.30 [47.71–77.62], ICA 56.27 ± 14.81, VA 65.71 ± 7.51; p = 0.655; CTA: MCA 58.97 ± 21.47, ICA 56.10 ± 15.51, VA 64.45 ± 9.52; p = 0.845). Similarly, the percentage of area stenosis was not significantly different between regions in either modality (DSA: MCA 84.20 [69.92–95.03], ICA 79.07 ± 10.65, VA 86.44 ± 5.39; p = 0.725; CTA: MCA 84.29 [69.82–93.34], ICA 77.74 ± 10.25, VA 87.37 ± 6.69; p = 0.607). In contrast, proximal vessel diameter showed significant differences across all territories for both imaging methods (DSA: MCA 2.19 ± 0.48 mm, ICA 3.26 ± 0.87 mm, VA 3.95 ± 0.89 mm; p < 0.01; CTA: MCA 2.15 ± 0.57 mm, ICA 3.26 ± 0.85 mm, VA 3.89 ± 0.67 mm; p < 0.01). A similar pattern was observed for the distal diameter, which was significantly larger in the ICA and VA compared to the MCA (DSA: MCA 2.00 ± 0.55 mm, ICA 2.78 ± 0.74 mm, VA 3.46 ± 0.49 mm; p < 0.01; CTA: MCA 2.04 ± 0.57 mm, ICA 2.73 ± 0.68 mm, VA 3.56 ± 0.54 mm; p < 0.01). Although the minimum lumen diameter tended to be larger in ICA and VA lesions, the differences did not reach statistical significance in either modality (DSA: MCA 0.93 ± 0.48 mm, ICA 1.35 ± 0.59 mm, VA 1.33 ± 0.62 mm; p = 0.105; CTA: MCA 0.94 ± 0.51 mm, ICA 1.39 ± 0.60 mm, VA 1.46 ± 0.57 mm; p = 0.071). Overall, both DSA and CTA consistently revealed that while stenosis severity (percent narrowing) remained comparable among vascular segments, vessel caliber—particularly proximal and distal diameters—differed significantly among MCA, ICA, and VA lesions. This suggests that anatomic location should be considered when evaluating lesion morphology or planning interventional strategies.

Table 2
www.frontiersin.org

Table 2. Morphological subgroup analysis of DSA and CTA.

3.3 Hemodynamic comparison of DSA and CTA

3.3.1 Hemodynamic assessment differential analysis

To comprehensively compare hemodynamic parameters obtained from DSA, CTA, and CMD, pairwise comparisons were performed for flow rate, pressure ratio (PR), wall shear stress ratio (WSSR), and High WSS area (Figure 3A). For flow rate, the median (IQR) values were 1.46 (0.88–1.99) ml/s for DSA and 1.44 (1.01–2.15) ml/s for CTA, with no significant difference (p = 0.116). For PR, the median values were 0.71 (0.42–0.89) for DSA, 0.74 (0.46–0.89) for CTA, and 0.73 (0.41–0.93) for CMD. Significant differences were observed between DSA and CTA (p = 0.041), whereas DSA vs. CMD (p = 0.068) and CTA vs. CMD (p = 0.452) showed no significant difference. For WSSR, the median values were 22.49 (6.11–36.10) for DSA, 18.02 (3.58–33.07) for CTA, and 21.43 (7.25–36.94) for CMD. Significant differences were noted between DSA and CTA (p < 0.01) and between CTA and CMD (p = 0.002), while no significant difference was found between DSA and CMD (p = 0.109). This suggests CMD-derived WSSR is closer to DSA estimates than CTA-derived values. For High WSS area, the values were 10.46 (4.72–14.98) for DSA, 10.15 (5.97–16.37) for CTA, and 11.72 (5.85–16.02) for CMD. Differences were statistically significant between DSA and CTA (p = 0.014) and between DSA and CMD (p = 0.019), but not between CTA and CMD (p = 0.735). Overall, hemodynamic measurements from DSA, CTA, and CMD demonstrated good consistency, although WSSR and high WSS area exhibited method-dependent variability.

Figure 3
Panel A features bar charts comparing flow rate, PR, WSSR, and high-WSS area across DSA, CTA, and CMD with significance marked by asterisks. Panel B shows scatter plots and Bland-Altman plots for PR and WSSR comparisons with mean difference lines. Panel C displays a color-coded bar chart of intraclass correlation coefficients for each parameter comparison, ranging from zero to one.

Figure 3. Comparative Analysis and Consistency Evaluation of Hemodynamic Parameters Derived from DSA, CTA, and CMD Models. (A) Dual Comparisons of Flow, Pressure Ratio (PR), Wall Shear Stress Ratio (WSSR), and High WSS Areas Among Digital Subtraction Angiography (DSA), Computed Tomography Angiography (CTA), and CMD. Box plots illustrate the median and interquartile range. (B) Bland-Altman Plots and Intraclass Correlation Coefficient (ICC) Analysis for Flow, PR, and WSSR Between Each Imaging Mode Pair (DSA vs. CTA, DSA vs. CMD, CTA vs. CMD). The blue dashed line represents the mean difference, while the red and green dashed lines indicate the Upper and Lower 95% Limits of Agreement (LoA), respectively. (C) ICC Analysis Results for Hemodynamic Parameters. Significance levels are marked as *p < 0.05, **p < 0.01, ***p < 0.001, ns: not significant.

3.3.2 Hemodynamic assessment consistency analysis

To evaluate the measurement consistency of hemodynamic parameters obtained by DSA, CTA, and CMD, we conducted intra-group correlation coefficient (ICC) analysis and Bland–Altman consistency assessment for four key indicators PR, WSSR, and high WSS area (Figure 3B). In terms of traffic, the results show that DSA and CTA have excellent consistency, with an ICC value of 0.954. The Bland–Altman plot shows the mean difference is −0.080 mL/s and the 95% consistency limit (LoA) range is −0.667 to 0.506 mL/s. For PR, all three methods showed strong commonality. The ICC values are 0.973 (DSA vs. CTA), 0.992 (DSA vs. CMD), and 0.964 (CTA vs. CMD), respectively. Bland–Altman analysis showed that the mean difference was very small, which was −0.018, −0.010 and +0.008, respectively, and the consistency boundary interval was narrow (such as the LoA of DSA and CMD was −0.072 to 0.052), indicating that the methods were highly consistent and had small deviations. In terms of WSSR, the ICC values also show excellent reliability: 0.983 (DSA vs. CTA), 0.989 (DSA vs. CMD), and 0.988 (CTA vs. CMD). The corresponding mean differences were 2.67, −0.77, and −3.45, with the consistency boundary ranging from about −6.08 to 11.41 for DSA and CTA to −11.99 to 5.09 for CTA and CMD. These results show that the WSSR estimates provided by CMD are more consistent with DSA than CTA. The intraclass correlation coefficient (ICC) demonstrated excellent agreement between measurement methods across all parameters, with values ranging from 0.902 to 0.996. The highest consistency was observed for diameter stenosis (ICC = 0.996) and area stenosis (ICC = 0.995), while minimum lumen diameter showed the lowest, yet still excellent, agreement (ICC = 0.902). All ICC values exceeded the 0.75 threshold for good reliability, and most surpassed the 0.90 threshold for excellent reliability (Figure 3C).

3.3.3 Hemodynamic subgroup analysis

To investigate the impact of the lesion site on hemodynamic characteristics, the study continued to categorize the stenosis based on its anatomical location into MCA (n = 31), VA (n = 3), and ICA (n = 6). Comparative analyses were conducted across these anatomical subgroups (Table 3). In terms of flow rate, both DSA and CTA data revealed significant differences among the three vascular regions (p values of 0.006 and 0.001, respectively), with the ICA lesion group exhibiting the highest mean flow rate (DSA: 2.86 ± 1.04 mL/s; CTA: 2.97 ± 1.09 mL/s), followed by the VA and MCA groups. Regarding PR, although CMD-derived PR tended to be lower in MCA lesions compared to VA and ICA, no statistically significant differences were observed among the three methods (DSA: p = 0.175; CTA: p = 0.205; CMD: p = 0.184). For WSSR, CTA identified significant intergroup differences (p = 0.040), with the VA lesion group showing the highest mean WSSR (36.06 ± 3.53), while the MCA and ICA groups exhibited more dispersed distributions. Although comparisons of WSSR based on DSA and CMD did not reach statistical significance (p > 0.05), they exhibited similar trends. Significant differences in high WSS regions were observed across anatomical locations in all three imaging modalities (DSA: p < 0.01; CTA: p < 0.01; CMD: p = 0.049). VA and ICA lesions consistently demonstrated larger high WSS areas (CTA: VA = 27.59 ± 8.03 mm2 vs. MCA = 9.59 ± 6.19 mm2), indicating more extensive local hemodynamic disturbances compared to MCA stenosis. In summary, these findings demonstrate that the stenosis location has a measurable influence on hemodynamic parameters, particularly flow magnitude and WSS distribution. Lesions in proximal arteries such as ICA and VA are often associated with higher flow rates and larger high WSS regions, which may reflect differences in vessel caliber, collateral circulation supply, and compensatory hemodynamic adjustments. As representatives, we present the hemodynamic analysis visualization of PR and WSSR using MCA in Figure 4.

Table 3
www.frontiersin.org

Table 3. Hemodynamic subgroup analysis of DSA and CTA.

Figure 4
Three pairs of 3D vascular models show simulations of pressure ratio (PR) and wall shear stress ratio (WSSR) in different imaging techniques: DSA, CTA, and CMD. Each pair includes a color scale corresponding to values on the side, illustrating variations in pressure and stress across the vascular structures.

Figure 4. Representative visualization of hemodynamic analysis in cases of MCA stenosis. Based on the results of patient-specific computational fluid dynamics (CFD) analysis of typical middle cerebral artery (MCA) lesions, the spatial distributions of pressure ratio (PR) and wall shear stress ratio (WSSR) are presented. These data are derived from digital subtraction angiography (DSA), computed tomography angiography (CTA), and their combined model (CMD). The color maps along the vessel surface indicate the intensity of each parameter, highlighting the specific differences in hemodynamic assessment among the various imaging modalities.

4 Discussion

This study systematically compared patient-specific CFD models constructed based on DSA, CTA, and CMD that integrates DSA blood flow data with CTA geometric morphology, in the assessment of morphological and hemodynamic characteristics of ICAS. The results demonstrated high consistency among the modalities in both anatomical and functional parameters, further validating the potential of CTA-based CFD as a non-invasive alternative to DSA for clinical evaluation.

The patient cohort included in this study predominantly featured intracranial arterial lesions and commonly exhibited traditional vascular risk factors such as hypertension and dyslipidemia, indicating that our sample population is representative within the context of ICAS research. DSA, long regarded as the “gold standard” for clinical diagnosis of vascular stenosis, has faced obstacles in its widespread adoption due to its high cost, invasive nature, and an estimated 1–2% incidence of procedure-related morbidity (31). Among the current non-invasive imaging methods, CTA has gained widespread adoption and application due to its high discriminatory power in identifying severe stenosis and occlusions (3235). At the morphological level, CTA and DSA exhibit highly concordant measurements of stenosis severity, reference diameters, and lesion lengths, with intraclass correlation coefficients (ICC > 0.90) demonstrating strong agreement. Additionally, the Bland–Altman plots indicate minimal measurement bias. These findings substantiate CTA as a reliable imaging modality for assessing the luminal morphology of ICAS, consistent with prior studies that highlight its high diagnostic accuracy for moderate to severe intracranial stenosis (9, 36, 37). However, subtle discrepancies have been observed in certain cases. CTA tends to slightly overestimate the distal reference diameter or the area of stenosis, a phenomenon potentially attributable to partial volume effects, calcification artifacts, or segmentation threshold variability, particularly in small-caliber vessels such as the MCA (35, 38, 39). Although these differences did not achieve statistical significance at the group level, they may influence therapeutic decision-making for critical lesions, thereby underscoring the critical importance of cautious interpretation of CTA measurements in clinical practice.

Compared to other hemodynamic assessment methodologies, the establishment of CFD necessitates the input of boundary blood flow information, which is often challenging to obtain under non-invasive conditions. The circulatory system within the human vascular network operates dynamically, influenced by a multitude of factors including instantaneous blood pressure, blood viscosity, collateral circulation, and vascular wall elasticity (40, 41). The assumptions regarding the boundary conditions and hemodynamic parameters within the CFD model significantly influence the accuracy of simulated blood flow. Wang et al. employed the TIMI frame count method, integrating DSA imaging with a three-dimensional model to compute the PR at the distal and proximal sites of intracranial stenosis, as well as the pressure wire measurements, achieving a strong correlation (r = 0.908, p < 0.001) (19). Liu et al. enhanced the consistency between CFD and pressure wire measurements by iteratively adjusting the inlet pressure boundary conditions after each CFD computation (22). However, in previous studies, CFD models based on CTA (9, 20, 21) or MRA (23) typically employed fixed flow rate information directly, which may introduce computational errors. Chen et al. integrated the stenotic area and blood flow velocity obtained from transcranial Doppler with MRA data and found that the WSSR exhibited a significant correlation with the degree of stenosis (25). This study utilized the TIMI method or neck vascular ultrasound to obtain specific flow rates at corresponding stenotic sites (42).

Previous studies have demonstrated that the morphological assessment of intracranial stenosis only partially reflects the physiological severity in patients with distal cerebral flow restriction (43, 44), whereas the PR quantifies the pressure gradient across the stenosis. In the present study, excellent concordance was observed in flow and PR measurements among DSA, CTA, and CFD simulations, with an intraclass correlation coefficient (ICC > 0.95). This finding substantiates the high reliability of CTA-based simulations when boundary conditions derived from DSA are incorporated. Previous studies have established that WSS, as a hemodynamic parameter, is associated with the degree of vascular stenosis and the pressure gradient. Elevated WSS (45, 46), on one hand, promotes the formation of atherosclerotic plaques, while on the other hand, it induces plaque rupture and thrombus formation. This study incorporated fixed quantitative WSSR and high WSS regions for comparative analysis. In terms of the WSSR, the results from CMD often exhibit closer alignment with DSA measurements compared to those obtained using CTA alone. Nonetheless, WSSR demonstrates the highest variability associated with imaging modalities, particularly in lesions of the MCA. This variability can be attributed to the complex geometry and flow dynamics near the proximal stenotic segment, especially around bifurcations, which complicates the definition of prestenotic WSS. As a ratio-based metric, WSSR is highly sensitive to minor variations in the denominator, thereby magnifying measurement errors. In contrast, high WSS regions, serving as a more stable and spatially integrated parameter, maintain consistency across different imaging modalities and vascular regions. Notably, previous studies have established a correlation between larger high WSS regions and plaque regression or vulnerability, reinforcing their significance as functional biomarkers of disease progression (24, 47).

Subgroup analysis demonstrates that the location of the lesion significantly influences both morphological and hemodynamic characteristics. Compared to MCA stenosis, lesions in the ICA and VA exhibit larger proximal and distal vessel diameters, higher blood flow rates, and more extensive regions of high WSS. These differences likely reflect the anatomical distinctions between proximal and distal cerebrovascular segments and the compensatory hemodynamic alterations that ensue (19, 21, 48). Interestingly, WSSR exhibited significant inter-regional variability exclusively in CTA-based models, suggesting that CTA-derived WSSR may be more susceptible to geometric distortions or segmentation variability. The feasibility of CMD as a hybrid approach warrants special mention. By integrating CTA’s high spatial resolution with the accurate flow dynamics derived from DSA, CMD offers a promising compromise between accessibility and precision. In this study, PR and WSSR derived from CMD demonstrated high concordance with DSA results, corroborating previous literature that integrating reliable flow data into non-invasive geometric models can significantly enhance the accuracy of CFD (21, 49).

Several limitations are inherent in this study: Initially, the small sample size necessitates the inclusion of a larger patient cohort to validate the accuracy of the findings, despite the preliminary explorations offering feasibility for subsequent in-depth investigations. Secondly, not all patients underwent all imaging examinations on the same day, implying that the luminal structure of the diseased vessels and the associated hemodynamic information may have varied temporally. Additionally, when employing measurements such as TIMI frame count or TCD to quantify patient-specific hemodynamic data, uncertainties arising from operator variability or imaging quality can introduce measurement errors. Lastly, akin to all previous CFD models, the simplified blood flow simulation model utilized in this study overlooks confounding factors such as inlet blood pulsatility and collateral circulation, which to some extent compromise the precision of vascular CFD simulation results. Future research should address these factors. It is crucial to carefully consider these limitations when interpreting the study outcomes.

5 Conclusion

This study demonstrates that the CTA-based CFD model exhibits a high degree of consistency with DSA in the anatomical and hemodynamic assessment of intracranial atherosclerotic stenosis. CTA model incorporating DSA hemodynamic information significantly enhances the accuracy of functional indices such as PR and WSSR. Although CTA alone may underestimate WSSR in the evaluation of middle cerebral artery lesions, it remains reliable in assessing the degree of stenosis and regions of high WSS. These findings support the use of CTA-based CFD technology as a feasible and non-invasive alternative to DSA, highlighting the potential of CMD in improving the precision of functional assessments.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors on request. Requests to access these datasets should be directed to Xulong Yin, eXhsMTU2NjIxOTA4ODZAMTYzLmNvbQ==.

Ethics statement

The studies involving humans were approved by the First Affiliated Hospital of Soochow University (Approval number: 2021422). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

RY: Resources, Writing – original draft. XY: Formal analysis, Software, Writing – original draft. GL: Data curation, Writing – original draft. JX: Supervision, Writing – original draft. QF: Methodology, Writing – original draft. HW: Funding acquisition, Writing – review & editing. BL: Conceptualization, Funding acquisition, Writing – original draft.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by Suzhou Science and Technology (SSD2024090), the Key Discipline Construction Project of Medicine in Jiangsu Province (No. 04999913) and BoXi Discipline-Class II Discipline-Neurology (No.07040120250011), Zhejiang Provincial “Pioneer” and “Leading Goose” Science and Technology Program (2024C03095).

Conflict of interest

GL and JX was employed by the ArteryFlow Technology Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Abbreviations

CFD, Computational fluid dynamics; DSA, Digital subtraction angiography; CTA, Computed tomographic angiography; DM, DSA based models; CMD, CTA based models through DSA blood flow boundary conditions; ICA, Internal carotid artery; MCA, Middle cerebral artery; TCD, Transcranial doppler; AS, Atherosclerotic stenosis; Ds, Smallest diameter; DS%, Percent diameter stenosis; WSS, Wall shear stress; MWSS, Maximum wall shear stress; WSSR, Translesional wall shear stress ratio High-WSS region The area of high wall shear stress; PR, Pressure ratio.

References

1. Wu, S, Wu, B, Liu, M, Chen, Z, Wang, W, Anderson, CS, et al. Stroke in China: advances and challenges in epidemiology, prevention, and management. Lancet Neurol. (2019) 18:394–405. doi: 10.1016/S1474-4422(18)30500-3

PubMed Abstract | Crossref Full Text | Google Scholar

2. Liu, L, Wang, D, Wong, KS, and Wang, Y. Stroke and stroke care in China: huge burden, significant workload, and a national priority. Stroke. (2011) 42:3651–4. doi: 10.1161/STROKEAHA.111.635755

PubMed Abstract | Crossref Full Text | Google Scholar

3. Gutierrez, J, Turan, TN, Hoh, BL, and Chimowitz, MI. Intracranial atherosclerotic stenosis: risk factors, diagnosis, and treatment. Lancet Neurol. (2022) 21:355–68. doi: 10.1016/S1474-4422(21)00376-8

PubMed Abstract | Crossref Full Text | Google Scholar

4. Chimowitz, MI, Kokkinos, J, Strong, J, Brown, MB, Levine, SR, Silliman, S, et al. The warfarin-aspirin symptomatic intracranial disease study. Neurology. (1995) 45:1488–93. doi: 10.1212/WNL.45.8.1488

PubMed Abstract | Crossref Full Text | Google Scholar

5. Chimowitz, MI, Lynn, MJ, Derdeyn, CP, Turan, TN, Fiorella, D, Lane, BF, et al. Stenting versus aggressive medical therapy for intracranial arterial stenosis. N Engl J Med. (2011) 365:993–1003. doi: 10.1056/NEJMoa1105335

PubMed Abstract | Crossref Full Text | Google Scholar

6. Kasner, SE, Chimowitz, MI, Lynn, MJ, Howlett-Smith, H, Stern, BJ, Hertzberg, VS, et al. Predictors of ischemic stroke in the territory of a symptomatic intracranial arterial stenosis. Circulation. (2006) 113:555–63. doi: 10.1161/CIRCULATIONAHA.105.578229

PubMed Abstract | Crossref Full Text | Google Scholar

7. Tao, L, Li, XQ, Hou, XW, Yang, BQ, Xia, C, Ntaios, G, et al. Intracranial atherosclerotic plaque as a potential cause of embolic stroke of undetermined source. J Am Coll Cardiol. (2021) 77:680–91. doi: 10.1016/j.jacc.2020.12.015

PubMed Abstract | Crossref Full Text | Google Scholar

8. Liebeskind, DS, Cotsonis, GA, Saver, JL, Lynn, MJ, Turan, TN, Cloft, HJ, et al. Collaterals dramatically alter stroke risk in intracranial atherosclerosis. Ann Neurol. (2011) 69:963–74. doi: 10.1002/ana.22354

PubMed Abstract | Crossref Full Text | Google Scholar

9. Leng, X, Lan, L, Ip, HL, Abrigo, J, Scalzo, F, Liu, H, et al. Hemodynamics and stroke risk in intracranial atherosclerotic disease. Ann Neurol. (2019) 85:752–64. doi: 10.1002/ana.25456

PubMed Abstract | Crossref Full Text | Google Scholar

10. Dakay, K, and Yaghi, S. Symptomatic intracranial atherosclerosis with impaired distal perfusion: a case study. Stroke. (2018) 49:e10–3. doi: 10.1161/STROKEAHA.117.019173

PubMed Abstract | Crossref Full Text | Google Scholar

11. Leng, X, Wong, KS, and Liebeskind, DS. Evaluating intracranial atherosclerosis rather than intracranial stenosis. Stroke. (2014) 45:645–51. doi: 10.1161/STROKEAHA.113.002491

PubMed Abstract | Crossref Full Text | Google Scholar

12. Pijls, NH, De Bruyne, B, Peels, K, Van Der Voort, PH, Bonnier, HJ, Bartunek, JKJJ, et al. Measurement of fractional flow reserve to assess the functional severity of coronary-artery stenoses. N Engl J Med. (1996) 334:1703–8.

Google Scholar

13. De Bruyne, B, Pijls, NH, Kalesan, B, Barbato, E, Tonino, PA, Piroth, Z, et al. Fractional flow reserve-guided PCI versus medical therapy in stable coronary disease. N Engl J Med. (2012) 367:991–1001. doi: 10.1056/NEJMoa1205361

PubMed Abstract | Crossref Full Text | Google Scholar

14. Koo, BK, Erglis, A, Doh, JH, Daniels, DV, Jegere, S, Kim, HS, et al. Diagnosis of ischemia-causing coronary stenoses by noninvasive fractional flow reserve computed from coronary computed tomographic angiograms. results from the prospective multicenter DISCOVER-FLOW (diagnosis of ischemia-causing Stenoses obtained via noninvasive fractional flow reserve) study. J Am Coll Cardiol. (2011) 58:1989–97. doi: 10.1016/j.jacc.2011.06.066

PubMed Abstract | Crossref Full Text | Google Scholar

15. Han, YF, Liu, WH, Chen, XL, Xiong, YY, Yin, Q, Xu, GL, et al. Severity assessment of intracranial large artery stenosis by pressure gradient measurements: a feasibility study. Catheter Cardiovasc Interv. (2016) 88:255–61. doi: 10.1002/ccd.26414

PubMed Abstract | Crossref Full Text | Google Scholar

16. Miao, Z, Liebeskind, DS, Lo, W, Liu, L, Pu, Y, Leng, X, et al. Fractional flow assessment for the evaluation of intracranial atherosclerosis: a feasibility study. Interv Neurol. (2016) 5:65–75. doi: 10.1159/000444333

PubMed Abstract | Crossref Full Text | Google Scholar

17. Yin, X, Yang, R, Li, Z, Wang, H, and Fang, Q. Integrating hemodynamic analysis with traditional imaging in intracranial atherosclerotic stenosis: current status and future perspectives. Front Neurol. (2025) 16:1589162. doi: 10.3389/fneur.2025.1589162

PubMed Abstract | Crossref Full Text | Google Scholar

18. Longevity, O. Retracted: hemodynamic and geometric risk factors for in-stent restenosis in patients with intracranial atherosclerotic stenosis. Oxidative Med Cell Longev. (2024) 2024:9863972. doi: 10.1155/2024/9863972

Crossref Full Text | Google Scholar

19. Wang, M, Leng, X, Mao, B, Zou, R, Lin, D, Gao, Y, et al. Functional evaluation of intracranial atherosclerotic stenosis by pressure ratio measurements. Heliyon. (2023) 9:e13527. doi: 10.1016/j.heliyon.2023.e13527

PubMed Abstract | Crossref Full Text | Google Scholar

20. Lan, L, Liu, H, Ip, V, Soo, Y, Abrigo, J, Fan, F, et al. Regional High Wall shear stress associated with stenosis regression in symptomatic intracranial atherosclerotic disease. Stroke. (2020) 51:3064–73. doi: 10.1161/STROKEAHA.120.030615

PubMed Abstract | Crossref Full Text | Google Scholar

21. Lan, L, Leng, X, Ip, V, Soo, Y, Abrigo, J, Liu, H, et al. Sustaining cerebral perfusion in intracranial atherosclerotic stenosis: the roles of antegrade residual flow and leptomeningeal collateral flow. J Cereb Blood Flow Metab. (2020) 40:126–34. doi: 10.1177/0271678X18805209

PubMed Abstract | Crossref Full Text | Google Scholar

22. Liu, J, Yan, Z, Pu, Y, Shiu, WS, Wu, J, Chen, R, et al. Functional assessment of cerebral artery stenosis: a pilot study based on computational fluid dynamics. J Cereb Blood Flow Metab. (2017) 37:2567–76. doi: 10.1177/0271678X16671321

PubMed Abstract | Crossref Full Text | Google Scholar

23. Zhang, D, Wu, X, Tang, J, Wang, P, Chen, GZ, and Yin, X. Hemodynamics is associated with vessel wall remodeling in patients with middle cerebral artery stenosis. Eur Radiol. (2021) 31:5234–42. doi: 10.1007/s00330-020-07607-w

PubMed Abstract | Crossref Full Text | Google Scholar

24. Dai, Y, Qian, Y, Zhang, M, Li, Y, Lv, P, Tang, X, et al. Associations between local haemodynamics and carotid intraplaque haemorrhage with different stenosis severities: a preliminary study based on MRI and CFD. J Clin Neurosci. (2019) 66:220–5. doi: 10.1016/j.jocn.2019.05.041

PubMed Abstract | Crossref Full Text | Google Scholar

25. Chen, Z, Qin, H, Liu, J, Wu, B, Cheng, Z, Jiang, Y, et al. Characteristics of wall shear stress and pressure of intracranial atherosclerosis analyzed by a computational fluid dynamics model: a pilot study. Front Neurol. (2019) 10:1372. doi: 10.3389/fneur.2019.01372

Crossref Full Text | Google Scholar

26. Antiga, L, Piccinelli, M, Botti, L, Ene-Iordache, B, Remuzzi, A, and Steinman, DA. An image-based modeling framework for patient-specific computational hemodynamics. Med Biol Eng Comput. (2008) 46:1097–112. doi: 10.1007/s11517-008-0420-1

PubMed Abstract | Crossref Full Text | Google Scholar

27. Gibson, CM, Cannon, CP, Daley, WL, Dodge, JT, Alexander, B Jr, Marble, SJ, et al. TIMI frame count: a quantitative method of assessing coronary artery flow. Circulation. (1996) 93:879–88.

Google Scholar

28. Wan, Y, Teng, X, Li, S, and Yang, Y. Application of transcranial Doppler in cerebrovascular diseases. Front Aging Neurosci. (2022) 14:1035086. doi: 10.3389/fnagi.2022.1035086

PubMed Abstract | Crossref Full Text | Google Scholar

29. Samuels, OB, Joseph, GJ, Lynn, MJ, Smith, HA, and Chimowitz, MI. A standardized method for measuring intracranial arterial stenosis. AJNR Am J Neuroradiol. (2000) 21:643–6.

PubMed Abstract | Google Scholar

30. Bartlett, ES, Symons, SP, and Fox, AJ. Correlation of carotid stenosis diameter and cross-sectional areas with CT angiography. AJNR Am J Neuroradiol. (2006) 27:638–42.

PubMed Abstract | Google Scholar

31. Wang, Z, Lu, M, Liu, W, Zheng, T, Li, D, Yu, W, et al. Assessment of carotid atherosclerotic disease using three-dimensional cardiovascular magnetic resonance vessel wall imaging: comparison with digital subtraction angiography. J Cardiovasc Magn Reson. (2020) 22:18. doi: 10.1186/s12968-020-0604-x

PubMed Abstract | Crossref Full Text | Google Scholar

32. Cumming, MJ, and Morrow, IM. Carotid artery stenosis: a prospective comparison of CT angiography and conventional angiography. AJR Am J Roentgenol. (1994) 163:517–23. doi: 10.2214/ajr.163.3.8079836

PubMed Abstract | Crossref Full Text | Google Scholar

33. Link, J, Brossmann, J, Grabener, M, Mueller-Huelsbeck, S, Steffens, JC, Brinkmann, G, et al. Spiral CT angiography and selective digital subtraction angiography of internal carotid artery stenosis. AJNR Am J Neuroradiol. (1996) 17:89–94.

Google Scholar

34. Link, J, Brossmann, J, Penselin, V, Glüer, CC, and Heller, M. Common carotid artery bifurcation: preliminary results of CT angiography and color-coded duplex sonography compared with digital subtraction angiography. AJR Am J Roentgenol. (1997) 168:361–5. doi: 10.2214/ajr.168.2.9016207

PubMed Abstract | Crossref Full Text | Google Scholar

35. Leclerc, X, Godefroy, O, Lucas, C, Benhaim, JF, Michel, TS, Leys, D, et al. Internal carotid arterial stenosis: CT angiography with volume rendering. Radiology. (1999) 210:673–82. doi: 10.1148/radiology.210.3.r99fe46673

PubMed Abstract | Crossref Full Text | Google Scholar

36. Leng, X, Scalzo, F, Ip, HL, Johnson, M, Fong, AK, Fan, FS, et al. Computational fluid dynamics modeling of symptomatic intracranial atherosclerosis may predict risk of stroke recurrence. PLoS One. (2014) 9:e97531. doi: 10.1371/journal.pone.0097531

PubMed Abstract | Crossref Full Text | Google Scholar

37. Tian, X, Fang, H, Lan, L, Ip, HL, Abrigo, J, Liu, H, et al. Risk stratification in symptomatic intracranial atherosclerotic disease with conventional vascular risk factors and cerebral haemodynamics. Stroke Vasc Neurol. (2023) 8:77–85. doi: 10.1136/svn-2022-001606

PubMed Abstract | Crossref Full Text | Google Scholar

38. Huang, J, Degnan, AJ, Liu, Q, Teng, Z, Yue, CS, Gillard, JH, et al. Comparison of NASCET and WASID criteria for the measurement of intracranial stenosis using digital subtraction and computed tomography angiography of the middle cerebral artery. J Neuroradiol. (2012) 39:342–5. doi: 10.1016/j.neurad.2011.11.005

PubMed Abstract | Crossref Full Text | Google Scholar

39. Silvennoinen, HM, Ikonen, S, Soinne, L, Railo, M, and Valanne, L. CT angiographic analysis of carotid artery stenosis: comparison of manual assessment, semiautomatic vessel analysis, and digital subtraction angiography. AJNR Am J Neuroradiol. (2007) 28:97–103.

PubMed Abstract | Google Scholar

40. Pavlin-Premrl, D, Boopathy, SR, Nemes, A, Mohammadzadeh, M, Monajemi, S, Ko, BS, et al. Computational fluid dynamics in intracranial atherosclerosis - lessons from cardiology: a review of CFD in intracranial atherosclerosis. J Stroke Cerebrovasc Dis. (2021) 30:106009. doi: 10.1016/j.jstrokecerebrovasdis.2021.106009

PubMed Abstract | Crossref Full Text | Google Scholar

41. Xu, P, Liu, X, Zhang, H, Ghista, D, Zhang, D, Shi, C, et al. Assessment of boundary conditions for CFD simulation in human carotid artery. Biomech Model Mechanobiol. (2018) 17:1581–97. doi: 10.1007/s10237-018-1045-4

PubMed Abstract | Crossref Full Text | Google Scholar

42. Dörfler, P, Puls, I, Schliesser, M, Mäurer, M, and Becker, G. Measurement of cerebral blood flow volume by extracranial sonography. J Cereb Blood Flow Metab. (2000) 20:269–71. doi: 10.1097/00004647-200002000-00007

PubMed Abstract | Crossref Full Text | Google Scholar

43. Zanaty, M, Rossen, JD, Roa, JA, Nakagawa, D, Hudson, JS, Kasab, SA, et al. Intracranial atherosclerosis: a disease of functional, not anatomic stenosis? How trans-stenotic pressure gradients can help guide treatment. Oper Neurosurg. (2020) 18:599–605. doi: 10.1093/ons/opz335

PubMed Abstract | Crossref Full Text | Google Scholar

44. Li, L, Yang, B, Dmytriw, AA, Wang, T, Luo, J, Li, Y, et al. Hemodynamic versus anatomic assessment of symptomatic atherosclerotic middle cerebral artery stenosis: the relationship between pressure wire Translesional gradient and angiographic lesion geometry. Front Neurol. (2021) 12:671778. doi: 10.3389/fneur.2021.671778

PubMed Abstract | Crossref Full Text | Google Scholar

45. Groen, HC, Gijsen, FJ, van der Lugt, A, Ferguson, MS, Hatsukami, TS, van dersteen, AF, et al. Plaque rupture in the carotid artery is localized at the high shear stress region: a case report. Stroke. (2007) 38:2379–81. doi: 10.1161/STROKEAHA.107.484766

Crossref Full Text | Google Scholar

46. Dong, J, Inthavong, K, and Tu, J. Image-based computational hemodynamics evaluation of atherosclerotic carotid bifurcation models. Comput Biol Med. (2013) 43:1353–62. doi: 10.1016/j.compbiomed.2013.06.013

PubMed Abstract | Crossref Full Text | Google Scholar

47. Stefanadis, C, Antoniou, CK, Tsiachris, D, and Pietri, P. Coronary atherosclerotic vulnerable plaque: current perspectives. J Am Heart Assoc. (2017) 6:543. doi: 10.1161/JAHA.117.005543

PubMed Abstract | Crossref Full Text | Google Scholar

48. Leng, X, Lan, L, Ip, HL, Fan, F, Ma, SH, Ma, K, et al. Translesional pressure gradient and leptomeningeal collateral status in symptomatic middle cerebral artery stenosis. Eur J Neurol. (2018) 25:404–10. doi: 10.1111/ene.13521

PubMed Abstract | Crossref Full Text | Google Scholar

49. Yin, X, Zhao, Y, Huang, F, Wang, H, and Fang, Q. Machine learning-based classification of anterior circulation cerebral infarction using computational fluid dynamics and CT perfusion metrics. Brain Sci. (2025) 15:399. doi: 10.3390/brainsci15040399

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: intracranial atherosclerotic stenosis, computational fluid dynamics, hemodynamics, computed tomography angiography, digital subtraction angiograph

Citation: Yang R, Yin X, Li G, Xiang J, Fang Q, Wang H and Li B (2025) CTA and DSA based computational fluid dynamics models for morphological and hemodynamic assessment of intracranial atherosclerotic stenosis. Front. Neurol. 16:1686189. doi: 10.3389/fneur.2025.1686189

Received: 15 August 2025; Accepted: 23 September 2025;
Published: 08 October 2025.

Edited by:

Gador Canton, University of Washington, United States

Reviewed by:

Tianqi Wang, University of Shanghai for Science and Technology, China
Patrick McGah, TerraPower LLC, United States

Copyright © 2025 Yang, Yin, Li, Xiang, Fang, Wang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Bo Li, bGlibzIwMDNzdXpob3VAMTYzLmNvbQ==; Hui Wang, ZHJfd2FuZ2h1aUBzaW5hLmNvbQ==

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.