A computational model for the dynamics of cerebrospinal fluid in the spinal subarachnoid space

Global models for the dynamics of coupled fluid compartments of the central nervous system (CNS) require simplified representations of the individual components which are both accurate and computationally efficient. This paper presents a one-dimensional model for computing the flow of cerebrospinal fluid (CSF) within the spinal subarachnoid space (SSAS) under the simplifying assumption that it consists of two coaxial tubes representing the spinal cord and the dura. A rigorous analysis of the first-order nonlinear system demonstrates that the system is elliptic-hyperbolic, and hence ill-posed, for some values of parameters, being hyperbolic otherwise. In addition, the system cannot be written in conservation-law form, and thus, an appropriate numerical approach is required, namely the path conservative approach. The designed computational algorithm is shown to be second-order accurate in both space and time, capable of handling strongly nonlinear discontinuities, and a method of coupling it with an unsteady inflow condition is presented. Such an approach is sufficiently rapid to be integrated into a global, closed-loop model for computing the dynamics of coupled fluid compartments of the CNS.