Problem
Derive a laminar flow profile for an incompressible fluid with density and viscosity flowing in a long pipe with a radius in steady conditions. Also, derive formulas for flow rate, maximal and average velocity!
Assumptions:
– the flow is laminar and steady,
– the fluid is incompressible,
– the flow is developed (i.e. the pipe is very long and we do not consider the section close to the inlet, )
– the pipe is very long so the flow is axisymmetric
The solution, in a couple of steps:
We tackle the problem in polar coordinate system. The origin is at the center of the pipe, with axis pointing towards the perimeter and the axis along the pipe. The problem is axisymmetric, so the axis is irrelevant. The flow velocity has two components , the component is the velocity along the pipe (along the axis) and the component is the velocity towards the perimeter of the pipe (along the axis).
The Navier-Stokes equations in this setting are: the continuity equation, the momentum transport equations in and in directions.
Consider the continuity equation:
The continuity equation for an incompressible fluid is:
The first term is zero, since the flow is developed, so there are no changes in the velocity along the pipe, i.e. . From the second term, we then see
As the flow can not penetrate the pipe wall, we know and so . The means that the constant in the above equation is zero, and finally we conclude that everywhere.
The momentum equation across the pipe
The steady state momentum equation in direction is:
Since this equations simplifies to
which means that the pressure is constant across the pipe and may only change along the pipe, i.e. .
The momentum equation along the pipe
The steady state momentum equation in direction is:
Taking into account and we simplify into
Since the left hand side of this equation depends solely on and the right hand side can only be a function of we conclude that the equation can only be satisfied if both sides are constant. This means that the pressure drop in the pipe is constant, i.e.
if is the pressure drop that occurs at a distance along the pipe. The minus sign comes from the fact that the pressure gradient is negative when pressure is decreasing in the direction of .
Solve the ordinary differential equation
Since we know, that the velocity depends only on , we can rewrite the steady state momentum equation in direction as a ordinary differential equation:
We integrate twice and introduce two integration constants: and as
At this point we consider boundary conditions to determine the integration constants. Due to axisymmetry, we know that the maximal velocity will be at the center of the pipe at . This must be true, as momentum diffuses from the fluid into the pipe wall and since it does that axisymmetrically it must be true that the velocity is largest in the center of the pipe. In this way, momentum diffuses from the centre in all directions towards the walls. The velocity has an extrema at the center of the pipe
To satisfy this, it follows that . The second boundary condition comes from the fact that there is no slip between the fluid and the pipe wall. The fluid at the pipe wall is stationary, i.e.
Since we already established that , we may work out that
The Hagen-Poiseuille flow profile
Using the expressions for the integration constant, we can work out the Hagen-Poiseuille flow profile
and the shear stress
The flow profile is parabolic with a maximum at the center of the pipe and velocity equal to zero at the pipe wall. The shear stress changes linearly from zero at the center of the pipe to the maximal value at the pipe wall.
The flow rate, the maximal velocity and the average velocity:
The maximal fluid velocity is at the center of the pipe, so:
The flow rate can be calculated by means of integration across the pipe.
Substituting and we obtain
The integral above evaluates to and we find an expression for the flow rate in a pipe:
We observe that the flow rate is linearly proportional to the pressure drop, scales inversely with fluid viscosity and is proportional to the radius to the power of four. This means, is one wants to pump the same amount of fluid per second in a pipe with half the diameter, one must provide 16 times larger pressure drop.
The average flow velocity is
It is interesting to note, that the maximal velocity is twice as large as the average velocity
The Darcy-Weisbach formula and the Darcy friction factor
Let us express the pressure drop in terms of the average velocity. Using the pipe diameter and introducing the flow Reynolds number as we may write:
A final rearrangement gives
Here, the factor is the Darcy friction factor. Expressing the pressure drop with head loss, i.e. , we can write
which is known as the Darcy-Weisbach formula.