Fluid dynamics in syringomyelia cavities: Effects of heart rate, CSF velocity, CSF velocity waveform and craniovertebral decompression.

Purpose How fluid moves during the cardiac cycle within a syrinx may affect its development. We measured syrinx fluid velocities before and after craniovertebral decompression in a patient and simulated syrinx fluid velocities for different heart rates, syrinx sizes and cerebrospinal fluid (CSF) flow velocities in a model of syringomyelia. Materials and methods With phase-contrast magnetic resonance we measured CSF and syrinx fluid velocities in a Chiari patient before and after craniovertebral decompression. With an idealized two-dimensional model of the subarachnoid space (SAS), cord and syrinx, we simulated fluid movement in the SAS and syrinx with the Navier-Stokes equations for different heart rates, inlet velocities and syrinx diameters. Results In the patient, fluid oscillated in the syrinx at 200 to 210 cycles per minute before and after craniovertebral decompression. Velocities peaked at 3.6 and 2.0 cm per second respectively in the SAS and the syrinx before surgery and at 2.7 and 1.5 cm per second after surgery. In the model, syrinx velocity varied between 0.91 and 12.70 cm per second. Increasing CSF inlet velocities from 1.56 to 4.69 cm per second increased peak syrinx fluid velocities in the syrinx by 151% to 299% for the three cycle rates. Increasing cycle rates from 60 to 120 cpm increased peak syrinx velocities by 160% to 312% for the three inlet velocities. Peak velocities changed inconsistently with syrinx size. Conclusions CSF velocity, heart rate and syrinx diameter affect syrinx fluid velocities, but not the frequency of syrinx fluid oscillation. Craniovertebral decompression decreases both CSF and syrinx fluid velocities.

Numerical study of intrathecal drug delivery to a permeable spinal cord: effect of catheter position and angle.

Intrathecal delivery is a procedure involving the release of therapeutic agents into the cerebrospinal fluid (CSF) hrough a catheter. It holds promise for treating high-impact central nervous system pathologies, for which systemic administration routes are ineffective. In this study we introduce a numerical model able to simultaneously account for solute transport in the fluid and in the spinal cord. Using a Discontinuous Galerkin method and a three-dimensional patient-specific geometry, we studied the effect of catheter position and angle on local spinal cord drug concentration. We considered twenty cardiac cycles to limit the computational cost of our approach, which resolves the physics both in space and time. We used clinically representative data for the drug injection speed and dose rate, and scaled drug diffusion/penetration properties to obtain observable effects during the considered simulation time. Based on our limited set of working parameters, lateral injection perpendicular to the cord turned out to be more effective than other configurations. Even if the adopted scaling does not allow for a direct clinical translation (a wider parametric assessment of the importance of CSF flow, geometry and diffusion properties is needed), it did not weaken our numerical approach, which can be used to systematically investigate multiple catheter, geometry and fluid/tissue properties configurations, thus paving the way for therapy control.

Multigrid block preconditioning for a coupled system of partial differential equations modeling the electrical activity in the heart.

The electrical activity of the heart may be modeled with a system of partial differential equations (PDEs) known as the bidomain model. Computer simulations based on these equations may become a helpful tool to understand the relationship between changes in the electrical field and various heart diseases. Because of the rapid variations in the electrical field, sufficiently accurate simulations require a fine-scale discretization of the equations. For realistic geometries this leads to a large number of grid points and consequently large linear systems to be solved for each time step. In this paper, we present a fully coupled discretization of the bidomain model, leading to a block structured linear system. We take advantage of the block structure to construct an efficient preconditioner for the linear system, by combining multigrid with an operator splitting technique.

Spinal fluid biomechanics and imaging: an update for neuroradiologists.

Flow imaging with cardiac-gated phase-contrast MR has applications in the management of neurologic disorders. Together with computational fluid dynamics, phase-contrast MR has advanced our understanding of spinal CSF flow. Phase-contrast MR is used to evaluate patients with Chiari I malformation who are candidates for surgical treatment. In theory, abnormal CSF flow resulting from the abnormal tonsil position causes syringomyelia and other neurologic signs and symptoms in patients with Chiari I. CSF flow imaging also has research applications in syringomyelia and spinal stenosis. To optimize MR acquisition and interpretation, neuroradiologists must have familiarity with healthy and pathologic patterns of CSF flow. The purpose of this review is to update concepts of CSF flow that are important for the practice of flow imaging in the spine.

Simulating CSF flow dynamics in the normal and the Chiari I subarachnoid space during rest and exertion.

