Civil Engineering Reference
In-Depth Information
Comparison with Program 8.2 will show essentially the same array declarations and
input parameters. Extra variables are the velocities
u
and
v
in the
x
and
y
directions,
ux
and
uy
respectively.
The problem chosen is the one-dimensional example shown in Figure 8.22 consisting of
a 56-m deep bed of fluid discretised by 40, 4-noded elements. There is a steady velocity in
∂f
v
c
y
f
=
∂
y
Surface
1
2
3
4
5
6
40 elements @ 1.4
= 56 m
77
78
80
79
y
82
81
bed
x
∂f
v
c
y
1.4 m
=
∂
y
nxe nye np_types
1 40 1
prop(cx,cy)
1.0e-6 0.49
etype(not needed)
x_coords y_coords
0.0 1.4
56.0 54.6 53.2 51.8 50.4 49.0 47.6 46.2 44.8 43.4
42.0 40.6 39.2 37.8 36.4 35.0 33.6 32.2 30.8 29.4
28.0 26.6 25.2 23.8 22.4 21.0 19.6 18.2 16.8 15.4
14.0 12.6 11.2 9.8 8.4 7.0 5.6 4.2 2.8 1.4
0.0
dtim nstep theta npri nres ntime ux uy
300.0 20 1.0 1 82 5 0.0 0.0135
nci
10
Figure 8.22
Mesh and data for Program 8.7 example
