Geoscience Reference
In-Depth Information
Given that the pressure in such a steady flow will be hydrostatic (
Equation 3.21
), the
pressure gradient can be written as:
0
@
1
A
¼
@
@
0
@
1
A
ð
z
ð
ð
0
0
z
@
p
x
¼
@
þ
gdz
gdz
gdz
@
@
x
x
ð
9
:
2
Þ
ð
0
@
@
þ
s
g
@
@
gz
@
@
x
þ
s
g
@
¼
x
gdz
x
¼
@
x
z
where, in the final step, the depth-uniform
x has been taken outside the integral
and r
s
is the density at the surface. Assuming (
Equation 3.40
) that the stress is related
to the velocity shear by t
x
¼
@
r/
@
N
z
0
@
u
z
where N
z
is the eddy viscosity, the force
@
balance can now be written as:
:
@
p
gz
@
@
x
þ
s
g
@
x
¼
@
t
x
@
z
¼
0
@
N
z
@
u
x
¼
ð
9
:
3
Þ
@
@
@
z
@
z
If N
z
is independent of depth and making the reasonable approximation that r
s
r
0
the reference density, this can be arranged to give:
@
2
u
g
N
z
g
N
z
z
2
¼
z
ð
9
:
4
Þ
@
where w
¼ @
/
@
x and x
¼
(1/r
0
)
@
r/
@
x. After integrating twice with respect to z,
we have:
2N
z
z
2
g
6N
z
z
3
g
u
ð
z
Þ¼
þ
Az
þ
B
ð
9
:
5
Þ
where A and B are constants of integration to be determined from the boundary
conditions: (i) at the surface (z
¼
0) the stress t
s
¼
0 and (ii) at the bottom (z
¼
h)
the velocity u
¼
0. Applying these conditions gives the constants as:
h
2
2N
z
h
3
6N
z
g
g
A
¼
0
;
B
¼
ð
9
:
6
Þ
so that the solution for u is:
2N
z
ð
g
6N
z
ð
g
z
2
h
2
z
3
h
3
u
ð
z
Þ¼
Þ
þ
Þ:
ð
9
:
7
Þ
Note that an alternative bottom boundary condition, setting the bed shear stress to
match the bottom stress given by a quadratic drag law t
b
¼
u
b
, yields similar
results which you can see in a numerical simulation at the website. In a steady state,
In a steady state, which you can see in a numerical simulation at the website, the
gradients w and x are not independent. They are related through the condition that the
net transport must match the freshwater input R
w
per unit width from the river, i.e.
k
b
r
j
u
b
j
Search WWH ::
Custom Search