Pressure Wave Propagation in Fluid-Filled Co-Axial Elastic Tubes Part 2: Mechanisms for the Pathogenesis of Syringomyelia

A computational fluid dynamics study to assess the impact of coughing on cerebrospinal fluid dynamics in Chiari type 1 malformation

Chiari type 1 malformation is a neurological disorder characterized by an obstruction of the cerebrospinal fluid (CSF) circulation between the brain (intracranial) and spinal cord (spinal) compartments. Actions such as coughing might evoke spinal cord complications in patients with Chiari type 1 malformation, but the underlying mechanisms are not well understood. More insight into the impact of the obstruction on local and overall CSF dynamics can help reveal these mechanisms. Therefore, our previously developed computational fluid dynamics framework was used to establish a subject-specific model of the intracranial and upper spinal CSF space of a healthy control. In this model, we emulated a single cough and introduced porous zones to model a posterior (OBS-1), mild (OBS-2), and severe posterior-anterior (OBS-3) obstruction. OBS-1 and OBS-2 induced minor changes to the overall CSF pressures, while OBS-3 caused significantly larger changes with a decoupling between the intracranial and spinal compartment. Coughing led to a peak in overall CSF pressure. During this peak, pressure differences between the lateral ventricles and the spinal compartment were locally amplified for all degrees of obstruction. These results emphasize the effects of coughing and indicate that severe levels of obstruction lead to distinct changes in intracranial pressure.

An analytic model for the flow induced in syringomyelia cavities

A simple two-dimensional fluid-structure-interaction problem, involving viscous oscillatory flow in a channel separated by an elastic membrane from a fluid-filled slender cavity, is analyzed to shed light on the flow dynamics pertaining to syringomyelia, a neurological disorder characterized by the appearance of a large tubular cavity (syrinx) within the spinal cord. The focus is on configurations in which the velocity induced in the cavity, representing the syrinx, is comparable to that found in the channel, representing the subarachnoid space surrounding the spinal cord, both flows being coupled through a linear elastic equation describing the membrane deformation. An asymptotic analysis for small stroke lengths leads to closed-form expressions for the leading-order oscillatory flow, and also for the stationary flow associated with the first-order corrections, the latter involving a steady distribution of transmembrane pressure. The magnitude of the induced flow is found to depend strongly on the frequency, with the result that for channel flow rates of non-sinusoidal waveform, as those found in the spinal canal, higher harmonics can dominate the sloshing motion in the cavity, in agreement with previous in vivo observations. Under some conditions, the cycle-averaged transmembrane pressure, also showing a marked dependence on the frequency, changes sign on increasing the cavity transverse dimension (i.e. orthogonal to the cord axis), underscoring the importance of cavity size in connection with the underlying hydrodynamics. The analytic results presented here can be instrumental in guiding future numerical investigations, needed to clarify the pathogenesis of syringomyelia cavities.

A one-dimensional model for the pulsating flow of cerebrospinal fluid in the spinal canal

The monitoring of intracranial pressure (ICP) fluctuations, which is needed in the context of a number of neurological diseases, requires the insertion of pressure sensors, an invasive procedure with considerable risk factors. Intracranial pressure fluctuations drive the wave-like pulsatile motion of cerebrospinal fluid (CSF) along the compliant spinal canal. Systematically derived simplified models relating the ICP fluctuations with the resulting CSF flow rate can be useful in enabling indirect evaluations of the former from non-invasive magnetic resonance imaging (MRI) measurements of the latter. As a preliminary step in enabling these predictive efforts, a model is developed here for the pulsating viscous motion of CSF in the spinal canal, assumed to be a linearly elastic compliant tube of slowly varying section, with a Darcy pressure-loss term included to model the fluid resistance introduced by the trabeculae, which are thin collagen-reinforced columns that form a web-like structure stretching across the spinal canal. Use of Fourier-series expansions enables predictions of CSF flow rate for realistic anharmonic ICP fluctuations. The flow rate predicted using a representative ICP waveform together with a realistic canal anatomy is seen to compare favourably with in vivo phase-contrast MRI measurements at multiple sections along the spinal canal. The results indicate that the proposed model, involving a limited number of parameters, can serve as a basis for future quantitative analyses targeting predictions of ICP temporal fluctuations based on MRI measurements of spinal-canal anatomy and CSF flow rate.

Slosh Simulation in a Computer Model of Canine Syringomyelia

The exact pathogenesis of syringomyelia is unknown. Epidural venous distention during raised intrathoracic pressure (Valsalva) may cause impulsive movement of fluid ("slosh") within the syrinx. Such a slosh mechanism is a proposed cause of syrinx dissection into spinal cord parenchyma resulting in craniocaudal propagation of the cavity. We sought to test the "slosh" hypothesis by epidural excitation of CSF pulse in a computer model of canine syringomyelia. Our previously developed canine syringomyelia computer model was modified to include an epidural pressure pulse. Simulations were run for: cord free of cavities; cord with small syringes at different locations; and cord with a syrinx that was progressively expanding caudally. If small syringes are present, there are peaks of stress at those locations. This effect is most pronounced at the locations at which syringes initially form. When a syrinx is expanding caudally, the peak stress is typically at the caudal end of the syrinx. However, when the syrinx reaches the lumbar region; the stress becomes moderate. The findings support the "slosh" hypothesis, suggesting that small cervical syringes may propagate caudally. However, when the syrinx is large, there is less focal stress, which may explain why a syrinx can rapidly expand but then remain unchanged in shape over years.

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.

