Journal List > J Korean Soc Radiol > v.77(3) > 1087833

Park and Ko: Computational Fluid Dynamics: Comparison of Prototype and Commercial Solutions for Intracranial Aneurysms with Different Entrance Length

Abstract

Purpose

The aim of this study is to evaluate the influence of entrance length into computational fluid dynamics (CFD), for distal internal carotid artery (ICA) aneurysms.

Materials and Methods

Two different solutions were used: Siemens prototype software (SW; prototype, not commercially available) and the package of commer-cial CFD SW (commercials). Twelve models of aneurysm at the distal ICA were ob-tained from three dimensional angiography. The original models had longer en-trance lengths from the cervical segment of ICA, whereas intracranial (IC) model lengths were shorter from the proximal cavernous ICA.

Results

The prototype showed faster flow than the commercial CFD, when comparing the maximum value on a velocity scale. Visualization of pathline within the aneu-rysm sac was seen well by the prototype, except one original model which showed a pathline leak. With IC models, both tools showed high correlation with the flow velocity.

Conclusion

We believe that the CFD might be influenced by different solvers and entrance length of distal ICA aneurysm models.

REFERENCES

1.Robertson AM., Watton PN. Computational fluid dynamics in aneurysm research: critical reflections, future directions. AJNR Am J Neuroradiol. 2012. 33:992–995.
crossref
2.Li MH., Chen SW., Li YD., Chen YC., Cheng YS., Hu DJ, et al. Prevalence of unruptured cerebral aneurysms in Chinese adults aged 35 to 75 years: a cross-sectional study. Ann In-tern Med. 2013. 159:514–521.
3.D'Urso PI., Lanzino G., Cloft HJ., Kallmes DF. Flow diversion for intracranial aneurysms: a review. Stroke. 2011. 42:2363–2368.
4.Moyle KR., Antiga L., Steinman DA. Inlet conditions for im-age-based CFD models of the carotid bifurcation: is it rea-sonable to assume fully developed flow? J Biomech Eng. 2006. 128:371–379.
crossref
5.Hoi Y., Wasserman BA., Lakatta EG., Steinman DA. Effect of common carotid artery inlet length on normal carotid bi-furcation hemodynamics. J Biomech Eng. 2010. 132:121008.
crossref
6.Hodis S., Kargar S., Kallmes DF., Dragomir-Daescu D. Artery length sensitivity in patient-specific cerebral aneurysm simulations. AJNR Am J Neuroradiol. 2015. 36:737–743.
crossref
7.Cebral JR., Castro MA., Appanaboyina S., Putman CM., Millan D., Frangi AF. Efficient pipeline for image-based patient-specif-ic analysis of cerebral aneurysm hemodynamics: technique and sensitivity. IEEE Trans Med Imaging. 2005. 24:457–467.
crossref
8.Hodis S., Kallmes DF., Dragomir-Daescu D. Adaptive grid gen-eration in a patient-specific cerebral aneurysm. Phys Rev E Stat Nonlin Soft Matter Phys. 2013. 88:052720.
crossref
9.Jansen IG., Schneiders JJ., Potters WV., van Ooij P., van den Berg R., van Bavel E, et al. Generalized versus patient-specif-ic inflow boundary conditions in computational fluid dy-namics simulations of cerebral aneurysmal hemodynamics. AJNR Am J Neuroradiol. 2014. 35:1543–1548.
crossref

Fig. 1.
Case 2. Segmentation and mesh generation. Prototype (A–D) and commercials (E, F). A. Virtual rendering shows small vessels. B. Model generation with surface triangular mesh after removing small vessels. C. Original model: model processing by adding cutting planes (arrows) and defining them as inlet or outlets. Note the in-equilateral triangles (yellow box) at inlet. D. Making IC model from original model. New inlet of IC model is made by adding cutting plane at the ascending cavernous segment of internal carotid artery. E. Imported model from prototype shows non-circular shape of outlet and in-equilateral triangles (arrowhead). F. After re-meshing, model has circular outlet (yellow, arrowheads) and even shaped triangles. IC = intracranial
jksr-77-133f1.tif
Fig. 2.
Case 3. Flowchart to make and determine the maximum value on velocity scale. A. With commercials, cutting line (arrow) was selected at outlet and is perpendicular to axis of the vessel. Cutting plan image (lower) has velocity scale and Cartesian coordinate system. B. With prototype, cutting plane (left lower) was made by moving axis (arrowheads) in other two multiplanar reconstruction window (upper left and upper right).
jksr-77-133f2.tif
Fig. 3.
Case 11. Capture images of pathline animation. A. Commercials with original model. B. Commercials with intracranial (IC) model. C. Prototype with original. D. Prototype with IC model. More pathlines are visible within the aneurysm sac on prototype, especially IC model. Clear vessel images overlap with the pathline on commercials. Non-segmented vessel image overlaps with the pathline on original model by prototype. No overlapped vessel image of IC one by prototype. The grading of pathline filling were 2 on commercial and 3 by prototype at each model.
jksr-77-133f3.tif
Fig. 4.
Case 7. 3D virtual rendering image (A) shows kissing artifact (arrowheads) between aneurysm sac and adjacent internal carotid artery (ICA). After removing the artifact by punching tool, 3D surface mesh (B) does not show kissing artifact. Leaking artifact (arrows) is noted on both original (C) and intracranial (IC) model (D). Indeed, pathline filling is poorer than the commercial (not shown) on original model contrast to IC model. Note that no leak between the proximal of ophthalmic artery and adjacent ICA (double arrowheads, A) after the punch out the kissing artifact.
jksr-77-133f4.tif
Table 1.
Comparison of Flow Velocities and Pathline Visualizations between Two Models and Solvers
Model l Size P-MV P-PF P-PL C-MV C-PF
1 05 × 3 Org 1884.64 1 0 992.24 0
IC 1502.74 2 0 833.42 0
2 04 × 3 Org 2140.04 2 0 1295.78 1
IC 1195.32 3 0 1000.5 1
3 05 × 3 Org 1607.15 2 0 1110.97 2
IC 1784.29 3 0 1015.63 2
4 06 × 4 Org 1791.64 1 0 797.31 1
IC 1643.12 2 0 763.27 1
5 09 × 6 Org 1229.66 3 1 512.74 1
IC 1101.45 3 0 520.14 1
6 04 × 4 Org 1881.96 3 0 1026.94 1
IC 1397.6 3 0 909.71 1
7 05 × 4 Org 2000.85 0 2 967.11 2
IC 1131.79 3 1 781.06 2
8 05 × 4 Org 1620.59 1 0 1372.07 1
IC 1483.02 2 0 1093.67 1
9 20 × 16 Org 1218.72 3 0 1495.75 2
IC 1610.02 3 0 1105.69 2
10 14 × 10 Org 1497.05 3 0 1169.4 2
IC 2419.36 2 0 1350.55 2
11 08 × 5 Org 1558.49 3 0 1211.42 2
IC 2310.74 3 0 1659.22 2
12 14 × 7 Org 1587.6 3 0 1169.4 1
IC 1847.96 3 0 1350.55 1

Grade 0 = none, grade 1 = < 1/3, grade 2 = 1/3–2/3, grade3 = > 2/3 of sac C = commercials, IC = intracranial model, MV = maximum value on veloci-ty scale (mm/sec), Org = original model, P = Siemens prototype, PF = path-line filling within aneurysm sac, PL = pathline leak which noted only on pro-totype, Size = size of aneurysm (mm)

TOOLS
Similar articles