Changeset 4720 for anuga_work/development/demos/island.py

Timestamp:

Sep 10, 2007, 6:32:46 PM (17 years ago)

Author:

ole

Message:

Played with the 'creep' problem. I am still not sure what is going on.

File:

• anuga_work/development/demos/island.py (modified) (5 diffs)

anuga_work/development/demos/island.py

@@ -3 +3 @@
 Island surrounded by water.
 This example investigates onshore 'creep'
+
+Creep is defined as water moving uphill over many timesteps.
+The initial adjustment at time = 0 is not 'creep'
 
 """
@@ -33 +36 @@
                                            'top': [2],
                                            'left': [3]},
-                          maximum_triangle_area = 5,
+                          maximum_triangle_area = 2,
                           filename = 'island.msh' ,
                           interior_regions=[ ([[50,25], [70,25], [70,75], [50,75]], 10.0)]
@@ -60 +63 @@
 #domain.alpha_balance = 10.0
 #domain.minimum_allowed_height = 1.0e-4 
-#domain.maximum_allowed_speed = 100.0 
+#domain.maximum_allowed_speed = 0 
 #domain.minimum_storable_height = 1.0e-4 
 
-domain.beta_h     = 0.0
-domain.tight_slope_limiters = 1
+#------------------------
+# Combinations with creep
+#------------------------
+
+#domain.maximum_allowed_speed = 0
+#domain.beta_h = 0.0
+#domain.tight_slope_limiters = 0
+
+#domain.maximum_allowed_speed = 0 
+#domain.beta_h = 0.0
+#domain.H0 = 0.01
+#domain.tight_slope_limiters = 1
+
+#---------------------------
+# Combinations with less creep ??
+#---------------------------
+
+domain.maximum_allowed_speed = 1
+domain.beta_h = 0.0
+domain.tight_slope_limiters = 0
+
+#domain.maximum_allowed_speed = 10 
+#domain.beta_h = 0.0
+#domain.H0 = 0.01
+#domain.tight_slope_limiters = 1
+
+
+
+
 
 #------------------------------------------------------------------------------
@@ -86 +116 @@
     return z
 
-stage_value = 15.0
+stage_value = 10.0
 #domain.set_quantity('friction', 0.1)  #Honky dory
 domain.set_quantity('friction', 0.01)     #Creep
@@ -110 +140 @@
 
 import time
-for t in domain.evolve(yieldstep = 1, finaltime = 5):
+for t in domain.evolve(yieldstep = 1, finaltime = 25):
     domain.write_time()
     #if allclose(t, 100):