Computer simulation of syringomyelia in dogs

Background

Syringomyelia is a pathological condition in which fluid-filled cavities (syringes) form and expand in the spinal cord. Syringomyelia is often linked with obstruction of the craniocervical junction and a Chiari malformation, which is similar in both humans and animals. Some brachycephalic toy breed dogs such as Cavalier King Charles Spaniels (CKCS) are particularly predisposed. The exact mechanism of the formation of syringomyelia is undetermined and consequently with the lack of clinical explanation, engineers and mathematicians have resorted to computer models to identify possible physical mechanisms that can lead to syringes. We developed a computer model of the spinal cavity of a CKCS suffering from a large syrinx. The model was excited at the cranial end to simulate the movement of the cerebrospinal fluid (CSF) and the spinal cord due to the shift of blood volume in the cranium related to the cardiac cycle. To simulate the normal condition, the movement was prescribed to the CSF. To simulate the pathological condition, the movement of CSF was blocked.

Results

For normal conditions the pressure in the SAS was approximately 400 Pa and the same applied to all stress components in the spinal cord. The stress was uniformly distributed along the length of the spinal cord. When the blockage between the cranial and spinal CSF spaces forced the cord to move with the cardiac cycle, shear and axial normal stresses in the cord increased significantly. The sites where the elevated stress was most pronounced coincided with the axial locations where the syringes typically form, but they were at the perimeter rather than in the central portion of the cord. This elevated stress originated from the bending of the cord at the locations where its curvature was high.

Conclusions

The results suggest that it is possible that repetitive stressing of the spinal cord caused by its exaggerated movement could be a cause for the formation of initial syringes. Further consideration of factors such as cord tethering and the difference in mechanical properties of white and grey matter is needed to fully explore this possibility.

Electronic supplementary material

The online version of this article (10.1186/s12917-018-1410-7) contains supplementary material, which is available to authorized users.

Longitudinal extension of myelomalacia by intramedullary and subdural hemorrhage in a canine model of spinal cord injury

BACKGROUND CONTEXT

In canine intervertebral disc (IVD) extrusion, a spontaneous animal model of spinal cord injury, hemorrhage is a consistent finding. In rodent models, hemorrhage might be involved in secondary tissue destruction by biochemical mechanisms.

PURPOSE

This study aimed to investigate a causal association between the extents of intramedullary, subdural and epidural hemorrhage and the severity of spinal cord damage following IVD extrusion in dogs.

STUDY DESIGN/SETTING

A retrospective study using histologic spinal cord sections from 83 dogs euthanized following IVD extrusion was carried out.

METHODS

The degree of hemorrhage (intramedullary, subdural, epidural), the degree of spinal cord damage in the epicenter (white and gray matter), and the longitudinal extent of myelomalacia were graded. Associations between the extent of hemorrhage and the degree of spinal cord damage were evaluated statistically.

RESULTS

Intramedullary and subdural hemorrhages were significantly associated with the degree of white (p<.001/ p=.004) and gray (both p<.001) matter damage, and with the longitudinal extension of myelomalacia (p<.001/p=.005). Intriguingly, accumulation of hemorrhagic cord debris inside or dorsal to a distended and ruptured central canal in segments distant to the epicenter of the lesion was observed exhibiting a wave-like pattern on longitudinal assessment. The occurrence of this debris accumulation was associated with high degrees of tissue destruction (all p<.001).

CONCLUSIONS

Tissue liquefaction and increased intramedullary pressure associated with hemorrhage are involved in the progression of spinal cord destruction in a canine model of spinal cord injury and ascending or descending myelomalacia. Functional and dynamic studies are needed to investigate this concept further.

Volume change theory for syringomyelia: A new perspective

Background:

The etiopathogenesis of syringomyelia is still an enigma. The authors present a novel theory based on fluid dynamics at the craniovertebral (CV) junction to explain the genesis of syringomyelia (SM). The changes in volume of spinal canal, spinal cord, central canal and spinal subarachnoid space (SSS) in relation to the posterior fossa have been analysed, specifically during postural movements of flexion and extension. The effect of fluctuations in volume of spinal canal and its contents associated with cerebrospinal fluid (CSF) flow dynamics at the CV junction have been postulated to cause the origin and propagation of the syringomyelia. The relevant literature on the subject has been reviewed and the author's theory has been discussed.

Conclusion:

Volume of spinal canal in flexion is always greater than that in extension. Flexion of spine causes narrowing of the ventral subarachnoid space (SAS) and widening of dorsal SAS while extension causes reverse changes leading to fluid movement in dorsal spinal SAS in flexion and ventral spinal SAS in extension. Cervical and lumbar spinal region with maximum bulk hence maximum area and volume undergo maximum deformation during postural changes. SSS CSF is the difference between the volume of spinal canal and spinal cord, varies in flexion and extension which is compensated by changes in posterior fossa (CSF) volume in normal circumstances. Blocked SAS at foramen magnum donot permit spinal SAS CSF exchange which during postural changes is compensated by cavitatory/cystic (syrinx) change at locations in cervical and lumbar spine with propensity for maximum deformation. Augmentation of posterior fossa volume by decompression helps by normalization of this CSF exchange dynamics but immobilizing the spinal movement theoretically will cease any dynamic volume changes thereby minimizing the destructive influence of the fluid exchange on the cord. Thus, this theory strengthens the rational of treating patients by either methodology.

