Svein Linge

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.

Numerical solution of the bidomain equations

Knowledge of cardiac electrophysiology is efficiently formulated in terms of mathematical models. However, most of these models are very complex and thus defeat direct mathematical reasoning founded on classical and analytical considerations. This is particularly so for the celebrated bidomain model that was developed almost 40 years ago for the concurrent analysis of extra- and intracellular electrical activity. Numerical simulations based on this model represent an indispensable tool for studying electrophysiology. However, complex mathematical models, steep gradients in the solutions and complicated geometries lead to extremely challenging computational problems. The greatest achievement in scientific computing over the past 50 years has been to enable the solving of linear systems of algebraic equations that arise from discretizations of partial differential equations in an optimal manner, i.e. such that the central processing unit (CPU) effort increases linearly with the number of computational nodes. Over the past decade, such optimal methods have been introduced in the simulation of electrophysiology. This development, together with the development of affordable parallel computers, has enabled the solution of the bidomain model combined with accurate cellular models, on geometries resembling a human heart. However, in spite of recent progress, the full potential of modern computational methods has yet to be exploited for the solution of the bidomain model. This paper reviews the development of numerical methods for solving the bidomain model. However, the field is huge and we thus restrict our focus to developments that have been made since the year 2000.

Characterization of Cyclic CSF Flow in the Foramen Magnum and Upper Cervical Spinal Canal with MR Flow Imaging and Computational Fluid Dynamics

SUMMARY: 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.

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.

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.

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.

Synchronizing computer simulations with measurement data for a case of atrial flutter.

Atrial flutter is a common supraventricular tachycardia that can be treated using radiofrequency catheter ablation, a procedure that is guided by electroanatomical mapping systems. In this paper, we propose an algorithm for incorporating mapping data into computer simulations of atrial electrical activity with the purpose of creating a more accurate map of electrical activation. The algorithm takes as input the extracellular potential values recorded at a number of sites throughout the atria and estimates the activation time for the entire atrial domain. We test the algorithm using synthetic mapping data and an anatomically detailed atrial geometry with an activation pattern typical of atrial flutter. The results show that the algorithm performs well with synthetic mapping data with information from relatively few mapping sites and in the presence of modeling and measurement error.

Effect of craniovertebral decompression on CSF dynamics in Chiari malformation type I studied with computational fluid dynamics: Laboratory investigation.

OBJECT The effect of craniovertebral decompression surgery on CSF flow dynamics in patients with Chiari malformation Type I (CM-I) has been incompletely characterized. The authors used computational fluid dynamics to calculate the effect of decompression surgery on CSF flow dynamics in the posterior fossa and upper cervical spinal canal. METHODS Oscillatory flow was simulated in idealized 3D models of the normal adult and the CM-I subarachnoid spaces (both previously described) and in 3 models of CM-I post-craniovertebral decompressions. The 3 postoperative models were created from the CM model by virtually modifying the CM model subarachnoid space to simulate surgical decompressions of different magnitudes. Velocities and pressures were computed with the Navier-Stokes equations in Star-CD for multiple cycles of CSF flow oscillating at 80 cycles/min. Pressure gradients and velocities were compared for 8 levels extending from the posterior fossa to the C3-4 level. Relative pressures and peak velocities were plotted by level from the posterior fossa to C3-4. The heterogeneity of flow velocity distribution around the spinal cord was compared between models. RESULTS Peak systolic velocities were generally lower in the postoperative models than in the preoperative CM model. With the 2 larger surgical defects, peak systolic velocities were brought closer to normal model velocities (equal values at C-3 and C-4) than with the smallest surgical defect. For the smallest defect, peak velocities were decreased, but not to levels in the normal model. In the postoperative models, heterogeneity in flow velocity distribution around the spinal cord increased from normal model levels as the degree of decompression increased. Pressures in the 5 models differed in magnitude and in pattern. Pressure gradients along the spinal canal in the normal and CM models were nonlinear, with steeper gradients below C3-4 than above. The CM model had a steeper pressure gradient than the normal model above C3-4 and the same gradient below. The postoperative models had lower pressure gradients than the CM model above C2-3. The most conservative decompression had lower pressure gradients than the normal model above C2-3. The two larger decompression defects had CSF pressure gradients below those in the normal model above C2-3. These 2 models had a less steep gradient above C-3 and a steeper gradient below. CONCLUSIONS In computer simulations, craniovertebral surgical defects generally diminished CSF velocities and CSF pressures.

On the frequency of automaticity during ischemia in simulations based on stochastic perturbations of the Luo-Rudy 1 model

AIMS To compute the effects of parameter perturbations for single ischemic cardiac cells, and to determine how perturbations influenced the tendency for the cells to undergo spontaneous depolarization (automaticity) during 20 min of acute ischemia. METHODS A modified Luo-Rudy 1 cell model was used. Since the range of biological variation and measurement errors is largely unknown, we conducted our study of the consequences of perturbations under the assumption that cell model parameters have a normal distribution with a 10% standard deviation. A total of 10000 random cell realizations were tested while varying important Luo-Rudy cell model parameters. Ischemia was modelled by deterministic functions chosen for the expected values of crucial ion concentrations and gating parameters as they developed with time, while realizing the respective parameter values from static normal distributions with a 10% standard deviation. RESULTS AND CONCLUSION It was found that the tendency towards automaticity did increase as the stochastic parameters were varied. In particular, cells with standard Luo-Rudy parameter values did not become automatic during ischemia, whereas a significant portion of the cells with randomized parameter values did. The relative importance of model parameter variations was also determined and a sodium m-gate activation parameter was identified as the most critical parameter. The frequency of arrhythmic events during acute ischemia is known to be bell-shaped, with a peak at around 7-8 min after the onset of ischemia. Our simulations display a similar peak in the frequency of automaticity.

Solving the heart mechanics equations with Newton and quasi Newton methods–a comparison

The non-linear elasticity equations of heart mechanics are solved while emulating the effects of a propagating activation wave. The dynamics of a 1 cm3 slab of active cardiac tissue was simulated as the electrical wave traversed the muscular heart wall transmurally. The regular Newton (Newton–Raphson) method was compared to two modified Newton approaches, and also to a third approach that delayed update only of some selected Jacobian elements. In addition, the impact of changing the time step (0.01, 0.1 and 1 ms) and the relative non-linear convergence tolerance (10−4, 10−3 and 10−2) was investigated. Updating the Jacobian only when slow convergence occured was by far the most efficient approach, giving time savings of 83–96%. For each of the four methods, CPU times were reduced by 48–90% when the time step was increased by a factor 10. Increasing the convergence tolerance by the same factor gave time savings of 3–71%. Different combinations of activation wave speed, stress rate and bulk modulus revealed that the fastest method became relatively even faster as stress rate and bulk modulus was decreased, while the activation speed had negligible influence in this respect.