Changeset 8159

Timestamp:

Mar 20, 2011, 10:56:33 PM (14 years ago)

Author:

steve

Message:

 

File:

• trunk/anuga_validation/validation_tests/dam_break_analytical.py (modified) (1 diff)

trunk/anuga_validation/validation_tests/dam_break_analytical.py

@@ -52 +52 @@
 
 
-
-
-
+def dam_break(x,t,h0=1.0,h1 = 10.0):
+    h2 = h2(h0,h1)
+    u2 = 2.0*(math.sqrt(g*h1) - math.sqrt(g*h2))
+    s = u2*h2/(h2 - h0)
+    c1 = math.sqrt(g*h1)
+    c0 = math.sqrt(g*h0)
+    c2 = math.sqrt(g*h2)
+    
+    if x < -t*c1:
+        h = h1
+        u = 0.0
+    elif x < t*(u2 - c2):
+        h = 1.0/g*(2.0/3.0 *c1 - 1.0/3.0*x/t)**2
+        u = 2.0/3.0*(c1 + x/t)
+    elif x < s*t:
+        h = h2
+        u = u2
+    else:
+        h = h0
+        u = u0
+        
+print dam_break(1.0,1.0)