Changeset 4773


Ignore:
Timestamp:
Oct 31, 2007, 4:56:02 PM (16 years ago)
Author:
duncan
Message:

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

Files:
3 added
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/uq_flume_2007/project.py

    r4750 r4773  
    66velocity_x_filename =  "velocity_x.csv"
    77velocity_y_filename =  "velocity_y.csv"
     8experiment_vel_filename =  "velocity_x_experiment.csv"
     9experiment_depth_filename =  "depth_experiment.csv"
    810
    911
  • anuga_validation/uq_flume_2007/run_dam.py

    r4750 r4773  
    11"""Script for running a dam break simulation of UQ's dam break tank.
     2
     3The simulation is based on the simulation Matt has presented and will
     4form part of the ANUGA validation paper.
     5
     6i want to compare the stage at .4m and the velocity at .45m.
    27
    38
     
    3338def main():
    3439     
    35     slopes = [0]
    36     frictions = [0.01] #can enter multiple frictions eg [0.00, 0.01]
    37     inital_depths = [0.2] #can enter multiple initial depths eg [0.1, 0.2]
    38     gate_positions = [0.75]
     40    slope= 0
     41    friction = 0.01
     42    inital_depth = 0.2
     43    gate_position = 0.75
    3944   
    40     scenario(0, 0.01, 0.02, 0.75)
     45    return scenario(slope, friction, inital_depth, gate_position)
    4146           
    4247
     
    4954    create_mesh.generate(project.mesh_filename,
    5055                         gate_position,
    51                          is_course=True) # this creates the mesh
    52                          #is_course=False) # this creates the mesh
     56                         #is_course=True) # this creates the mesh
     57                         is_course=False) # this creates the mesh
    5358
    5459    head,tail = path.split(project.mesh_filename)
     
    101106        print 'finished'
    102107       
    103     points = [[0.4,0.2]]
     108    points = [[0.4,0.2], [0.45,0.2]]
    104109    #-------------------------------------------------------------------------
    105110    # Calculate gauge info
     
    108113    interpolate_sww2csv( project.basename +".sww",
    109114                        points,
    110                         project.depth_filename + '.csv',
    111                         project.velocity_x_filename + '.csv',
    112                         project.velocity_y_filename + '.csv')
     115                        project.depth_filename ,
     116                        project.velocity_x_filename,
     117                        project.velocity_y_filename )
    113118
    114119
  • anuga_work/publications/anuga_2007/anuga_validation.tex

    r4768 r4773  
    483483This successful replication of the tsunami wave tank simulation on a
    484484complex 3D beach is a positive first step in validating the \ANUGA{}
    485 modelling capability. Subsequent validation will be conducted as
    486 additional datasets become available.
    487 
     485modelling capability.
     486
     487
     488\begin{figure}[htbp]
     489\centerline{\includegraphics[width=4in]{uq-flume-depth.eps}}
     490\caption{Comparison of wave tank and \ANUGA{} water height at .4 m
     491  from the gate}\label{fig:uq-flume-depth}
     492\end{figure}
     493
     494\begin{figure}[htbp]
     495\centerline{\includegraphics[width=4in]{uq-flume-velocity.eps}}
     496\caption{Comparison of wave tank and \ANUGA{} water velocity at .45 m
     497  from the gate}\label{fig:uq-flume-depth}
     498\end{figure}
     499
     500Flume data from the University of Queensland has also been used for
     501validating \ANUGA{}.  The experiments were carried out at the Gordon
     502McKay Hydraulics Laboratory at St Lucia.  The Flume was set up for
     503Dam-break experiments, having a water reservior at one end.  The flume
     504was glass-sided, 3m long, 0.4m in wide, and 0.4m deep (check
     505depth), with a PVC bottom. The reservoir in the flume was 0.75m long.
     506For this experiment the reservoir water was 0.2m deep. At time zero
     507the reservoir gate is opened and the water flows into the other side
     508of the flume.
    488509
    489510\section{Conclusions}
Note: See TracChangeset for help on using the changeset viewer.