r sqrt(z^2+y^2) "radial coordinate of vessel"
u u_max*(1-r^2/R_vessel^2) "Hagen-Poiseuille velocity"
epsl epsl_blood*max(1-phils,1e-3) "Plasma volume fraction"