Changeset 7996

Timestamp:

Sep 2, 2010, 5:39:00 PM (14 years ago)

Author:

habili

Message:

Confirming if the correct Culvert Velocity is being used, that is the resultant of Delta_Total_Energy

Location:

trunk/anuga_core/source/anuga/structures

Files:

• boyd_box_operator.py (modified) (7 diffs)
• structure_operator.py (modified) (3 diffs)

trunk/anuga_core/source/anuga/structures/boyd_box_operator.py

@@ -28 +28 @@
                  use_momentum_jet=True,
                  use_velocity_head=True,
+                 description=None,
                  verbose=False):
                      
@@ -39 +40 @@
                                                        manning,
                                                        enquiry_gap,
+                                                       description,
                                                        verbose)            
         
@@ -120 +122 @@
 
         # Stats
-        self.discharge  = outflow_extra_height*self.outflow.get_area()/timestep
-        self.velocity = self.discharge/outlet_depth/self.width
+        self.discharge  = Q#outflow_extra_height*self.outflow.get_area()/timestep
+        self.velocity = barrel_speed#self.discharge/outlet_depth/self.width
 
         new_outflow_height = self.outflow.get_average_height() + outflow_extra_height
@@ -185 +187 @@
 
             if self.use_velocity_head :
-                driving_energy = self.inflow.get_enquiry_specific_energy()
+                self.driving_energy = self.inflow.get_enquiry_specific_energy()
             else:
-                driving_energy = self.inflow.get_enquiry_height()
+                self.driving_energy = self.inflow.get_enquiry_height()
 
             height = self.culvert_height
             width = self.culvert_width
             flow_width = self.culvert_width
-
-            Q_inlet_unsubmerged = 0.540*g**0.5*width*driving_energy**1.50 # Flow based on Inlet Ctrl Inlet Unsubmerged
-            Q_inlet_submerged = 0.702*g**0.5*width*height**0.89*driving_energy**0.61  # Flow based on Inlet Ctrl Inlet Submerged
+            # intially assume the culvert flow is controlled by the inlet
+            # check unsubmerged and submerged condition and use Min Q
+            # but ensure the correct flow area and wetted perimeter are used
+            Q_inlet_unsubmerged = 0.544*g**0.5*width*self.driving_energy**1.50 # Flow based on Inlet Ctrl Inlet Unsubmerged
+            Q_inlet_submerged = 0.702*g**0.5*width*height**0.89*self.driving_energy**0.61  # Flow based on Inlet Ctrl Inlet Submerged
 
             # FIXME(Ole): Are these functions really for inlet control?
@@ -202 +206 @@
                 if dcrit > height:
                     dcrit = height
-                flow_area = width*dcrit
+                    flow_area = width*dcrit
+                    perimeter= 2.0*(width+dcrit)
+                else: # dcrit < height
+                    flow_area = width*dcrit
+                    perimeter= 2.0*dcrit+width
                 outlet_culvert_depth = dcrit
                 case = 'Inlet unsubmerged Box Acts as Weir'
-            else:
+            else: # Inlet Submerged but check internal culvert flow depth
                 Q = Q_inlet_submerged
-                flow_area = width*height
-                outlet_culvert_depth = height
+                dcrit = (Q**2/g/width**2)**0.333333
+                if dcrit > height:
+                    dcrit = height
+                    flow_area = width*dcrit
+                    perimeter= 2.0*(width+dcrit)
+                else: # dcrit < height
+                    flow_area = width*dcrit
+                    perimeter= 2.0*dcrit+width
+                outlet_culvert_depth = dcrit
                 case = 'Inlet submerged Box Acts as Orifice'
 
             dcrit = (Q**2/g/width**2)**0.333333
-
+            # May not need this .... check if same is done above
             outlet_culvert_depth = dcrit
             if outlet_culvert_depth > height:
@@ -223 +238 @@
                 perimeter = width+2*outlet_culvert_depth
                 case = 'INLET CTRL Culvert is open channel flow we will for now assume critical depth'
-
-            if self.delta_total_energy < driving_energy:
+            # Initial Estimate of Flow for Outlet Control using energy slope 
+            #( may need to include Culvert Bed Slope Comparison)
+            hyd_rad = flow_area/perimeter
+            culvert_velocity = math.sqrt(self.delta_total_energy/((self.sum_loss/2/g)+(self.manning**2*self.culvert_length)/hyd_rad**1.33333))
+            Q_outlet_tailwater = flow_area * culvert_velocity
+            
+            
+            if self.delta_total_energy < self.driving_energy:
                 # Calculate flows for outlet control
 
@@ -252 +273 @@
                     log_to_file(self.log_filename, s)
 
-                # Outlet control velocity using tail water
+                # Final Outlet control velocity using tail water
                 culvert_velocity = math.sqrt(self.delta_total_energy/((self.sum_loss/2/g)+(self.manning**2*self.culvert_length)/hyd_rad**1.33333))
                 Q_outlet_tailwater = flow_area * culvert_velocity

trunk/anuga_core/source/anuga/structures/structure_operator.py

@@ -25 +25 @@
                  manning,
                  enquiry_gap,
+                 description,
                  verbose):
         
@@ -42 +43 @@
         self.manning = manning
         self.enquiry_gap = enquiry_gap
+        self.description = description
         self.verbose = verbose
 
         self.discharge = 0.0
         self.velocity = 0.0
+        self.delta_total_energy = 0.0
+        self.driving_energy = 0.0
         
         self.__create_exchange_polygons()
@@ -156 +160 @@
 
         message = '---------------------------\n'
-        message += 'Structure report:\n'
+        message += 'Structure report for structure %s:\n' % self.description
         message += '--------------------------\n'
         message += 'Discharge [m^3/s]: %.2f\n' % self.discharge
         message += 'Velocity  [m/s]: %.2f\n' % self.velocity
-#        message += 'Total boundary outflow [m^3/s]: %.2f\n' % total_boundary_outflow
+        message += 'Inlet Driving Energy %.2f\n' % self.driving_energy
+        message += 'delta total energy %.2f\n' % self.delta_total_energy
 #        message += 'Net boundary flow by tags [m^3/s]\n'
 #        for tag in boundary_flows: