^{1}

^{*}

^{1}

The present study involves estimation of open channel flow parameters having different bed materials invoking data of Gradual Varied Flow (GVF). Use of GVF data facilitates estimation of flow parameters. The necessary data base was generated by conducting laboratory. In the present study, the efficacy of the Genetic Algorithm (GA) optimization technique is assessed in estimation of open channel flow parameters from the collected experimental data. Computer codes are developed to obtain optimal flow parameters Optimization Technique. Applicability, adequacy and robustness of the developed code are tested using sets of theoretical data generated by experimental work. A simulation model was developed to compute GVF depths at preselected discrete sections for given downstream head and discharge rate. This model is linked to an optimizer to estimate optimal value of decision variables. The proposed model is employed to a set of laboratory data for three bed materials. Application of proposed model reveals that optimal value of fitting parameter ranges from 1.42 to 1.48 as the material gets finer and optimal decision variable ranges from 0.015 to 0.024. The optimal estimates of Manning’s n of three different bed conditions of experimental channel appear to be higher than the corresponding reported/Strickler’s estimates.

Parameter identification techniques have been widely used in the field of hydrology, meteorology, and oceanography. [

The issue of parameter identification based on the optimal control theories in oceanography can be traced from the early work of [

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. [

Genetic programming (GP―an extension of genetic algorithms to the domain of computer programs [

One of the very few studies which dealt with the identification of compound channel flow parameters is the one by [

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 [

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.

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 (d_{50} = 8 mm and d_{50} = 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.

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 d_{50} = 25 mm, d_{50} = 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

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.

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 ford_{50} = 25 mm, d_{50} = 8 mm and lined concrete respectively.

Genetic Algorithms& Model Description

Optimization techniques were successfully used by [

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

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 y ( x i , Q k , H i ) 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 M^{th} 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.

Governing differential equation used for simulation of GVF is given as:

d y d x = S o − S f 1 − Q 2 T g A 3 (1)

In this equation d y d x is change in depth y with distance x; S_{f} is energy slope and T is top width. S_{f} can be calculated by using Manning’s formula as:

S f = n c 2 Q 2 A 2 R 4 / 3 (2)

where n c 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:

y i + 1 = y i − β Δ x (3)

where, y i + 1 and y i is depth of water level at i + 1^{th} and i^{th} section respectively, Δ x is the distance between them and β is the average slope which is given as follow:

β = d y d x | y i + d y d x | y i + 1 2 (4)

where, ( d y d x | y i ) and ( d y d x | y i + 1 ) are the change in the depth of flow with channel distance x at i^{th} and i+1^{th} section. Equation (3) can be further elaborated using previously mentioned equation as:

y i + 1 = y i − ( S o − n c i 2 Q 2 T 2 y i 10 / 3 1 − Q 2 g T 2 y i 3 S o − n c i 2 Q 2 + T 2 y i + 1 10 / 3 1 − Q 2 g T 2 y i + 1 3 ) Δ x 2 (5)

where, n c i and n c i + 1 are the composite roughness of i^{th} and i + 1^{th} section. An iterative procedure is adopted for the computation of y i + 1 . In this procedure, y i + 1 l + 1 is calculated where l is the number of iteration as:

y i + 1 l + 1 = y i − β Δ x (6)

And the iteration ends when it met the converging criterion, which is given as:

| y i + 1 l + 1 − y i + 1 l | < ∈ (7)

where, ∈ is a constant term. Thus, using the above mentioned approach y i + 1 is computed for each discrete step up to M^{th} step and this leads to the simulation of GVF profiles.

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

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

Simulator

Multi roughness channels are not uncommon in field application of open channel flow hydraulics. [

As mentioned below (For composite channel roughness look

Discharge rates (m^{3}/s) | 9.50 × 10^{−3 } | 1.25 × 10^{−3} | 1.45 × 10^{−3} |
---|---|---|---|

Downstream depths (m) | 0.25 m | 0.30 m | 0.35 m |

Bed materials (d_{50} in mm) | d_{50} = 25 mm | d_{50} = 8 mm | Lined concrete |

0.2 | 0.7 | 1.2 | 1.7 | 2.2 | 2.7 | 3.7 | 4.7 | 5.7 | 6.7 |
---|---|---|---|---|---|---|---|---|---|

7.7 | 8.7 | 9.7 | 10.7 | 12.7 | 14.7 | 16.7 | 18.7 | 20.7 | 22.7 |

n c = ( n 1 ∝ ∗ B + n 2 ∝ ∗ y + n 3 ∝ ∗ y ) 1 / ∝ ( B + 2 y ) 1 / ∝ (8)

where, n c is the composite Manning’s n, n 1 , n 2 , and n 3 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, n c 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. d_{50} = 25 mm, d_{50} = 8 mm and lined concrete as bed materials.

Decision Variables:

( n i , i = 1 , 2 , 3 ) ; and ∝ .

Objective Function:

Min Z = ∑ l 3 ∑ k 3 ∑ i M w i [ y ( x i , Q k , H l ) − y ︷ i k l ] 2 (9)

where, y ( x i , Q k , H l ) and y ︷ i k l are simulated and experimentally measured depth at i^{th} discrete section, k^{th} discharge rate and l^{th} downstream head respectively; M is a subset of the locations where the observed depth is larger than 1.01 x normal depth; w i is the weight assigned to the mismatch at i^{th} location. In the present study the weights are assigned to index the length discretized by the discrete sections. Thus ( w i ) is defined as follows:

w i = ( x i + 1 − x i − 1 ) 2 (10)

Constraint:

1) Following six constraints were assigned to impose upper and lower limits of the segment roughness coefficients (nmax_{i} and nmin_{i}, i = 1, 2, 3).

n m a x i ≥ n m i n i , i = 1 , 2 , 3 (11)

The adopted values of the limits are given in

2) Following three constraints were assigned to ensure realistic relative roughness of the three roughness coefficients.

n 1 ≥ n 2 ≥ n 3 (12)

3) Following constraints was assigned to impose upper and limits of fitting parameters ( ∝ ).

2 ≥ ∝ ≥ 1 (13)

n_{1} | n_{2} | n_{3} | |
---|---|---|---|

nmax_{i} | 0.1 | 0.1 | 0.1 |

nmin_{i} | 0.001 | 0.001 | 0.001 |

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. d_{50} = 25 mm, d_{50} = 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

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).

Bed materials | n_{1} | n_{2} | n_{3} | ∝ | Min Z (m^{2}) |
---|---|---|---|---|---|

d_{50} = 25 mm | 0.032 | 0.016 | 0.018 | 1.42 | 1.16 × 10^{−4 } |

d_{50} = 8 mm | 0.028 | 0.016 | 0.018 | 1.46 | 1.62 × 10^{−4 } |

Lined concrete | 0.027 | 0.015 | 0.017 | 1.48 | 1.09 × 10^{−4 } |

The observed reproduction of GVF profiles for Q = 9.601 × 10^{−3} m^{3}/s and d_{50} = 25 mm is flotted as shown in

The Optimal reproduction of GVF profiles for Q = 9.601 × 10^{−3} m^{3}/s and d_{50} = 25 mm is flotted as shown in

The observed reproduction of GVF profiles for Q = 10.233 × 10^{−3} m^{3}/s and d_{50} = 8 mm is flotted as shown in

The Optimal reproduction of GVF profiles for Q = 10.233 × 10^{−3} m^{3}/s and d_{50} = 8 mm is flotted as shown in

The observed reproduction of GVF profiles for Q = 10.314 × 10^{−3} m^{3}/s and lined concrete is flotted as shown in

The Optimal reproduction of GVF profiles for Q = 10.314 × 10^{−3} m^{3}/s and lined concrete is flotted as shown in

Estimated parameters

The bed roughness ( n 1 ) varies from 0.027 to 0.034 as bed material/condition changes from lined concrete to gravel (d_{50} = 25 mm). The corresponding reported/Strickler’s estimates are given in

The roughness coefficient of glass and GI sheet sides as optimized for various runs are presented in

The estimated roughness coefficients satisfy the known inequality ( n 2 < n 3 ) 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.

Bed material/condition | Reported/Strickler’s Estimation | Optimal estimates |
---|---|---|

d_{50} = 20 mm | 0.0247 | 0.032 |

d_{50} = 6 mm | 0.0202 | 0.028 |

Lined concrete | 0.013 - 0.015 | 0.027 |

Side d_{50 } | d_{50} = 25 mm | d_{50} = 8 mm | Lined concrete | Tabulated values |
---|---|---|---|---|

Glass | 0.016 | 0.016 | 0.015 | 0.010 |

GI sheet | 0.018 | 0.018 | 0.017 | 0.012 |

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.

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 (d_{50} = 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

n c = ( ∑ i = 1 N n i ∝ P i ) 1 / ∝ ( ∑ i = 1 N P i ) 1 / ∝

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.

The authors declare no conflicts of interest regarding the publication of this paper.

Gadissa, E. and Teshome, A. (2018) Estimation of Open Channel Flow Parameters by Using Genetic Algorithm. Open Journal of Optimization, 7, 51-64. https://doi.org/10.4236/ojop.2018.73003