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