A non-hydrostatic two-layer staggered scheme for transient waves due to anti-symmetric seabed thrust
The development of transient waves generated by bottom motion is studied numericallyin this work. A non-hydrostatic numerical scheme, based on solving the two dimensionalEuler equations using two-layer approximation for the vertical direction, is implemented.The dispersion relation of this scheme is shown to agree with the analytical dispersionrelation over a wide range of kd0, where k denotes the wave number and d0 the characteristic water depth. To ensure that a good balance between non-linearity and dispersion is accommodated by the scheme, the propagation of a solitary wave (undisturbed in shape) was simulated. Our next focus was on the simulation of transient waves generatedby bottom motion. After conducting a benchmark test against Hammack's experimentalresults for downward bottom motion, an anti-symmetric bottom thrust was considered.The resulting transient waves developed different behavior depending on the water depth.Finally, to mimic the December 2004 tsunami, a seabed motion was generated over Acehbathymetry. This simulation showed that a package of wave trains developed and propa-gated towards the Aceh coast, and exhibited inter alia the feature of shoreline withdrawaloften observed.