def test_boyd_6(self):
"""test_boyd_6
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g=9.81
culvert_slope=1 # Downward
inlet_depth=1.50
inlet_velocity= 4.0
outlet_depth=0.8
outlet_velocity=4.0
culvert_length=10.0
culvert_width=3.60
culvert_height=1.20
culvert_type='box'
manning=0.013
sum_loss=1.5
inlet_specific_energy=inlet_depth + 0.5*inlet_velocity**2/g
z_in = 10.0
z_out = 10.0-culvert_length*culvert_slope/100
E_in = z_in+inlet_depth + 0.5*inlet_velocity**2/g
E_out = z_out+outlet_depth + 0.5*outlet_velocity**2/g
delta_total_energy = E_in-E_out
Q, v, d = boyd_generalised_culvert_model(inlet_depth,
outlet_depth,
inlet_velocity,
outlet_velocity,
inlet_specific_energy,
delta_total_energy,
g,
culvert_length,
culvert_width,
culvert_height,
culvert_type,
manning,
sum_loss)
#print ('%s,%.3f'%('SPEC_E = ',inlet_specific_energy))
#print ('%s,%.3f'%('Delta E = ',delta_total_energy))
#print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d))
#print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',13.546, 3.136, 1.20))
assert num.allclose(Q, 13.546, rtol=1.0e-2) #inflow
assert num.allclose(v, 3.136, rtol=1.0e-2) #outflow velocity
assert num.allclose(d, 1.20, rtol=1.0e-2) #depth at outlet used to calc v
def test_boyd_5(self):
"""test_boyd_5
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g=9.81
culvert_slope=10 # Downward
inlet_depth=1.50
inlet_velocity= 1.0
outlet_depth=2.5
outlet_velocity=0.5
culvert_length=10.0
culvert_width=0.0
culvert_height=1.20
culvert_type='circle'
manning=0.013
sum_loss=1.5
inlet_specific_energy=inlet_depth + old_div(0.5*inlet_velocity**2,g)
z_in = 10.0
z_out = 10.0-old_div(culvert_length*culvert_slope,100)
E_in = z_in+inlet_depth + old_div(0.5*inlet_velocity**2,g)
E_out = z_out+outlet_depth + old_div(0.5*outlet_velocity**2,g)
delta_total_energy = E_in-E_out
Q, v, d = boyd_generalised_culvert_model(inlet_depth,
outlet_depth,
inlet_velocity,
outlet_velocity,
inlet_specific_energy,
delta_total_energy,
g,
culvert_length,
culvert_width,
culvert_height,
culvert_type,
manning,
sum_loss)
#print ('%s,%.3f'%('SPEC_E = ',inlet_specific_energy))
#print ('%s,%.3f'%('Delta E = ',delta_total_energy))
#print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST05 Q-v-d',Q,v,d))
#print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',0.759, 0.671, 1.20))
assert num.allclose(Q, 0.759, rtol=1.0e-2) #inflow
assert num.allclose(v, 0.671, rtol=1.0e-2) #outflow velocity
assert num.allclose(d, 1.20, rtol=1.0e-2) #depth at outlet used to calc v
def test_boyd_1(self):
"""test_boyd_1
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g=9.81
inlet_depth=0.150
outlet_depth=0.15
inlet_velocity=1.00
outlet_velocity=0.5
culvert_length=10.0
culvert_width=3.6
culvert_height=1.20
culvert_type='box'
manning=0.013
sum_loss=1.5
inlet_specific_energy=inlet_depth + 0.5*inlet_velocity**2/g
culvert_slope=1 # % Downward
z_in = 10.0
z_out = -culvert_length*culvert_slope/100
E_in = z_in+inlet_depth + 0.5*inlet_velocity**2/g
E_out = z_out+outlet_depth + 0.5*outlet_velocity**2/g
delta_total_energy = E_in-E_out
inlet_specific_energy=inlet_depth + 0.5*inlet_velocity**2/g
Q, v, d = boyd_generalised_culvert_model(inlet_depth,
outlet_depth,
inlet_velocity,
outlet_velocity,
inlet_specific_energy,
delta_total_energy,
g,
culvert_length,
culvert_width,
culvert_height,
culvert_type,
manning,
sum_loss)
#print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST01 Q-v-d',Q,v,d))
#print('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.5526, 1.146, 0.1339))
assert num.allclose(Q, 0.5526, rtol=1.0e-1) #inflow
assert num.allclose(v, 1.146, rtol=1.0e-1) #outflow velocity
assert num.allclose(d, 0.1339, rtol=1.0e-1) #depth at outlet used to calc v
def test_boyd_1(self):
"""test_boyd_1
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g=9.81
inlet_depth=0.150
outlet_depth=0.15
inlet_velocity=1.00
outlet_velocity=0.5
culvert_length=10.0
culvert_width=0.0
culvert_height=1.20
culvert_type='circle'
manning=0.013
sum_loss=1.5
inlet_specific_energy=inlet_depth + old_div(0.5*inlet_velocity**2,g)
culvert_slope=10.0 # % Downward
z_in = 10.0
z_out = old_div(-culvert_length*culvert_slope,100)
E_in = z_in+inlet_depth + old_div(0.5*inlet_velocity**2,g)
E_out = z_out+outlet_depth + old_div(0.5*outlet_velocity**2,g)
delta_total_energy = E_in-E_out
inlet_specific_energy=inlet_depth + old_div(0.5*inlet_velocity**2,g)
Q, v, d = boyd_generalised_culvert_model(inlet_depth,
outlet_depth,
inlet_velocity,
outlet_velocity,
inlet_specific_energy,
delta_total_energy,
g,
culvert_length,
culvert_width,
culvert_height,
culvert_type,
manning,
sum_loss)
#print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST01 Q-v-d',Q,v,d))
#print('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.113, 0.297, 0.443))
assert num.allclose(Q, 0.113, rtol=1.0e-1) #inflow
assert num.allclose(v, 0.297, rtol=1.0e-1) #outflow velocity
assert num.allclose(d, 0.443, rtol=1.0e-1) #depth at outlet used to calc v
def test_boyd_2(self):
"""test_boyd_2
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g=9.81
culvert_slope=10 # Downward
inlet_depth=0.500
outlet_depth=0.700
inlet_velocity=1.0
outlet_velocity=0.50
culvert_length=10.0
culvert_width=0.0
culvert_height=1.20
culvert_width=0.0
culvert_type='circle'
manning=0.013
sum_loss=1.5
inlet_specific_energy=inlet_depth + 0.5*inlet_velocity**2/g
z_in = 0.0
z_out = -culvert_length*culvert_slope/100
E_in = z_in+inlet_depth + 0.5*inlet_velocity**2/g
E_out = z_out+outlet_depth + 0.5*outlet_velocity**2/g
delta_total_energy = E_in-E_out
Q, v, d = boyd_generalised_culvert_model(inlet_depth,
outlet_depth,
inlet_velocity,
outlet_velocity,
inlet_specific_energy,
delta_total_energy,
g,
culvert_length,
culvert_width,
culvert_height,
culvert_type,
manning,
sum_loss)
#print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST02 Q-v-d',Q,v,d))
#print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.585, 0.873, 0.687))
assert num.allclose(Q, 0.585, rtol=1.0e-1) #inflow
assert num.allclose(v, 0.873, rtol=1.0e-1) #outflow velocity
assert num.allclose(d, 0.687, rtol=1.0e-1) #depth at outlet used to calc v
def Xtest_boyd_00(self):
"""test_boyd_00
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g=9.81
culvert_slope=0.1 # Downward
inlet_depth=0.2
outlet_depth=0.0
inlet_velocity=0.0,
outlet_velocity=0.0,
culvert_length=4.0
culvert_width=1.2
culvert_height=0.75
culvert_type='box'
manning=0.013
sum_loss=0.0
inlet_specific_energy=inlet_depth #+0.5*v**2/g
z_in = 0.0
z_out = -culvert_length*culvert_slope/100
E_in = z_in+inlet_depth # +
E_out = z_out+outlet_depth # +
delta_total_energy = E_in-E_out
Q, v, d = boyd_generalised_culvert_model(inlet_depth,
outlet_depth,
inlet_velocity,
outlet_velocity,
inlet_specific_energy,
delta_total_energy,
g,
culvert_length,
culvert_width,
culvert_height,
culvert_type,
manning,
sum_loss)
#print Q, v, d
assert num.allclose(Q, 0.185, rtol=1.0e-3)
def test_boyd_0(self):
"""test_boyd_0
This tests the Boyd routine with data obtained from ??? by Petar Milevski
This test is the only one that passed in late February 2009
"""
g=9.81
culvert_slope=0.1 # Downward
inlet_depth=2.0
outlet_depth=0.0
inlet_velocity=0.0,
outlet_velocity=0.0,
culvert_length=4.0
culvert_width=1.2
culvert_height=0.75
culvert_type='box'
manning=0.013
sum_loss=0.0
inlet_specific_energy=inlet_depth #+0.5*v**2/g
z_in = 0.0
z_out = -culvert_length*culvert_slope/100
E_in = z_in+inlet_depth # +
E_out = z_out+outlet_depth # +
delta_total_energy = E_in-E_out
Q, v, d = boyd_generalised_culvert_model(inlet_depth,
outlet_depth,
inlet_velocity,
outlet_velocity,
inlet_specific_energy,
delta_total_energy,
g,
culvert_length,
culvert_width,
culvert_height,
culvert_type,
manning,
sum_loss)
#print Q, v, d
assert num.allclose(Q, 3.118, rtol=1.0e-3)
def Xtest_boyd_2(self):
"""test_boyd_2
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g=9.81
culvert_slope=0.01 # Downward
inlet_depth=1.135
outlet_depth=0.0
culvert_length=4.0
culvert_width=0.75
culvert_height=0.75
culvert_type='pipe'
manning=0.013
sum_loss=1.5
inlet_specific_energy=inlet_depth #+0.5*v**2/g
z_in = 0.0
z_out = -culvert_length*culvert_slope/100
E_in = z_in+inlet_depth #+ 0.5*v**2/g
E_out = z_out+outlet_depth #+ 0.5*v**2/g
delta_total_energy = E_in-E_out
Q, v, d = boyd_generalised_culvert_model(inlet_depth,
outlet_depth,
inlet_velocity,
outlet_velocity,
inlet_specific_energy,
delta_total_energy,
g,
culvert_length,
culvert_width,
culvert_height,
culvert_type,
manning,
sum_loss)
print Q, v, d
assert num.allclose(Q, 1.00, rtol=1.0e-2) #inflow
assert num.allclose(v, 2.59, rtol=1.0e-2) #outflow velocity
assert num.allclose(d, 0.563, rtol=1.0e-2) #depth at outlet used to calc v
def test_boyd_6(self):
"""test_boyd_6
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g = 9.81
culvert_slope = 10 # Downward
inlet_depth = 1.50
inlet_velocity = 4.0
outlet_depth = 0.80
outlet_velocity = 4.0
culvert_length = 10.0
culvert_width = 3.60
culvert_height = 1.20
culvert_type = 'box'
manning = 0.013
sum_loss = 1.5
inlet_specific_energy = inlet_depth + 0.5 * inlet_velocity**2 / g
z_in = 10.0
z_out = 10.0 - culvert_length * culvert_slope / 100
E_in = z_in + inlet_depth + 0.5 * inlet_velocity**2 / g
E_out = z_out + outlet_depth + 0.5 * outlet_velocity**2 / g
delta_total_energy = E_in - E_out
Q, v, d = boyd_generalised_culvert_model(
inlet_depth, outlet_depth, inlet_velocity, outlet_velocity,
inlet_specific_energy, delta_total_energy, g, culvert_length,
culvert_width, culvert_height, culvert_type, manning, sum_loss)
#print ('%s,%.3f'%('SPEC_E = ',inlet_specific_energy))
#print ('%s,%.3f'%('Delta E = ',delta_total_energy))
#print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d))
#print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',15.537, 3.597, 1.20))
assert num.allclose(Q, 15.537, rtol=1.0e-2) #inflow
assert num.allclose(v, 3.597, rtol=1.0e-2) #outflow velocity
assert num.allclose(d, 1.20,
rtol=1.0e-2) #depth at outlet used to calc v
def Xtest_boyd_10(self):
"""test_boyd_9
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g = 9.81
culvert_slope = 0.1 # Downward
inlet_depth = 1.901
outlet_depth = 1.5
culvert_length = 4.0
culvert_width = 0.75
culvert_height = 0.75
culvert_type = 'pipe'
manning = 0.013
sum_loss = 1.5
inlet_specific_energy = inlet_depth #+0.5*v**2/g
z_in = 0.0
z_out = -culvert_length * culvert_slope / 100
E_in = z_in + inlet_depth #+ 0.5*v**2/g
E_out = z_out + outlet_depth #+ 0.5*v**2/g
delta_total_energy = E_in - E_out
Q, v, d = boyd_generalised_culvert_model(inlet_depth, outlet_depth,
inlet_specific_energy,
delta_total_energy, g,
culvert_length, culvert_width,
culvert_height, culvert_type,
manning, sum_loss)
print Q, v, d
assert num.allclose(Q, 1.00, rtol=1.0e-2) #inflow
assert num.allclose(v, 2.204, rtol=1.0e-2) #outflow velocity
assert num.allclose(d, 0.76,
rtol=1.0e-2) #depth at outlet used to calc v
def Xtest_boyd_00(self):
"""test_boyd_00
This tests the Boyd routine with data obtained from ??? by Petar Milevski
"""
# FIXME(Ole): This test fails (20 Feb 2009)
g = 9.81
culvert_slope = 0.1 # Downward
inlet_depth = 0.2
outlet_depth = 0.0
inlet_velocity = 0.0,
outlet_velocity = 0.0,
culvert_length = 4.0
culvert_width = 1.2
culvert_height = 0.75
culvert_type = 'box'
manning = 0.013
sum_loss = 0.0
inlet_specific_energy = inlet_depth #+0.5*v**2/g
z_in = 0.0
z_out = -culvert_length * culvert_slope / 100
E_in = z_in + inlet_depth # +
E_out = z_out + outlet_depth # +
delta_total_energy = E_in - E_out
Q, v, d = boyd_generalised_culvert_model(
inlet_depth, outlet_depth, inlet_velocity, outlet_velocity,
inlet_specific_energy, delta_total_energy, g, culvert_length,
culvert_width, culvert_height, culvert_type, manning, sum_loss)
#print Q, v, d
assert num.allclose(Q, 0.185, rtol=1.0e-3)
def test_boyd_0(self):
"""test_boyd_0
This tests the Boyd routine with data obtained from ??? by Petar Milevski
This test is the only one that passed in late February 2009
"""
g = 9.81
culvert_slope = 0.1 # Downward
inlet_depth = 2.0
outlet_depth = 0.0
inlet_velocity = 0.0,
outlet_velocity = 0.0,
culvert_length = 4.0
culvert_width = 1.2
culvert_height = 0.75
culvert_type = 'box'
manning = 0.013
sum_loss = 0.0
inlet_specific_energy = inlet_depth #+0.5*v**2/g
z_in = 0.0
z_out = -culvert_length * culvert_slope / 100
E_in = z_in + inlet_depth # +
E_out = z_out + outlet_depth # +
delta_total_energy = E_in - E_out
Q, v, d = boyd_generalised_culvert_model(
inlet_depth, outlet_depth, inlet_velocity, outlet_velocity,
inlet_specific_energy, delta_total_energy, g, culvert_length,
culvert_width, culvert_height, culvert_type, manning, sum_loss)
#print Q, v, d
assert num.allclose(Q, 3.118, rtol=1.0e-3)
