Estimation of Open Channel Flow Parameters by Using Genetic Algorithm ()
1. Introduction
Parameter identification techniques have been widely used in the field of hydrology, meteorology, and oceanography. [1] used the adjoint equation method to identify a profile of Manning’s n in an idealized trapezoidal open channel. In Addition to this, [1] used Lagrangian multipliers and a least square error criterion to estimate roughness coefficients. Research involving the GMS equation traditionally focuses on the determination of the roughness coefficient, (n), under different flow regimes [2] . Optimization techniques were successfully used by [3] , to identify parameters for regular prismatic channels having simple cross-sections.
The issue of parameter identification based on the optimal control theories in oceanography can be traced from the early work of [4] , carried out early an adjoint parameter identification for bottom drag coefficient in a tidal channel. An effective methodology has been proposed to evaluate optimal design of cross sectional area of a channel having composite roughness using Manning’s roughness equation. [5] [6] estimated the bottom friction and water depth in a two-dimensional tidal flow. More recently, [7] used the quasi-Newton method to identify Manning’s roughness coefficients in shallow water flows. Nevertheless, the above studies considered only the case of in-bank flow. Therefore, there is a need to extend the method to out-bank flow, where flood plain roughness will obviously have to be considered.
Genetic algorithms are computationally simple yet powerful search algorithms that seek to produce mathematically the mechanics of natural selection and natural genetics, according to the biological processes of survival and adaptation. [8] [9] identified a constant Manning’s n in an open channel flow with a movable bed. [10] identified the friction parameter in 1D open channel considering the selection of performance function and effect of uncertainty in observed data. And also [10] used a nonlinear least square technique with three types of objective function and identified open channel friction parameters by a modified Gauss-Newton method.
Genetic programming (GP―an extension of genetic algorithms to the domain of computer programs [11] ), a technique generated from the seminal work of numerous researchers in the 1970s and 1980s, generates possible solutions that fit Manning and Albert Strickler. A superior algorithm was proposed for the tree type network which involves the segmentation of channel network into small parts followed by their individual solution using forth order Runge-Kutta method and linking the solution of smaller units to yield the solution of entire channel network by applying shooting method [12] .
One of the very few studies which dealt with the identification of compound channel flow parameters is the one by [13] . In this study, roughness coefficients in the main channel and flood plains were identified as two different parameters using an automatic optimization method. The model was applied to Duong River in Vietnam, where roughness coefficients of the main channel and the flood plain were presented as different constant values as well as polynomial functions of stage. [14] solved the inverse problem of identifying the roughness coefficient in a channel network using the sequential quadratic programming algorithm.
Research involving the GMS equation traditionally focuses on the determination of the roughness coefficient, (n), under different flow regimes and/or for different riverbed materials as even the presence of biological soil crusts can affect the surface roughness, runoff and erodibility of the channel [15] [16] , estimated flood discharges using the Levenberg?Marquardt minimization algorithm.
Therefore, the objectives his study are to: 1) identify open channel flow parameters by using Genetic Algorithm optimization Technique, 2) generate and monitor gradually varied flow profiles corresponding to different bed materials, discharge and ponded depths, 3) Invoking the observed data of the GVF profiles and the linked simulation optimization approach to estimates Manning’s n corresponding to different channel bed materials in the experimental channel, and 4) minimize errors by using optimization methods.
2. Methods and Materials
2.1. Methods
This study was carried out to identify open channel flow parameters by using Genetic Algorithm optimization technique. Manning’s roughness coefficient and other parameters are estimated for different bed materials used (d50 = 8 mm and d50 = 25 mm grain size and Lined concrete bed materials). Also, GVF flow profile is identified. Crank-Nicolson method is used to solve the governing differential equation. Parameter optimization technique is used to find the optimal value of coefficient roughness for three different bed materials. Estimation of roughness coefficient is based on Manning’s equation for estimation of manning roughness coefficient and corresponding manning roughness parameters. This estimation invokes the data of observed GVF profiles and such accounts for different bed materials with the flow depth. Experimental works is done to several sets of data monitored in Hydraulics Laboratory of Civil Engineering Department.
2.2. Materials
Flume
A rectangular tilting flume of length 30 m, width 0.205 m and height 0.50 m was used. The bed of the flume was made up of lined concrete and the other two sides were made up of glass and GI sheet. Discharge was released through an inlet pipe of 0.010 m diameter into the flume. The entrance of the channel was provided with flow suppressors in order to make the flow stable. In order to maintain desired depth of water at the downstream of the channel, a tail gate was fitted at the end of the channel. Water discharging from the tail gate, passed to the sump which was circulated again through a 15 hp centrifugal pump for further experimentation.
Experimental Procedures
The experiments were conducting by adopting the following steps as mentioned below:
Slope Measurement
All the sets of experiment were performed on a particular slope of the channel. The slope was measured by using two steel containers connected with a long rubber tube. Both the containers were placed on the channel bed separated by the rubber tube along the length of the channel. One of the containers placed at higher elevation was filled with water and simultaneously care was taken to remove air bubble from the connecting tube. They are left undisturbed for sufficient amount of time around 24 hours. Then the water levels were measured. Then, slope of the channel was computed.
Sieve Analysis
Sieve analysis was performed to determine the particle size of the material used to create artificial bed roughness. Results of sieve analysis were plotted to investigate the particle size of the bed material used in the present study. Experiments were conducted on two different bed materials. First on one rough bed condition having gravel as a bed particle size d50 = 25 mm, d50 = 8 mm and then on the smooth condition having lined concrete as bed material. The sieve analysis of aggregate is done by a set of standard sieves. The procedure of determination is presented below. Total mass of sample is 1kg (Look Figure 1 and Figure 2).
Calibration of orifice meter
Orifice meter was provided in the inlet pipe for the measurement of discharge. Orifice plate was made up of GI sheet having diameter of 0.06 m and the diameter of inlet pipe was 0.10 m. Ultrasonic flow meter was used for the calibration of coefficient of discharge of orifice meter. Different discharges were noted corresponding to varying head. This result was plotted and the best fitted line was used.
Figure 1. Gradation curve for d50 = 25 mm.
Figure 2. Gradation curve for d50 = 8 mm.
Measurement of water surface profiles
1) Water was released into the rectangular flume by opening the valve of inlet pipe.
2) The desired depth of flow was maintained at the downstream end by operating sluice gate provided at the end of the channel. The depth of water was measured using pointer gauge.
3) After a while when the flow become steady in the channel and the desired depth was maintained at the downstream end, the water surface profile was being measured.
4) Starting from the maintained depth at the downstream end (0.00 m), the water surface profile is measured towards upstream at ten (21) discrete locations that are 0.00 m, 0.20 m, 0.70 m, 1.20 m, 1.70 m, 2.20 m, 2.70 m, 3.70 m, 4.70 m, 5.70 m, 6.70 m, 7.70 m, 8.70 m, 9.70 m, 10.70 m, 12.70 m, 14.70 m, 16.70 m, 18.70 m, 20.70 m and 22.70 m.
5) The above mentioned steps were repeated for three different downstream depths, Discharges rates and bed roughness are mentioned in Tables.
Collection of data
The data obtained for experimental measured water surface profiles corresponding to different bed materials is presented in Tables ford50 = 25 mm, d50 = 8 mm and lined concrete respectively.
Genetic Algorithms& Model Description
Optimization techniques were successfully used by [17] , to identify parameters for regular prismatic channels having simple cross-sections. These researchers used the same optimization algorithm (the so-called “Influence Coefficient” Algorithm) which, mathematically, is closely related to both quasi linearization and the gradient method.
Genetic Algorithm combines survival of the fittest among string structures (used to represent the specific entity to be optimized) with a structured, yet randomized information exchange to form a search algorithm. In every generation, a new set of artificial entities (or strings) is created using bits and pieces of the fittest of the old; an occasional new part is tried for good measure. While randomized, Gas are no simple random walk. They efficiently exploit historical information to speculate on new search points, with expected new performance. GAs, in essence, are heuristic (nonexact), probabilistic (stochastic), combinatorial (discrete), search based optimization technique and the continuing price or performance improvements of computational systems have made them attractive for various types of optimization problems.
In particular, genetic algorithms work very well on mixed (continuous and discrete), combinatorial problems. They are less susceptible to getting “stuck” at local optima than gradient search or other traditional optimization methods; rather they search for global optima. The most important feature of GA is their robust nature and the balance between efficiency and efficacy necessary for survival in many different environments. GAs is theoretically and empirically proven to provide a robust search in complex search space. Usually, a GA for solving a particular problem has five major components, which are succinctly described in the succeeding sections. A general flowchart of GA is shown in Figure 3.
Simulation Model
The optimization problem posed in the preceding section is solved by employing the linked optimization problem. This approach would require development of a model for simulation of GVF depths at preselected discrete sections for given downstream head and discharge rate. Subsequently this simulation model is linked to an optimizer for addressing the optimization problem. Effectively the simulation model would provide the vector of computed depths
appearing in the objective function. The details of the simulation model in the following sections.
Discretization of reach
In the simulation model the entire channel reach is discretized into M small space steps such that depth of water level at Mth step is greater than 1.01 × normal depth.
Governing differential equation
Some advanced numerical techniques like standard forth order Runge-Kutta method, Kutta-Merson method etc. are also used to solve the governing differential equation of the flow.
Figure 3. Flowchart of genetic algorithm.
Governing differential equation used for simulation of GVF is given as:
(1)
In this equation
is change in depth y with distance x; Sf is energy slope and T is top width. Sf can be calculated by using Manning’s formula as:
(2)
where
is composite roughness coefficient given by Equation (8).
Simulation strategy
Crank-Nicolson method is used to solve the governing differential equation mentioned in above section. In this method, depth of water level at next space step is calculated as:
(3)
where,
and
is depth of water level at i + 1th and ith section respectively,
is the distance between them and
is the average slope which is given as follow:
(4)
where,
and
are the change in the depth of flow with channel distance x at ith and i+1th section. Equation (3) can be further elaborated using previously mentioned equation as:
(5)
where,
and
are the composite roughness of ith and i + 1th section. An iterative procedure is adopted for the computation of
. In this procedure,
is calculated where l is the number of iteration as:
(6)
And the iteration ends when it met the converging criterion, which is given as:
(7)
where,
is a constant term. Thus, using the above mentioned approach
is computed for each discrete step up to Mth step and this leads to the simulation of GVF profiles.
2.3. Model Application
Data Base
As discussed, the laboratory channel is rectangular in section with roughening bed material lying on the bed width and the two sides made up of glass and GI sheet. Thus there are three segments of wetted perimeter i.e. bed and two sides one of glass and other of GI sheet. The experiments were performed for three types of bed conditions and several sets of discharge rates downstream head as enumerated in Table 1.
The depths were measured at 20 locations upstream of the control section. The co-ordinates of the locations (measured upstream of the control section) are shown in Table 2.
Simulator
Multi roughness channels are not uncommon in field application of open channel flow hydraulics. [18] , used singular value decomposition to calibrate the Manning’s roughness in one-dimensional (1D) Saint Venant equations. Due to different roughness of wetted perimeter the overall roughness of the channel is given by composite roughness. Composite roughness comprises of individual roughness effect of channel cross section. Seventeen different equations based on several assumptions along with six different techniques to sub divide the channel cross section were given by numerous investigators [19] . The credibility of these equations would be assessed by employing experimental data.
As mentioned below (For composite channel roughness look Figure 4), the experimental channel consists of three types of wetted perimeter; accordingly following equation is used in the simulator for computing the composite roughness
:
Table 1. Data used for experimental measurement of water surface profiles.
Table 2. Co-ordinates of the locations of observed GVF profiles (measured in meter).
(8)
where,
is the composite Manning’s n,
,
, and
are value of Manning’s n for bed and sides respectively. B is bed width and y is the depth of flow. Since composite roughness depends on the depth of the flow, which is not constant in the present scenario. Therefore,
is computed at each section of the water surface flow profile. The value of
is take n as 0.001 in Equation (11).
Optimization
The following problem was solved three times corresponding to different bed conditions i.e. d50 = 25 mm, d50 = 8 mm and lined concrete as bed materials.
Decision Variables:
; and
.
Objective Function:
(9)
where,
and
are simulated and experimentally measured depth at ith discrete section, kth discharge rate and lth downstream head respectively; M is a subset of the locations where the observed depth is larger than 1.01 x normal depth;
is the weight assigned to the mismatch at ith location. In the present study the weights are assigned to index the length discretized by the discrete sections. Thus (
) is defined as follows:
(10)
Constraint:
1) Following six constraints were assigned to impose upper and lower limits of the segment roughness coefficients (nmaxi and nmini, i = 1, 2, 3).
(11)
The adopted values of the limits are given in Table 3.
2) Following three constraints were assigned to ensure realistic relative roughness of the three roughness coefficients.
(12)
3) Following constraints was assigned to impose upper and limits of fitting parameters (
).
(13)
Table 3. Upper and lower limits of roughness coefficients.
Since the reported value of
1.5, a range of 1 to 2 was prescribed.
Linked simulation optimization approach is used to estimate the optimal values of the parameters for three bed conditions i.e. d50 = 25 mm, d50 = 8 mm and lined concrete as bed materials and their corresponding GVF profiles were simulated.
Optimal values
Optimal values of decision variables and their corresponding minimized objective function value for different bed materials are mentioned in Table 4.
Optimal reproduction of GVF profiles
Computed GVF profiles corresponding to the optimal parameter values and the variation of composite roughness are in the following figures. The profile is plotted for three different bed materials corresponding to discharge rates and water depth (as shown in Figures 5-10).
Figure 5. Observed reproduction of GVF profiles (Q = 9.601 × 10−3 m3/s and d50 = 25 mm).
Figure 6. Optimal reproduction of GVF profiles (Q = 9.601 × 10−3 m3/s and d50 = 25 mm).
Table 4. Optimal values of decision variables and objective function.
Figure 7. Observed reproduction of GVF profiles (Q = 10.233 × 10−3 m3/sand d50 = 8 mm).
Figure 8. Optimal reproduction of GVF profiles (Q = 10.233 × 10−3 m3/s and d50 = 8 mm).
Figure 9. Observed reproduction of GVF profiles (Q = 10.314 × 10−3 m3/s and lined concrete).
Figure 10. Optimal reproduction of GVF profiles (Q = 10.314 × 10−3 m3/s and lined concrete).
The observed reproduction of GVF profiles for Q = 9.601 × 10−3 m3/s and d50 = 25 mm is flotted as shown in Figure 5.
The Optimal reproduction of GVF profiles for Q = 9.601 × 10−3 m3/s and d50 = 25 mm is flotted as shown in Figure 6.
The observed reproduction of GVF profiles for Q = 10.233 × 10−3 m3/s and d50 = 8 mm is flotted as shown in Figure 7.
The Optimal reproduction of GVF profiles for Q = 10.233 × 10−3 m3/s and d50 = 8 mm is flotted as shown in Figure 8.
The observed reproduction of GVF profiles for Q = 10.314 × 10−3 m3/s and lined concrete is flotted as shown in Figure 9.
The Optimal reproduction of GVF profiles for Q = 10.314 × 10−3 m3/s and lined concrete is flotted as shown in Figure 10.
Estimated parameters
The bed roughness (
) varies from 0.027 to 0.034 as bed material/condition changes from lined concrete to gravel (d50 = 25 mm). The corresponding reported/Strickler’s estimates are given in Table 5; by using strickler’s equation. It may be seen that optimal roughness estimates are higher than Strickler’s estimates.
The roughness coefficient of glass and GI sheet sides as optimized for various runs are presented in Table 6.
The estimated roughness coefficients satisfy the known inequality (
) and are higher than the tabulated values. This establishes the credibility of the proposed model.
The optimal value of
(fitting parameter) ranges from 1.42 to 1.48, which differs from the reported value i.e. 1.5. The optimal value of
increases as the bed materials get finer.
Reproduction of observed profile
Computed GVF profiles corresponding to the optimal parameter values match quite well with corresponding observed profiles.
Table 5. Reported/Strickler’s estimated optimal estimates for bed materials.
Table 6. Reported/Strickler’s estimates and optimal estimates for sides.
Variability of composite roughness
It can be observed that composite roughness reduces with increase in flow depth. Apparently because of increase in weight age of side resistance, the value of composite roughness increase.
3. Conclusions
This study was carried out to identify open channel flow parameters. Manning’s roughness coefficient and other parameters are estimated for different bed materials used (d50 = 25 mm grain size, 8 mm grain size particles and Lined concrete bed materials). Also, based on the estimated value of Manning roughness coefficient and flow depths, GVF flow profile is identified.
An optimization method is applied to identify the parameters based on Manning formula for estimation of manning roughness coefficient and corresponding manning roughness parameters. This estimation invokes the data of observed GVF profiles and such accounts for different bed materials with the flow depth.
Experimental works are done to several sets of data monitored in Hydraulics Laboratory of Civil Engineering Department. The application led to the following conclusions.
1) The GVF profile computed on the basis of estimated parameters matches quite closely with the corresponding observed profiles.
2) Strickler’s formula underestimates the roughness due to the bed material.
3) The following commonly used formula is calibrated for Manning coefficient estimation
4) The currently documented value of
is 1.5. However, the present work reveals that it varies from 1.42 to 1.48. The value of
generally decreases as the bed material gets coarser.