1 | ANUGA CONVERGENCE STUDY USING TRUE-SCALE VERSION OF THE |
---|
2 | OKUSHIRI ISLAND TSUNAMI WAVETANK EXPERIMENT |
---|
3 | |
---|
4 | This directory currently contains code to scale-up to true scale, the 1:400 |
---|
5 | wave tank simulation of the 1993 Okushiri island tsunami as described at |
---|
6 | the Third International Conference on Long Wave Runup: |
---|
7 | http://www.cee.cornell.edu/longwave/index.cfm?page=benchmark&problem=2. |
---|
8 | Once this up-scaling has been completed and verified, then the files will be used for conducting |
---|
9 | a convergence study in ANUGA. This component is still in development. |
---|
10 | |
---|
11 | Data files available in this directory are |
---|
12 | |
---|
13 | okushiri_truescale_bathymetry.txt: |
---|
14 | The true-scale digital elevation model |
---|
15 | |
---|
16 | okushiri_truescale_input.txt: |
---|
17 | The true-scale timeseries applied at the western boundary |
---|
18 | |
---|
19 | okushiri_output_truescale_ch5-7-9.txt: |
---|
20 | Experimental data measured at three gauge locations in the original wavetank experiment |
---|
21 | which has been up-scaled to true-scale |
---|
22 | |
---|
23 | |
---|
24 | The ANUGA scripts to run are |
---|
25 | |
---|
26 | project_truescale.py: |
---|
27 | This script contains project filenames and is called in the create and run scripts below. |
---|
28 | |
---|
29 | create_okushiri_truescale.py: |
---|
30 | This script will convert the text files to native |
---|
31 | ANUGA netcdf formats and also create a suitable triangular mesh. |
---|
32 | |
---|
33 | run_okushiri_truescale.py: |
---|
34 | This script will run a numerical simulation based on the the bathymetry |
---|
35 | and the given boundary condition and store the model output in an ANUGA |
---|
36 | sww file which can be viewed using animate, or further interrogated by ANUGA. |
---|
37 | |
---|
38 | compare_timeseries.py |
---|
39 | This script will extract timeseries from the sww file and plot them |
---|
40 | together with the experimental data provided. Numerical similarity |
---|
41 | measures will also be computed. |
---|
42 | |
---|
43 | |
---|
44 | Methodology for true-scale transformation: |
---|
45 | |
---|
46 | Positions (ie bathymetry, polygon definitions, gauge locations) were derived by |
---|
47 | carrying out a scalar multiplication of all x, y, z values in the original files by 400. |
---|
48 | |
---|
49 | Input waveform: Truescale input wave with period T' is assumed to have a wavelength |
---|
50 | and amplitude 400 times the original 1:400 waveform. Given the relationship: |
---|
51 | |
---|
52 | T = wavelength / sqrt (g * h) |
---|
53 | |
---|
54 | Then T' = 400 * wavelength / sqrt (g * 400 * h) |
---|
55 | = 20 * wavelength / sqrt (g * h) |
---|
56 | = 20 * T |
---|
57 | |
---|
58 | Therefore, time (s) is multiplied by 20 and water surface (m) (ie amplitude) is multiplied by 400. |
---|
59 | |
---|