Changeset 2379

Timestamp:

Feb 13, 2006, 12:12:04 PM (19 years ago)

Author:

nicholas

Message:

 

Location:

development/momentum_sink

Files:

• build.tsh (added)
• building_generator.py (added)
• building_generator.pyc (added)
• compare_sww.py (modified) (1 diff)
• create_buildings.py (modified) (2 diffs)
• create_buildings.pyc (added)
• friction_block.py (modified) (2 diffs)
• friction_block.pyc (added)
• number (added)
• number.py (added)
• project.pyc (added)
• project_friction.py (added)
• project_friction.pyc (added)
• rotate_buildings.py (added)
• rotate_buildings.pyc (added)
• run_buildings.py (modified) (3 diffs)
• run_friction.py (modified) (6 diffs)
• run_rotate.py (added)
• test.tsh (added)

development/momentum_sink/compare_sww.py

@@ -11 +11 @@
 gauge_locations = [[10,10],[20,20]]
 
-quantities = ['stage']
+quantities = ['stage','elevation']
 gauges = file_function(swwfile,
                   quantities = quantities,

development/momentum_sink/create_buildings.py

@@ -1 +1 @@
 
-
-
+import Numeric
+import math
 # add inundation dir to your pythonpath
 from pmesh.mesh import Mesh
@@ -18 +18 @@
 
     # Boundary of problem
-    WidtH = 200 # width of boudary in metres
+    WidtH = 500 # width of boudary in metres
     #W = WidtH/8
     #L = W
-    outer_polygon = [[0,0],[2*WidtH,0],[2*WidtH,WidtH],[0,WidtH]]
+    outer_polygon = [[0,0],[1.5*WidtH,0],[1.5*WidtH,WidtH],[0,WidtH]]
     print outer_polygon
     m.add_region_from_polygon(outer_polygon, tags={'wall':[0,1,2], 'wave':[3]})
 
     # inner polygons => building boundaries
-    depth = 15 # depth of building side to oncoming wave
+    depth = 18 # depth of building side to oncoming wave
     wh = depth/2 
-    breadth = 15 # breadth of building
+    breadth = 18 # breadth of building
     lh = breadth/2 
+      
     print "building footprint"
     print depth * breadth , "m^2"
     block = 625
     BL = block**0.5
-    print BL
-    #[12.5,37.5,62.5,87.5,112.5,137.5,162.5,187.5]
-    D = WidtH/2 + BL/2
-    B = BL/2
-    for D in [112.5,137.5,162.5,187.5]:
-        for B in [12.5,37.5,62.5,87.5,112.5,137.5,162.5,187.5]:
+    ForDep = (0.2*WidtH) + (BL/2)
+    RearDep = 1.2*WidtH 
+    porosity = breadth/BL
+    print porosity, " Building porosity"
+    
+    Breadths = Numeric.arrayrange( (BL/2), WidtH, (BL))
+    print Breadths, "Breadths"
+    Depths = Numeric.arrayrange( ForDep, RearDep, BL )
+    print Depths, "Depths"
+
+    for D in Depths:
+        Breadths = Breadths + BL/2
+        for B in Breadths:
             polygon = [[D-wh,B-lh],[D+wh,B-lh],[D+wh,B+lh],[D-wh,B+lh]]
             m.add_hole_from_polygon(polygon, tags={'wall':[0,1,2,3]})# Adds holes with reflective boundaries.        
+        
+
             #print polygon
-            
-            print B
-
-        print D
-
     m.generate_mesh(maximum_triangle_area=maximum_triangle_area)
     triangle_count = m.get_triangle_count()

development/momentum_sink/friction_block.py

@@ -18 +18 @@
 
     # Boundary of problem
-    WidtH = 202 # width of boudary in metres
+    WidtH = 200 # width of boudary in metres
     #W = WidtH/8
     #L = W
@@ -28 +28 @@
     #Boundary
     dict = {}
-    dict['points'] = [[0,0],[0.5*WidtH,0],[1.5*WidtH,0],[2*WidtH,0],[2*WidtH,WidtH] \
-                      ,[1.5*WidtH,WidtH],[0.5*WidtH,WidtH],[0,WidtH]]
+    dict['points'] = [[0,0],[0.2*WidtH,0],[1.2*WidtH,0],[3.5*WidtH,0],[3.5*WidtH,WidtH] \
+                      ,[1.2*WidtH,WidtH],[0.2*WidtH,WidtH],[0,WidtH]]
     
     dict['segments'] = [[0,1], [1,2], [2,3], [3,4],

development/momentum_sink/run_buildings.py

@@ -25 +25 @@
 t0 = time.time()
 
-meshname, triagle_count = cache(create_mesh,(10),
+meshname, triagle_count = cache(create_mesh,(100),
                                 {'mesh_file':meshname,
                                  'triangles_in_name':True}
@@ -55 +55 @@
 
 #Setup Initial Conditions
-domain.set_quantity('friction', 0)
+domain.set_quantity('friction', 0.01)
 domain.set_quantity('stage', 0)
 
@@ -65 +65 @@
 Br = Reflective_boundary(domain)
 Bt = Transmissive_boundary(domain)
-Bd = Dirichlet_boundary([6,0,0])
+Bd = Dirichlet_boundary([10,0,0])
 Bw = Time_boundary(domain=domain,
                    f=lambda t: [(60<t<660)*4, 0, 0])
 
 domain.set_boundary( {'wall': Br,'wave': Bd, 'back': Bt, 'exterior':Bd} ) 
-domain.set_region(Set_region('mound', 'friction', 3, location='unique vertices'))
+
 #Evolve
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 40):
+for t in domain.evolve(yieldstep = 1, finaltime = 20):
     domain.write_time()      
 

development/momentum_sink/run_friction.py

@@ -6 +6 @@
 import time
 
-import project
+import project_friction
 
 from pyvolution.pmesh2domain import pmesh_to_domain_instance 
@@ -21 +21 @@
 from pmesh.mesh import importMeshFromFile
 
-meshname = project.meshname
-outputname = project.outputname
+meshname = project_friction.meshname
+outputname = project_friction.outputname
 
 t0 = time.time()
@@ -47 +47 @@
 
 
-domain.set_name(project.basename + '_%d' %triagle_count)
-domain.set_datadir(project.outputdir)
+domain.set_name(project_friction.basename + '_%d' %triagle_count)
+domain.set_datadir(project_friction.outputdir)
 domain.store = True
 domain.quantities_to_be_stored = ['stage']
@@ -56 +56 @@
 
 #Setup Initial Conditions
-domain.set_quantity('friction', 0)
+domain.set_quantity('friction', 0.01)
 domain.set_quantity('stage', 0)
-domain.set_region(Set_region('mound', 'friction', 700, location='unique vertices'))
+domain.set_region(Set_region('mound', 'friction', 0.8)) #, location='unique vertices'))
 #Setup Boundary Conditions
 print domain.get_boundary_tags()
@@ -66 +66 @@
 Br = Reflective_boundary(domain)
 Bt = Transmissive_boundary(domain)
-Bd = Dirichlet_boundary([6,0,0])
+Bd = Dirichlet_boundary([36,0,0])
 Bw = Time_boundary(domain=domain,
                    f=lambda t: [(60<t<660)*4, 0, 0])
@@ -75 +75 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 40):
+for t in domain.evolve(yieldstep = 1, finaltime = 80):
     domain.write_time()