Spontaneous Improvement of Chiari I Malformation and Syringomyelia in a Patient With Cystic Fibrosis

BACKGROUND AND IMPORTANCE:

Syringomyelia is highly associated with Chiari I malformation, but the pathophysiologic mechanism of syrinx formation and its relation to downward cerebellar tonsillar displacement remains elusive. Cough, Valsalva maneuver, and other physiological strains transiently exacerbate the clinical symptoms of these conditions, exert profound effects on the flow dynamics across the craniospinal junction, and are thought to play an important role in the pathogenesis of syringomyelia.

CLINICAL PRESENTATION:

We report the case of a patient with cystic fibrosis who presented during an exacerbation of bronchiectasis and was found to have a Chiari I malformation with associated syringomyelia. Eight months later, when the patient had returned to baseline pulmonary status, repeat imaging showed interval improvement in both the size of the syrinx and descent of cerebellar tonsils.

CONCLUSION:

This rare case of spontaneous improvement of syringomyelia and Chiari I malformation attributable to relief from chronic cough offers interesting insight into the mechanism of these disorders.

Poro-elastic modeling of Syringomyelia - a systematic study of the effects of pia mater, central canal, median fissure, white and gray matter on pressure wave propagation and fluid movement within the cervical spinal cord

Syringomyelia, fluid-filled cavities within the spinal cord, occurs frequently in association with a Chiari I malformation and produces some of its most severe neurological symptoms. The exact mechanism causing syringomyelia remains unknown. Since syringomyelia occurs frequently in association with obstructed cerebrospinal fluid (CSF) flow, it has been hypothesized that syrinx formation is mechanically driven. In this study we model the spinal cord tissue either as a poro-elastic medium or as a solid linear elastic medium, and simulate the propagation of pressure waves through an anatomically plausible 3D geometry, with boundary conditions based on in vivo CSF pressure measurements. Then various anatomic and tissue properties are modified, resulting in a total of 11 variations of the model that are compared. The results show that an open segment of the central canal and a stiff pia (relative to the cord) both increase the radial pressure gradients and enhance interstitial fluid flow in the central canal. The anterior median fissure, anisotropic permeability of the white matter, and Poisson ratio play minor roles.

Effect of the central canal in the spinal cord on fluid movement within the cord

Computational studies are used to demonstrate the effect of oscillating CSF flow on pressures within the spinal cord. We tested the hypothesis that the central canal in the spinal cord affects spinal cord pressure gradients resulting from oscillatory CSF flow. Two computational models of the spinal cord were created with the same dimensions. Model 1 had a homogeneous porous structure. Model 2 had the same structure with the addition of a central fluid filled space, representing the central canal of the cord. We simulated oscillatory flow in the fluid space using standard computational fluid dynamics tools. For all phases of the CSF flow cycle and for specific projections through the model we calculated pressure gradients and fluid movement in the cord models. Pressures in the models varied through the flow cycle. Model 1 had linearly varying pressure along its long axis that varied with the cycle and had no pressure gradients across the cord. Model 2 had nonlinear varying pressure along its long axis varying with the time in the cycle and had transient centrifugal and centripetal pressure gradients with a central fluid space. The radial pressures varied linearly with distance from the fluid space. Centrifugal and centripetal pressure gradients resulted in radially directed fluid flow in the cord. The central canal within the spinal cord alters the pressure fields occurring during oscillatory CSF flow and creates centrifugal and centripetal fluid flux in the cord.

Comparison of 4D phase-contrast MRI flow measurements to computational fluid dynamics simulations of cerebrospinal fluid motion in the cervical spine

Cerebrospinal fluid (CSF) dynamics in the cervical spinal subarachnoid space (SSS) have been thought to be important to help diagnose and assess craniospinal disorders such as Chiari I malformation (CM). In this study we obtained time-resolved three directional velocity encoded phase-contrast MRI (4D PC MRI) in three healthy volunteers and four CM patients and compared the 4D PC MRI measurements to subject-specific 3D computational fluid dynamics (CFD) simulations. The CFD simulations considered the geometry to be rigid-walled and did not include small anatomical structures such as nerve roots, denticulate ligaments and arachnoid trabeculae. Results were compared at nine axial planes along the cervical SSS in terms of peak CSF velocities in both the cranial and caudal direction and visual interpretation of thru-plane velocity profiles. 4D PC MRI peak CSF velocities were consistently greater than the CFD peak velocities and these differences were more pronounced in CM patients than in healthy subjects. In the upper cervical SSS of CM patients the 4D PC MRI quantified stronger fluid jets than the CFD. Visual interpretation of the 4D PC MRI thru-plane velocity profiles showed greater pulsatile movement of CSF in the anterior SSS in comparison to the posterior and reduction in local CSF velocities near nerve roots. CFD velocity profiles were relatively uniform around the spinal cord for all subjects. This study represents the first comparison of 4D PC MRI measurements to CFD of CSF flow in the cervical SSS. The results highlight the utility of 4D PC MRI for evaluation of complex CSF dynamics and the need for improvement of CFD methodology. Future studies are needed to investigate whether integration of fine anatomical structures and gross motion of the brain and/or spinal cord into the computational model will lead to a better agreement between the two techniques.

Syrinx fluid transport: modeling pressure-wave-induced flux across the spinal pial membrane

