Changeset 4362

Timestamp:

Apr 5, 2007, 4:13:27 PM (17 years ago)

Author:

ole

Message:

Work on convergence of solitary wave

Location:

anuga_validation/solitary_waves

Files:

• convergence_results.txt (modified) (1 diff)
• solitary_wave_runup.py (modified) (5 diffs)

anuga_validation/solitary_waves/convergence_results.txt

@@ -1 +1 @@
 CONVERGENCE RESULTS
-        
-        
-        
-max_area = 1.0  
+
+
+
+**********************************************
+Coastline elevation = -0.00 at (x,y)=(300.05, 85.17)
+Max coastline elevation = 4.24 at (88.25, 85.17)
+Run up distance = 211.81
+Max area in runup zone = 2.00
+**********************************************  
+
 **********************************************
 Coastline elevation = -0.00 at (x,y)=(300.05, 99.61)
 Max coastline elevation = 3.91 at (104.56, 99.61)
 Run up distance = 195.49
+Max area in runup zone = 1.00
+**********************************************  
+        
 **********************************************
+Coastline elevation = -0.00 at (x,y)=(300.05, 0.55)
+Max coastline elevation = 3.39 at (130.32, 0.55)
+Run up distance = 169.74
+Max area in runup zone = 0.50
+**********************************************
+
+**********************************************
+Coastline elevation = -0.00 at (x,y)=(300.05, 11.05)
+Max coastline elevation = 2.99 at (150.74, 11.05)
+Run up distance = 149.31
+Max area in runup zone = 0.25
+**********************************************
+
+

anuga_validation/solitary_waves/solitary_wave_runup.py

@@ -3 +3 @@
 Water driven up a linear slope and time varying boundary,
 similar to a beach environment.
+
+For the case (H=0.5m, tan beta=1:50, and still water depth 20m),
+the analytical solution for runup height (R) should be 3.97997m.
 
 """
@@ -27 +30 @@
 
 
+
 #------------------------------------------------------------------------------
 # Model constants
+
+max_area = 0.25 # Maximal triangle area in runup zone (interior_polygon)
+              # This is the variable that we use to measure convergence.
+              # Range in [1.0/8; 8.0] 
+
+
+
 
 slope = -0.02       # 1:50 Slope, reaches h=20m 1000m from western bndry,
@@ -53 +64 @@
 #------------------------------------------------------------------------------
 
-# Structured mesh
-dx = 40.           # Resolution: Length of subdivisions on x axis (length)
-dy = 40.          # Resolution: Length of subdivisions on y axis (width)
-
 length = east-west
 width = north-south
-##points, vertices, boundary = rectangular_cross(length/dx, width/dy,
-##                                               len1=length, len2=width,
-##                                               origin = (west, south)) 
-##
-##domain = Domain(points, vertices, boundary) # Create domain
-
 
 # Unstructured mesh
@@ -70 +71 @@
 interior_polygon = [[400,north],[west+10,north],
                     [west+10,south],[400,south]]
-res = 1.
-simulation_name = 'runup_convergence' + str(res)
+
+simulation_name = 'runup_convergence' + str(max_area)
 meshname = simulation_name + '.msh'
 create_mesh_from_regions(polygon,
                          boundary_tags={'top': [0], 'left': [1],
                                         'bottom': [2], 'right': [3]},
-                         maximum_triangle_area=100.,#dx*dy/4.,
+                         maximum_triangle_area=100.,
                          filename=meshname,
-                         #interior_regions=[[interior_polygon,dx*dy/32.]])
-                         interior_regions=[[interior_polygon,res]])
+                         interior_regions=[[interior_polygon, max_area]])
 
 
 
 domain = Domain(meshname, use_cache=True, verbose = True)
-domain.set_minimum_storable_height(0.01)
-
+#domain.set_minimum_storable_height(0.01)
+#domain.set_minimum_allowed_height(0.01)
+domain.beta_h = 0.0
+domain.limit2007 = 1
 domain.set_name(simulation_name)
 
@@ -175 +177 @@
 print 'Max coastline elevation = %.2f at (%.2f, %.2f)' %(w_max, x_max, y_max)
 print 'Run up distance = %.2f' %sqrt( (x_max-x0)**2 + (y_max-y0)**2 )
+print 'Max area in runup zone = %.2f' %max_area 
 print '**********************************************'