Consider two stationary plane walls with a constant volume flow rate is injected/sucked at the point of intersection of plane walls and let the angle subtended by two walls be. Take the cylindrical coordinate system with representing point of intersection and the centerline and are the corresponding velocity components. The resulting flow is two-dimensional if the plates are infinitely long in the axial direction, or the plates are longer but finite, if one were neglect edge effects and for the same reason the flow can be assumed to be entirely radial i.e.,. Then the continuity equation and the incompressible Navier–Stokes equations reduce to The boundary conditions are no-slip condition at both walls and the third condition is derived from the fact that the volume flux injected/sucked at the point of intersection is constant across a surface at any radius.
Formulation
The first equation tells that is just function of, the function is defined as Different authors defines the function differently, for example, Landau defines the function with a factor. But following Whitham, Rosenhead the momentum equation becomes Now letting the and momentum equations reduce to and substituting this into the previous equation results in Multiplying by and integrating once, where are constants to be determined from the boundary conditions. The above equation can be re-written conveniently with three other constants as roots of a cubic polynomial, with only two constants being arbitrary, the third constant is always obtained from other two because sum of the roots is. The boundary conditions reduce to where is the corresponding Reynolds number. The solution can be expressed in terms of elliptic functions. For convergent flow, the solution exists for all, but for the divergent flow, the solution exists only for a particular range of.
The equation takes the same form as an undamped nonlinear oscillator one can pretend that is time, is displacement and is velocity of a particle with unit mass, then the equation represents the energy equation with zero total energy, then it is easy to see that the potential energy is where in motion. Since the particle starts at for and ends at for, there are two cases to be considered.
First case are complex conjugates and. The particle starts at with finite positive velocity and attains where its velocity is and acceleration is and returns to at final time. The particle motion represents pure outflow motion because and also it is symmetric about.
Second case , all constants are real. The motion from to to represents a pure symmetric outflow as in the previous case. And the motion to to with for all time represents a pure symmetric inflow. But also, the particle may oscillate between, representing both inflow and outflow regions and the flow is no longer need to symmetric about.
The rich structure of this dynamical interpretation can be found in Rosenhead.
Pure outflow
For pure outflow, since at, integration of governing equation gives and the boundary conditions becomes The equations can be simplified by standard transformations given for example in Jeffreys.
The limiting condition is obtained by noting that pure outflow is impossible when, which implies from the governing equation. Thus beyond this critical conditions, no solution exists. The critical angle is given by where where is the complete elliptic integral of the first kind. For large values of, the critical angle becomes. The corresponding critical Reynolds number or volume flux is given by where is the complete elliptic integral of the second kind. For large values of, the critical Reynolds number or volume flux becomes.
Pure inflow
For pure inflow, the implicit solution is given by and the boundary conditions becomes Pure inflow is possible only when all constants are real and the solution is given by where is the complete elliptic integral of the first kind.
Limiting form
As Reynolds number increases, the flow tends to become uniform, except forboundary layers near the walls. Since is large and is given, it is clear from the solution that must be large, therefore. But when,, the solution becomes It is clear that everywhere except in the boundary layer of thickness. The volume flux is so that and the boundary layers have classical thickness.
