Research Article | Open Access
Farshad Ehsani, Seyed Ghorban Hosseini, Hossein Soury, "Modeling of Unsteady Flow through the Canals by Semiexact Method", Modelling and Simulation in Engineering, vol. 2014, Article ID 495715, 12 pages, 2014. https://doi.org/10.1155/2014/495715
Modeling of Unsteady Flow through the Canals by Semiexact Method
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.
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  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 . 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 , 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.
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.
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  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 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 .