Changeset 8353 for trunk/anuga_work/development/gareth/tests/channel_floodplain/channel_floodplain1.py

Timestamp:

Mar 7, 2012, 9:49:06 PM (13 years ago)

Author:

davies

Message:

Updates to balanced_dev and associated tests

File:

• trunk/anuga_work/development/gareth/tests/channel_floodplain/channel_floodplain1.py (modified) (6 diffs)

trunk/anuga_work/development/gareth/tests/channel_floodplain/channel_floodplain1.py

@@ -10 +10 @@
 import numpy
 from anuga.structures.inlet_operator import Inlet_operator
-
+from anuga import *
+#from swb_domain import domain
 #from anuga import *
 #from balanced_basic import *
@@ -23 +24 @@
 floodplain_width = 14.0 # Model domain width
 floodplain_slope = 1./300.
-chan_initial_depth = 0.8 # Initial depth of water in the channel
+chan_initial_depth = 0.65 # Initial depth of water in the channel
 chan_bankfull_depth = 1.0 # Bankfull depth of the channel
 chan_width = 10.0 # Bankfull width of the channel
 bankwidth = 2. # Width of the bank regions -- note that these protrude into the channel
 man_n=0.03 # Manning's n
-l0 = 2.009 # Length scale associated with triangle side length in channel (min_triangle area = 0.5*l0^2)
+l0 = 1.000 # Length scale associated with triangle side length in channel (min_triangle area = 0.5*l0^2)
 
 assert chan_width < floodplain_width, \
@@ -81 +82 @@
 
 
-domain.set_name('channel_floodplain1_test') # Output name
-domain.extrapolate_velocity_second_order=True
+domain.set_name('channel_floodplain1_bal_dev_lowvisc') # Output name
+#domain.set_store_vertices_uniquely(True)
+#domain.use_edge_limiter=False
+#domain.extrapolate_velocity_second_order=False
 #------------------------------------------------------------------------------
 #
@@ -139 +142 @@
 
 # Define inlet operator 
+flow_in_yval=100.0
 if True:
-    line1 = [ [floodplain_width/2. - chan_width/2., 0.],\
-              [floodplain_width/2. + chan_width/2., 0.] \
+    line1 = [ [floodplain_width/2. - chan_width/2., flow_in_yval],\
+              [floodplain_width/2. + chan_width/2., flow_in_yval] \
               ]
     Qin = 0.5*(floodplain_slope*(chan_width*chan_initial_depth)**2.*man_n**(-2.)\
@@ -191 +195 @@
 Br = anuga.Reflective_boundary(domain) # Solid reflective wall
 Bt = anuga.Transmissive_boundary(domain) # Transmissive boundary
-#Bout_sub = anuga.Dirichlet_boundary( \
-#        [-floodplain_length*floodplain_slope - chan_bankfull_depth + \
-#        chan_initial_depth, 0., 0.]) #An outflow boundary for subcritical steady flow
+
+
+Bout_sub = anuga.Dirichlet_boundary( \
+        [-floodplain_length*floodplain_slope - chan_bankfull_depth + \
+        chan_initial_depth, 0., 0.]) #An outflow boundary for subcritical steady flow
 
 def outflow_stage_boundary(t):
@@ -226 +232 @@
 #------------------------------------------------------------------------------
 
-for t in domain.evolve(yieldstep=1.0, finaltime=800.0):
+for t in domain.evolve(yieldstep=2.0, finaltime=3200.0):
     print domain.timestepping_statistics()
+    xx=domain.quantities['ymomentum'].centroid_values
+    dd=(domain.quantities['stage'].centroid_values - domain.quantities['elevation'].centroid_values)
+    dd=dd*(dd>0.)
+
+    tmp = xx/(dd+1.0e-06)*(dd>0.0)
+    print tmp.max(), tmp.argmax(), tmp.min(),  tmp.argmin()
+
+    # Compute flow through cross-section -- check that the inflow boundary condition is doing its job
+    # This also provides another useful steady-state check
+    if( numpy.floor(t/100.) == t/100. ):
+        print '#### COMPUTING FLOW THROUGH CROSS-SECTIONS########'
+        s1 = domain.get_flow_through_cross_section([[0., floodplain_length-300.0], [floodplain_width, floodplain_length-300.0]])
+        s2 = domain.get_flow_through_cross_section([[0., floodplain_length-1.0], [floodplain_width, floodplain_length-1.0]])
+        
+        print 'Cross sectional flows: ', s1, s2 
+        print '$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$'
+
+    #vv = domain.get_flow_through_cross_section
     #print domain.quantities['ymomentum'].get_integral()
     #print 'Qin = ', Qin