Changeset 8683

Timestamp:

Jan 30, 2013, 7:34:10 PM (11 years ago)

Author:

steve

Message:

Pressure now calculated via 1/2 g h^{2 \rho_w}

Location:

trunk/anuga_core/validation_tests/Tests/Experimental_data/dam_break_Yeh_Petroff

Files:

• numerical_Yeh_Petroff.py (modified) (2 diffs)
• plot_results_force.py (modified) (1 diff)

trunk/anuga_core/validation_tests/Tests/Experimental_data/dam_break_Yeh_Petroff/numerical_Yeh_Petroff.py

@@ -55 +55 @@
 from anuga.utilities.argparsing import parse_standard_args
 alg, cfl = parse_standard_args()
-domain.set_flow_algorithm(alg)
+
+alg = '2_0_limited'
+cfl = 0.9
+print alg, cfl
+domain.set_flow_algorithm('2_0_limited')
 domain.set_CFL(cfl)
 domain.set_minimum_allowed_height(0.01)
@@ -87 +91 @@
 domain.set_quantity('stage', stage)
 domain.set_quantity('elevation',elevation)
-domain.set_quantity('friction', 0.0)
+domain.set_quantity('friction', 0.03)
 
 #-----------------------------------------------------------------------------

trunk/anuga_core/validation_tests/Tests/Experimental_data/dam_break_Yeh_Petroff/plot_results_force.py

@@ -113 +113 @@
 for i in range(N):
     for k in range(len(indicesL)):
-        Sh[i] += p2.stage[i,indicesL[k]]
-        Sh[i] -= p2.stage[i,indicesR[k]]
-force = 9.8*Sh
+        Sh[i] += p2.stage[i,indicesL[k]]**2
+        Sh[i] -= p2.stage[i,indicesR[k]]**2
+
+fudge_factor = 1.0
+length = 0.12
+g = 9.8
+rho_w = 1024
+average_pressure = 0.5*g*Sh/len(indicesL)*rho_w
+
+force = fudge_factor*average_pressure*length