Changes between Initial Version and Version 1 of InternalBoundaries

Timestamp:

Jun 9, 2010, 3:40:25 PM (15 years ago)

Author:

James Hudson

Comment:

--

InternalBoundaries

@@ +1 @@
+= Buildings, Bridges and Dams =
+
+A new feature in ANUGA 1.2 lets you create create arbitrary geometry to represent solid structures. These take the form of holes within the mesh, and can be given arbitrary boundary properties, just like the edges of the domain.
+
+The demo code can be found within demos/buildings.py, part of which is reproduced here:
+
+{{{
+length = 50
+width = 10
+resolution = 0.15 # make this number smaller to make the simulation more accurate
+
+# Create the "world" as a long, skinny channel
+boundary_poly = poly_from_box(0, length, 0, width)
+
+# Place 3 buildings downstream
+building_polys = [  poly_from_box(10, 15, 2.5, 7.5),        # upstream box
+                    poly_from_box(22.5, 27.5, 1.5, 6.5),    # middle box
+                    poly_from_box(35, 40, 3.5, 8.5)]        # downstream box
+
+# create a domain mesh, with 3 building holes in it
+domain = anuga.create_domain_from_regions(boundary_poly,
+                                            boundary_tags={'left': [0],
+                                                   'bottom': [1],
+                                                   'right': [2],
+                                                   'top': [3]},
+                                            maximum_triangle_area = resolution,
+                                            mesh_filename = 'building.msh',
+                                            interior_holes = building_polys,
+                                            use_cache=True, # to speed it up
+                                            verbose=True)   # log output on
+}}}
+
+Note the new parameter to create_domain_from_regions:
+{{{
+interior_holes = building_polys
+}}}
+This lets you specify a list of polygons to be cut out of the mesh.
+
+The internal boundaries have their own boundary tag:
+{{{exterior}}}
+which can be set just like any other boundary.
+