BACKGROUND AND PURPOSE: CSF fluid dynamics in healthy subjects and patients with Chiari I have been characterized during rest with phase-contrast MR imaging and CFD. CSF flow velocities and pressures in the nonresting state have not been adequately characterized. We used computer simulations to study CSF dynamics during increased heart rates in the normal and Chiari I subarachnoid space. MATERIALS AND METHODS: Cyclic CSF flow was simulated for multiple cycles in idealized 3D models of the subarachnoid space for normal and Chiari I malformation subarachnoid spaces, with flow cycles corresponding to 80 or 120 heart beats per minute. Flow velocities and pressures were computed by the Navier-Stokes equations. Synchronous bidirectional flow and flow patterns were displayed in Star-CD and inspected visually. Peak velocities and pressure differences in the 2 models were compared for the 2-cycle frequencies. RESULTS: Elevating the cycle rate from 80 to 120 cpm increased peak superior-inferior pressure gradients (top-bottom) by just 0.01% in the normal model and 2% in the Chiari model. Corresponding average pressure gradients increased by 92% and 100%, respectively. In addition, in both models, the range of synchronous bidirectional flow velocities increased. Systolic velocities had smaller increases with faster cycling. For each cycle rate, peak and average pressure gradients in the Chiari model were greater than in the normal model by 11%-16%. CONCLUSIONS: Raising the cycle rate from 80 to 120 cpm increased superior-inferior average pressure gradients and the range of synchronous bidirectional flow velocities in the normal and Chiari I models.

Estimation of CSF flow resistance in the upper cervical spine.

Chiari I patients have increased CSF velocities in the foramen magnum due hypothetically to increased pressure gradients or reduced flow resistance. We calculated flow resistance in the cervical spinal canal in a group of subjects with and without the Chiari malformation. Eight subjects including healthy volunteers and Chiari I patients were studied. From 3D high resolution MR images of the cervical spine mathematical models of the subarachnoid spaces were created by means of standard programs for segmentation and discretization. Oscillatory flow through the subarachnoid space was simulated. Cross-sectional area of the subarachnoid space was computed at each level from C1 through C4 and the length of this spinal canal segment was measured. Peak caudad CSF flow velocity at each level was plotted against cross-section area. CSF volumetric flux and resistance were calculated for each subject. The correlation between velocity and resistance was calculated. In all subjects, peak velocities increased progressively from C1 to C4 by 0.6 to 0.7 cm/s per level. Spinal canal areas diminished from C1 to C5 in each subject at a rate of -0.25 to -0.29 cm(2) per level. Resistance averaged 4.3 pascal/ml/s in the eight subjects; 3.8 pascal/ml/s in patients with tonsilar herniation and 6.0 pascal/ml/s in volunteers. Velocity correlated inversely with resistance (R(2) = 0.6). CSF velocities correlated inversely with the flow resistance in the upper cervical spinal canal. Resistance tends to be lower in Chiari I patients than in healthy volunteers.

CSF pressure and velocity in obstructions of the subarachnoid spaces.

According to some theories, obstruction of CSF flow produces a pressure drop in the subarachnoid space in accordance with the Bernoulli theorem that explains the development of syringomyelia below the obstruction. However, Bernoulli's principle applies to inviscid stationary flow unlike CSF flow. Therefore, we performed a series of computational experiments to investigate the relationship between pressure drop, flow velocities, and obstructions under physiologic conditions. We created geometric models with dimensions approximating the spinal subarachnoid space with varying degrees of obstruction. Pressures and velocities for constant and oscillatory flow of a viscid fluid were calculated with the Navier-Stokes equations. Pressure and velocity along the length of the models were also calculated by the Bernoulli equation and compared with the results from the Navier-Stokes equations. In the models, fluid velocities and pressure gradients were approximately inversely proportional to the percentage of the channel that remained open. Pressure gradients increased minimally with 35% obstruction and with factors 1.4, 2.2 and 5.0 respectively with 60, 75 and 85% obstruction. Bernoulli's law underestimated pressure changes by at least a factor 2 and predicted a pressure increase downstream of the obstruction, which does not occur. For oscillatory flow the phase difference between pressure maxima and velocity maxima changed with the degree of obstruction. Inertia and viscosity which are not factored into the Bernoulli equation affect CSF flow. Obstruction of CSF flow in the cervical spinal canal increases pressure gradients and velocities and decreases the phase lag between pressure and velocity.

Patient-specific 3D simulation of cyclic CSF flow at the craniocervical region.

BACKGROUND AND PURPOSE: Flow simulations in patient-specific models of the subarachnoid space characterize CSF flow in more detail than MR flow imaging. We extended previous simulation studies by including cyclic CSF flow and patient-specific models in multiple patients with Chiari I. We compared simulation results with MR flow measurements. MATERIALS AND METHODS: Volumetric high resolution image sets acquired in 7 patients with Chiari I, 3 patients who had previous craniovertebral decompression, and 3 controls were segmented and converted to mathematical models of the subarachnoid space. CSF flow velocities and pressures were calculated with high spatial and temporal resolution during simulated oscillatory flow in each model with the Navier-Stokes equations. Pressures, velocities, and bidirectional flow were compared in the groups (with Student t test). Peak velocities in the simulations were compared with peak velocities measured in vivo with PCMR. RESULTS: Flow visualization for patients and volunteers demonstrated nonuniform reversing patterns resembling those observed with PCMR. Velocities in the 13 subjects were greater between C2 and C5 than in the foramen magnum. Chiari patients had significantly greater peak systolic and diastolic velocities, synchronous bidirectional flow, and pressure gradients than controls. Peak velocities measured in PCMR correlated significantly (P = .003; regression analysis) despite differences between them. CONCLUSIONS: In simulations of CSF, patients with Chiari I had significantly greater peak systolic and diastolic velocities, synchronous bidirectional flow, and pressure gradients than controls.

