سرعت جریان در کانال مرکب پیچان تحت تأثیر ضریب خمیدگی

نوع مقاله : پژوهشی اصیل (کامل)

نویسندگان
1 دانشجوی دکتری مهندسی عمران- آب و سازه های هیدرولیکی، دانشکده فنی، دانشگاه ارومیه
2 دانشیار مهندسی عمران- هیدرولیک و مکانیک مهندسی رودخانه، دانشکده فنی، دانشگاه ارومیه.
3 استادیار گروه مهندسی آب، دانشکده کشاورزی، دانشگاه زنجان.
چکیده
رودخانه­های طبیعی برای تنظیم شیب خط انرژی، به­ندرت در مسیری مستقیم جریان دارند و معمولاً دارای مسیری انحنادار هستند که به آن­ها پیچان­رود گفته می­شود. پس از پیدایش رودخانه­های پیچان، با گذشت زمان و حرکت جانبی پیچ­ها، میزان پیشروی قوس برونی افزایش یافته و مقدار ضریب خمیدگی یا سینوسی در پلان زیاد می شود. درکانال­های پیچان، میزان انحنا با عدد بی بعدی به نام ضریب خمیدگی تعریف می­شود که نسبت طول مسیر انحنادارکانال اصلی به طول مستقیم سیلاب دشت می­باشد. با افزایش مقدار ضریب خمیدگی، شیب، سرعت جریان و ظرفیت انتقال دبی رودخانه کاهش می یابد. در نتیجه خطر سیلابی شدن به طور قابل توجهی افزایش یافته و در هنگام بروز سیلاب، سطح آب از مقطع اصلی رودخانه فراتر رفته و وارد دشت های سیلابی آن می شود. در این حالت با توجه به اندرکنش بین جریان با سرعت بالا در کانال اصلی و جریان کندتر در سیلاب­دشت­ها و تبادل اندازه حرکت بین این دو ناحیه، مشخصات جریان دائماً در حال تغییر می­کند. در این تحقیق، مشخصات هیدرولیکی جریان شامل توزیع عرضی سرعت متوسط عمقی، پروفیل سطح آب، توزیع سرعت در طول، نسبت سرعت عرضی به طولی در محور مرکزی کانال اصلی، میزان سرعت متوسط و دبی عبوری از کانال اصلی، بردارهای سرعت در مقاطع مختلف عرضی و توزیع برآیند سرعت در سطوح مختلف افقی با توجه به تغییر ضریب خمیدگی برای شش نوع کانال با مقادیر سینوسی مختلف به صورت عددی مورد بررسی قرارگرفته است. برای بررسی اثر ضریب خمیدگی کانال مرکب پیچان روی مشخصات هیدرولیکی جریان، از نرم­افزار FLOW3D استفاده شده است و مدل آشفتگی این نرم­افزار نیز طوری انتخاب گردیدکه انطباق بهتری با داده های آزمایشگاهی داشته باشد. برای این منظور، از دو مدل آشفتگی پرکاربرد RNG و k-ε استفاده شد و عملکرد این دو مدل در شبیه­سازی مشخصه هیدرولیکی جریان مورد بررسی قرار گرفت و مشخص گردید که مدل آشفتگی RNG دارای دقت بیشتری می­باشد. در ادامه، این مدل به عنوان مدل آشفتگی نهایی برای شبیه سازی عددی معرفی گردید. در نهایت با بررسی کانال مرکب پیچان تحت تأثیر ضریب خمیدگی مختلف مشخص شد که با افزایش ضریب خمیدگی کانال از مقدار 1 به 1.641، مقدار سرعت متوسط مقطع کانال اصلی به طور میانگین 54% و میزان دبی عبوری از کانال اصلی به طور متوسط 38% کاهش می یابد. با افزایش ضریب خمیدگی، مقدار حداکثر سرعت متوسط عمقی از 0.55 به 0.38 متر برثانیه کاهش و مقدار حداکثر ارتفاع سطح آب از 0.305 به 0.332 متر در قوس بیرونی مقطع CS1 افزایش می یابد. افزایش مقدار ضریب خمیدگی باعث می شود نسبت سرعت عرضی به سرعت طولی در محور مرکزی کانال اصلی افزایش یافته به­طوری­که مقدار آن در بحرانی­ترین حالت از صفر به 0.4 می رسد. با افزایش ضریب خمیدگی، مقدار سرعت طولی حداکثر به سمت سیلاب­دشت سمت راست (خم داخلی) حرکت کرده و از مقدار آن کاسته می شود. به طوری­که با افزایش ضریب خمیدگی از 1 به 1.641، مقدار حداکثر سرعت طولی از 0.55 به 0.42 متر برثانیه کاهش یافته و موقعیت آن از مرکز کانال اصلی به سمت خم داخلی در بالای عمق لبریز کانال اصلی حرکت می­کند.

