Geology Reference
In-Depth Information
This derivation extends the classic derivation (Ashley and Landahl, 1965) for
the streamfunction < to problems with non-zero /. If we now differentiate
Equation
7.
3.15c with respect to y and Equation
7.
3.15d with respect to x, and
eliminate the velocity potential, we obtain the governing equation
<
xx
+ <
yy
= 0 (7.3.15e)
This takes the same form as Equation
7
.3.15a so that the same solution
algorithm applies. The boundary conditions used are obtained from Equations
7
.3.15c and
7
.3.15d since ) is already known. The normal Neumann derivative
<
y
is applied along the horizontal upper and lower box boundaries, while the
normal derivative <
x
is applied along the vertical left and right boundaries.
How is < used to trace streamlines? Note that streamline slope dy/dx is
kinematically equal to the velocity ratio )
y
/)
x
, that is,
dy/dx = )
y
/)
x
= - <
x
/<
y
+ /y/<
y
(7.3.15f)
Then the total differential satisfies
d< = <
x
dx + <
y
dy = /y dx (7.3.15g)
When annular convergences are small, / can be neglected so that d< = 0. Thus,
< is constant along a streamline. In this limit, streamlines are easily constructed
by using a contour plotter for the converged numerical field <(x,y).
Figure 7.12.
Periodic cascade of upstream and downstream lobes.
Force
Figure 7.13.
Periodic boundary value problem flow domain.
Search WWH ::
Custom Search