Changeset 8523

Timestamp:

Aug 17, 2012, 1:19:37 AM (12 years ago)

Author:

davies

Message:

Adding 'anuga_tsunami' directory for version suited to tsunami

Location:

trunk

Files:

• anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py (modified) (1 diff)
• anuga_core/validation_tests/Tests/Experimental/Isolated_Building/results.tex (modified) (1 diff)
• anuga_work/development/gareth/experimental/anuga_tsunami (added)
• anuga_work/development/gareth/experimental/anuga_tsunami/__init__.py (added)
• anuga_work/development/gareth/experimental/anuga_tsunami/compile_all.py (added)
• anuga_work/development/gareth/experimental/anuga_tsunami/runup.py (added)
• anuga_work/development/gareth/experimental/anuga_tsunami/swb2_boundary_conditions.py (added)
• anuga_work/development/gareth/experimental/anuga_tsunami/swb2_domain.py (added)
• anuga_work/development/gareth/experimental/anuga_tsunami/swb2_domain_ext.c (added)
• anuga_work/development/gareth/tests/urban_flow/plotme.py (modified) (1 diff)

trunk/anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py

@@ -149 +149 @@
 barrier()
 t0 = time.time()
-
 for t in domain.evolve(yieldstep = yieldstep, finaltime = finaltime):
     if myid == 0:

trunk/anuga_core/validation_tests/Tests/Experimental/Isolated_Building/results.tex

@@ -1 +1 @@
 
-\section{Isolated Building}
+\section{Isolated Obstacle}
 
-Will show a comparison of \anuga{} to the isolated building test.
+Herein \anuga{} is tested against the isolated obstacle test, which involves a dam-break flow past an isolated building. This benchmark experiment is reported in:
+
+Soares-Frazao, S. and Zech, Y. (2007) Experimental study of dam-break flow against an isolated obstacle. Journal of Hydraulic Resaerch 45(Extra Issue):27-36
+
+The data is provided in a DVD-ROM associated with the above JHR issue.
+
+Note that many models have been applied to this problem, and their performance is reviewed in Soares Frazao, S., Noel, B., Spinewine, B., and Zech, Y. (YEAR) The isolated building test case: results from the IMPACT benchmark. At the time of writing, the performance of \anuga{} on this problem is qualitatively similar to the performance of most other models -- i.e. it simulates the main aspects of the flow, but often not the details. It seems likely that these experimental flows may well partly violate the shallow water approximations, and to the extent that this is true, we would not expect perfect agreement.
 
 \subsection{Results}
+For details of the experiment which are required to run the model, consult Soares-Frazao and Zech (2007). A zoom of the model geometry around the building and reservoir mouth is shown in Figure~\ref{urbanflow_geom}, along with measurement stations. 
+\begin{figure}
+\includegraphics{XXX}
+\caption{Zoom of the flume planform geometry near the isolated building. Red areas are obstacles which are never submerged by the flow. See Soares-Frazo and Zech (2007) for further details on the flume geometry.}
+\label{urbanflow_geom}
+\end{figure}
 
-
-Will have a bunch of figures.
-
+Water was released from the reservoir, and allowed to flow past the building. Figures XX - XX2 compare the stage and velocity at each of these measurement points with the modelled result. 
+\begin{figure}
+\includegraphics{XXX}
+\caption{Comparison betweeen modelled and observed stage, x-velocity and y-velocity at sites G1-G6.}
+\label{urbanflow_geom}
+\end{figure}
 
 

trunk/anuga_work/development/gareth/tests/urban_flow/plotme.py

@@ -54 +54 @@
         pyplot.savefig(fig_filename)
         pyplot.clf()
+
+
+### Plot of experimental set-up
+pyplot.figure(figsize=(8,3))
+pyplot.scatter(p.x, p.y, c=p.elev>0.2,edgecolors='none', s=5.0, marker='8')
+pyplot.gca().set_aspect('equal')
+pyplot.xlim((4,15))
+pyplot.title('Zoom of the flume geometry')
+#pyplot.colorbar()
+# Add gauge locations
+for i in range(Gauge_xy.shape[0]):
+    gauge_name='G' + str(i+1)
+    pyplot.scatter(Gauge_xy[i,0], Gauge_xy[i,1],marker='x',color='white')
+    pyplot.text(Gauge_xy[i,0], Gauge_xy[i,1], s=gauge_name, fontweight='black', color='white', ha='left')
+
+pyplot.text(7.0, 1.7, 'Reservoir \n Outflow',color='white', size=8)
+pyplot.quiver(7,1.5, 1.0, 0.0,color='white')
+pyplot.text(4.5, 1.0, 'Reservoir',color='white')
+
+pyplot.text(11.7, 1.5, 'Isolated \n Building',color='white',size=8)
+#pyplot.quiver(12.5,1.5, -1.0, 1.0, color='white')
+
+fig_filename=fig_dir+'Model_domain.png'
+pyplot.savefig(fig_filename, orientation='landscape')
+pyplot.clf()