کلیدواژه‌ها

موضوعات


عنوان مقاله English

Flow Velocity in Meandering Compound Channel under the Influence of Sinusoidal Change

نویسندگان English

mohammad naghavi 1
mir ali mohammadi 2
Ghorban Mahtabi 3
1 Ph.D. Candidate in Water & Hydraulic Structures, Department of Civil Engineering, Faculty of Eng., Urmia University
2 Associate Professor in Water & Hydraulic Structures, Department of Civil Engineering, Faculty of Eng.
3 Assistant Professor in Water Engineering, Department of Water Engineering, Faculty of Agric., University of Zanjan, Zanjan, Iran.
چکیده English

Natural rivers are rarely in direct flow because of regulating the energy grade-line, and usually have a curved path to which it is referred to as "meandering channels". After the appearance of meandering rivers, with the passage of time and the lateral movement of the meanders, the external bending progression and the sinusoidal or curvature is increased. In the meandering channels, the curvature of the meandering sections with a dimensionless number can be defined as the sinusoidal which is the ratio of meander length of the main channel to the floodplain length. By increasing sinusoidal slope number, flow velocity and river discharge capacity decrease. As a result, the risk of flood has increased significantly and during floods the water level exceeds to the main river boundary and enters to the floodplains. In this case, due to the interaction between higher velocities in the main channel and the slower velocities in the floodplains and the exchange of the momentum between these two regions, the flow profile is constantly changing. In this research, the hydraulic characteristics of the flow including the depth-averaged velocity, the free surface water profile, longitudinal velocity distributions, ratios of transverse to longitudinal velocities in the central axis of the main channel and the mean velocity and flow rate of the main channel along the meandering compound channel have been investigated numerically, regarding the change in the sinusoidal ratio for Six types of channels with different sinusoidal ratios. In order to investigate the effect of the sinusoidal ratio in meandering compound channel on the hydraulic characteristics of the flow, the FLOW3D software is used, the software was also chosen so that the turbulence model with experimental data have better compliance. For this purpose, two RNG and k-ε turbulence models were used and the performance of these two models was investigated in simulating the important hydraulic characteristics of the flow, such as the flow velocity, and it was determined that the RNG turbulence model has a better accuracy than the k- ε turbulence model. In the following, this model was introduced as the final turbulence model for numerical simulation. Numerical simulation results show that by increasing the sinusoidal ratio of channel from 1 to 1.641, the mean velocity of the main channel section is decreased by 54% on average and the flow rate of the main channel decreases by the average of 38%. Also, by increasing the sinusoidal ratio, the maximum depth-averaged velocity decreases from 0.55 m/s to 0.38 m/s, and the maximum free surface height of the water rises from 0.305 m to 0.332 m in the outer bend of the CS1 cross section. Increasing the sinusoidal ratio causes the ratio of the transverse velocity to be increased longitudinally in the central axis of the main channel, so that its value in the most critical state reaches from zero to 0.4. As the sinusoidal ratio increases, the maximum length velocity moves towards the right side flood plain (internal bend) and decreases its value so that by increasing the sinusoidal ratio from 1 to 1.641, the maximum longitudinal velocity 0.55 m/s to 0.42 m/s and its position moves from the center of the main channel to the inner bend over the depth of the main channel overflow.

کلیدواژه‌ها English

