Geoscience Reference
In-Depth Information
The above is not a satisfactory equation for finding c because we have an
unknown constant p t . To eliminate p t , we must find an independent equation
in which p t appears. We consider the horizontal equation of motion again, but
this time we consider it everywhere—not just at z ¼ 0. The inviscid, steady-state,
two-dimensional equation of motion for this case is
p 0
u
@
u
=@
x þ w
@
u
=@
z þ 0 @
=@
x ¼ 0
ð 3
:
21 Þ
which is similar to (3.18), except that the vertical advection term is retained.
Combining (3.21) with the Boussinesq two-dimensional equation of continuity
@
u
=@
x þ@
w
=@
z ¼ 0
ð 3
:
22 Þ
we find that
u 2
p 0
@
=@
x þ@=@
z ð uw Þþ 0 @
=@
x ¼ 0
ð 3
:
23 Þ
We now integrate (3.23) over the x-z-plane from x ¼1 to x ¼1 and from
z ¼ 0toz ¼ H (i.e., over the entire domain), making use of (3.16) and (3.17), and
eliminate u using the mass conservation relation (3.13), to find that
c 2 Hh
2 g h 2
1
H h Þ¼ 0 p t H
ð 1 0 Þ= 0
ð 3
:
24 Þ
Eliminating p t from (3.24) and (3.20), we find that
c 2
¼ g h ½ð 1 0 Þ= 0 f 2 ½ð H h Þ=ð H þ h Þ½ð H
1
2 h Þ=
H g
ð 3
:
25 Þ
Alas, it is not obvious whether or not (3.25) will work for all values of h
H,sowe
must find yet another constraint to nail down c. To do so, we now make use of
the Boussinesq vorticity (the y-component) equation (2.51) in flux form, subject to
the conditions of being steady state and two dimensional (
=
@=@
y ¼ 0), and to the
continuity equation (3.22):
u
Þ=@
x þ@ð w
Þ=@
z ¼@
B
=@
x
ð 3
:
26 Þ
where
x, the y-component of the vorticity vector. Integrating
(3.26) over the entire domain and using the upper and lower kinematic boundary
conditions w ¼ 0atz ¼ 0 and z ¼ H, and assuming that w ¼ 0atx ¼1
(altogether assuming that w ¼ 0 along all the edges of the domain), it is found
that
¼ @
u
=@
z @
w
=@
2 u 2 H ;1 ¼ g h ½ð 1 0 Þ= 0
1
ð 3
:
27 Þ
We need to find an expression now for c in terms of u H ;1 and eliminate the
latter. To do so, we use of the conservation of mass relation (3.13), noting that
u u H ;1 and find that
c 2
2
28 Þ
We now have two independent equations for c 2 . Setting (3.28) and (3.25) equal to
each other, we obtain the following relationship between H and h:
¼ g h ½ð 1 0 Þ= 0 2 ½ð H h Þ=
H
ð 3
:
2
H 2
1
ð H h Þ
=
¼½ð H h Þ=ð H þ h Þ½ð H
2 h Þ=
H
ð 3
:
29 Þ
This equation has solutions h ¼ H
=
2, h ¼ H,andH ¼1 . The solution h ¼ H is
Search WWH ::




Custom Search