Research Papers

Characterization of the Discrepancies Between Four-Dimensional Phase-Contrast Magnetic Resonance Imaging and In-Silico Simulations of Cerebrospinal Fluid Dynamics

[+] Author and Article Information
Soroush Heidari Pahlavian

Conquer Chiari Research Center,
Department of Mechanical Engineering,
The University of Akron,
Akron, OH 44325
e-mail: sh113@uakron.edu

Alexander C. Bunck

Department of Radiology,
University Hospital of Cologne,
Cologne 50923, Germany
e-mail: alexander.bunck@uk-koeln.de

Francis Loth

Conquer Chiari Research Center,
Department of Mechanical Engineering,
The University of Akron,
Akron, OH 44325
e-mail: loth@uakron.edu

R. Shane Tubbs

Children’s of Alabama,
Birmingham AL 35233
Department of Anatomical Sciences,
St. George’s University,
St. George, Grenada
Centre for Anatomy and Human Identification,
Dundee University,
Dundee DD1 4HN, UK
e-mail: Shane.Tubbs@childrensal.org

Theresia Yiallourou

Laboratory of Hemodynamics and Cardiovascular Technology,
Lausanne 1015, Switzerland
e-mail: thiresia.giallourou@epfl.ch

Jan Robert Kroeger

Department of Radiology,
University Hospital of Cologne,
Cologne 50923, Germany
e-mail: janrobertkroeger@gmail.com

Walter Heindel

Department of Clinical Radiology,
University Hospital of Muenster,
Muenster 48149, Germany
e-mail: heindel@uni-muenster.de

Bryn A. Martin

Conquer Chiari Research Center,
Department of Mechanical Engineering,
The University of Akron,
Akron, OH 44325
e-mail: director@chiari-research.org

Manuscript received August 18, 2014; final manuscript received January 21, 2015; published online February 20, 2015. Assoc. Editor: Tim David.

J Biomech Eng 137(5), 051002 (May 01, 2015) (8 pages) Paper No: BIO-14-1406; doi: 10.1115/1.4029699 History: Received August 18, 2014; Revised January 21, 2015; Online February 20, 2015