Syrinxes are fluid-filled cavities of the spinal cord that characterize syringomyelia, a disease involving neurological damage. Their formation and expansion is poorly understood, which has hindered successful treatment. Syrinx cavities are hydraulically connected with the spinal subarachnoid space (SSS) enveloping the spinal cord via the cord interstitium and the network of perivascular spaces (PVSs), which surround blood vessels penetrating the pial membrane that is adherent to the cord surface. Since the spinal canal supports pressure wave propagation, it has been hypothesized that wave-induced fluid exchange across the pial membrane may play a role in syrinx filling. To investigate this conjecture a pair of one-dimensional (1-d) analytical models were developed from classical elastic tube theory coupled with Darcy's law for either perivascular or interstitial flow. The results show that transpial flux serves as a mechanism for damping pressure waves by alleviating hoop stress in the pial membrane. The timescale ratio over which viscous and inertial forces compete was explicitly determined, which predicts that dilated PVS, SSS flow obstructions, and a stiffer and thicker pial membrane-all associated with syringomyelia-will increase transpial flux and retard wave travel. It was also revealed that the propagation of a pressure wave is aided by a less-permeable pial membrane and, in contrast, by a more-permeable spinal cord. This is the first modeling of the spinal canal to include both pressure-wave propagation along the spinal axis and a pathway for fluid to enter and leave the cord, which provides an analytical foundation from which to approach the full poroelastic problem.

A one-dimensional model of the spinal cerebrospinal-fluid compartment

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.

The presence of arachnoiditis affects the characteristics of CSF flow in the spinal subarachnoid space: a modelling study

Syringomyelia is a neurological disorder characterised by high pressure fluid-filled cysts within the spinal cord. As syringomyelia is associated with abnormalities of the central nervous system that obstruct cerebrospinal fluid (CSF) flow, it is thought that changes in CSF dynamics play an important role in its pathogenesis. Using three-dimensional computational models of the spinal subarachnoid space (SAS), this study aims to determine SAS obstructions, such as arachnoiditis, change in CSF dynamics in the SAS. The geometry of the SAS was reconstructed from a series of MRI images. CSF is modelled as an incompressible Newtonian fluid with a dynamic viscosity of 1 mPa s. Three computational models simulated CSF flow in either the unobstructed SAS, or with the SAS obstructed by a porous region simulating dorsal or circumferential arachnoiditis. The permeability of this porous obstruction was varied for the model with dorsal arachnoiditis. The results show that arachnoiditis increases flow resistance in the SAS and this is accompanied by a modest increase in magnitude and/or shift in timing (with respect to the cardiac cycle) of the CSF pressure drop across the region of arachnoiditis. This study suggests that syrinx formation may be related to a change in temporal CSF pulse pressure dynamics.

A coupled hydrodynamic model of the cardiovascular and cerebrospinal fluid system

Coupling of the cardiovascular and cerebrospinal fluid (CSF) system is considered to be important to understand the pathophysiology of cerebrovascular and craniospinal disease and intrathecal drug delivery. A coupled cardiovascular and CSF system model was designed to examine the relation of spinal cord (SC) blood flow (SCBF) and CSF pulsations along the spinal subarachnoid space (SSS). A one-dimensional (1-D) cardiovascular tree model was constructed including a simplified SC arterial network. Connection between the cardiovascular and CSF system was accomplished by a transfer function based on in vivo measurements of CSF and cerebral blood flow. A 1-D tube model of the SSS was constructed based on in vivo measurements in the literature. Pressure and flow throughout the cardiovascular and CSF system were determined for different values of craniospinal compliance. SCBF results indicated that the cervical, thoracic, and lumbar SC each had a signature waveform shape. The cerebral blood flow to CSF transfer function reproduced an in vivo-like CSF flow waveform. The 1-D tube model of the SSS resulted in a distribution of CSF pressure and flow and a wave speed that were similar to those in vivo. The SCBF to CSF pulse delay was found to vary a great degree along the spine depending on craniospinal compliance and vascular anatomy. The properties and anatomy of the SC arterial network and SSS were found to have an important impact on pressure and flow and perivascular fluid movement to the SC. Overall, the coupled model provides predictions about the flow and pressure environment in the SC and SSS. More detailed measurements are needed to fully validate the model.

Spinal subarachnoid space pressure measurements in an in vitro spinal stenosis model: implications on syringomyelia theories

Full explanation for the pathogenesis of syringomyelia (SM), a neuropathology characterized by the formation of a cystic cavity (syrinx) in the spinal cord (SC), has not yet been provided. It has been hypothesized that abnormal cerebrospinal fluid (CSF) pressure, caused by subarachnoid space (SAS) flow blockage (stenosis), is an underlying cause of syrinx formation and subsequent pain in the patient. However, paucity in detailed in vivo pressure data has made theoretical explanations for the syrinx difficult to reconcile. In order to understand the complex pressure environment, four simplified in vitro models were constructed to have anatomical similarities with post-traumatic SM and Chiari malformation related SM. Experimental geometry and properties were based on in vivo data and incorporated pertinent elements such as a realistic CSF flow waveform, spinal stenosis, syrinx, flexible SC, and flexible spinal column. The presence of a spinal stenosis in the SAS caused peak-to-peak cerebrospinal fluid CSF pressure fluctuations to increase rostral to the stenosis. Pressure with both stenosis and syrinx present was complex. Overall, the interaction of the syrinx and stenosis resulted in a diastolic valve mechanism and rostral tensioning of the SC. In all experiments, the blockage was shown to increase and dissociate SAS pressure, while the axial pressure distribution in the syrinx remained uniform. These results highlight the importance of the properties of the SC and spinal SAS, such as compliance and permeability, and provide data for comparison with computational models. Further research examining the influence of stenosis size and location, and the importance of tissue properties, is warranted.

