Developing Semi-Analytical solutions for Saint-Venant Equations in the Uniform Flow Region

آب و خاک

Article 2, Volume 29, Issue 6, January 2015, Pages 1427-1437   

Document Type: Research Article

DOI: 10.22067/jsw.v29i6.32422

Authors

M.M. Heidari^{* 1}; S. Kouchakzadeh²

¹Razi University, Kermanshah

²University of Tehran

Abstract

Introduction: Unsteady flow in irrigation systems is the result of operations in response to changes in water demand that affect the hydraulic performance networks. The increased hydraulic performance needed to recognize unsteady flow and quantify the factors affecting it. Unsteady flow in open channels is governed by the fully dynamic Saint Venant equation, which express the principles of conservation of mass and momentum. Unsteady flow in open channels can be classified into two types: routing and operation-type problems. In the routing problems, The Saint Venant equations are solved to get the discharge and water level in the time series. Also, they are used in the operation problem to compute the inflow at the upstream section of the channel according to the prescribed downstream flow hydrographs. The Saint Venant equation has no analytical solution and in the majority cases of such methods use numerical integration of continuity and momentum equations, and are characterized by complicated numerical procedures that are not always convenient for carrying out practical engineering calculations. Therefore, approximate methods deserve attention since they would allow the solution of dynamic problems in analytical form with enough exactness. There are effective methods for automatic controller synthesis in control theory that provide the required performance optimization. It is therefore important to get simplified models of irrigation canals for control design. It would be even more interesting to have linear models that explicitly depend on physical parameters. Such models would allow one to, handle the dynamics of the system with fewer parameters, understand the impact of physical parameters on the dynamics, and facilitate the development a systematic design method. Many analytical models have been proposed in the literature, Most of them have been obtained in the frequency domain by applying Laplace transform to linearized Saint-Venant equations. The got transcendental function can then be simplified using various methods to get a model expressed as a rational function of s (the Laplace variable), possibly including a time delay. It is therefore important to develop simple analytical models able to accurately reproduce the dynamic behavior of the system in realistic conditions.
Materials and Methods: Changes in water demand can create transient flow in irrigation networks. The Saint Venant equations are the equations governing open channel flow when unsteady flow propagates. In this research, the finite volume method using the time splitting scheme was employed to develop a computer code for solving the one dimensional unsteady flow equations. Considering stationary regime and small variations around it, the Saint-Venant equations around initial condition was linearized.
The Laplace transform is applied to the linearized saint venant equations, leading to an ordinary differential equation in the space variable x and parameterized by the Laplace variable s. The integration of this equation lead to a transfer matrix, and gives the discharge Q*(x, s) at any location with respect for the upstream discharge. This matrix is coupled with the downstream boundary condition and developed an equation that solved using Simpson integration algorithm. It should be noted numerical solution of developed equation is easier than solving fully dynamic saint venant and is less sensitive to the spatial step and the researcher simply writing code.
Results and Discussion: Froud Number (F), variation of inflow discharge (ΔQ/Q), and dimensionless parameter of KF2 in which K is the kinematic flow number, are effective factors on accuracy of developed equation. In order to determine the effect of the factors on accuracy of presenting formula, several simulations were performed using numerical model. The presented formula and numerical model were compared for 10, 20 and 30 percent discharge variation and error calculated, the maximum error increases with increasing ΔQ/Q.
To assess the importance of Froud Number and KF2, also several simulations were carried out, the results showed that the maximum error in the development equation for various Froud Number and KF2>1, is less than 3.8 percent.
Conclusion: Using Laplace transform to the saint venant equations and with respect to upstream and downstream boundary a formula for routing discharge presented. Investigation of the applicability range of presenting formula and cognitive effective factors on accuracy is necessary. So, the finite volume method using the time splitting scheme was employed to develop a computer code for solving the one dimensional unsteady flow equation. Then some tests of unsteady flow were simulated and verified the equations. The results showed that the maximum error increases with decreasing KF2 and increasing the rate of sudden changes of discharge. The maximum error in the presented formula for all tests with KF2>1, less than 3.8 percent.

Keywords

Irrigation systems; Laplace Transform; Linearized Saint-Venant Equations; Unsteady flow

References

1- Chaudry M.H. 2008. Open Channel Flow. Springer, NewYork.

2- Cunge J.A., Holly F.M., and Vervey J.A. 1980. Practical Aspect of Computational River Engineering. Pitman Advanced Publishing Program, London.

3- Fan P., and Li J.C. 2006. Diffusive wave solutions for open channel flows with uniform and concentrated lateral inflow. Advances in Water Resources, 29: 1000-1019.

4- Fennema R.J., and Chaudry M.H. 1987. Simulation of one dimensional dam break flows. Journal of Hydraulic Research, 25(1):41-51.

5- Ghobadian R., and Gholzari S. 2012. Mathematical model to simulate unsteady flow in irrigation networks with internal boundary conditions (a case study, Darband irrigation network). Journal of Water and Soil, 26(1):150-161. (in Persian with English abstract)

6- Heidari M.M., Kouchakzadeh S., Darvishi E., and Hoorfar A.H. 2012. Determination of applicability range of Travel time and response time formulas in Irrigation networks. Iranian Journal of Irrigation and Drainage, 5:120-132. (in Persian with English abstract)

7- Islam A., Raghuwanshi N.S., and Singh R. 2008. Development and application of hydraulic simulation model for irrigation canal network. Journal of Irrigation and Drainage Engineering, 134(1): 49-59.

8- Jeffrey A. 2002. Advanced Engineering Mathematics. Harcourt Press, Newcastle.

9- Kazezyilmaz-Alhan C.M., and Medina M.A. 2007. Kinematic and diffusion waves: analytical and numerical solutions to overland and channel flow. Journal of Hydraulic Engineering, 133(2): 217-228.

10- Kouchakzadeh S., Monem M.J., and Kasbdouz S. 1999. Determination of optimal water distribution policy in an irrigation network (case study: Qurichay network). Iranian Journal of Agriculture Science, 30(2): 369-378. (in Persian with English abstract)

11- Litrico X., and Fromion V. 2004. Frequency modeling of open channel Flow. Journal of Hydraulic Engineering, 130(8): 806-815.

12- Litrico X., Pomet J.B., and Guinot V. 2010. Simplified nonlinear modeling of river flow routing. Advances in Water Resources, 33: 1015-1023.

13- Namin M.M., Bohluly A., Alahdin S., Ranjbaran L., and Mohammadzadeh M. 2007. Flow simulation in river networks by time-splitting method and comparison with mike11 model. 7th International River Engineering Conference Shahid Chamran University, Ahwaz, Iran. (in Persian with English abstract)

14- Onizuka K., and Odai S.N. 1998. Burgers’ equation model for unstedy flow in open channels. Journal of Hydraulic Engineering, 124(5): 509-513.

15- Preissmann A. 1960. Propagation des intumescences dans les canaux et rivières. p. 433-442. First Congress of the French Association for Computation, AFCAL Grenoble, France.

16- Robinson J.C. 2004. An Introduction to Ordinary Differential Equations. Cambridge University Press, Cambridge.

17- Wang G.T., and Chen S. 2003. A semianalytical solution of the saint-venant equations for channel flood routing. Water Resources Research, 39(4): 1-10.