UZH-Logo

Three-dimensional computational modeling of subject-specific cerebrospinal fluid flow in the subarachnoid space


Gupta, S; Soellinger, M; Boesiger, P; Poulikakos, D; Kurtcuoglu, V (2009). Three-dimensional computational modeling of subject-specific cerebrospinal fluid flow in the subarachnoid space. Journal of Biomechanical Engineering, 131(2):021010.1-021010.11.

Abstract

This study aims at investigating three-dimensional subject-specific cerebrospinal fluid (CSF) dynamics in the inferior cranial space, the superior spinal subarachnoid space (SAS), and the fourth cerebral ventricle using a combination of a finite-volume computational fluid dynamics (CFD) approach and magnetic resonance imaging (MRI) experiments. An anatomically accurate 3D model of the entire SAS of a healthy volunteer was reconstructed from high resolution T2 weighted MRI data. Subject-specific pulsatile velocity boundary conditions were imposed at planes in the pontine cistern, cerebellomedullary cistern, and in the spinal subarachnoid space. Velocimetric MRI was used to measure the velocity field at these boundaries. A constant pressure boundary condition was imposed at the interface between the aqueduct of Sylvius and the fourth ventricle. The morphology of the SAS with its complex trabecula structures was taken into account through a novel porous media model with anisotropic permeability. The governing equations were solved using finite-volume CFD. We observed a total pressure variation from -42 Pa to 40 Pa within one cardiac cycle in the investigated domain. Maximum CSF velocities of about 15 cm/s occurred in the inferior section of the aqueduct, 14 cm/s in the left foramen of Luschka, and 9 cm/s in the foramen of Magendie. Flow velocities in the right foramen of Luschka were found to be significantly lower than in the left, indicating three-dimensional brain asymmetries. The flow in the cerebellomedullary cistern was found to be relatively diffusive with a peak Reynolds number (Re) = 72, while the flow in the pontine cistern was primarily convective with a peak Re =386. The net volumetric flow rate in the spinal canal was found to be negligible despite CSF oscillation with substantial amplitude with a maximum volumetric flow rate of 109 ml/min. The observed transient flow patterns indicate a compliant behavior of the cranial subarachnoid space. Still, the estimated deformations were small owing to the large parenchymal surface. We have integrated anatomic and velocimetric MRI data with computational fluid dynamics incorporating the porous SAS morphology for the subject-specific reconstruction of cerebrospinal fluid flow in the subarachnoid space. This model can be used as a basis for the development of computational tools, e.g., for the optimization of intrathecal drug delivery and computer-aided evaluation of cerebral pathologies such as syrinx development in syringomelia.

This study aims at investigating three-dimensional subject-specific cerebrospinal fluid (CSF) dynamics in the inferior cranial space, the superior spinal subarachnoid space (SAS), and the fourth cerebral ventricle using a combination of a finite-volume computational fluid dynamics (CFD) approach and magnetic resonance imaging (MRI) experiments. An anatomically accurate 3D model of the entire SAS of a healthy volunteer was reconstructed from high resolution T2 weighted MRI data. Subject-specific pulsatile velocity boundary conditions were imposed at planes in the pontine cistern, cerebellomedullary cistern, and in the spinal subarachnoid space. Velocimetric MRI was used to measure the velocity field at these boundaries. A constant pressure boundary condition was imposed at the interface between the aqueduct of Sylvius and the fourth ventricle. The morphology of the SAS with its complex trabecula structures was taken into account through a novel porous media model with anisotropic permeability. The governing equations were solved using finite-volume CFD. We observed a total pressure variation from -42 Pa to 40 Pa within one cardiac cycle in the investigated domain. Maximum CSF velocities of about 15 cm/s occurred in the inferior section of the aqueduct, 14 cm/s in the left foramen of Luschka, and 9 cm/s in the foramen of Magendie. Flow velocities in the right foramen of Luschka were found to be significantly lower than in the left, indicating three-dimensional brain asymmetries. The flow in the cerebellomedullary cistern was found to be relatively diffusive with a peak Reynolds number (Re) = 72, while the flow in the pontine cistern was primarily convective with a peak Re =386. The net volumetric flow rate in the spinal canal was found to be negligible despite CSF oscillation with substantial amplitude with a maximum volumetric flow rate of 109 ml/min. The observed transient flow patterns indicate a compliant behavior of the cranial subarachnoid space. Still, the estimated deformations were small owing to the large parenchymal surface. We have integrated anatomic and velocimetric MRI data with computational fluid dynamics incorporating the porous SAS morphology for the subject-specific reconstruction of cerebrospinal fluid flow in the subarachnoid space. This model can be used as a basis for the development of computational tools, e.g., for the optimization of intrathecal drug delivery and computer-aided evaluation of cerebral pathologies such as syrinx development in syringomelia.

Citations

25 citations in Web of Science®
39 citations in Scopus®
Google Scholar™

Altmetrics

Downloads

1 download since deposited on 18 Nov 2009
0 downloads since 12 months
Detailed statistics

Additional indexing

Item Type:Journal Article, refereed, original work
Communities & Collections:04 Faculty of Medicine > Institute of Biomedical Engineering
Dewey Decimal Classification:170 Ethics
610 Medicine & health
Language:English
Date:2009
Deposited On:18 Nov 2009 12:56
Last Modified:05 Apr 2016 13:33
Publisher:American Society of Mechanical Engineers
ISSN:0148-0731
Publisher DOI:10.1115/1.3005171
PubMed ID:19102569
Permanent URL: http://doi.org/10.5167/uzh-24014

Download

[img]
Filetype: PDF - Registered users only
Size: 2MB
View at publisher

TrendTerms

TrendTerms displays relevant terms of the abstract of this publication and related documents on a map. The terms and their relations were extracted from ZORA using word statistics. Their timelines are taken from ZORA as well. The bubble size of a term is proportional to the number of documents where the term occurs. Red, orange, yellow and green colors are used for terms that occur in the current document; red indicates high interlinkedness of a term with other terms, orange, yellow and green decreasing interlinkedness. Blue is used for terms that have a relation with the terms in this document, but occur in other documents.
You can navigate and zoom the map. Mouse-hovering a term displays its timeline, clicking it yields the associated documents.

Author Collaborations