Computerized biomechanical simulation of cerebrospinal fluid hydrodynamics: Challenges and opportunities

Feasibility of assessing non-invasive intracranial compliance using FSI simulation-based and MR elastography-based brain stiffness

Intracranial compliance (ICC) refers to the change in intracranial volume per unit change in intracranial pressure (ICP). Magnetic resonance elastography (MRE) quantifies brain stiffness by measuring the shear modulus. Our objective is to investigate the relationship between ICC and brain stiffness through fluid-structure interaction (FSI) simulation, and to explore the feasibility of using MRE to assess ICC based on brain stiffness. This is invaluable due to the clinical importance of ICC, as well as the fast and non-invasive nature of the MRE procedure. We employed FSI simulation in hydrocephalus patients with aqueductal stenosis to non-invasively calculate ICP which is the basis of the calculation of ICC and FSI-based brain stiffness. The FSI simulated parameters used have been validated with experimental data. Our results showed that there is no relationship between FSI simulated-based brain stiffness and ICC in hydrocephalus patients. However, MRE-based brain stiffness may be sensitive to changes in intracranial fluid dynamic parameters such as cerebral perfusion pressure (CPP), cerebral blood flow (CBF), and ICP, as well as to mechano-vascular changes in the brain, which are determining parameters in ICC assessment. Although optimism has been found regarding the assessment of ICC using MRE-based brain stiffness, especially for acute-onset brain disorders, further studies are necessary to clarify their direct relationship.

Why Intracranial Compliance Is Not Utilized as a Common Practical Tool in Clinical Practice

Intracranial compliance (ICC) holds significant potential in neuromonitoring, serving as a diagnostic tool and contributing to the evaluation of treatment outcomes. Despite its comprehensive concept, which allows consideration of changes in both volume and intracranial pressure (ICP), ICC monitoring has not yet established itself as a standard component of medical care, unlike ICP monitoring. This review highlighted that the first challenge is the assessment of ICC values, because of the invasive nature of direct measurement, the time-consuming aspect of non-invasive calculation through computer simulations, and the inability to quantify ICC values in estimation methods. Addressing these challenges is crucial, and the development of a rapid, non-invasive computer simulation method could alleviate obstacles in quantifying ICC. Additionally, this review indicated the second challenge in the clinical application of ICC, which involves the dynamic and time-dependent nature of ICC. This was considered by introducing the concept of time elapsed (TE) in measuring the changes in volume or ICP in the ICC equation (volume change/ICP change). The choice of TE, whether short or long, directly influences the ICC values that must be considered in the clinical application of the ICC. Compensatory responses of the brain exhibit non-monotonic and variable changes in long TE assessments for certain disorders, contrasting with the mono-exponential pattern observed in short TE assessments. Furthermore, the recovery behavior of the brain undergoes changes during the treatment process of various brain disorders when exposed to short and long TE conditions. The review also highlighted differences in ICC values across brain disorders with various strain rates and loading durations on the brain, further emphasizing the dynamic nature of ICC for clinical application. The insight provided in this review may prove valuable to professionals in neurocritical care, neurology, and neurosurgery for standardizing ICC monitoring in practical application related to the diagnosis and evaluation of treatment outcomes in brain disorders.

A mathematical framework for the dynamic interaction of pulsatile blood, brain, and cerebrospinal fluid

BACKGROUND

Shedding light on less-known aspects of intracranial fluid dynamics may be helpful to understand the hydrocephalus mechanism. The present study suggests a mathematical framework based on in vivo inputs to compare the dynamic interaction of pulsatile blood, brain, and cerebrospinal fluid (CSF) between the healthy subject and the hydrocephalus patient.

METHOD

The input data for the mathematical formulations was pulsatile blood velocity, which was measured using cine PC-MRI. Tube law was used to transfer the created deformation by blood pulsation in the vessel circumference to the brain domain. The pulsatile deformation of brain tissue with respect to time was calculated and considered to be inlet velocity in the CSF domain. The governing equations in all three domains were continuity, Navier-Stokes, and concentration. We used Darcy law with defined permeability and diffusivity values to define the material properties in the brain.

RESULTS

We validated the preciseness of the CSF velocity and pressure through the mathematical formulations with cine PC-MRI velocity, experimental ICP, and FSI simulated velocity and pressure. We used the analysis of dimensionless numbers including Reynolds, Womersley, Hartmann, and Peclet to evaluate the characteristics of the intracranial fluid flow. In the mid-systole phase of a cardiac cycle, CSF velocity had the maximum value and CSF pressure had the minimum value. The maximum and amplitude of CSF pressure, as well as CSF stroke volume, were calculated and compared between the healthy subject and the hydrocephalus patient.

