1 | This directory contains files for the prototype study for stochasticity conducted by Suresh Kumar of ACFR and Ole Nielsen of GA. |
---|
2 | |
---|
3 | |
---|
4 | |
---|
5 | Project outline: |
---|
6 | |
---|
7 | The case study is based on a validation conducted against a wavetank experiment emulating the 1993 Okushiri Island tsunami. |
---|
8 | |
---|
9 | |
---|
10 | 1)Investigate independent, identically distributed errors in bathymetry. Decide a standard deviation of the bathymetry errors. |
---|
11 | |
---|
12 | |
---|
13 | The important files are: |
---|
14 | project.py: Controls number of samples etc |
---|
15 | create_realisatons.py: Makes a set of perturbed bathymetries |
---|
16 | run_model.py: Run the simulation for each realisation and output three timeseries for each |
---|
17 | |
---|
18 | |
---|
19 | |
---|
20 | |
---|
21 | (2)Establish that the Monte-Carlo estimates have converged |
---|
22 | |
---|
23 | Brute force validation where we run say 1000 simulations and then |
---|
24 | 2000, and evaluate whether the mean solution at the 3 monitoring |
---|
25 | stations has significantly changed. |
---|
26 | |
---|
27 | |
---|
28 | |
---|
29 | (3)Measures of performance at the three channels where measurements |
---|
30 | are available |
---|
31 | |
---|
32 | (a) 1 standard deviation from the mean (time averaged) |
---|
33 | |
---|
34 | (b) Windowed RMS error measure |
---|
35 | |
---|
36 | (c) Raw RMS error |
---|
37 | |
---|
38 | (d)Visual inspection |
---|
39 | |
---|
40 | |
---|
41 | Milestone: Evaluate the effects of the stochasticity of bathymetry on |
---|
42 | the computed solution using the measures of performance |
---|
43 | listed here to decide on continuance of this line of work. |
---|
44 | |
---|
45 | |
---|
46 | |
---|
47 | Second mini-project (Pending success of earlier project): |
---|
48 | |
---|
49 | Introduce spatially correlated errors in bathymetry and perform |
---|
50 | Monte-Carlo simulations with sparse, locally correlated deviates. |
---|
51 | |
---|
52 | |
---|
53 | Research Questions: |
---|
54 | |
---|
55 | * Which resolution is necessary? |
---|
56 | * What is the effect of noise in the friction terms? |
---|
57 | * Effect of Boundary conditions, source terms and initial condition. |
---|
58 | * Can the whole process be done stochastically - |
---|
59 | perhaps a hybrid where CFD kick starts a stochastic |
---|
60 | predictive process (like ITO) |
---|
61 | |
---|
62 | |
---|
63 | |
---|
64 | |
---|
65 | |
---|