Changeset 2014 for inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py

Timestamp:

Nov 9, 2005, 4:46:14 PM (19 years ago)

Author:

steve

Message:

 

File:

• inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py (modified) (6 diffs)

inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py

@@ -29 +29 @@
 Q = 9.985/1000.0
 wh0 = Q/(2.0*pi*0.1)
-h0 = bed[2] + 0.005
-wh1 = -Q/(2.0*pi*0.5)
-h1 = 0.562
+stage0 = bed[2] + 0.005
+wh1 = -Q/(2.0*pi*0.503)
+stage1 = 0.562
 Manning = 0.009
 
@@ -90 +90 @@
 
 #---------------------------------
-# Function for initial water depth
+# Function for initial water elevation
+# (stage)
 def level(x,y):
     z = bed_z(x,y)
     n = x.shape[0]
-    w = 0*x
+    stage = 0*x
     for i in range(n):
-        w[i] = h0
-        h = w[i] - z[i]
-    return w
+        stage[i] = stage0
+    return stage
 
 
-def outflow_height(t):
-    return [h1, 0 , 0]
+#def outflow_stage(t):
+#    return [stage1, 0 , 0]
 
 
@@ -109 +109 @@
 #---------------------------
 # Set up boundary conditions
-DD_BC_INNER = Dirichlet_Discharge_boundary(domain, h0, wh0)
-DD_BC_OUTER = Dirichlet_Discharge_boundary(domain, h1, wh1)
+DD_BC_INNER = Dirichlet_Discharge_boundary(domain, stage0, wh0)
+DD_BC_OUTER = Dirichlet_Discharge_boundary(domain, stage1, wh1)
 
 domain.set_boundary({'inner': DD_BC_INNER, 'outer': DD_BC_OUTER})
@@ -129 +129 @@
 
 
-# from realtime_visualisation_new import Visualiser
-# vxmom = Visualiser(domain,title='xmomentum',scale_z=10.0)
-# vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
+#domain.initialise_visualiser()
+#    #domain.visualiser.coloring['stage'] = True
+#domain.visualiser.scale_z['stage'] = 2.0
+#domain.visualiser.scale_z['elevation'] = 0.05
+#
+#
+#from realtime_visualisation_new import Visualiser
+##vxmom = Visualiser(domain,title='xmomentum',scale_z=10.0)
+##vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
 
-w = domain.quantities['stage']
+stage = domain.quantities['stage']
 
 #----------
@@ -142 +148 @@
 
 t0 = time.time()
-for t in domain.evolve(yieldstep = 1.0, finaltime = 1000.):
+for t in domain.evolve(yieldstep = .1, finaltime = 3.0):
     domain.write_time()
 
@@ -148 +154 @@
     radial_momentum = domain.create_quantity_from_expression(exp)
 
-    outer_stage = w.get_values(location='centroids', indices=[typical_outer[0]])
-    outer_radial_mom = radial_momentum.get_values(location='centroids', indices=[typical_outer[0]])
-    
-    inner_stage = w.get_values(location='centroids', indices=[typical_inner[0]])
-    inner_radial_mom = radial_momentum.get_values(location='centroids', indices=[typical_inner[0]])
+    print 'outer stage      ', stage.get_values(location='vertices',
+                                            indices=[typical_outer[0]])
+    print '      radial mom ', \
+          radial_momentum.get_values(location='centroids',
+                                     indices=[typical_outer[0]])
 
-    print inner_stage
-    print inner_radial_mom
-    print outer_stage
-    print outer_radial_mom
+    print 'inner stage      ', stage.get_values(location='centroids',
+                                            indices=[typical_inner[0]])
+    print '      radial mom ', \
+          radial_momentum.get_values(location='centroids',
+                                     indices=[typical_inner[0]])
 
 #    f.write('time = %25.15e wall clock time %g \n' % (domain.time, time.time()))
     f.write('%10.3f %25.15e %25.15e %25.15e %25.15e \n' % (domain.time, inner_stage, inner_radial_mom, outer_stage, outer_radial_mom))
 
+    f.write('time = %25.15e wall clock time %g \n' % (domain.time, time.time()))
+    f.write('%g \n' % stage.get_values(location='centroids',
+                                   indices=[typical_outer[0]])[0])
 
 f.close()