Pulsatile Flow in a Circular Tube

The velocity distribution, , is computed numerically using NDSolve for a pulsatile pressure-driven 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 (i.e. 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 long-time 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 quasi-steady 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.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


The oscillating pressure-driven flow in a tube obeys the following equations:
, ,
, ,
, ,
, , , and ,
where , , and are the dimensionless velocity, radial position, and time.
The variable can be considered as a Reynolds number, since it appears as the ratio of inertial forces to viscous forces. There are two characteristic times for this problem: , the period of the imposed pressure gradient, and , the time for diffusion of momentum across the tube.
[1] L. G. Leal, Laminar Flow and Convective Transport Processes, Boston: Butterworth–Heinemann, 1992.
[2] L. G. Leal, Advanced Transport Phenomena: Fluid Mechanics and Convective Transport Processes, Cambridge: Cambridge University Press, 2007.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.