Modelling and Simulation in Engineering

Volume 2014 (2014), Article ID 495715, 12 pages

http://dx.doi.org/10.1155/2014/495715

Research Article

Modeling of Unsteady Flow through the Canals by Semiexact Method

Farshad Ehsani,¹ Seyed Ghorban Hosseini,² and Hossein Soury³

¹Department of Mechanical Engineering, Yasouj Branch, Islamic Azad University, Yasouj, Iran

²Department of Chemistry, Malek Ashtar University of Technology, P.O. Box 16765-3454, Tehran, Iran

³Department of Mechanical Engineering, Tarbiat Modares University, Tehran, Iran

Received 23 July 2013; Accepted 13 November 2013; Published 20 February 2014

Academic Editor: Ligang Wu

Copyright © 2014 Farshad Ehsani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The study of free-surface and pressurized water flows in channels has many interesting application, one of the most important being the modeling of the phenomena in the area of natural water systems (rivers, estuaries) as well as in that of man-made systems (canals, pipes). For the development of major river engineering projects, such as flood prevention and flood control, there is an essential need to have an instrument that be able to model and predict the consequences of any possible phenomenon on the environment and in particular the new hydraulic characteristics of the system. The basic equations expressing hydraulic principles were formulated in the 19th century by Barre de Saint Venant and Valentin Joseph Boussinesq. The original hydraulic model of the Saint Venant equations is written in the form of a system of two partial differential equations and it is derived under the assumption that the flow is one-dimensional, the cross-sectional velocity is uniform, the streamline curvature is small and the pressure distribution is hydrostatic. The St. Venant equations must be solved with continuity equation at the same time. Until now no analytical solution for Saint Venant equations is presented. In this paper the Saint Venant equations and continuity equation are solved with homotopy perturbation method (HPM) and comparison by explicit forward finite difference method (FDM). For decreasing the present error between HPM and FDM, the st.venant equations and continuity equation are solved by HAM. The homotopy analysis method (HAM) contains the auxiliary parameter ħ that allows us to adjust and control the convergence region of solution series. The study has highlighted the efficiency and capability of HAM in solving Saint Venant equations and modeling of unsteady flow through the rectangular canal that is the goal of this paper and other kinds of canals.

1. Introduction

Hydraulics has a long tradition of providing a scientific basis for engineering applications [1, 2]. Firstly, conceptual models were designed starting from empirical relations obtained from field observations or model scale experiments. Lately, mathematics started playing an important role not only to describe the properties of these relations but also to formulate analytical solutions of particular model situations in order to capture the essential features of those phenomena. Actually, the research and the applications in the field of computational fluid hydraulics and fluid dynamics evolved with the advent of electronic computers. The first applications in computational hydraulics concerned programming analytical formulae rather than deriving generic numerical schemes and techniques based on physical principles like conservation laws for mass and momentum. Later developments extended the research and the applications in this field towards simulating complicated flow phenomena in arbitrarily shaped geometries.

The basic equations expressing hydraulic principles were formulated in the 19th century by Barre de Saint Venant and Valentin Joseph Boussinesq. The original hydraulic model of the Saint Venant equations [3] is written in the form of a system of two partial differential equations and it is derived under the assumption that the flow is one-dimensional, the cross-sectional velocity is uniform, the streamline curvature is small, and the pressure distribution is hydrostatic [4]. One-dimensional flows do not actually exist in nature, but the equations remain valid provided that the flow is approximately one-dimensional: as pointed out by Steffler and Jin [5], they are inappropriate to analyze free surface flow problems with horizontal length scales close to flow depth.

The goal of paper is to present a semiexact solution of Saint Venant equations. These partial differential equations are solved simultaneously with continuity equation by homotopy perturbation method [6–17]. Then they are solved by homotopy analysis method [18–22]. In the end both of above methods are compared with the finite difference method as a numerical method that it is shown the capability and suitability of HAM in solving Saint Venant equations.

2. Symbols

The following symbols are used in this derivation:

  = The cross-sectional area of the section  = depth of flow at the section  = elevation of surface above a datum at the section  = mean velocity at the section  = discharge at the section  = width of the top of the section  = position of the section measured from the upstream end  = time  = acceleration due to gravity  = mass density of the fluid.

Other symbols are defined in the text at the point when they are introduced.

3. The Derivation of the Continuity Equation

Consider a short length of channel and assuming that there is no lateral inflow, then

This has the partial derivative since is changing with both and time . Now the volume of water between Sections 1 and 2 is increasing as a rate of

where is the top width. As cross-sectional area then this is equivalent to

The terms are equal in magnitude but of opposite sign, so

The continuity equation is

So the Continuity equation for rectangular canals becomes

4. The Derivation of Dynamic or Momentum Equation

By applying Newton’s 2nd law to our elemental length of channel

Also

By using partial differential theory, the value of acceleration at the direction of water motion becomes

