LittleInch thanks for your reply, you have right, politeness is a must, will strive to improve.
I deleted my previeus post, prefered a shorter message.
Very interesting what you image heppen: "flow ... keeps changing cross section area". It's like a wave-flow regime, even if the feeded flowrate is costant.
In the first range of time there is a very dynamic step with generation of waves, and the time needed to have the stady state could be very long. But when achieved, do you think it could be imagened a regular decreasing level profile?
Bimr, the Chezy/Manning's equation [Vavg = (1/n) x R^2/3 x S^1/2] is not directly applicable because the slope of pipe S is 0%, for this reason i'm looking for a deepening, without bother finite difference method. Circular geometry is a complication.