Changeset 5568

Timestamp:

Jul 24, 2008, 4:05:04 PM (15 years ago)

Author:

sexton

Message:

incorporated definitions from web and Synolakis paper

File:

• anuga_work/development/landslide_4/benchmark_problem4.py (modified) (5 diffs)

anuga_work/development/landslide_4/benchmark_problem4.py

@@ -30 +30 @@
 #interior_regions = [[project_slide.poly_syd1, local_res],
 #                    [project_slide.poly_coast, coast_res]]
-length = 10.
-width = 10.
-depth = 10.
-bounding_polygon = [[0,width],[0,0],[length,0],[length,width]]
+length = 104.
+width = 3.7
+
+bounding_polygon = [[0,0],[length,0],[length,width],[0,width]]
 meshname = 'slide.msh'
+remainder_res = 1.
+
 from caching import cache
 _ = cache(create_mesh_from_regions,
           bounding_polygon,
-           {'boundary_tags': {'top': [0], 'left': [1], 
-                              'bottom': [2], 'right': [3]},
+           {'boundary_tags': {'bottom': [0], 'right': [1], 
+                              'top': [2], 'left': [3]},
            'maximum_triangle_area': remainder_res,
            'filename': meshname},
@@ -63 +65 @@
 #-------------------------------------------------------------------------------
 from math import tan
+depth = 4.6
+b = 9.144 # from Synolakis paper
 def topography(x,y):
-    z = -depth + tan(0.5)*x
+    z = []
+    for xi in x:
+        if xi < b: 
+            z.append(-xi/2.)
+        else:
+            z.append(-depth)
+        #print 'hello', xi,b, z
     return z
 domain.set_quantity('stage', 0.0)
@@ -78 +88 @@
 #-------------------------------------------------------------------------------
 from smf import slide_tsunami
+from math import atan, degrees
 
-tsunami_source = slide_tsunami(length=2.,
-                               width=1.,
-                               depth=1.,
-                               slope=45,
-                               thickness=0.1, 
-                               x0=length/2, 
+a = .4572 # from Synolakis paper
+tsunami_source = slide_tsunami(length=0.91,
+                               width=0.61,
+                               depth=.1, # delta from Synolakis paper
+                               slope=degrees(atan(0.5)),
+                               thickness=a, 
+                               x0=0.75, 
                                y0=width/2, 
                                alpha=0.,
-                                         #dx=0.001,
-                               domain=domain)
-                               #verbose=True)
+                               #dx=0.001,
+                               domain=domain,
+                               verbose=True)
 
 #-------------------------------------------------------------------------------                                 
@@ -97 +109 @@
 
 Br = Reflective_boundary(domain)
-Bd = Dirichlet_boundary([tide,0,0])
+Bd = Dirichlet_boundary([0,0,0])
 
-domain.set_boundary( {'e0': Bd,  'e1': Bd, 'e2': Bd, 'e3': Bd} )
+domain.set_boundary( {'top': Bd,  'bottom': Bd, 'right': Bd, 'left': Bd} )
 
 
@@ -111 +123 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 10): 
+for t in domain.evolve(yieldstep = 1, finaltime = 500): 
     domain.write_time()
     stagestep = domain.get_quantity('stage')