Changeset 5085 for anuga_work/debug

Timestamp:

Feb 27, 2008, 10:12:20 AM (17 years ago)

Author:

ole

Message:

Further simplification of run_dam debugging.

File:

• anuga_work/debug/bad_time_step/run_dam.py (modified) (7 diffs)

anuga_work/debug/bad_time_step/run_dam.py

@@ -1 +1 @@
 """
+Script for debugging potential problem with the boundary
+condition: transmissive_momentum_set_stage
 
-Script for running a breaking wave simulation of Jon Hinwoods wave tank.
-Note: this is based on the frinction_ua_flume_2006 structure.
+When the wave reflects back towards the boundary excessive values
+for momentum are generated leading to degenerate timesteps.
 
-
-Duncan Gray, GA - 2007 
+Ole Nielsen and Duncan Gray, GA - 2008 
 """
 
@@ -23 +24 @@
 
 from anuga.abstract_2d_finite_volumes.util import file_function
-#from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
+from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 
@@ -49 +50 @@
 #-------------------------------------------------------------------------
 friction=0.01
-maximum_triangle_area=0.01
-#maximum_triangle_area=0.1 #OK
 
 basename = 'debug_boundary'
-mesh_filename = basename + '.msh'
 
 boundary_file = 'boundary.tms'
@@ -64 +62 @@
 #-------------------------------------------------------------------------
 
-# Original unstructured mesh ()
-#create_mesh.generate(mesh_filename,
-#                     maximum_triangle_area=maximum_triangle_area)
-
-
-# Simple mesh with same dimensions
-
 xright  = 19.0
 ybottom = 0
 ytop    = 0.45
 xleft = 0.0
-#xslope = 4.0 
 
-#Outline 
-point_sw = [xleft, ybottom]
-point_se = [xright, ybottom]
-point_nw = [xleft, ytop]    
-point_ne = [xright, ytop]
-
-
-create_mesh_from_regions([point_sw,
-                          point_se,
-                          point_ne,
-                          point_nw],
-                         {'wall': [0,1,2], 'wave': [3]},
-                         maximum_triangle_area=maximum_triangle_area,
-                         filename=mesh_filename)
-
-
-#TODO
-#points, elements, boundary = rectangular_cross(m, n, len1=1.0, len2=1.0, origin = (0.0, 0.0)):
+#Very simple mesh
+points, elements, boundary = rectangular_cross(200, 5,
+                                               len1=xright,
+                                               len2=ytop,
+                                               origin = (0.0, 0.0))
 
 
@@ -101 +78 @@
 # Setup computational domain
 #-------------------------------------------------------------------------
-domain = Domain(mesh_filename, use_cache = False, verbose = True)
+#domain = Domain(mesh_filename, use_cache = False, verbose = True)
+domain = Domain(points, elements, boundary, verbose = True)
 print domain.statistics()
 
@@ -125 +103 @@
 Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function)
 Br = Reflective_boundary(domain)
-domain.set_boundary( {'wall': Br, 'wave': Bts} )
+domain.set_boundary( {'top': Br, 'right': Br, 'bottom':Br, 'left': Bts} )
 
 
@@ -134 +112 @@
 for t in domain.evolve(0.1, 15):
     print domain.timestepping_statistics(track_speeds=False)
-    print domain.boundary_statistics(tags=['wave'])
+    print domain.boundary_statistics(tags=['left'])
     
 print 'finished'