28.01.2013 Views

Zienkiewicz O.C., Taylor R.L. Vol. 3. The finite - tiera.ru

Zienkiewicz O.C., Taylor R.L. Vol. 3. The finite - tiera.ru

Zienkiewicz O.C., Taylor R.L. Vol. 3. The finite - tiera.ru

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

8.3 Dif®culties in modelling surface waves<br />

<strong>The</strong> main defects of the simple surface-wave model described above are the following:<br />

1. inaccuracy when the wave height becomes large. <strong>The</strong> equations are no longer valid<br />

when becomes large, and for very large , the waves will break, which introduces<br />

energy loss.<br />

2. lack of modelling of bed friction. This will be discussed below.<br />

<strong>3.</strong> lack of modelling of separation at re-entrant corners. At re-entrant corners there is<br />

a singularity in the velocity of the form 1= r<br />

p . <strong>The</strong> velocities become large, and<br />

physically the viscous e€ects, neglected above, become important. <strong>The</strong>y cause<br />

retardation, ¯ow separation and eddies. This e€ect can only be modelled in an<br />

approximate way.<br />

Now the response can be determined for a given excitation frequency, as discussed in<br />

Chapter 17 of <strong>Vol</strong>ume 1.<br />

8.4 Bed friction and other effects<br />

<strong>The</strong> Che zy bed friction term is non-linear and if it is included in its original form it<br />

makes the equations very di cult to solve. <strong>The</strong> usual procedure is to assume that<br />

its main e€ect is to damp the system, by absorbing energy, and to introduce a<br />

linear term, which in one period absorbs the same amount of energy as the Che zy<br />

term. <strong>The</strong> linearized bed friction version of Eq. (8.2) is<br />

r T …Hr †‡ !2<br />

g<br />

ÿ i!M ˆ 0 or<br />

@<br />

@x i<br />

H @<br />

@x i<br />

‡ !2<br />

g<br />

ÿ i!M ˆ 0 …8:11†<br />

where M is a linearized bed friction coe cient, which can be written as<br />

M ˆ 8u max=3 C 2 H, C is the Che zy constant and u max is the maximum velocity at<br />

the bed at that point. In general the results for will now be complex, and iteration<br />

has to be used, since M depends upon the unknown u max. From the ®nite element<br />

point of view, there is no longer any need to separate the `sti€ness' and `mass'<br />

matrices. Instead, Eq. (8.11) is weighted using the element shape function and the<br />

entire complex element matrix is formed. <strong>The</strong> matrix right-hand side arises from<br />

whatever exciting forces are present. <strong>The</strong> re-entrant corner e€ect and wave-absorbing<br />

walls and permeable breakwaters can also be modelled in a similar way, as both of<br />

these introduce a damping e€ect, due to viscous dissipation. <strong>The</strong> method is explained<br />

in reference 15, where an example showing ¯ow through a perforated wall in an<br />

o€shore st<strong>ru</strong>cture is solved.<br />

8.5 <strong>The</strong> short-wave problem<br />

<strong>The</strong> short-wave problem 245<br />

Short-wave di€raction problems are those in which the wavelength is much smaller<br />

than any of the dimensions of the problem. Such problems arise in surface waves<br />

on water, acoustics and pressure waves, electromagnetic waves and elastic waves.<br />

<strong>The</strong> methods described in this chapter will solve the problems, but the requirement

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!