The purpose of the present study was to compare subject-specific magnetic resonance imaging (MRI)-based computational fluid dynamics (CFD) simulations with time-resolved three-directional (3D) velocity-encoded phase-contrast MRI (4D PCMRI) measurements of the cerebrospinal fluid (CSF) velocity field in the cervical spinal subarachnoid space (SSS). Three-dimensional models of the cervical SSS were constructed based on MRI image segmentation and anatomical measurements for a healthy subject and patient with Chiari I malformation. CFD was used to simulate the CSF motion and compared to the 4D PCMRI measurements. Four-dimensional PCMRI measurements had much greater CSF velocities compared to CFD simulations (1.4 to 5.6× greater). Four-dimensional PCMRI and CFD both showed anterior and anterolateral dominance of CSF velocities, although this flow feature was more pronounced in 4D PCMRI measurements compared to CFD. CSF flow jets were present near the nerve rootlets and denticulate ligaments (NRDL) in the CFD simulation. Flow jets were visible in the 4D PCMRI measurements, although they were not clearly attributable to nerve rootlets. Inclusion of spinal cord NRDL in the cervical SSS does not fully explain the differences between velocities obtained from 4D PCMRI measurements and CFD simulations.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Bunck, A. C., Kroeger, J. R., Juettner, A., Brentrup, A., Fiedler, B., Crelier, G. R., Martin, B. A., Heindel, W., Maintz, D., Schwindt, W., and Niederstadt, T., 2012, “Magnetic Resonance 4D Flow Analysis of Cerebrospinal Fluid Dynamics in Chiari I Malformation With and Without Syringomyelia,” Eur. Radiol., 22(9), pp. 1860–1870. [CrossRef] [PubMed]
Clarke, E. C., Stoodley, M. A., and Bilston, L. E., 2013, “Changes in Temporal Flow Characteristics of CSF in Chiari Malformation Type I With and Without Syringomyelia: Implications for Theory of Syrinx Development,” J. Neurosurg., 118(5), pp. 1135–1140. [CrossRef] [PubMed]
Martin, B. A., Labuda, R., Royston, T. J., Oshinski, J. N., Iskandar, B., and Loth, F., 2010, “Spinal Subarachnoid Space Pressure Measurements in an In Vitro Spinal Stenosis Model: Implications on Syringomyelia Theories,” ASME J. Biomech. Eng., 132(11), p. 111007. [CrossRef]
Bradley, W. G., Jr., Scalzo, D., Queralt, J., Nitz, W. N., Atkinson, D. J., and Wong, P., 1996, “Normal-Pressure Hydrocephalus: Evaluation With Cerebrospinal Fluid Flow Measurements at MR Imaging,” Radiology, 198(2), pp. 523–529. [CrossRef] [PubMed]
Hsu, Y., Hettiarachchi, H. D., Zhu, D. C., and Linninger, A. A., 2012, “The Frequency and Magnitude of Cerebrospinal Fluid Pulsations Influence Intrathecal Drug Distribution: Key Factors for Interpatient Variability,” Anesth. Analg., 115(4), pp. 879–879. [CrossRef]
Papisov, M. I., Belov, V. V., and Gannon, K. S., 2013, “Physiology of the Intrathecal Bolus: The Leptomeningeal Route for Macromolecule and Particle Delivery to CNS,” Mol. Pharm., 10(5), pp. 1522–1532. [CrossRef] [PubMed]
Haughton, V. M., Korosec, F. R., Medow, J. E., Dolar, M. T., and Iskandar, B. J., 2003, “Peak Systolic and Diastolic CSF Velocity in the Foramen Magnum in Adult Patients With Chiari I Malformations and in Normal Control Participants,” AJNR. Am. J. Neuroradiol., 24(2), pp. 169–176. [PubMed]
Iskandar, B. J., Quigley, M., and Haughton, V. M., 2004, “Foramen Magnum Cerebrospinal Fluid Flow Characteristics in Children With Chiari I Malformation Before and After Craniocervical Decompression,” J. Neurosurg., 101(2 Suppl), pp. 169–178. [CrossRef] [PubMed]
Quigley, M. F., Iskandar, B., Quigley, M. E., Nicosia, M., and Haughton, V., 2004, “Cerebrospinal Fluid Flow in Foramen Magnum: Temporal and Spatial Patterns at MR Imaging in Volunteers and in Patients With Chiari I Malformation,” Radiology, 232(1), pp. 229–236. [CrossRef] [PubMed]
Krueger, K. D., Haughton, V. M., and Hetzel, S., 2010, “Peak CSF Velocities in Patients With Symptomatic and Asymptomatic Chiari I Malformation,” AJNR. Am. J. Neuroradiol., 31(10), pp. 1837–1841. [CrossRef] [PubMed]
Bunck, A., Kroeger, J., Juettner, A., Brentrup, A., Fiedler, B., Crelier, G., Martin, B., Heindel, W., Maintz, D., Schwindt, W., and Niederstadt, T., 2012, “Magnetic Resonance 4D Flow Analysis of Cerebrospinal Fluid Dynamics in Chiari I Malformation With and Without Syringomyelia,” Eur. Radiol., 22(9), pp. 1860–1870. [CrossRef] [PubMed]
Bunck, A. C., Kroger, J. R., Juttner, A., Brentrup, A., Fiedler, B., Schaarschmidt, F., Crelier, G. R., Schwindt, W., Heindel, W., Niederstadt, T., and Maintz, D., 2011, “Magnetic Resonance 4D Flow Characteristics of Cerebrospinal Fluid at the Craniocervical Junction and the Cervical Spinal Canal,” Eur. Radiol., 21(8), pp. 1788–1796. [CrossRef] [PubMed]
Santini, F., Wetzel, S. G., Bock, J., Markl, M., and Scheffler, K., 2009, “Time-Resolved Three-Dimensional (3D) Phase-Contrast (PC) Balanced Steady-State Free Precession (BSSFP),” Magn. Reson. Med., 62(4), pp. 966–974. [CrossRef] [PubMed]
Sweetman, B., and Linninger, A. A., 2011, “Cerebrospinal Fluid Flow Dynamics in the Central Nervous System,” Ann. Biomed. Eng., 39(1), pp. 484–496. [CrossRef] [PubMed]
Matsumae, M., Hirayama, A., Atsumi, H., Yatsushiro, S., and Kuroda, K., 2013, “Velocity and Pressure Gradients of Cerebrospinal Fluid Assessed With Magnetic Resonance Imaging,” J. Neurosurg., 120(1), pp. 218–227. [CrossRef] [PubMed]
Stadlbauer, A., Salomonowitz, E., Brenneis, C., Ungersbock, K., van der Riet, W., Buchfelder, M., and Ganslandt, O., 2012, “Magnetic Resonance Velocity Mapping of 3D Cerebrospinal Fluid Flow Dynamics in Hydrocephalus: Preliminary Results,” Eur. Radiol., 22(1), pp. 232–242. [CrossRef] [PubMed]
Linge, S. O., Haughton, V., Lovgren, A. E., Mardal, K. A., Helgeland, A., and Langtangen, H. P., 2011, “Effect of Tonsillar Herniation on Cyclic CSF Flow Studied With Computational Flow Analysis,” AJNR. Am. J. Neuroradiol., 32(8), pp. 1474–1481. [CrossRef] [PubMed]
Linge, S. O., Haughton, V., Lovgren, A. E., Mardal, K. A., and Langtangen, H. P., 2010, “CSF Flow Dynamics at the Craniovertebral Junction Studied With an Idealized Model of the Subarachnoid Space and Computational Flow Analysis,” AJNR. Am. J. Neuroradiol., 31(1), pp. 185–192. [CrossRef] [PubMed]
Heidari Pahlavian, S., Yiallourou, T., Tubbs, R. S., Bunck, A. C., Loth, F., Goodin, M., Raisee, M., and Martin, B. A., 2014, “The Impact of Spinal Cord Nerve Roots and Denticulate Ligaments on Cerebrospinal Fluid Dynamics in the Cervical Spine,” PLoS ONE, 9(4), p. e91888. [CrossRef] [PubMed]
Shaffer, N., Martin, B., and Loth, F., 2011, “Cerebrospinal Fluid Hydrodynamics in Type I Chiari Malformation,” Neurol. Res., 33(3), pp. 247–260. [CrossRef] [PubMed]
Yiallourou, T. I., Kröger, J. R., Stergiopulos, N., Maintz, D., Martin, B. A., and Bunck, A. C., 2012, “Comparison of 4D Phase-Contrast MRI Flow Measurements to Computational Fluid Dynamics Simulations of Cerebrospinal Fluid Motion in the Cervical Spine,” PLoS ONE, 7(12), p. e52284. [CrossRef] [PubMed]
Alleyne, C. H., Jr., Cawley, C. M., Barrow, D. L., and Bonner, G. D., 1998, “Microsurgical Anatomy of the Dorsal Cervical Nerve Roots and the Cervical Dorsal Root Ganglion/Ventral Root Complexes,” Surg. Neurol., 50(3), pp. 213–218. [CrossRef] [PubMed]
Lang, J., 1993, Clinical Anatomy of the Cervical Spine, G. Thieme Verlag, University of Wurzburg, Germany.
Tubbs, R. S., Salter, G., Grabb, P. A., and Oakes, W. J., 2001, “The Denticulate Ligament: Anatomy and Functional Significance,” J. Neurosurg., 94(2 Suppl), pp. 271–275. [CrossRef] [PubMed]
Bloomfield, I. G., Johnston, I. H., and Bilston, L. E., 1998, “Effects of Proteins, Blood Cells and Glucose on the Viscosity of Cerebrospinal Fluid,” Pediatr. Neurosurg, 28(5), pp. 246–251. [CrossRef] [PubMed]
Cheng, S., Stoodley, M. A., Wong, J., Hemley, S., Fletcher, D. F., and Bilston, L. E., 2012, “The Presence of Arachnoiditis Affects the Characteristics of CSF Flow in the Spinal Subarachnoid Space: A Modelling Study,” J. Biomech., 45(7), pp. 1186–1191. [CrossRef] [PubMed]
Loth, F., Yardimci, M. A., and Alperin, N., 2001, “Hydrodynamic Modeling of Cerebrospinal Fluid Motion Within the Spinal Cavity,” ASME J. Biomech. Eng., 123(1), pp. 71–79. [CrossRef]
Rutkowska, G., Haughton, V., Linge, S., and Mardal, K. A., 2012, “Patient-Specific 3D Simulation of Cyclic CSF Flow at the Craniocervical Region,” AJNR. Am. J. Neuroradiol., 33(9), pp. 1756–1762. [CrossRef] [PubMed]
Clarke, E. C., Fletcher, D. F., Stoodley, M. A., and Bilston, L. E., 2013, “Computational Fluid Dynamics Modelling of Cerebrospinal Fluid Pressure in Chiari Malformation and Syringomyelia,” J. Biomech., 46(11), pp. 1801–1809. [CrossRef] [PubMed]
Bunck, A. C., Kroeger, J. R., Juettner, A., Brentrup, A., Fiedler, B., Crelier, G. R., Martin, B. A., Heindel, W., Maintz, D., Schwindt, W., and Niederstadt, T., 2012, “Magnetic Resonance 4D Flow Analysis of Cerebrospinal Fluid Dynamics in Chiari I Malformation With and Without Syringomyelia,” Eur. Radiol., 22(9), pp. 1860–1870. [CrossRef] [PubMed]
Malloy, J. J., and Low, F. N., 1976, “Scanning Electron Microscopy of the Subarachnoid Space in the Dog. IV. Subarachnoid Macrophages,” J. Comp. Neurol., 167(3), pp. 257–283. [CrossRef] [PubMed]
Key, E., and Retzius, M., 1876, Studien in der Anatomie des Nervensystems und des Bindegewebes, Samson & Wallin, Stockholm, Sweden.
Nauta, H. J., Dolan, E., and Yasargil, M. G., 1983, “Microsurgical Anatomy of Spinal Subarachnoid Space,” Surg. Neurol., 19(5), pp. 431–437. [CrossRef] [PubMed]
Oldfield, E. H., Muraszko, K., Shawker, T. H., and Patronas, N. J., 1994, “Pathophysiology of Syringomyelia Associated With Chiari I Malformation of the Cerebellar Tonsils. Implications for Diagnosis and Treatment,” J. Neurosurg., 80(1), pp. 3–15. [CrossRef] [PubMed]
Mikulis, D. J., Wood, M. L., Zerdoner, O. A., and Poncelet, B. P., 1994, “Oscillatory Motion of the Normal Cervical Spinal Cord,” Radiology, 192(1), pp. 117–121. [CrossRef] [PubMed]
Cousins, J., and Haughton, V., 2009, “Motion of the Cerebellar Tonsils in the Foramen Magnum During the Cardiac Cycle,” AJNR. Am. J. Neuroradiol., 30(8), pp. 1587–1588. [CrossRef] [PubMed]
Hofmann, E., Warmuth-Metz, M., Bendszus, M., and Solymosi, L., 2000, “Phase-Contrast MR Imaging of the Cervical CSF and Spinal Cord: Volumetric Motion Analysis in Patients With Chiari I Malformation,” AJNR. Am. J. Neuroradiol., 21(1), pp. 151–158. [PubMed]
Milhorat, T. H., Chou, M. W., Trinidad, E. M., Kula, R. W., Mandell, M., Wolpert, C., and Speer, M. C., 1999, “Chiari I Malformation Redefined: Clinical and Radiographic Findings for 364 Symptomatic Patients,” Neurosurgery, 44(5), pp. 1005–1017. [CrossRef] [PubMed]
Kalata, W., Martin, B. A., Oshinski, J. N., Jerosch-Herold, M., Royston, T. J., and Loth, F., 2009, “MR Measurement of Cerebrospinal Fluid Velocity Wave Speed in the Spinal Canal,” IEEE Trans. Biomed. Eng., 56(6), pp. 1765–1768. [CrossRef] [PubMed]
Cheng, S., Fletcher, D., Hemley, S., Stoodley, M., and Bilston, L., 2014, “Effects of Fluid Structure Interaction in a Three Dimensional Model of the Spinal Subarachnoid Space,” J. Biomech., 47(11), pp. 2826–2830. [CrossRef] [PubMed]
Sigmund, E. E., Suero, G. A., Hu, C., McGorty, K., Sodickson, D. K., Wiggins, G. C., and Helpern, J. A., 2012, “High-Resolution Human Cervical Spinal Cord Imaging at 7 T,” NMR Biomed., 25(7), pp. 891–899. [CrossRef] [PubMed]
Elkins, C. J., and Alley, M. T., 2007, “Magnetic Resonance Velocimetry: Applications of Magnetic Resonance Imaging in the Measurement of Fluid Motion,” Exp. Fluids, 43(6), pp. 823–858. [CrossRef]