CONCLUSION

The present in vivo-based mathematical framework has the potential to gain insight into the less-known points in the physiological function of intracranial fluid dynamics and the hydrocephalus mechanism.

Long-term recovery behavior of brain tissue in hydrocephalus patients after shunting

The unpredictable complexities in hydrocephalus shunt outcomes may be related to the recovery behavior of brain tissue after shunting. The simulated cerebrospinal fluid (CSF) velocity and intracranial pressure (ICP) over 15 months after shunting were validated by experimental data. The mean strain and creep of the brain had notable changes after shunting and their trends were monotonic. The highest stiffness of the hydrocephalic brain was in the first consolidation phase (between pre-shunting to 1 month after shunting). The viscous component overcame and damped the input load in the third consolidation phase (after the fifteenth month) and changes in brain volume were stopped. The long-intracranial elastance (long-IE) changed oscillatory after shunting and there was not a linear relationship between long-IE and ICP. We showed the long-term effect of the viscous component on brain recovery behavior of hydrocephalic brain. The results shed light on the brain recovery mechanism after shunting and the mechanisms for shunt failure.

Hydrodynamic comparison of shunt and endoscopic third ventriculostomy in adult hydrocephalus using in vitro models and fluid-structure interaction simulation

OBJECTIVE

The comparison of the efficiency of shunt placement and endoscopic third ventriculostomy (ETV) in treating of adult hydrocephalus patients with various intensities and different obstruction intensities in the aqueduct of Sylvius (AS).

METHODS

In vitro models with separated ventricles were simulated and implemented for modeling shunt and ETV surgeries in one healthy subject and hydrocephalus patients with various intensities, as well as three different obstruction intensities in AS and under two cerebrospinal fluid (CSF) dynamic conditions. The fluid-structure interaction simulation was also carried out to validate in vitro results.

RESULTS

The efficiency of both methods in reducing the maximum CSF pressure in the subarachnoid space (MCPS) decreased by an increase in the patient's intensities. Contrary to shunting, the efficiency of ETV in reducing MCPS demonstrated a decline (8.3-16.4%) by an increase in obstruction levels in AS. Based on the findings, shunt efficiency in decreasing MCPS in patients with low intensity was more remarkable compared to ETV. However, ETV was more efficient than shunt in the patient with intracranial hypertension. Further, shunt placement and ETV led to a significant reduction in the amplitude of CSF pressure in the SAS (ACPS) in patients with sneezing, coughing, Valsalva maneuver, and exercising effects in contrast to other patients. Moreover, ACPS reduction was not related to the intensity of the disease in both treatment methods. In contrast to shunt, an increase in the obstruction level in AS led to a reduction in ACPS in ETV in both CSF dynamic conditions.

CONCLUSIONS

The noises from irregular disorders increased the discharging of CSF after shunt placement, and activities such as sneezing, coughing, Valsalva maneuvers, and exercising increased the risk of shunt overdrainage by 10.4~47.8%, especially in the patient with intracranial hypertension. Based on the proposed in vitro ETV and shunt models, an increase of head compliance was higher in ETV compared to the shunt. Eventually, an increase in the obstruction level of AS after ETV led to a decline in head compliance in contrast to shunt.

Boundary conditions investigation to improve computer simulation of cerebrospinal fluid dynamics in hydrocephalus patients

Three-D head geometrical models of eight healthy subjects and 11 hydrocephalus patients were built using their CINE phase-contrast MRI data and used for computer simulations under three different inlet/outlet boundary conditions (BCs). The maximum cerebrospinal fluid (CSF) pressure and the ventricular system volume were more effective and accurate than the other parameters in evaluating the patients' conditions. In constant CSF pressure, the computational patient models were 18.5% more sensitive to CSF volume changes in the ventricular system under BC "C". Pulsatile CSF flow rate diagrams were used for inlet and outlet BCs of BC "C". BC "C" was suggested to evaluate the intracranial compliance of the hydrocephalus patients. The results suggested using the computational fluid dynamic (CFD) method and the fully coupled fluid-structure interaction (FSI) method for the CSF dynamic analysis in patients with external and internal hydrocephalus, respectively.