CSF Flow in Chiari I and Syringomyelia from the Perspective of Computational Fluid Dynamics.

Phase contrast MR in patients with the Chiari I malformation demonstrates abnormal CSF flow in the foramen magnum and upper cervical spinal canal, related to abnormal pressure gradients. The purpose of this study was to analyze the role of CSF pressure in the pathogenesis of syringomyelia, with computational models. The spinal cord was modeled as a cylindrical poro-elastic structure with homogenous and isotropic permeability. The permeability was then made heterogeneous and anisotropic to represent the different properties of the central canal, gray and white matter. Fluid with a defined pressure, varying both in time and space, was prescribed in the SAS. Simulations were performed to quantify deformations and fluid movement within the cord. In the simulations with uniform permeability fluid moved into the cord in regions of higher pressure and out of the cord in regions of lower pressure. With permeability differences simulating gray and white matter the pattern was more complex, but similar. Adding the central spinal canal, fluid moved into the cord as in the previous case. However, preferential flow along the central canal hindered fluid from flowing back into the SAS. Pressure gradients in the SAS produce movement of fluid in the spinal cord. Assuming different relative permeability in gray matter, white matter and the central spinal canal, abnormal CSF gradients lead to accumulation of fluid within and adjacent to the spinal cord central canal.

Effect of tonsillar herniation on cyclic CSF flow studied with computational flow analysis.

BACKGROUND AND PURPOSE: The Chiari I malformation, characterized by tonsils extending below the foramen magnum, has increased CSF velocities compared with those in healthy subjects. Measuring the effect of tonsillar herniation on CSF flow in humans is confounded by interindividual variation. The goal of this study was to determine the effect of herniated tonsils on flow velocity and pressure dynamics by using 3D computational models. MATERIALS AND METHODS: A previously described 3D mathematic model of the normal subarachnoid space was modified by extending the tonsils inferiorly. The chamber created was compared with the anatomy of the subarachnoid space. Pressures and velocities were calculated by CFA methods for sinusoidal flow of a Newtonian fluid. Results were displayed as 2D color-coded plots and 3D animations. Pressure gradients and flow velocities were compared with those in the normal model. Velocity distributions were also compared with those in clinical images of CSF flow. RESULTS: The model represented grossly the subarachnoid space of a patient with Chiari I malformation. Fluid flow patterns in the Chiari model were complex, with jets in some locations and stagnant flow in others. Flow jets, synchronous bidirectional flow, and pressure gradients were greater in the Chiari model than in the normal model. The distribution of flow velocities in the model corresponded well with those observed in clinical images of CSF flow in patients with Chiari I. CONCLUSIONS: Tonsillar herniation per se increases the pressure gradients and the complexity of flow patterns associated with oscillatory CSF flow.

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.

CSF flow dynamics at the craniovertebral junction studied with an idealized model of the subarachnoid space and computational flow analysis.

BACKGROUND AND PURPOSE: How CSF flow varies with the anatomy of the subarachnoid space has not been sufficiently well studied. The goal of this study was to develop an idealized 3D computational model of the subarachnoid space and then to use this model to study the detailed spatiotemporal effects of anatomic variations on CSF pressures and velocities. MATERIALS AND METHODS: We created a geometric model with a computer-assisted design program. The model contained a central structure for the brain and spinal cord axis and a second surrounding structure for the peripheral borders of the subarachnoid space. Model dimensions were adjusted to capture the main characteristics of the normal human posterior fossa and cervical spinal anatomy. CSF flow was modeled as water with a sinusoidal flow pattern in time. Velocities and pressures during craniocaudal and caudocranial flow were calculated with computational fluid dynamics (CFD) software. Simulated flow was compared with published phase-contrast MR imaging measurements of CSF flow in healthy human subjects. RESULTS: The model contained geometric characteristics of the posterior fossa and spinal canal. Flow velocities varied with the time in the cycle and location in space. Flow velocities had spatial variations that resembled those in healthy human subjects. Reynolds numbers were moderate, showing a laminar flow regime. Pressure varied uniformly along the long axis of the model during craniocaudal and caudocranial flow. CONCLUSIONS: In an idealized geometric approximation of the human subarachnoid space, CSF velocities and pressures can be studied in spatiotemporal detail with mathematic models.