source: anuga_work/development/convergence_okushiri_2008/README.txt @ 5604

Last change on this file since 5604 was 5366, checked in by Leharne, 15 years ago
File size: 2.3 KB
Line 
1ANUGA CONVERGENCE STUDY USING TRUE-SCALE VERSION OF THE
2OKUSHIRI ISLAND TSUNAMI WAVETANK EXPERIMENT
3
4This directory currently contains code to scale-up to true scale, the 1:400
5wave tank simulation of the 1993 Okushiri island tsunami as described at
6the Third International Conference on Long Wave Runup:
7http://www.cee.cornell.edu/longwave/index.cfm?page=benchmark&problem=2.
8Once this up-scaling has been completed and verified, then the files will be used for conducting
9a convergence study in ANUGA. This component is still in development.
10
11Data files available in this directory are
12
13okushiri_truescale_bathymetry.txt:
14  The true-scale digital elevation model
15
16okushiri_truescale_input.txt:
17  The true-scale timeseries applied at the western boundary   
18
19okushiri_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
24The ANUGA scripts to run are
25
26project_truescale.py:
27  This script contains project filenames and is called in the create and run scripts below.
28
29create_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 
33run_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     
38compare_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
44Methodology for true-scale transformation:
45
46Positions (ie bathymetry, polygon definitions, gauge locations) were derived by
47carrying out a scalar multiplication of all x, y, z values in the original files by 400.
48
49Input waveform: Truescale input wave with period T' is assumed to have a wavelength
50and amplitude 400 times the original 1:400 waveform. Given the relationship:
51
52T = wavelength / sqrt (g * h)
53
54Then T' = 400 * wavelength / sqrt (g * 400 * h)
55           = 20 * wavelength / sqrt (g * h)
56           = 20 * T
57
58Therefore, time (s) is multiplied by 20 and water surface (m) (ie amplitude) is multiplied by 400.
59
Note: See TracBrowser for help on using the repository browser.