Cirovic Srdjan, Kim Minsuok
The Centre for Biomedical Engineering, University of Surrey, Guildford, United Kingdom.
J Biomech Eng. 2012 Feb;134(2):021005. doi: 10.1115/1.4005853.
Modeling of the cerebrospinal fluid (CSF) system in the spine is strongly motivated by the need to understand the origins of pathological conditions such as the emergence and growth of fluid-filled cysts in the spinal cord. In this study, a one-dimensional (1D) approximation for the flow in elastic conduits was used to formulate a model of the spinal CSF compartment. The modeling was based around a coaxial geometry in which the inner elastic cylinder represented the spinal cord, middle elastic tube represented the dura, and the outermost tube represented the vertebral column. The fluid-filled annuli between the cord and dura, and the dura and vertebral column, represented the subarachnoid and epidural spaces, respectively. The system of governing equations was constructed by applying a 1D form of mass and momentum conservation to all segments of the model. The developed 1D model was used to simulate CSF pulse excited by pressure disturbances in the subarachnoid and epidural spaces. The results were compared to those obtained from an equivalent two-dimensional finite element (FE) model which was implemented using a commercial software package. The analysis of linearized governing equations revealed the existence of three types of waves, of which the two slower waves can be clearly related to the wave modes identified in previous similar studies. The third, much faster, wave emanates directly from the vertebral column and has little effect on the deformation of the spinal cord. The results obtained from the 1D model and its FE counterpart were found to be in good general agreement even when sharp spatial gradients of the spinal cord stiffness were included; both models predicted large radial displacements of the cord at the location of an initial cyst. This study suggests that 1D modeling, which is computationally inexpensive and amenable to coupling with the models of the cranial CSF system, should be a useful approach for the analysis of some aspects of the CSF dynamics in the spine. The simulation of the CSF pulse excited by a pressure disturbance in the epidural space, points to the possibility that regions of the spinal cord with abnormally low stiffness may be prone to experiencing large strains due to coughing and sneezing.
对脊柱中脑脊液(CSF)系统进行建模,其强烈动机源于需要了解诸如脊髓中充满液体的囊肿的出现和生长等病理状况的起源。在本研究中,采用弹性管道中流动的一维(1D)近似来构建脊髓CSF腔室的模型。建模基于同轴几何结构,其中内部弹性圆柱体代表脊髓,中间弹性管代表硬脑膜,最外层管代表脊柱。脊髓与硬脑膜之间以及硬脑膜与脊柱之间充满液体的环形区域分别代表蛛网膜下腔和硬膜外间隙。通过将一维形式的质量和动量守恒应用于模型的所有部分来构建控制方程组。所开发的一维模型用于模拟蛛网膜下腔和硬膜外间隙中压力扰动激发的CSF脉冲。将结果与使用商业软件包实现的等效二维有限元(FE)模型得到的结果进行比较。对线性化控制方程的分析揭示了三种类型波的存在,其中两种较慢的波可以明显与先前类似研究中确定的波模式相关。第三种快得多的波直接源自脊柱,对脊髓的变形影响很小。即使包括脊髓刚度的急剧空间梯度,一维模型及其有限元对应模型得到的结果在总体上也吻合良好;两个模型都预测在初始囊肿位置脊髓有较大的径向位移。本研究表明,一维建模计算成本低且易于与颅CSF系统模型耦合,对于分析脊柱中CSF动力学的某些方面应是一种有用的方法。硬膜外间隙中压力扰动激发的CSF脉冲模拟表明,脊髓中刚度异常低的区域可能因咳嗽和打喷嚏而容易经历大应变。