Characterization of cyclic CSF flow in the foramen magnum and upper cervical spinal canal with MR flow imaging and computational fluid dynamics

CSF flow has been shown to exhibit complex patterns in MR images in both healthy subjects and in patients with Chiari I. Abnormal CSF flow oscillations, according to prevailing opinion, cause syringomyelia and other clinical manifestations that affect some patients with the Chiari I malformation. For this article, we reviewed the literature on PC MR of CSF flow, collected the published CFD studies relevant to CSF flow, and performed flow simulations. PC MR creates cine and still images of CSF flow and measurements of flow velocities. CFD, a technique used to compute flow and pressure in liquid systems, simulates the CSF flow patterns that occur in a specific geometry or anatomy of the SAS and a specific volume of flow. Published PC MR studies show greater peak CSF velocities and more complex flow patterns in patients with Chiari I than in healthy subjects, with synchronous bidirectional flow one of the characteristic markers of pathologic flow. In mathematic models of the SAS created from high-resolution MR images, CFD displays complex CSF flow patterns similar to those shown in PC MR in patients. CFD shows that the pressure and flow patterns vary from level to level in the upper spinal canal and differ between patients with Chiari and healthy volunteers. In models in which elasticity and motion are incorporated, CFD displays CSF pressure waves in the SAS. PC MR and CFD studies to date demonstrate significant alterations of CSF flow and pressure patterns in patients with Chiari I. CSF flow has nonlaminar complex spatial and temporal variations and associated pressure waves and pressure gradients. Additional simulations of CSF flow supplemented by PC MR will lead to better measures for distinguishing pathologic flow abnormalities that cause syringomyelia, headaches, and other clinical manifestations in Chiari I malformations.

Evaluation by Fluid/Structure-Interaction Spinal-Cord Simulation of the Effects of Subarachnoid-Space Stenosis on an Adjacent Syrinx

A finite-element numerical model was constructed of the spinal cord, pia mater, filum terminale, cerebrospinal fluid in the spinal subarachnoid space (SSS), and dura mater. The cord was hollowed out by a thoracic syrinx of length 140 mm, and the SSS included a stenosis of length 30 mm opposite this syrinx. The stenosis severity was varied from 0% to 90% by area. Pressure pulse excitation was applied to the model either at the cranial end of the SSS, simulating the effect of cranial arterial pulsation, or externally to the abdominal dura mater, simulating the effect of cough. A very short pulse was used to examine wave propagation; a pulse emulating cardiac systole was used to examine the effects of fluid displacement. Additionally, repetitive sinusoidal excitation was applied cranially. Bulk fluid flow past the stenosis gave rise to prominent longitudinal pressure dissociation (“suck”) in the SSS adjacent to the syrinx. However, this did not proportionally increase the longitudinal motion of fluid in the syrinx. The inertia of the fluid in the SSS, together with the compliance of this space, gave a resonance capable of being excited constructively or destructively by cardiac or coughing impulses. The main effect of mild stenosis was to lower the frequency of this resonance; severe stenosis damped out to-and-fro motions after the end of the applied excitation. Syrinx fluid motion indicated the fluid momentum and thus the pressure developed when the fluid was stopped by the end of the syrinx; however, the tearing stress in the local cord material depended also on the instantaneous local SSS pressure and was therefore not well predicted by syrinx fluid motion. Stenosis was also shown to give rise to a one-way valve effect causing raised SSS pressure caudally and slight average cord displacement cranially. The investigation showed that previous qualitative predictions of the effects of suck neglected factors that reduced the extent of the resulting syrinx fluid motion and of the cord tearing stress, which ultimately determines whether the syrinx lengthens.

The influence of coughing on cerebrospinal fluid pressure in an in vitro syringomyelia model with spinal subarachnoid space stenosis

Background

The influence of coughing, on the biomechanical environment in the spinal subarachnoid space (SAS) in the presence of a cerebrospinal fluid flow stenosis, is thought to be an important etiological factor in craniospinal disorders, including syringomyelia (SM), Chiari I malformation, and hydrocephalus. The aim of this study was to investigate SAS and syrinx pressures during simulated coughing using in vitro models and to provide information for the understanding of the craniospinal fluid system dynamics to help develop better computational models.

Methods

Four in vitro models were constructed to be simplified representations of: 1) non-communicating SM with spinal SAS stenosis; 2) non-communicating SM due to spinal SAS stenosis with a distensible spinal column; 3) non-communicating SM post surgical removal of a spinal SAS stenosis; and 4) a spinal SAS stenosis due to spinal trauma. All of the models had a flexible spinal cord. To simulate coughing conditions, an abrupt CSF pressure pulse (~ 5 ms) was imposed at the caudal end of the spinal SAS by a computer-controlled pump. Pressure measurements were obtained at 4 cm intervals along the spinal SAS and syrinx using catheter tip transducers.

Results

Pressure measurements during a simulated cough, showed that removal of the stenosis was a key factor in reducing pressure gradients in the spinal SAS. The presence of a stenosis resulted in a caudocranial pressure drop in the SAS, whereas pressure within the syrinx cavity varied little caudocranially. A stenosis in the SAS caused the syrinx to balloon outward at the rostral end and be compressed at the caudal end. A >90% SAS stenosis did not result in a significant Venturi effect. Increasing compliance of the spinal column reduced forces acting on the spinal cord. The presence of a syrinx in the cord when there was a stenosis in the SAS, reduced pressure forces in the SAS. Longitudinal pressure dissociation acted to suck fluid and tissue caudocranially in the SAS with a stenosis.

