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
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