Numerical model of CSF flow and transport in the spinal canal

Among other things, the files allow the calculation of solute transport with densities different from that of CSF, for a relevant Richardson number of order one. This model replicates the numerical simulations performed in the project using canonical configurations and also forms the basis for simulations in more realistic anatomical and physiological configurations.

The model presents a constant eccentricity between the dura and pia membranes equal to beta=0.5, computed as the distance between axes divided by the mean characteristic radial gap between the membranes. The simulation is implemented for a Womersley number equal to 3, computed as the square root of the mean characteristic radial gap power to two times the angular frequency divided by the kinematic viscosity, and a solute of molecular such that the Schmidt number is equal to 1000.