Changeset 4773

Timestamp:

Oct 31, 2007, 4:56:02 PM (17 years ago)

Author:

duncan

Message:

Working on the anuga-2007 paper. Adding results from the UQ flume.

Files:

• anuga_validation/uq_flume_2007/plot.py (added)
• anuga_validation/uq_flume_2007/project.py (modified) (1 diff)
• anuga_validation/uq_flume_2007/run_dam.py (modified) (5 diffs)
• anuga_work/publications/anuga_2007/anuga_validation.tex (modified) (1 diff)
• anuga_work/publications/anuga_2007/uq_flume_depth.eps (added)
• anuga_work/publications/anuga_2007/uq_flume_velocity.eps (added)

anuga_validation/uq_flume_2007/project.py

@@ -6 +6 @@
 velocity_x_filename =  "velocity_x.csv"
 velocity_y_filename =  "velocity_y.csv"
+experiment_vel_filename =  "velocity_x_experiment.csv"
+experiment_depth_filename =  "depth_experiment.csv"
 
 

anuga_validation/uq_flume_2007/run_dam.py

@@ -1 +1 @@
 """Script for running a dam break simulation of UQ's dam break tank.
+
+The simulation is based on the simulation Matt has presented and will
+form part of the ANUGA validation paper.
+
+i want to compare the stage at .4m and the velocity at .45m.
 
 
@@ -33 +38 @@
 def main():
      
-    slopes = [0]
-    frictions = [0.01] #can enter multiple frictions eg [0.00, 0.01]
-    inital_depths = [0.2] #can enter multiple initial depths eg [0.1, 0.2]
-    gate_positions = [0.75]
+    slope= 0
+    friction = 0.01 
+    inital_depth = 0.2
+    gate_position = 0.75
    
-    scenario(0, 0.01, 0.02, 0.75)
+    return scenario(slope, friction, inital_depth, gate_position)
             
 
@@ -49 +54 @@
     create_mesh.generate(project.mesh_filename,
                          gate_position,
-                         is_course=True) # this creates the mesh
-                         #is_course=False) # this creates the mesh
+                         #is_course=True) # this creates the mesh
+                         is_course=False) # this creates the mesh
 
     head,tail = path.split(project.mesh_filename)
@@ -101 +106 @@
         print 'finished'
         
-    points = [[0.4,0.2]]
+    points = [[0.4,0.2], [0.45,0.2]]
     #-------------------------------------------------------------------------
     # Calculate gauge info
@@ -108 +113 @@
     interpolate_sww2csv( project.basename +".sww",
                         points,
-                        project.depth_filename + '.csv',
-                        project.velocity_x_filename + '.csv',
-                        project.velocity_y_filename + '.csv')
+                        project.depth_filename ,
+                        project.velocity_x_filename,
+                        project.velocity_y_filename )
 
 

anuga_work/publications/anuga_2007/anuga_validation.tex

@@ -483 +483 @@
 This successful replication of the tsunami wave tank simulation on a
 complex 3D beach is a positive first step in validating the \ANUGA{}
-modelling capability. Subsequent validation will be conducted as
-additional datasets become available.
-
+modelling capability. 
+
+
+\begin{figure}[htbp]
+\centerline{\includegraphics[width=4in]{uq-flume-depth.eps}}
+\caption{Comparison of wave tank and \ANUGA{} water height at .4 m
+  from the gate}\label{fig:uq-flume-depth}
+\end{figure}
+
+\begin{figure}[htbp]
+\centerline{\includegraphics[width=4in]{uq-flume-velocity.eps}}
+\caption{Comparison of wave tank and \ANUGA{} water velocity at .45 m
+  from the gate}\label{fig:uq-flume-depth}
+\end{figure}
+
+Flume data from the University of Queensland has also been used for
+validating \ANUGA{}.  The experiments were carried out at the Gordon
+McKay Hydraulics Laboratory at St Lucia.  The Flume was set up for
+Dam-break experiments, having a water reservior at one end.  The flume
+was glass-sided, 3m long, 0.4m in wide, and 0.4m deep (check
+depth), with a PVC bottom. The reservoir in the flume was 0.75m long.
+For this experiment the reservoir water was 0.2m deep. At time zero
+the reservoir gate is opened and the water flows into the other side
+of the flume.
 
 \section{Conclusions}