Hello!
I have tested the new tube fonction with the 039_bridge test case.
I am a bit suprised by the result shown.
Here a the main parameters (there are 2 identical tubes):
- S1=8,66
- S2=7,10
- RD1=RD2=1
- HAU=1,5
- LAR=2,5
- Ce=0,5
- Cs=1,0
- L=0,2
The computed discharge through each tube is ca 25 m3/s.
This result is obtained by the following formula (buse.f):
Q = SECT * SQRT( 2.D0*GRAV*(S1-S2)/(CE1+L) )
I would have expected the last one to be used instead, valid when both inlet and outlet are submerged (which is the case here):
Q = SECT * SQRT( 2.D0*GRAV*(S1-S2)/(L+CS2+CE1) ), which here only gives a discharge of ca 16 m3/s using the computed (S1-S2).
When looking at the first equation, not using CS2 coefficient with the given (S1-S2) implies an equivalent orifice coefficient of 1,19, which seems strange (V cannot be > (2gh)^0,5), whereas it is 0,77 in the last one, which is logical.
Being not familiar with programming, I am not sure how water surface profiles and each formula are linked.
A last question: why didn't you use the same submerged weir formula than the one used in loinoy.f?
Thanks in advance,
PL