By applying Newton’s 2nd law to our elemental length of channel and assuming the hydrostatics pressure

The First term of left side of (10) showed the difference between the hydrostatics forces and second term of left side showed the friction due of water motion.

From (10), it can be deduced that

Substituting (8) at (11) gives

Equation (12) can be rearranged to the below form

That

  = hydraulic radius  = energy slope  = acceleration slope  = friction slope.

By using (7)

By Adding (14) and (12)

This is the Saint Venant equation.

Then by using the homotopy perturbation method the simultaneous partial differential equations (15) and (6) are solved and obtained the semiexact solution for deep and velocity of flow.

5. Basic Idea of Homotopy Perturbation Method

The homotopy perturbation method is combination of the classical perturbation technique and homotopy technique. We start with the following nonlinear differential equation to explain the basic idea of the HPM:

Subject to boundary condition

where is a general differential operator, is a boundary operator, is a known analytical function, is the boundary of domain , and denotes differentiation along the normal drawn outwards from . The operator can, generally speaking, be divided into two parts: a linear part and a nonlinear part .

Equation (17) therefore can be rewritten as follows:

In case the nonlinear equation (16) has no “small parameter”, we can construct the following homotopy:

where is the parameter of homotopy. According to the homotopy perturbation method, the approximation solution of (19) can be expressed as a series of the power of ; that is,

For St. Venant equation (15) and continuity equation (6) we obtain

These equations are solved with each other and step by step.

For the continuity equation (6)

For the momentum equation (15)

Changing from to , consider as follows:

The best approximation for the solution is

Comparison of the expressions with the same powers of the parameter gives the following:

The above partial differential equations must be supplemented by conditions ensuring a uniqueness of the solution. For above equations we assume the following conditions in the example.

Consider the following Saint Venant and continuity equations with initial value problem:

For an example to show the capability of HPM and HAM the boundary and initial conditions are considered as below:

is the order of in (24). By assuming the above initial conditions a solution for equations system is as follows:

Important Point. We assume that because after 8iteration the sum of homotopy perturbation sentences converged

6. Basic Idea of Homotopy Analysis Method (HAM)

In this section we employ the homotopy analysis method to the discussed problem. To show the basic idea, let us consider the following differential equation:

where is a nonlinear operator, denotes independent variable, and is an unknown function, respectively. For Simplicity, we ignore all boundary or initial conditions, which can be treated in the similar way. By means of generalizing the traditional homotopy method, Liao [18] constructs the so-called zero-order deformation equation

where is the embedding parameter, is a nonzero auxiliary parameter, is an auxiliary function, is an auxiliary linear operator, (τ) is an initial guess of (τ), and ) is an unknown function, respectively; it is important that one has great freedom to choose auxiliary things in HAM. Obviously, when and , it holds

respectively. Thus as increases from to 1, the solution varies from the initial guess to the solution . Expanding in the Taylor series with respect to , one has

If the auxiliary linear operator, the initial guess, the auxiliary parameter , and the auxiliary function must be chosen so properly, the series (33) converges at , and one has

which must be one of solutions of original nonlinear equation, as proved by Liao; as and , (31) becomes

This is mostly used in HPM, whereas the solution can be obtained directly without using the Taylor series. According to (36), the governing equation can be deduced from the zero-order deformation equation (31). The vector is defined as

Differentiating (31) times with respect to the embedding parameter , and then setting and finally dividing them by we will have the so-called th-order deformation equation as

where

It should be emphasized that for is governed by the linear equation (37) with the linear boundary conditions coming from the original problem, which can be easily solved using the symbolic computation software.

7. Application of Homotopy Analysis Method

For Saint Venant equation (15) and continuity equation (6) these obtained

These equations are solved with each other and step by step.

For the continuity equation (6)

and for the momentum equation (15)

In the following, we apply HAM to solve St. Venant and continuity equations in the canals. To obey both the rule of solution expression and the rule of the coefficient periodicity, the corresponding auxiliary function can be determined uniquely. Then

8. Following the Homotopy Analysis Method

Consider

By using first order explicit finite difference method (FDM), (15) and (6) are solved simultaneously and the capability of HPM and HAM is compared with FDM (finite difference method) as a numerical method for solving Saint Venant equations.

9. Result and Conclusion

In Figures 1, 2, 3, 4, 5, and 6, two-dimensional plot for the comparison of HPM, HAM, and FDM for the solutions (, ) and (, ) for different values of , , and is shown. It could be observe a good accuracy between these methods. In the end we have shown , and , by three-dimensional plots in Figures 7, 8, 9, and 10 for , and observed that the results converged after 8 sentences of HPM and there is no need to write the term of homotopy perturbation method. Then it’s shown that HAM is converged after just 5 sentences of and .

Figure 1: The comparison of HPM, HAM, and FDM for the solution for different values of and ; .

Figure 2: The comparison of HPM, HAM, and FDM for the solution for different values of and ; .

