Changeset 1876

Timestamp:

Oct 6, 2005, 3:48:58 PM (19 years ago)

Author:

steve

Message:

fixed problem with BC for circular hydraulic jump.

File:

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

inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py

@@ -29 +29 @@
 wh0 = Q/(2.0*pi*0.1)
 h0 = bed[2] + 0.005
+wh1 = -Q/(2.0*pi*0.5)
 h1 = 0.562
 Manning = 0.009
@@ -44 +45 @@
 # Set a default tagging
 
+
 for id, face in domain.boundary:
     domain.boundary[(id,face)] = 'outer'
     point = domain.get_vertex_coordinate(id,(face+1)%3)
     radius2 = point[0]*point[0] + point[1]*point[1]
+    typical_outer = (id,face)
     if radius2 < 0.1:
         domain.boundary[(id,face)] = 'inner'
+        typical_inner = (id,face)
 
 
@@ -104 +108 @@
 #---------------------------
 # Set up boundary conditions
-DD_BC = Dirichlet_Discharge_boundary(domain, h0, wh0)
-TM_BC = Transmissive_Momentum_Set_Stage_boundary(domain, outflow_height)
-D_BC =  Dirichlet_boundary([h1, 0, 0])
+DD_BC_INNER = Dirichlet_Discharge_boundary(domain, h0, wh0)
+DD_BC_OUTER = Dirichlet_Discharge_boundary(domain, h1, wh1)
 
-domain.set_boundary({'inner': DD_BC, 'outer': D_BC})
+domain.set_boundary({'inner': DD_BC_INNER, 'outer': DD_BC_OUTER})
 
 #------------------
 # Order of accuracy
-domain.default_order = 2
-domain.CFL = 1.0
+domain.default_order = 1
+domain.CFL = 0.75
 #domain.beta_w = 0.5
 #domain.beta_h = 0.2
@@ -128 +131 @@
 
 from realtime_visualisation_new import Visualiser
-vxmom = Visualiser(domain,title='xmomentum',scale_z=10.0)
-vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
+#vxmom = Visualiser(domain,title='xmomentum',scale_z=10.0)
+#vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
 
 
@@ -136 +139 @@
 import time
 
-print 'PROBLEM WITH INSTABILITY AFTER t'
 t0 = time.time()
 for t in domain.evolve(yieldstep = .01, finaltime = 10):
     domain.write_time()
-    vxmom.update_quantity('xmomentum')
-    vymom.update_quantity('ymomentum')
+    #vxmom.update_quantity('xmomentum')
+    #vymom.update_quantity('ymomentum')
+    print 'outer stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_outer[0]])
+    print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2 +
+                                   domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2)
+    print 'inner stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_inner[0]])
+    print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2 +
+                                   domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2)
+
 
 print 'That took %.2f seconds' %(time.time()-t0)
-