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

• Published in: Computer methods and programs in biomedicine Vol. 231; p. 107209
• Main Authors: Gholampour, Seifollah, Balasundaram, Hemalatha, Thiyagarajan, Padmavathi, Droessler, Julie
• Format: Journal Article
• Language: English
• Published: Ireland Elsevier B.V 01.04.2023
• Subjects: Blood dynamics; Brain; Brain - physiology; Cerebrospinal fluid dynamics; Humans; Hydrocephalus; Hydrocephalus - blood; Hydrocephalus - cerebrospinal fluid; Hydrocephalus - diagnostic imaging; Hydrodynamics; Intracranial pressure; Magnetic Resonance Imaging; Mathematical framework; Pulsatile interaction; Stroke - blood; Stroke - cerebrospinal fluid; Brain; Cerebrospinal fluid dynamics; Blood dynamics; Hydrocephalus; Pulsatile interaction; Intracranial pressure; Mathematical framework
• ISSN: 0169-2607; 1872-7565; 1872-7565
• DOI: 10.1016/j.cmpb.2022.107209

•Reynolds and Hartmann in blood flow had no relationship with velocity change.•Peclet number in CSF and blood flow changes directly with velocity/concentration.•Womersley and Peclet in brain fluid change inversely with velocity/concentration.•The phase lag between CSF velocity and pressure diagrams is 57º. 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. 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. 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. 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.