# Correction on Chang et al. equation factors of suspended load estimation (case study: Central Alborz Rivers of Iran).

IntroductionYang [5]: Rivers always have erosion and sediment transport, thus in river hydraulic and its morphology, survey of sediment transport capacity of stream and sediment transport mechanism are of high importance. Sediment of rivers transported in two general forms: suspended load and bed load and in most of natural rivers, sediment mostly transported as suspended load.

Shafai Bajestan and Ostad Asgari [3]: because in most of the rivers there is not any sedimentography and hydrometry station for sediment evaluation and also because of difficulty in sediment measurement, and need to experts and costly professional equipments, researchers offered various equation for simpler and better measurement of sediment in rivers and channels, in the recent years.

Some of these methods evaluate just suspended load, some do only bed load and some others evaluate collection of them, called bed materials. The majority of these methods are formed based on laboratory works and always their application is in question.

Shafai Bajestan [2]: there are different equation for suspended load determination such as Lane and Kalinske, Einstein, Brooks, Chang, Simons and Richardson, Bagnold, Toffaleti, Samaga and so on.

Any of these equations have various factors for determining suspended sediment amount and in most of the cases researchers offer graphs and integrals for simplifying these factors among which are Einstein, Toffaleti and Chang et al.

The goal of this study is checking the authenticity of [I.sub.1] and [I.sub.2] factors in Chang et al. equation and also offering the best integral and graph for [I.sub.1] and [I.sub.2] factors on this equation.

Materials and methods

In this study for correction of Chang et al. factors, data of Central Alborz Rivers are used.

Chang, Simons and Richardson (1965):

By applying the velocity distribution (du/d[xi] = [U.sub.*]/K[xi][square root of 1 - [xi]]) over the depth, Chang et al. obtained the following expression for the mass transfer coefficient:

[[epsilon].sub.s] = [beta]KD[xi][U.sub.*] [(l - [xi]).sup.1/2] Eq.1