Conclusions

Pressures in the spinal SAS during a simulated cough in vitro had similar peak, transmural, and longitudinal pressures to in vivo measurements reported in the literature. The pressure wave velocities and pressure gradients during coughing (longitudinal pressure dissociation and transmural pressure) were impacted by alterations in geometry, compliance, and the presence of a syrinx and/or stenosis.

The Pathogenesis of Syringomyelia: A Re-Evaluation of the Elastic-Jump Hypothesis

Syringomyelia is a disease in which fluid-filled cavities, called syrinxes, form in the spinal cord causing progressive loss of sensory and motor functions. Invasive monitoring of pressure waves in the spinal subarachnoid space implicates a hydrodynamic origin. Poor treatment outcomes have led to myriad hypotheses for its pathogenesis, which unfortunately are often based on small numbers of patients due to the relative rarity of the disease. However, only recently have models begun to appear based on the principles of mechanics. One such model is the mathematically rigorous work of Carpenter and colleagues (2003, “Pressure Wave Propagation in Fluid-Filled Co-Axial Elastic Tubes Part 1: Basic Theory,” ASME J. Biomech. Eng., 125(6), pp. 852–856; 2003, “Pressure Wave Propagation in Fluid-Filled Co-Axial Elastic Tubes Part 2: Mechanisms for the Pathogenesis of Syringomyelia,” ASME J. Biomech. Eng., 125(6), pp. 857–863). They suggested that a pressure wave due to a cough or sneeze could form a shocklike elastic jump, which when incident at a stenosis, such as a hindbrain tonsil, would generate a transient region of high pressure within the spinal cord and lead to fluid accumulation. The salient physiological parameters of this model were reviewed from the literature and the assumptions and predictions re-evaluated from a mechanical standpoint. It was found that, while the spinal geometry does allow for elastic jumps to occur, their effects are likely to be weak and subsumed by the small amount of viscous damping present in the subarachnoid space. Furthermore, the polarity of the pressure differential set up by cough-type impulses opposes the tenets of the elastic-jump hypothesis. The analysis presented here does not support the elastic-jump hypothesis or any theory reliant on cough-based pressure impulses as a mechanism for the pathogenesis of syringomyelia.

A Coaxial Tube Model of the Cerebrospinal Fluid Pulse Propagation in the Spinal Column

The dynamics of the movement of the cerebrospinal fluid (CSF) may play an important role in the genesis of pathological neurological conditions such as syringomyelia, which is characterized by the presence of a cyst (syrinx) in the spinal cord. In order to provide sound theoretical grounds for the hypotheses that attribute the formation and growth of the syrinx to impediments to the normal movement of the CSF, it is necessary to understand various modes through which CSF pulse in the spinal column propagates. Analytical models of small-amplitude wave propagation in fluid-filled coaxial tubes, where the outer tube represents dura, the inner tube represents the spinal cord, and the fluid is the CSF, have been used to that end. However, so far, the tendency was to model one of the two tubes as rigid and to neglect the effect of finite thickness of the tube walls. The aim of this study is to extend the analysis in order to address these two potentially important issues. To that end, classical linear small-amplitude analysis of wave propagation was applied to a system consisting of coaxial tubes of finite thickness filled with inviscid incompressible fluid. General solutions to the governing equations for the case of harmonic waves in the long wave limit were replaced with the boundary conditions to yield the characteristic (dispersion) equation for the system. The four roots of the characteristic equation correspond to four modes of wave propagation, of which the first three are associated with significant motion of the CSF. For the normal range of parameters the speeds of the four modes are c1=13m∕s, c2=14.7m∕s, c3=30.3m∕s, and c4=124.5m∕s, which are well within the range of values previously reported in experimental and theoretical studies. The modes with the highest and the lowest speeds of propagation can be attributed to the dura and the spinal cord, respectively, whereas the remaining two modes involve some degree of coupling between the two. When the thickness of the spinal cord, is reduced below its normal value, the first mode becomes dominant in terms of the movement of the CSF, and its speed drops significantly. This suggests that the syrinx may be characterized by an abnormally low speed of the CSF pulse.

A unifying hypothesis for hydrocephalus, Chiari malformation, syringomyelia, anencephaly and spina bifida

