Although cylindrical coordinates are useful, transforming equations from Cartesian to cylindrical coordinates often results in some 'hairy' expressions. For our purposes today, we'll need the definitions for the velocity potential function and the streamfunction in this new coordinate system:
Velocity Potential: vr = d/dr(PHI), vtheta = d/dtheta(PHI/r), vz = d/dz(PHI)
Streamfunction: vr = d/dtheta(PSI/r), vtheta = -d/dr(PSI)
Next, we'll apply these relations to a few potential flow scenarios.
Source and Sink
A source/sink is a flow radially outward/inward from/to a point at some volumetric flowrate, m. We'll define the velocity components for this type of flow as:
For a source type of flow, m > 0; for a sink, m < 0. We should also note that at r = 0, the velocity becomes infinite, which is physically impossible. "Wait," you may be asking yourself, "If this flow is physically impossible, why do we care?" The answer to this question is that you'll have to wait till next week to find out. (Don't you just love cliff-hangers?) Let's figure out the velocity potential and stream functions that govern a source (realizing that if we simply make m negative, we'll also have the functions that govern a sink):
It seems odd that we could consider a vortex to be irrotational; however, we must remember that rotation refers to the orientation of a fluid element and not the path followed by the element. An irrotational vortex is called a free vortex, while a rotational vortex is called a forced vortex. (Recall that the velocity potential function is only valid for irrotational flows.) A combined vortex is one with a forced vortex at its core and a free vortex outside its core.
The property of circulation is often associated with vortex motion. Circulation, GAMMA, is defined as the line integral of the tangential component of the velocity taken around a closed curve in the flow field. The velocity potential and streamfunctions for free vortices are often written in terms of the circulation of the flow*.
For a free vortex, we will define the following velocity components:
where K is the magnitude of the velocity.
Cool, huh? Remember that if we plot constant values of PHI and PSI, we will obtain a flow net that helps us visualize the flow. Next week, we'll do some fun things with the results from today and last week.
*PHI = (0.5*GAMMA*theta)/pi and PSI = -(0.5*GAMMA*ln(r))/pi