The velocity distribution is computed numerically using NDSolve for a pulsatile pressuredriven flow in a tube This model considerably simplifies the actual flow through veins and arteries When the kinematic viscosity is large and momentum diffusion is fast ie the Strouhal number is small as shown in the first snapshot the velocity at the center of the tube shown in red can instantaneously adjust to changing pressure shown in blue The velocity is then in phase with the pressure gradient On the other hand when the kinematic viscosity is small and momentum diffusion is slow then is large see the second snapshot and the velocity at the center of the tube lags the pressure gradient by 90° The longtime analytical solution given in Details and shown in green on the plot for matches the numerical solution very well after a short transient period This derivation can be found in 1–2 One can see that when everything is in phase then the flow behaves locally as a quasisteady Poisueille flow in which the green and blue curves coincide the green curve being the parabolic solution and the blue curve taking inertial effects into account At a high Strouhal number inertial effects become significant and the velocity field is no longer locally parabolic green and blue curves are distinct Both limiting cases are shown in the last two snapshots