This work is a modified version of the Casey Holter Memorial prize essay presented to the Society for Research into Hydrocephalus and Spina Bifida, June 29th 2007, Heidelberg, Germany. It describes the origin and consequences of the Chiari malformation, and proposes that hydrocephalus is caused by inadequate central nervous system (CNS) venous drainage. A new hypothesis regarding the pathogenesis, anencephaly and spina bifida is described.Any volume increase in the central nervous system can increase venous pressure. This occurs because veins are compressible and a CNS volume increase may result in reduced venous blood flow. This has the potential to cause progressive increase in cerebrospinal fluid (CSF) volume. Venous insufficiency may be caused by any disease that reduces space for venous volume. The flow of CSF has a beneficial effect on venous drainage. In health it moderates central nervous system pressure by moving between the head and spine. Conversely, obstruction to CSF flow causes localised pressure increases, which have an adverse effect on venous drainage.The Chiari malformation is associated with hindbrain herniation, which may be caused by low spinal pressure relative to cranial pressure. In these instances, there are hindbrain-related symptoms caused by cerebellar and brainstem compression. When spinal injury occurs as a result of a Chiari malformation, the primary pathology is posterior fossa hypoplasia, resulting in raised spinal pressure. The small posterior fossa prevents the flow of CSF from the spine to the head as blood enters the central nervous system during movement. Consequently, intermittent increases in spinal pressure caused by movement, result in injury to the spinal cord. It is proposed that posterior fossa hypoplasia, which has origins in fetal life, causes syringomyelia after birth and leads to damage to the spinal cord in spina bifida. It is proposed that hydrocephalus may occur as a result of posterior fossa hypoplasia, where raised pressure occurs as a result of obstruction to flow of CSF from the head to the spine, and cerebral injury with raised pressure occurs in anencephaly by this mechanism.The current view of dysraphism is that low central nervous system pressure and exposure to amniotic fluid, damage the central nervous system. The hypothesis proposed in this essay supports the view that spina bifida is a manifestation of progressive hydrocephalus in the fetus. It is proposed that that mesodermal growth insufficiency influences both neural tube closure and central nervous system pressure, leading to dysraphism.

Imaging diagnosis--degenerative intraspinal cyst associated with an intervertebral disc

A 6-year-old neutered male Rottweiler had chronic episodic signs of thoracolumbar pain and inability to stand that did not improve after rest and nonsteroidal antiinflammatory medication. In magnetic resonance images, an extradural mass was identified between the spinal cord and the T13/L1 intervertebral disc; the mass was compressing the cord in a dorsoventral direction. The lesion had a low signal intensity rim that enhanced slightly after gadolinium administration, and contents with similar signal intensity to cerebrospinal fluid. The lesion was removed surgically by severing its attachments to the annulus fibrosus. The histologic diagnosis was degenerative intraspinal cyst. The origin of such cysts, and their relationship to synovial and ganglion cysts is discussed.

Ascending central canal dilation and progressive ependymal disruption in a contusion model of rodent chronic spinal cord injury

BACKGROUND

Chronic spinal cord injury (SCI) can lead to an insidious decline in motor and sensory function in individuals even years after the initial injury and is accompanied by a slow and progressive cytoarchitectural destruction. At present, no pathological mechanisms satisfactorily explain the ongoing degeneration.

METHODS

Adult female Sprague-Dawley rats were anesthetized laminectomized at T10 and received spinal cord contusion injuries with a force of 250 kilodynes using an Infinite Horizon Impactor. Animals were randomly distributed into 5 groups and killed 1 (n = 4), 28 (n = 4), 120 (n = 4), 450 (n = 5), or 540 (n = 5) days after injury. Morphometric and immunohistochemical studies were then performed on 1 mm block sections, 6 mm cranial and 6 mm caudal to the lesion epicenter. The SPSS 11.5 t test was used to determine differences between quantitative measures.

RESULTS

Here, we document the first report of an ascending central canal dilation and progressive ependymal disruption cranial to the epicenter of injury in a contusion model of chronic SCI, which was characterized by extensive dural fibrosis and intraparenchymal cystic cavitation. Expansion of the central canal lumen beyond a critical diameter corresponded with ependymal cell ciliary loss, an empirically predictable thinning of the ependymal region, and a decrease in cell proliferation in the ependymal region. Large, aneurysmal dilations of the central canal were accompanied by disruptions in the ependymal layer, periependymal edema and gliosis, and destruction of the adjacent neuropil.

CONCLUSION

Cells of the ependymal region play an important role in CSF homeostasis, cellular signaling and wound repair in the spinal cord. The possible effects of this ascending pathology on ependymal function are discussed. Our studies suggest central canal dilation and ependymal region disruption as steps in the pathogenesis of chronic SCI, identify central canal dilation as a marker of chronic SCI and provide novel targets for therapeutic intervention.

Focal spinal arachnoiditis increases subarachnoid space pressure: a computational study

BACKGROUND

Enlarging fluid filled cystic cavitations form within the spinal cord in up to 28% of spinal cord injured patients. These post-traumatic syrinxes can cause neurological deterioration and current treatment results are unsatisfactory. Localized scar tissue (arachnoiditis) within the subarachnoid space at the level of injury has been suggested to be involved in the pathogenesis of syrinx formation. This study tests the hypothesis that pressure pulses in the subarachnoid space are accentuated adjacent to regions of arachnoiditis, which may drive fluid into the spinal cord and contribute to syrinx formation.

METHODS

An axisymmetric, cylindrical computational fluid dynamics model was developed to represent the subarachnoid space under normal physiological conditions and in the presence of arachnoiditis. Cerebrospinal fluid flow into the model was estimated from magnetic resonance imaging flow studies. Arachnoiditis was modelled as a porous obstruction in the subarachnoid space.

FINDINGS

Peak fluid pressures were higher above the obstruction than in the absence of obstruction. The peak pressures were strongly dependent on the permeability of the obstruction.

INTERPRETATION

Elevations in subarachnoid space pressures due to arachnoiditis may facilitate fluid flow into the spinal cord, enhancing syrinx formation. This suggests that it may be worthwhile to investigate strategies that inhibit arachnoiditis or minimize systolic pressure peaks for treating or preventing syringomyelia.

Unraveling the riddle of syringomyelia