Grahic Jump Location
Fig. 1

Modeling process for the cervical SSS in the healthy case. (a) MR image segmentation to reconstruct the SSS. (b) Addition of idealized NRDL. (c) CSF flow simulation using CFD.

Grahic Jump Location
Fig. 2

In vivo CSF flow waveforms measured by 4D PCMRI for the healthy subject (top) and CMI patient (bottom) along the spine from C1 to the 1st thoracic vertebrae (th). Note, the magnitude of CSF flow amplitude changes along the cervical spine; a change likely due to SSS compliance. For the present study, the waveform with the greatest amplitude was used as the CFD flow boundary condition. This waveform was located at C3 for the healthy volunteer and at C2M for the CMI patient (C2P = C2 spinous process, C2M = middle of C2 vertebrae).

Grahic Jump Location
Fig. 3

Comparison of peak systolic (negative values) and peak diastolic (positive values) CSF velocities along the cervical spine between 4D PCMRI and CFD for the healthy volunteer (top) and CMI patient (bottom). Peak systolic and peak diastolic velocities were defined as the maximum magnitude of through-plane velocities that occurred anywhere within the designated cross sections at the time points corresponding to peak systole and peak diastole, respectively.

Grahic Jump Location
Fig. 4

Distribution of through-plane peak systolic velocity magnitude at different locations along the cervical spine plotted (a) for the healthy case and (b) for the CMI patient. In each set of contours, the left and right column represents the results obtained from the CFD simulation and 4D PCMRI measurements, respectively. Note: CFD results are shown in two columns; one with a velocity range identical to 4D PCMRI and another with a smaller velocity range to show detail in the CFD velocity profiles.

Grahic Jump Location
Fig. 6

Formation of concentrated velocity regions (jets) between nerve rootlets (NR) and denticulate ligaments (DL). The anterior and posterior side of the geometry is represented by “A” and “P,” respectively. Note that due to the variable insertion lines of the idealized DLs into the spinal cord, these structures were not completely extended inferiorly between the rootlets at some levels, such as the depicted level.

Grahic Jump Location
Fig. 5

Distribution of midsagittal peak systolic velocity magnitude for the healthy subject obtained from CFD (left) and 4D PCMRI (right). Note that the magnitude images used to obtain the CFD geometry (left) had higher resolution than those obtained by 4D PCMRI (right).



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In