Compound Channel
Meandering
Sinusoidal
Depth-Averaged Velocity
Flow3D
[1] Toebes G.H. & Sooky A. A. 1967 Hydraulics of meandering rivers with flood plains. A. S. C. E. J. of Waterways and Harbors, 33, 213–236.
[2] Knight D.W. & Demetriou J.D.1983 Flood plain and main channel flow interaction. Journal of Hydraulic Engineering,109 (8), 1073–1092.
[3] Kiely G.1990 Overbank Flow in Meandering Channels the Important Mechanisms. International Conference on River Flood Hydraulics, 207–217.
[4] Ervine D.A., Willetts B.B., Sellin R.H.J. & Lorena M.1993 Factors Affecting Conveyance in Meandering Compound Flows. Journal of Hydraulic Engineering, ASCE, 119(12), 1383-1399.
[5] Shiono K. & Muto Y.1998 Complex flow mechanisms in compound meandering channels with overbank flow. Journal of Fluid Mechanics, 376, 221–261.
[6] Shiono K., Muto Y., Knight D. W. & Hyde, A. F. L. 1999 Energy losses due to secondary flow and turbulence in meandering channels with over bank flow. Journal of Hydraulic Research, 37(5), 641–664.
[7] Shiono K., Al-Romaih J.S. & Knight D.W. 1999 Stage-discharge assessment in compound meandering channels. Journal of Hydraulic Engineering, 125 (1), 66–77.
[8] Spooner J. & Shiono. K. 2003 Compound Meandering Channels with overbank flow. Water and Maritime Engineering,156(3), 225-233.
[9] Wormleaton P.R., Sellin R.H.J., Bryant T., Loveless J.H., Hey R.D. & Catmur S.E. 2004 Flow structures in a two-stage channel with a mobile bed. Journal of Hydraulic Research, 42(2), 145–162.
[10] De Marchis M. & Napoli E. 2008 The effect of geometrical parameters on the discharge capacity of meandering compound channels. Advances in Water Resources, 31,1662–1673.
[11] Rameshwaran P., Sun X., Shiono K., Chandler J.H. & Sellin, R. H. J. 2008 The Modelling of Compound Channel Flow: Physical Model of River Blackwater, Wallingford. UK.
[12] Shiono K., Spooner J., Chan TL., Rameshwaran P. & Chandler J.H. 2008 Flow characteristics in meandering channels with non-mobile and mobile beds for overbank flows. Journal of Hydraulic Research, 46(1), 113–32.
[13] Liu C., Wright N., Liu X. & Yang K. 2014 An analytical model for lateral depth-averaged velocity distributions along a meander in curved compound channels. Advances in Water Resources, 74, 26–43.
[14] Mera I., Francat M.J., Anta J. & Peña E. 2014 Turbulence anisotropy in a compound meandering channel with different submergence conditions. Advances in Water Resources.
[15] Liu C., Shan Y., Liu X. & Yang K. 2016 Method for assessing stage-discharge in meandering compound channels. ICE-Water Manage,169(1), 17–29.
[16] Liu C., Shan Y., Liu X., Yang K. & Liu H. 2016 The effect of floodplain grass on the flow characteristics of meandering compound channels. Journal of Hydrology, 542, 1-17.
[17] Shan Y., Liu X., Yang K. & Liu C. 2017 Analytical model for stage-discharge estimation in meandering compound channels with submerged flexible vegetation. Advances in Water Resources, 108, 170–183.
[18] Saeed-Golanik S. & Mohammadi, M. 2017 A Numerical Consideration of Hydraulic Aspedts of Flow in Meandering Prismatic Compound Channels by Flow3D Model. Proceedings, 5th International Congress on: Civil Eng., Architecture and Urban Development; 26-28 December, Shahid Beheshti University, Tehran, Iran. (in Farsi).
[19] Shan Y., Huang S., Liu C., Guo Y. & Yang K. 2018 Prediction of the depth-averaged two-dimensional flow direction along a meander in compound channels. Journal of Hydrology, 565, 318–330.
[20] Moncho-Esteve I., García-Villalba M., Muto Y., Shiono K. & Palau-Salvador G. 2018 A numerical study of the complex flow structure in a compound meandering channel. Advances in Water Resources, 116, 95–116.
[21]Muste M. & Ettema R. 2004 Scale effects in flume experiments on flow around a spur dike in flatbed channel. Journal of Hydraulic Eng, 130, 635- 646.
[22]González A.E. & Bombardelli F.A. 2005 Two phase flow theoretical and numerical models for hydraulic jumps including air entrainment. In Proc. XXXI IAHR Congress. Seoul. Korea.
[23]Acharya A. 2001 Experimental study and numerical simulation of flow and sediment transport around a series of spur dikes. PhD Thesis.In the Graduate College.the university of Arizona.