Changeset 1834

Timestamp:

Sep 15, 2005, 5:06:08 PM (18 years ago)

Author:

steve

Message:

Added inflow BC Dirichlet_Discharge_boundary

Location:

inundation

Files:

• analytical solutions/Analytical_solution_circular_hydraulic_jump.py (added)
• analytical solutions/circular.pt (added)
• pyvolution/shallow_water.py (modified) (1 diff)

inundation/pyvolution/shallow_water.py

@@ -1009 +1009 @@
 
 
+class Dirichlet_Discharge_boundary(Boundary):
+    """Returns same momentum conserved quantities as
+    those present in its neighbour volume. Sets stage
+
+    Underlying domain must be specified when boundary is instantiated
+    """
+
+    def __init__(self, domain = None, h0=None, wh0=None):
+        Boundary.__init__(self)
+
+        if domain is None:
+            msg = 'Domain must be specified for this type boundary'
+            raise msg
+
+        if h0 is None:
+            h0 = 0.0
+
+        if wh0 is None:
+            wh0 = 0.0
+
+        self.domain   = domain
+        self.h0  = h0
+        self.wh0 = wh0
+
+    def __repr__(self):
+        return 'Dirichlet_Discharge_boundary(%s)' %self.domain
+
+    def evaluate(self, vol_id, edge_id):
+        """Dirichlet_Discharge_boundary boundaries return the edge momentum
+        values of the volume they serve.
+        """
+
+        normal = self.domain.get_normal(vol_id,edge_id)
+        q = [self.h0, -self.wh0*normal[0], -self.wh0*normal[1]]
+        return q
+
+
+        #FIXME: Consider this (taken from File_boundary) to allow
+        #spatial variation
+        #if vol_id is not None and edge_id is not None:
+        #    i = self.boundary_indices[ vol_id, edge_id ]
+        #    return self.F(t, point_id = i)
+        #else:
+        #    return self.F(t)