Figure 3: The comparison of HPM, HAM, and FDM for the solution for different values of and ; .

Figure 4: The comparison of HPM, HAM, and FDM for the solution for different values of and ; .

Figure 5: The comparison of HPM, HAM, and FDM for the solution for different values of and ; .

Figure 6: The comparison of HPM, HAM, and FDM for the solution for different values of and ; .

(a)

(b)

Figure 7: Three-dimensional plot for the solution of Saint Venant equations via homotopy analysis method obtained (a) = + + + + (b) = + + + + + .

(a)

(b)

Figure 8: Three-dimensional plot for the solution of Saint Venant equations via homotopy analysis method obtained (a) = + + + + (b) = + + + + + .

(a)

(b)

Figure 9: Three-dimensional plot for the solution of Saint Venant equations via homotopy perturbation method obtained (a) = + + + + + + + (b) = + + + + + + + + .

(a)

(b)

Figure 10: Three-dimensional plot for the solution of Saint Venant equations via homotopy analysis method obtained (a) = + + + + + + + (b) = + + + + + + + + .

In this paper, HPM has been successfully applied to finding the solutions of Saint Venant equations in rectangular canals. The obtained solution is compared with finite difference method. The homotopy perturbation method had a little difference with finite difference method, so we solve Saint Venant equations by homotopy analysis method and by changing , we could control error and observed that the accuracy of HAM is more than HPM for solution of Saint Venant equations. All the figures show that the results of the homotopy analysis method are in approximate agreement with FDM.

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.

References

A. T. Ippen, Estuary and Coastline Hydrodynamics, McGraw-Hill, 1966.

H. Rouse, Engineering Hydraulics, John Wiley & Sons, New York, NY, USA, 1950.

B. de Saint Venant, “Théorie du mouvement non-permanent des eaux avec application aux crues des riviéres atàl’introdution des marées dans leur lit,” Comptes Rendus de l'Académie des Sciences, vol. 73, pp. 147–154, 1871.

B. C. Yen, “Open-channel flow equations revisited,” ASCE Engineering Mechanics Division, vol. 99, no. 5, pp. 979–1009, 1973.

P. M. Steffler and Y.-C. Jin, “Depth averaged and moment equations for moderately shallow free surface flow,” Journal of Hydraulic Research, vol. 31, no. 1, pp. 5–17, 1993.

J.-H. He, “Coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.

J. H. He, “Application of topological technology to construction of a perturbation system for a strongly nonlinear equation,” Journal of the Juliusz Schauder Center, vol. 20, pp. 77–83, 2002.

J.-H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003.

J.-H. He, “The homotopy perturbation method for nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004.

J. Biazar and H. Ghazvini, “Exact solutions for non-linear Schrödinger equations by He's homotopy perturbation method,” Physics Letters A, vol. 366, no. 1-2, pp. 79–84, 2007.

J. Biazar, R. Ansari, K. Hosseini, and P. Gholamin, “Solution of the linear and non-linear schrödinger equations using homotopy perturbation and Adomian Decomposition methods,” International Mathematical Forum, vol. 38, pp. 1891–1897, 2008.

J. Biazar and H. Aminikhah, “Study of convergence of homotopy perturbation method for systems of partial differential equations,” Computers and Mathematics with Applications, vol. 58, no. 11-12, pp. 2221–2230, 2009.

B. Ganjavi, H. Mohammadi, D. D. Ganji, and A. Barari, “Homotopy perturbation method and variational iteration method for solving Zakharov-Kuznetsov equation,” American Journal of Applied Sciences, vol. 5, no. 7, pp. 811–817, 2008.

M. A. Abdou and A. A. Soliman, “Variational iteration method for solving Burger's and coupled Burger's equations,” Journal of Computational and Applied Mathematics, vol. 181, no. 2, pp. 245–251, 2005.

S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, 1992.

F. Ehsani and F. Ehsani, A Comparison of Homotopy Perturbation Method (HPM), Adomian’s Decomposition Method (ADM) and Homotopy Analysis Method (HAM) in Solving Gardner Equation, 2012

F. Ehsani, F. Ehsani, A. Hadi, and N. Hadi, Analytical Solution of Phi-Four Equation, 2013.

S. J. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, Chapman and Hall/CRC Press, Boca Raton, Fla, USA, 2003.

S.-J. Liao, “On the analytic solution of magnetohydrodynamic flows of non-Newtonian fluids over a stretching sheet,” Journal of Fluid Mechanics, no. 488, pp. 189–212, 2003.

S. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004.

S. Liao, “A new branch of solutions of boundary-layer flows over an impermeable stretched plate,” International Journal of Heat and Mass Transfer, vol. 48, no. 12, pp. 2529–2539, 2005.

D. Q. Kern and D. A. Kraus, Extended Surface Heat Transfer, McGraw-Hill, New York, NY, USA, 1972.