29.07.2013 Views

Numerical modeling of waves for a tsunami early warning system

Numerical modeling of waves for a tsunami early warning system

Numerical modeling of waves for a tsunami early warning system

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

and<br />

<strong>Numerical</strong> <strong>modeling</strong> <strong>of</strong> <strong>waves</strong> <strong>for</strong> a <strong>tsunami</strong> <strong>early</strong> <strong>warning</strong> <strong>system</strong><br />

(ccg) x<br />

ax = <br />

<br />

(ccg) u<br />

ux +(ccg) v<br />

vx [auux + avvx] =<br />

= u2 x (ccga) u +(ccg) u<br />

uxvxav +(ccg) v<br />

uxvxau + v2 <br />

x (ccga) v =<br />

= au u2 <br />

x (ccg) u + uxvx (ccg) v + av v2 <br />

x (ccg) v + uxvx (ccg) u<br />

(ccg) y ay = <br />

<br />

[auuy + avvy] =<br />

(ccg) u<br />

uy +(ccg) v<br />

vy<br />

= u2 y (ccga) u +(ccg) u<br />

uyvyav +(ccg) v<br />

uyvyau + v2 <br />

y (ccga) v =<br />

= au u2 <br />

y (ccg) u + uyvy (ccg) v + av v2 <br />

y (ccg) v + uyvy (ccg) u<br />

Now we are ready to write equation (B.6)<br />

ccg [auuu2 x + avvv2 x +2auvuxvx + au (ux) x + av <br />

(vx) x ]+<br />

+ccg auuu2 y + avvv2 y +2auvuyvy + au (uy) y + av (vy) y + k2a <br />

+<br />

<br />

+au u2 <br />

x (ccg) u + uxvx (ccg) v + av v2 <br />

x (ccg) v + uxvx (ccg) u +<br />

<br />

+au u2 <br />

y (ccg) u + uyvy (ccg) v + av v2 <br />

y (ccg) v + uyvy (ccg) u =0<br />

which can be rewritten as follows:<br />

<br />

ccgauu u2 x + u2 <br />

y + ccgavv v2 x + v2 <br />

y<br />

+ccg (au∇2u + av∇2v)+ccgk2 <br />

a+<br />

+au (ccg) u u2 x + u2 <br />

y<br />

+av (ccg) v<br />

<br />

v 2 x + v 2 y<br />

+ au (ccg) v (uxvx + uyvy)+<br />

+2ccgauv (uxvx + uyvy)+<br />

+ av (ccg) u (uxvx + uyvy) =0<br />

(B.10)<br />

(B.11)<br />

(B.12)<br />

(B.13)<br />

the first and the sixth terms <strong>of</strong> equation (B.13) can be assembled as the<br />

second and the seventh, yielding to<br />

<br />

u2 x + u2 <br />

y (ccgau) u + <br />

v2 x + v2 <br />

y<br />

(uxvx + uyvy) <br />

<br />

au (ccg) v + av (ccg) u +2ccgauv<br />

(ccgav) v + ccgk 2 a+<br />

+ ccg (au∇ 2 u + av∇ 2 v)=0<br />

(B.14)<br />

while the derivatives <strong>of</strong> a above are taken with respect to the mapped<br />

coordinates (u, v), the coefficients <strong>of</strong> the derivatives still involve the<br />

Università degli Studi di Roma Tre - DSIC 115

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

Saved successfully!

Ooh no, something went wrong!