Where, [xi] = y / D and ([U.sub.*] = [(gds).sup.1/2]. Substituting Eq.1 into Eq.2

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Eq.2

Where, [C.sub.a] and C =concentration by weight of sediment in (a) and (y) distances above bed respectively and it obtained from Rouse equation.

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Eq.3

In this relation, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

[Z.sub.2] = 2W/([beta]U,K) and [[xi].sub.a] = a/d.

Then the expression for the suspended load becomes:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Eq.4

Where, [I.sub.1] and [I.sub.2] are integrals that can be evaluated using figs.1 and 2. The transport [q.sub.sw] is measured in weight per volume of water-sediment mixture.

If [q.sub.sw] is expressed in weight per second per unit channel width and [C.sub.a] is concentration by weight, then:

[q.sub.sw] = [gamma]D[C.sub.a][[VI.sub.1] - 2[U.sub.*]/k [I.sub.2]) Eq.5

[FIGURE 1 OMITTED]

Examining Chang et al. equation:

For estimation of suspended load in a river by Chang et al. equation, we need to extract [I.sub.1] and [I.sub.2] factors by use of relative graph or integral and substituting factors in equation.

Kazemi [1]: when data are low, determining the numeral amount of I1 and I2 factors by corresponding Graphs is simply possible, but in this study which is, applied for estimation of suspended load by Chang method in rivers of Central Alborz of Iran, for estimation of suspended load in Jajrod river, the stream simulation was done in 17 cross section of this river and 30-year annual average discharge was used. Totally about 1020 data for selection of [I.sub.1] and [I.sub.2] was obtained. Also in Taleghan river simulation was done in 25 cross section and 30 years annual average discharge was used. Totally about 1500 data for selection Of [I.sub.1] and [I.sub.2] was obtained. Finding these data by graphs need a lot of time and thereupon amount of error will increase.

[FIGURE 2 OMITTED]

Hence, for solving this problem it was suggested that numeral solution of integrals of these factors be used instead of graphs. Therefore, a program was written by FORTRAN language and use of Yang reference to obtain numeral amount of integrals by entering required parameters of [I.sub.1] and [I.sub.2] factors. After writing FORTRAN program, for assurance about the program accuracy, some data were entered in the program and at the same time the answer of data was determined by use of graphs. Comparison between output of FORTRAN program and graphs showed that there is an incompatibility between them. After assurance about authenticity of written program, error Possibility in offered graphs and integrals in different references were examined. For this reason, first Yang reference and then older references of this equation (Simons and Senturk reference) were examined.

Results:

Survey of two mentioned references about Chang et al. equation showed that there are some errors in graphs and integrals of this equation.

Examining graphs and integrals of [I.sub.1] factor in Simons and Senturk reference:

Examining [I.sub.1] factor in Simons and Senturk reference showed that integral of this factor (Eq.6) is correspondent with output of Fortran program but checking of graphs of this factor (fig.3) showed that [Z.sub.2] value in this graph that started from 0.1 is correct until 0.6 but from 0.6 to end, [Z.sub.2] value has error and every of this value is presented 10 times its actual value.

Examining [I.sub.1] factor in Yang reference showed that the integral of this factor has some errors (Eq.7) but the presented graph in this reference is correct (fig.4).

[I.sub.1] factor in Simons and Senturk reference:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Eq.6

[I.sub.1] factor in Yang reference:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Eq.7

Examining graphs and integrals of I2 factor in Simons and Senturk reference:

For checking [I.sub.2] factor a comparison was also done between mentioned references and written program by FORTRAN. Scrutiny showed that the considered Integrals in both of references have some error (eq.8 and 9), also checking graphs showed that graph of Simons and Senturk reference (fig.5) has error and graph of Yang reference (fig.6) is correct.

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

[I.sub.2] integral in Simons and Senturk reference

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Eq.8

[I.sub.2] integral in Yang reference

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Eq.9

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

Selection of best graph and integral for [I.sub.1] and [I.sub.2] factors:

After checking mentioned references and comparing them with written program by FORTRAN it was clear that integral of Simons and Senturk reference (Eq.6) and graph of Yang reference (Eq.4) offer best answer for [I.sub.1] factor.

For [I.sub.2] factor examines showed that graph of Yang book (fig.6) offer the best answer but none of relative integrals have accordance with FORTRAN program.

Therefore, all of status were first written for integral of [I.sub.2], and then by FORTRAN program different status of integral of [I.sub.2] factor were obtained. Comparing the results of FORTRAN program with graphs data showed that the following equation (Eq.10) has the best answer for [I.sub.2] factor.

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Eq.10

Conclusion:

Solving some of the sediment transport equations like those of Einstein, Chang, Simons Richardson and Toffaleti for estimation of sediment transport amount, needs extracting many of relative Factors on these equations by use of graphs and diagrams of these equations or solving difficult integrals of these equations.

Obtaining different factors by use of graph and diagram just is possible when data are a few is possible and when data are a lot, integrals of these graphs must be solved. Studying different references shows error in graphs and integrals of sediment transport equations and in many cases a sediment transport equation is offered in different forms in various references and it causes error in sediment estimation by use of these equations.

Studying Chang et al. equation in different references showed that integral and graph of [I.sub.1] and [I.sub.2] factors of this equation in various references have some errors. These errors are corrected by use of a program written in FORTRAN language and the best graph and integral for [I.sub.1] and [I.sub.2] factors of Chang et al. equation were obtained.

Results of this study clearly show the possibility of error in different equation for sediment transport estimation. The result of error on these equations leads to error on actual sediment amount.

Based on result of this study, we suggest being careful in using the sediment transport estimation and it will be better to use the main reference of these equations. We also suggest using different computer programs which are capable of solving difficult equations, simulation of diagrams and graphs of sediment transport equations to be sure about the authenticity of various equations.

Notation:

[I.sub.1] and [I.sub.2]: parameters in Chang's transport equation

[U.sub.*]: shear velocity

K: constant

[xi]: Relative depth = y/D

[[epsilon].sub.s] : momentum diffusion coefficient for sediment

[beta] : Constant

D: average flow depth

g: gravitational

d: sediment particle diameter

s: water surface or energy slope; or slope

C: sediment concentration

[C.sub.a]: sediment concentration

[omega]: fall velocity of sediment particle

a: thicknesses of bed layers

u: local velocities

[gamma]: Specific weight of water

k: Von Karman -Prandtl universal constant (= 0.4)

V: average flow velocity

References

[1.] Kazemi, Y., 2008. Estimation of Bedload to Suspend Load ratio in Central Alborz River (Jajrod and Taleghan River), M.Sc. Thesis, Tehran University.

[2.] Shafai Bajestan, M., 2008. Basic Theory and Practice of Hydraulics of Sediment Transport. (Second Edition). Shahid Chamran University.

[3.] Shafai Bajestan, M. and M. Ostad Asgari, 2000. A Mathematical Model to Evaluate the Bed and Total Load by the Modified Einstein Procedure. Journal of Sciences and Technology of Agriculture and Natural Resources, 4(2).

[4.] Simons, D.B. and F. Senturk, 1977. Sediment Transport Technology, 807 pp. Water Recourse. Publ., Highlands Ranch.

[5.] Yang, C.T., 1993. Sediment Transport: Theory and Practice. The McGraw-Hill companies, Inc.

(1) A. Salajegheh, (2) Y. Kazemi, (3) N. Rostami, (4) M.M. Heidari

(1) Associated Professor, Natural Resources Faculty, University of Tehran, Iran

(2) M.Sc. of Watershed Management, University of Tehran, Iran

(3) Ph.D. candidate of Watershed Management sciences and engineering, University of Tehran, Iran

(4) Ph.D. student of Irrigation and Reclamation engineering Department, University of Tehran, Iran

Corresponding Author

N. Rostami, Ph.D. candidate of Watershed Management sciences and engineering, University of Tehran, Iran

Tel: +98 9183405107, E-mail: noredin_rostami@ut.ac.ir

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Original Article |
---|---|

Author: | Salajegheh, A.; Kazemi, Y.; Rostami, N.; Heidari, M.M. |

Publication: | Advances in Environmental Biology |

Article Type: | Report |

Geographic Code: | 7IRAN |

Date: | Feb 1, 2012 |

Words: | 1839 |

Previous Article: | Lead and zinc levels in Scomberomorus guttatus, Scomberomorus commerson and Otolithes ruber from Hendijan, Iran. |

Next Article: | The Effects of varying nitrogen doses on yield and some yield components of black cumin (Nigella Sativa L.). |

Topics: |