18
Land-Ocean Systems in the Siberian Arctic: Dynamics and History
't~ , 't; - components of the wind stress; 't~ ,'t~ - components of the bottom stress; r - the
density of sea water.
Boundary conditions on the solid boundary are:
U = V =0
(2)
And on the open sea boundary:
~=O
(3)
The bottom stress is parametrized in terms of a quadratic low:
(4)
K
where K =
w 2 ; Kw - the coefficient of bottom stress, takes the value 0.0026.
(H+~)
Wind stress on the sea surface was computed from:
(5)
where ra - density of air; W - the wind velocity; CS = (1.76+0.81W)1O- 3 .
Solution
Using algorithm of finite element's, we divide the horizontal space domain W into N arbitrary
triangle elements. We can write for each element:
6
6
A
6
U = L Uj'Pj; V = L Vj'Pj; ~ = L~j'Pj
(6)
i=l
i=l
i=l
where Ui, Vi and Xi - values of U, V and X respectively in i-points; ji - quadratic basis function
(Connor and Brebbia, 1979)
For the whole domain:
M
U= :L,.UjFj
(7)
j;l
M
V= :L,.VjFj
(8)
i=l
(9)
where M - the number of the points all over the domain; Fi - functions, defined as: