Simulation of the relations of the Stage-Discharge, erosion and sedimentation of the Yalfan River in Ekbatan Dam using by FLUENT numerical model
Document Type : Original Article
Authors
1
Department of environment, faculty of Natural resources and environment ,Malayer University ,Malayer,Iran
2
Faculty of Natural Resources and Environment, Malayer University
10.22034/gmpj.2025.543684.1579
Abstract
Introduction
Currently, using discharge-Eschle curves and measuring discharge in natural rivers is generally difficult and error-prone. Therefore, the use of new models plays an effective role in reducing the time and cost associated with hydrometry (Ahmad et al. 2021). Due to the importance of this phenomenon and because the use of physical models requires a large space, high cost, and long time to conduct experiments, many river engineering problems are investigated with mathematical models (Azizi et al., 2019). For this reason, to better understand these phenomena, numerical and computational fluid dynamics models should be used to simulate flow, which are more accurate than other simulation methods (Ildoromi et al. 2020).
Bagheri (2023) evaluated the methods for calculating the discharge-Eschle curve based on bedform resistance relations at the Ablo Nekarud hydrometric station in Mazandaran and showed that the Einstein-Barbarossa method has the lowest absolute error compared to other methods. Mohammadi et al. (2023) showed in estimating the discharge-echel curve in natural rivers with unsteady flow that with a new method based on constant velocity meters, the discharge-echel curve of unsteady flow and depth and discharge hydrographs can be tracked.Kim et al. (2023) evaluated the sedimentation and erosion of riverbanks and channels in the Mekong Delta using the MK-BEHI index and showed that differences in bank stress and riverbank structures cause changes in the erosion rate and riverbank load.Boustani et al. (2023) studied the distribution of bed shear stress using numerical simulation and compared the results of Ansys Fluent and Flow-3D, showing that the numerical results of the Fluent model were closer to the experimental data compared to the Flow 3D model.
Methodology
The Yalfan River is located in the Ekbatan Dam watershed in the northern part of Alvand, 10 kilometers southeast of Hamedan city. In this study, open channel flow in a one-dimensional stable has been investigated.manner in channels with conventional rectangular and trapezoidal cross-sections in the Yalfan River of Ekbatan Dam, Hamedan, Given that the shape of the sides of the channel cross-section has an important effect on the discharge-Ecshle relationships, and the distribution of runoff in channels with different cross-section shapes and surfaces has an important effect on causing erosion and instability of the riverbed.Therefore, numerical values of runoff advance in rectangular and trapezoidal sections were calculated at different times and compared with experimental results for different turbulence models and percentage of error of each turbulence model was calculated and compared with each other.In this research, the second-order discretization method of the Fluent model was used to separate the load and the effect of rectangular and trapezoidal sections to investigate the propagation of runoff flow in scour and edge erosion, and the Reynolds stress numbers were used to determine the movement threshold and to analyze the movement pattern and settling of sediment particles in different flow tubes in a laboratory glass channel, using Fluent software.
Results and Discussion
The results showed that in a rectangular cross section, runoff distribution and pressure are distributed almost equally across the channel surface, and the flow rates do not show significant changes.In a trapezoidal channel with a closed cross-section, discharge - Eschle's measurements showed that the spread of runoff with higher velocity and pressure at the bottom of the channel causes increased erosion at the bottom.The results of the study in a trapezoidal channel with an open cross-section, the values of the discharge - Eschle and the histograms of pressure and diffusion velocity, showed that the pressure and diffusion velocity are higher on the sides of the channel and the erosion rate is higher on the sides.The results of the assessment of the discharge-Ecshle relationships showed that the more closed the trapezoidal channel opening is, the higher the pressure at the bottom and the lower the pressure at the channel wall, and as a result, the greater the erosion rate at the bed.Conversely, the more the trapezoid opening is opened, the greater the runoff pressure on the sides of the bed and, consequently, the greater the erosion rate.The results showed that due to the change in the relationship between the discharge-Ecshle and the flow velocity and pressure at different sections, the rate of erosion and sedimentation changes drastically.The study of runoff distribution showed that the channel shape has an important effect on creating riverbed instability and causes meandering and shifting of the river course over time, which should be considered in river engineering works.
Conclusion
Comparison of the output results of the Fluent hydraulic model shows that this model has good ability to estimate changes in river cross sections.The results show that changes in the channel sections of the Yalfan River along the longitudinal profile of the river cause changes in bank erosion and bank erosion.So that at the beginning of the flow with a trapezoidal cross-section, due to changes in flow propagation speed and pressure, erosion and foaming are greater than at the sides.While in rectangular sections in the middle and end of the stream, lateral erosion and bed sedimentation are observed, which must be carefully considered in the location and design of appropriate river engineering control and protection structures.The results showed that modeling the runoff propagation phenomenon in the Fluent environment has acceptable accuracy, but generating the geometry file and solution field network in Gambit software for study using topographic maps requires a lot of time.Therefore, it is recommended that other software be used to create geometry and modeling for modeling in large watersheds, and that the Fluent model be used for modeling in small watersheds that require high accuracy.
Keywords: Erosion, Cross Sections, Model, Flow Velocity, Fluent
Currently, using discharge-Eschle curves and measuring discharge in natural rivers is generally difficult and error-prone. Therefore, the use of new models plays an effective role in reducing the time and cost associated with hydrometry (Ahmad et al. 2021). Due to the importance of this phenomenon and because the use of physical models requires a large space, high cost, and long time to conduct experiments, many river engineering problems are investigated with mathematical models (Azizi et al., 2019). For this reason, to better understand these phenomena, numerical and computational fluid dynamics models should be used to simulate flow, which are more accurate than other simulation methods (Ildoromi et al. 2020).
Bagheri (2023) evaluated the methods for calculating the discharge-Eschle curve based on bedform resistance relations at the Ablo Nekarud hydrometric station in Mazandaran and showed that the Einstein-Barbarossa method has the lowest absolute error compared to other methods. Mohammadi et al. (2023) showed in estimating the discharge-echel curve in natural rivers with unsteady flow that with a new method based on constant velocity meters, the discharge-echel curve of unsteady flow and depth and discharge hydrographs can be tracked.Kim et al. (2023) evaluated the sedimentation and erosion of riverbanks and channels in the Mekong Delta using the MK-BEHI index and showed that differences in bank stress and riverbank structures cause changes in the erosion rate and riverbank load.Boustani et al. (2023) studied the distribution of bed shear stress using numerical simulation and compared the results of Ansys Fluent and Flow-3D, showing that the numerical results of the Fluent model were closer to the experimental data compared to the Flow 3D model.
Methodology
The Yalfan River is located in the Ekbatan Dam watershed in the northern part of Alvand, 10 kilometers southeast of Hamedan city. In this study, open channel flow in a one-dimensional stable has been investigated.manner in channels with conventional rectangular and trapezoidal cross-sections in the Yalfan River of Ekbatan Dam, Hamedan, Given that the shape of the sides of the channel cross-section has an important effect on the discharge-Ecshle relationships, and the distribution of runoff in channels with different cross-section shapes and surfaces has an important effect on causing erosion and instability of the riverbed.Therefore, numerical values of runoff advance in rectangular and trapezoidal sections were calculated at different times and compared with experimental results for different turbulence models and percentage of error of each turbulence model was calculated and compared with each other.In this research, the second-order discretization method of the Fluent model was used to separate the load and the effect of rectangular and trapezoidal sections to investigate the propagation of runoff flow in scour and edge erosion, and the Reynolds stress numbers were used to determine the movement threshold and to analyze the movement pattern and settling of sediment particles in different flow tubes in a laboratory glass channel, using Fluent software.
Results and Discussion
The results showed that in a rectangular cross section, runoff distribution and pressure are distributed almost equally across the channel surface, and the flow rates do not show significant changes.In a trapezoidal channel with a closed cross-section, discharge - Eschle's measurements showed that the spread of runoff with higher velocity and pressure at the bottom of the channel causes increased erosion at the bottom.The results of the study in a trapezoidal channel with an open cross-section, the values of the discharge - Eschle and the histograms of pressure and diffusion velocity, showed that the pressure and diffusion velocity are higher on the sides of the channel and the erosion rate is higher on the sides.The results of the assessment of the discharge-Ecshle relationships showed that the more closed the trapezoidal channel opening is, the higher the pressure at the bottom and the lower the pressure at the channel wall, and as a result, the greater the erosion rate at the bed.Conversely, the more the trapezoid opening is opened, the greater the runoff pressure on the sides of the bed and, consequently, the greater the erosion rate.The results showed that due to the change in the relationship between the discharge-Ecshle and the flow velocity and pressure at different sections, the rate of erosion and sedimentation changes drastically.The study of runoff distribution showed that the channel shape has an important effect on creating riverbed instability and causes meandering and shifting of the river course over time, which should be considered in river engineering works.
Conclusion
Comparison of the output results of the Fluent hydraulic model shows that this model has good ability to estimate changes in river cross sections.The results show that changes in the channel sections of the Yalfan River along the longitudinal profile of the river cause changes in bank erosion and bank erosion.So that at the beginning of the flow with a trapezoidal cross-section, due to changes in flow propagation speed and pressure, erosion and foaming are greater than at the sides.While in rectangular sections in the middle and end of the stream, lateral erosion and bed sedimentation are observed, which must be carefully considered in the location and design of appropriate river engineering control and protection structures.The results showed that modeling the runoff propagation phenomenon in the Fluent environment has acceptable accuracy, but generating the geometry file and solution field network in Gambit software for study using topographic maps requires a lot of time.Therefore, it is recommended that other software be used to create geometry and modeling for modeling in large watersheds, and that the Fluent model be used for modeling in small watersheds that require high accuracy.
Keywords: Erosion, Cross Sections, Model, Flow Velocity, Fluent
Keywords
ایلدرمی، ع.، شیخی، آ. .، 1395. بررسی تغییرات مورفولوژی رودخانه و نقش آن در فرسایش و رسوبگذاری با استفاده از HEC-RAS (مطالعه موردی: رودخانه خرم آباد-دوآب ویسان)، پژوهشهای ژئومورفولوژی کمی، سال 5، شماره 3، صص 146-163.
ایلدرمی، ع، ر، محمدپناه مقدم، م .،1400. شبیهسازی تأثیر تغییر مورفولوژی محل تلاقی شبکه زهکشی بر الگوی فرسایش و رسوبگذاری رودخانه سیمینه رود همدان با استفاده از مدل عددی فلوئنت، نشریه پژوهشهای ژئومورفولوژی کمی، دوره 9، شماره 4، صص 112- 129.
ایلدرمی، ع، ر.، 1401. شبیهسازی الگوریتم الگوی حرکت ذرات رسوبی رودخانه سد اکباتان با استفاده از مدل عددی فلوئنت، نشریه پژوهشهای ژئومورفولوژی کمی، دوره 11، شماره 2، صص 14-31.
بایزیدی، م.، کرمی، ن.، 1396. پیشبینی روند رسوبگذاری در رودخانه قره سو با استفاده از مدل GSTARS 3، مجله محیط زیست و مهندسی آب، دوره 3، شماره 1، صص 66-80.
باقری، ع .، 1402. ارزیابی روشهای محاسبه منحنی دبی-اشل بر اساس روابط مقاومت فرم بستر در ایستگاه هیدرومتری آبلو نکارود مازندران، نشریه آبیاری و زهکشی ایران، دوره 17، شماره 3، صص 493-503.
بوستانی، آ.، خداشناس، س، ر.، اسماعیلی، ک.، 1402. بررسی توزیع تنش برشی بستر در کانال های مرکب مستقیم با استفاده از شبیهسازی عددی و مقایسه نتایج Ansys Fluent و Flow-3D، نشریه آبیاری و زهکشی ایران، دوره 17، شماره 2.
جباری، ا.، قبادیان، ر.، خان محمدی، ف.، 1401. شبیهسازی تغییرات بستری پیوندگاه رودخانههای گاماسیاب و قره سو در سیلاب نشریه پژوهشهای ژئومورفولوژی کمی، دوره 11، شماره 2، صص 71-87.
جمالی، ع، ر.، امین نژاد، ب.، ظهیری، ع، ر.، 1402. استخراج رابطه دبی-اشل مقاطع مرکب چند سیلاب دشتی با مدلهای یک بعدی و دوبعدی، مجله پژوهشهای حفاظت آب و خاک، دوره 30، شماره 4، صص 103-123.
خانجانی زوراب، م.، ایلدرمی، ع، ر.، نوری ح.، 1399. بررسی انتقال رسوب و تغییرات بستر رودخانه یلفان سداکباتان همدان، نشریه پژوهشهای ژئومورفولوژی کمی، دوره 9، شماره 3، صص 50-69 .
رستمی، م.، امینی، ر.، ظهیری، ع، ر.، 1403. مدل سهبعدی توسعه داده شده برای استخراج رابطه دبی-اشل در مقاطع مرکب چندسیلاب دشتی مستقیم، نشریه هیدرولیک، دوره 19، شماره 2 ، صص 119-138.
عزیزی، ش.، ایلدرمی، ع، ر.، نوری، ح.، 1398. تأثیر تغییر مقاطع عرضی بر طغیان و ظرفیت انتقال رودخانه آبشینه همدان، نشریه پژوهشهای ژئومورفولوژِی کمی، سال8، شماره 2، صص 209-189.
علی نوری، پ.، بی جن خان ، م.، 1402.بررسی میزان تعمیم پذیری نتایج حاصل از مدل های آزمایشگاهی به رودخانههای طبیعی برای تخمین منحنی دبی– اشل، نشریه مهندسی منابع، دوره16، شماره 57، صص31-46.
گودرزی، ا.، اله یونسی، ح.، شاهینژاد، ب.، ترابی، ح.، 1401. مدلسازی عددی انتقال رسوب و فرسایش کنار رودخانهای در شرایط سیلابی مطالعه موردی حوضه آبریز کشکان، نشریه هیدروژئومورفولوژی، دوره 9، شماره 31 ، صص 137-157.
صالحی، م.، فرهادی، ع، ر.، خانی تملیه، ذ.، احمدپور، ظ.، 1397. مدلسازی عددی امواج ناشی از شکست سد سنگی-ملاتی در زیرحوضه نوشان آذربایجان غربی، نشریه محیط زیست و توسعه، دوره 4، شماره 2، صص 90-100.
ظهیری، ع.، قلی نژاد، ج.، دهقانی، ا.، 1398.محاسبه ظرفیت انتقال رسوب در رودخانهها به کمک مدل ریاضی شبه دوبعدی، پژوهشنامه مدیریت حوزه آبخیز، سال 10، شماره 19، صص 142- 153.
محمدی، م.، فغفور مغربی، م.، 1402. تخمین منحنی دبی- اشل در رودخانههای طبیعی با جریان غیردائمی با استفاده از روش مبتنی بر کنتورهای هم سرعت، نشریه هیدرولیک ، دوره 18، شماره 2، صص 1-20.
مظفری، ح.، اصغری سراسکانرود، ص.، اسفندیاری درآباد، ف.، 1401. مدلسازی تغییرات فرسایش و رسوب رودخانه سجاسرود قبل و بعد از ساخت سد گلابر با روش GCD، نشریه پژوهشهای ژئومورفولوژی کمی، دوره 11، شماره 3، صص 57-8
محمدزاده وطن چی، س.، فغفور مغربی، م.، 1401. تحلیل عدم قطعیت منحنیهای دبی- اشل در رودخانهها، نشریه مهندسی عمران و محیط زیست، دوره 52، شماره 106،صص 157- 167.
نبیزاده ولوکلائی، م.، ظهیری، ع، ر.، د دهقانی، ا.، مفتاح هلقی، م.، 1402. واسنجی هیدرولیک جریان و رسوب مدل ریاضی شبهدوبعدی شیونو و نایت در رودخانه زارمرود، نشریه آب و آبیاری ایران، دوره 13، شماره 4، صص 66-84.
Ahmad, D., Afzal, M., and Ishaq, M. (2024). Impacts of riverbank erosion and flooding on communities along the Indus River, Pakistan. Natural Hazards, 120(1), 131-152. DOI: 10.1007/s11069-023-06216-2
Ahmadi, A., and Maghrebi, M.F. (2021). A robust approach for rating curves estimation in open channels using isovel contours. International Journal of River Basin Management, 19(3), 281-296.
FLUENT Inc. (2006). FLUENT 6.3. User s and Tutorial Guide Manual, Lebanon, USA.
Ghosh, D., and Saha, S. (2024). Identifying river bank erosion potential zones through geo-spatial and binary logistic regression modeling approach: a case study of river Ganga in Malda district (India). Modeling Earth Systems and Environment, 10(1), 81-98 .doi.org/10.1007/s40808-023-01740-3
Huang, P. C. (2024). Estimation of riverbank erosion by combining channel morphological models with AI techniques. Geomatics, Natural Hazards and Risk, 15(1), 2359983. doi.org/10.1080/19475705.2024.2359983
Kim, T. T., Ngoc, P., Nga, T. N. Q., Nguyet, N. T. T., Truong, H. N., Diem, P. T. M., and Bay, N. T. (2023). Modifying BEHI (Bank Erosion Hazard Index) to map and assess the levels of potential riverbank erosion of highly human impacted rivers: a case study for Vietnamese Mekong River system. Environmental Earth Sciences, 82(23), 554. DOI:10.1007/s12665-023-11249-8
Lai, Y. G., and Wu, K. (2019). A Three-Dimensional Flow and Sediment Transport Model for Free-Surface Open Channel Flows on Unstructured Flexible Meshes. Fluids, 4(1), 18. /doi.org/10.3390/fluids4010018
Salawu, O. G., and Leke, D. (2024). Assessment of Changes in Channel Width and Braiding: The Implications on Flooding Upstream of Shiroro Reservoir in river Kaduna, Kaduna State, Nigeria. Journal of Applied Sciences and Environmental Management, 28(1), 179-186. DOI:10.4314/jasem.v28i1.20
Zhao, Q., Misra, S., Svendsen, I., and Kirby, j. (2004). Numerical study of turbulent hydraulic jump. ASCE Engineering Mechanics Conference, University of Delaware Newmark.