The pathophysiology of syringomyelia development is not fully understood. Current prevailing theories suggest that increased pulse pressure in the subarachnoid space forces cerebrospinal fluid (CSF) through the spinal cord into the syrinx. It is generally accepted that the syrinx consists of CSF. The here-proposed intramedullary pulse pressure theory instead suggests that syringomyelia is caused by increased pulse pressure in the spinal cord and that the syrinx consists of extracellular fluid. A new principle is introduced implying that the distending force in the production of syringomyelia is a relative increase in pulse pressure in the spinal cord compared to that in the nearby subarachnoid space. The formation of a syrinx then occurs by the accumulation of extracellular fluid in the distended cord. A previously unrecognized mechanism for syrinx formation, the Bernoulli theorem, is also described. The Bernoulli theorem or the Venturi effect states that the regional increase in fluid velocity in a narrowed flow channel decreases fluid pressure. In Chiari I malformations, the systolic CSF pulse pressure and downward motion of the cerebellar tonsils are significantly increased. This leads to increased spinal CSF velocities and, as a consequence of the Bernoulli theorem, decreased fluid pressure in narrow regions of the spinal CSF pathways. The resulting relatively low CSF pressure in the narrowed CSF pathway causes a suction effect on the spinal cord that distends the cord during each systole. Syringomyelia develops by the accumulation of extracellular fluid in the distended cord. In posttraumatic syringomyelia, the downwards directed systolic CSF pulse pressure is transmitted and reflected into the spinal cord below and above the traumatic subarachnoid blockage, respectively. The ensuing increase in intramedullary pulse pressure distends the spinal cord and causes syringomyelia on both sides of the blockage. The here-proposed concept has the potential to unravel the riddle of syringomyelia and affords explanations to previously unanswered clinical and theoretical problems with syringomyelia. It also explains why syringomyelia associated with Chiari I malformations may develop in any part of the spinal cord including the medullary conus. Syringomyelia thus preferentially develops where the systolic CSF flow causes a suction effect on the spinal cord, i.e., at or immediately caudal to physiological or pathological encroachments of the spinal subarachnoid space.

Syringomyelia Hydrodynamics: An In Vitro Study Based on In Vivo Measurements

A simplified in vitro model of the spinal canal, based on in vivo magnetic resonance imaging, was used to examine the hydrodynamics of the human spinal cord and subarachnoid space with syringomyelia. In vivo magnetic resonance imaging (MRI) measurements of subarachnoid (SAS) geometry and cerebrospinal fluid velocity were acquired in a patient with syringomyelia and used to aid in the in vitro model design and experiment. The in vitro model contained a fluid-filled coaxial elastic tube to represent a syrinx. A computer controlled pulsatile pump was used to subject the in vitro model to a CSF flow waveform representative of that measured in vivo. Fluid velocity was measured at three axial locations within the in vitro model using the same MRI scanner as the patient study. Pressure and syrinx wall motion measurements were conducted external to the MR scanner using the same model and flow input. Transducers measured unsteady pressure both in the SAS and intra-syrinx at four axial locations in the model. A laser Doppler vibrometer recorded the syrinx wall motion at 18 axial locations and three polar positions. Results indicated that the peak-to-peak amplitude of the SAS flow waveform in vivo was approximately tenfold that of the syrinx and in phase (SAS∼5.2±0.6ml∕s,syrinx∼0.5±0.3ml∕s). The in vitro flow waveform approximated the in vivo peak-to-peak magnitude (SAS∼4.6±0.2ml∕s,syrinx∼0.4±0.3ml∕s). Peak-to-peak in vitro pressure variation in both the SAS and syrinx was approximately 6 mm Hg. Syrinx pressure waveform lead the SAS pressure waveform by approximately 40 ms. Syrinx pressure was found to be less than the SAS for ∼200ms during the 860-ms flow cycle. Unsteady pulse wave velocity in the syrinx was computed to be a maximum of ∼25m∕s. LDV measurements indicated that spinal cord wall motion was nonaxisymmetric with a maximum displacement of ∼140μm, which is below the resolution limit of MRI. Agreement between in vivo and in vitro MR measurements demonstrates that the hydrodynamics in the fluid filled coaxial elastic tube system are similar to those present in a single patient with syringomyelia. The presented in vitro study of spinal cord wall motion, and complex unsteady pressure and flow environment within the syrinx and SAS, provides insight into the complex biomechanical forces present in syringomyelia.

The Origins of Syringomyelia: Numerical Models of Fluid/Structure Interactions in the Spinal Cord

A two-dimensional axi-symmetric numerical model is constructed of the spinal cord, consisting of elastic cord tissue surrounded by aqueous cerebrospinal fluid, in turn surrounded by elastic dura. The geometric and elastic parameters are simplified but of realistic order, compared with existing measurements. A distal reflecting site models scar tissue formed by earlier trauma to the cord, which is commonly associated with syrinx formation. Transients equivalent to both arterial pulsation and percussive coughing are used to excite wave propagation. Propagation is investigated in this model and one with a central canal down the middle of the cord tissue, and in further idealized versions of it, including a model with no cord, one with a rigid cord, one with a rigid dura, and a double-length untapered variant of the rigid-dura model. Analytical predictions for axial and radial wave-speeds in these different situations are compared with, and used to explain, the numerical outcomes. We find that the anatomic circumstances of the spinal cerebrospinal fluid cavity probably do not allow for significant wave steepening phenomena. The results indicate that wave propagation in the real cord is set by the elastic properties of both the cord tissue and the confining dura mater, fat, and bone. The central canal does not influence the wave propagation significantly.