Changeset 5300

Timestamp:

May 9, 2008, 3:08:19 PM (17 years ago)

Author:

steve

Message:

Working version of wave.py and dam_break.py

Location:

anuga_validation/convergence_study

Files:

• dam_break.py (modified) (4 diffs)
• wave.py (modified) (2 diffs)

anuga_validation/convergence_study/dam_break.py

@@ -41 +41 @@
 # Setup domain
 #------------------------------------------------------------------------------
-dx = 1000.
+dx = 200.
 dy = dx
 L = 100000.
-W = 50*dx
+W = 10*dx
 
 # structured mesh
@@ -64 +64 @@
 domain.use_edge_limiter = True
 domain.tight_slope_limiters = True
+domain.use_centroid_velocities = False
 
 domain.CFL = 1.0
@@ -80 +81 @@
 
 h0 = 10.0
-h1 = 0.1
+h1 = 0.0
 
 def height(x,y):
@@ -110 +111 @@
 vis = RealtimeVisualiser(domain)
 vis.render_quantity_height("stage", zScale =h0*500, dynamic=True)
-vis.colour_height_quantity('stage', (1.0, 0.5, 0.5))
+vis.colour_height_quantity('stage', (0.0, 0.5, 1.0))
 vis.start()
 #===============================================================================

anuga_validation/convergence_study/wave.py

@@ -42 +42 @@
 # Setup domain
 #------------------------------------------------------------------------------
-dx = 100.
+dx = 1000.
 dy = dx
 L = 100000.
@@ -64 +64 @@
 
 domain.use_edge_limiter = True
-domain.tight_slope_limiters = True
+domain.tight_slope_limiters = False
+domain.use_centroid_velocities = False
 
 domain.CFL = 1.0
 
-domain.beta_w      = 0.6
+domain.beta_w      = 1.0
 domain.beta_w_dry  = 0.0
-domain.beta_uh     = 0.6
+domain.beta_uh     = 1.0
 domain.beta_uh_dry = 0.0
-domain.beta_vh     = 0.6
+domain.beta_vh     = 1.0
 domain.beta_vh_dry = 0.0
 domain.beta_h    = 0.0
 
-interactive_visualisation = False
+interactive_visualisation = True