تحلیل عددی هیدرولیک جریان در سیستم کالورت-سرریز تاج‌دایره‌ای همراه با دریچه

حیدری, دنیا; ایزدی نیا, الهام

doi:https://doi.org/10.30486/TSWS.2024.1191046

رقم المقالة : 140308291191046 زيارة : 65 الصفحة: 114 - 126

https://doi.org/10.30486/TSWS.2024.1191046

نوع المخطوط: ابحاث

تحلیل عددی هیدرولیک جریان در سیستم کالورت-سرریز تاج‌دایره‌ای همراه با دریچه

الموضوعات :

دنیا حیدری ¹ , الهام ایزدی نیا ²

1 - دانشجوی کارشناسی ارشد،گروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه علم و صنعت ایران، تهران، ایران.
2 - استادیار گروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه شهید اشرفی اصفهانی، اصفهان، ایران.

تاريخ الإرسال : 13 الجمعة , جمادى الأولى, 1446 تاريخ التأكيد : 08 الإثنين , جمادى الثانية, 1446 تاريخ الإصدار : 23 الثلاثاء , جمادى الثانية, 1446

الکلمات المفتاحية: سرریز دایره‌ای, دینامیک سیالات محاسباتی(CFD), کالورت- سرریز ترکیبی, هیدرولیک جریان, شبیه‌سازی عددی,

ملخص المقالة :

سرریزها نقشی حیاتی در پروژه‌های مهندسی آب دارند و به مدیریت جریان‌های سطحی کمک می‌کنند. یکی از انواع پرکاربرد این سازه‌ها، سرریز با تاج‌دایره‌ای است که به دلیل طراحی ساده و توانایی بالا در عبور جریان، در بسیاری از پروژه‌ها مورد استفاده قرار می‌گیرد. در این پژوهش، هدف اصلی بررسی ویژگی‌های هیدرولیکی جریان در سیستم کالورت- سرریز با تاج‌دایره‌ای و تحلیل تأثیر پارامترهای کلیدی مانند شعاع تاج و بازشدگی دریچه بر ضریب دبی است. به این منظور، پروفیل‌های سرعت، سطح آب، و سایر پارامترهای جریان به عنوان متغیرهای اصلی تحلیل شده‌اند. مدل‌های مختلف سیستم کالورت- سرریز با تاج‌دایره‌ای در نرم‌افزار گمبیت طراحی و شبکه‌بندی شده و سپس به نرم‌افزار فلوئنت منتقل شده‌اند تا شبیه‌سازی‌های دقیق دینامیک سیالات محاسباتی (CFD) صورت گیرد. در این مطالعه، ضریب دبی تاج سرریز، دریچه، و کل سازه سرریز برای مقادیر مختلف بار آبی و پارامترهای متغیر محاسبه و مقایسه شده است. نتایج این پژوهش نشان می‌دهد که ضریب دبی تاج سرریز در حضور دریچه بیشتر از حالت بدون دریچه بوده و با کاهش شعاع تاج و افزایش بازشدگی دریچه، ضریب دبی افزایش می‌یابد. همچنین، مقایسه نتایج شبیه‌سازی‌های عددی با داده‌های آزمایشگاهی معتبر نشان‌دهنده مطابقت قابل‌قبول و دقت بالای مدل‌سازی عددی است. این تحقیق تلاش دارد تا با ارائه داده‌های دقیق و تحلیل‌های کاربردی، راهکارهایی برای بهبود طراحی و عملکرد سازه‌های آبی ترکیبی ارائه دهد

المصادر:

Bagheri, S., & Heydarpour, M. (2010). Application of Euler equations in determining the discharge coefficient of circular crest spillways. Iranian Journal of Water and Soil Research, 40(1)67-73. (In Persian)
Beigi Poor, G., & Bina, M. (2005). Hydraulic flow in circular and cylindrical crested weirs. 5th Iranian Hydraulic Conference, Shahid Bahonar University of Kerman, Kerman, Iran, 1-8.
Bhajantri, M.R., Eldho, T.I., & Deolalikar, P.B. (2006). Hydrodynamic modelling of flow over a spillway using a two-dimensional finite volume-based numerical model. Sadhana, 31(6), 743-754. https://doi.org/10.1007/BF02716893
Ghobadian, R., Farmani Fard, M., & Karami, A. (2012). Effect of crest roughness and body on discharge coefficient and hydraulic conditions of flow in circular crested weirs using a physical model. Journal of Soil and Water Scienc-es, 22(1), 1-8.
Izadinia, E., & Heidari, D. (2023). Flow simulation in a combined culvert-circular crown spillway structure, 22nd Iranian Hydraulic Conference, Maragha. https://civilica.com/doc/1963858
Izadinia, E., & Heidari, D. (2022). Flow simulation in a combined structure of culvert and circular crown spillway, 21st National Conference of Iranian Hydraulics, Ahvaz. https://civilica.com/doc/1694896
Negm, A.M., Al-Brahim, A.M., & Alhamid, A.A. (2002). Combined-free flow over weirs and below gates. Journal of Hydraulic Research, 40(3), 359-365. https://doi.org/10.1080/00221680209499950
Ramamurthy, A.S., & Vo, N.D. (1993). Characteristics of circular crested weir. Journal of Hydraulic Engineering, 119(9), 1055–1062.
Samani, J., & Mazaheri, M. (2009). Combined flow over weir and under gate. Journal of Hydraulic Engineering, 135(3), 224–227. https://doi.org/10.1061/(ASCE)0733-9429(2009)135:3(224)
Schmocker, L., Halldórsdóttir, B., & Hager, W. (2011). Effect of weir face angles on circular-crested weir flow. Jour-nal of Hydraulic Engineering, 137(6), 637–643. https://doi.org/10.1061/(ASCE)HY.1943-7900.0000346
Severi, A., Masoudian, M., Kordi, E., & Roettcher, K. (2014). Discharge coefficient of combined-free over-under flow on a cylindrical weir-gate. ISH Journal of Hydraulic Engineering. 42-52. https://doi.org/10.1080/09715010.2014.939503

نص كامل:

Technical Strategies in Water Systems

https://sanad.iau.ir/journal/tsws

ISSN (Online): 2981-1449

Summer 2024: Vol 2, Issue 2, 114-126

https://doi.org/10.30486/TSWS.2024.1191046

$C:\Users\s_mohseni\Desktop\vol 1-3\Backup_of_iraji logo rahbord.jpg$

Research Article

Numerical analysis of flow hydraulics in a culvert-weir system with a circular crest and gate

Donya Heidari 1, Elham Izadinia2*

1 Master's Student, Department of Civil Engineering, Faculty of Engineering and Technology, Iran University of Science and Technology, Tehran, Iran.

2 Assistant Professor, Department of Civil Engineering, Faculty of Engineering and Technology, Shahid Ashrafi Esfahani University, Isfahan, Iran.

Corresponding Author email: elham.izadinia@gmail.com

Received: 16 Oct 2024

Accepted: 09 Dec 2024

Published: 24 Dec 2024

Abstract

Weirs are essential components in water engineering projects, effectively managing surface flows. Among these structures, circular-crested weirs are particularly popular due to their simple design and high flow capacity. This study primarily focuses on investigating the hydraulic characteristics of flow in a culvert-weir system with a circular crest and analyzing the effects of critical parameters such as crest radius and gate opening on the discharge coefficient. For this purpose, velocity profiles, water surface levels, and other flow parameters were analyzed as the main variables. Various models of the culvert-weir system with circular crests were designed and meshed using Gambit software, then imported to Fluent software for detailed computational fluid dynamics (CFD) simulations. The discharge coefficients of the weir crest, gate, and the entire structure were calculated and compared under different water load conditions and independent parameters. The results indicate that the discharge coefficient of the weir crest is higher when a gate is present compared to a gate-free condition. Furthermore, the discharge coefficient increases as the crest radius decreases and the gate opening expands. Additionally, the comparison between numerical simulation results and reliable experimental data shows a high level of accuracy and satisfactory agreement, validating the numerical modeling approach. This study provides precise data and practical analyses to improve the design and performance of combined hydraulic structures.

Keywords: Circular Weir, Computational Fluid Dynamics (CFD), Culvert-Weir Combination, Flow Hydraulics, Numerical Simulation

مقاله پژوهشی

تحلیل عددی هیدرولیک جریان در سیستم کالورت-سرریز تاج‌دایره‌ای همراه با دریچه

دنیا حیدری1، الهام ایزدینیا2*

1. دانشجوی کارشناسی ارشد،گروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه علم و صنعت ایران، تهران، ایران.

2. استادیار گروه مهندسی عمران، دانشکده فنی و مهندسی، دانشگاه شهید اشرفی اصفهانی، اصفهان، ایران.

ایمیل نویسنده مسئول: elham.izadinia@gmail.com

چاپ: 04/10/1403

پذیرش: 19/09/1403

دریافت: 25/7/1403

چکیده

واژه‌های کلیدی: سرریز دایره‌ای، دینامیک سیالات محاسباتی(CFD)، کالورت- سرریز ترکیبی، هیدرولیک جریان، شبیه‌سازی عددی

1- مقدمه

سرریزهای تاج‌دایره‌ای، که به عنوان یکی از انواع سرریزهای لبه‌باریک شناخته می‌شوند، از جمله سازه‌های هیدرولیکی پرکاربرد در پروژه‌های مهندسی آب به شمار می‌روند. این نوع سرریزها در اواخر قرن نوزدهم و اوایل قرن بیستم، قبل از معرفی سرریزهای اوجی شکل، به طور وسیعی استفاده می‌شدند و به دلیل ویژگی‌های خاص خود در بسیاری از کاربردهای مدرن حضور دارند. در طول قرن نوزدهم، با افزایش تقاضا برای ظرفیت‌های بالاتر در تخلیه آب از سازه‌ها، طراحی سرریزهای تاج‌دایره‌ای مورد توجه قرار گرفت و این تلاش‌ها به شکوفایی این نوع سرریزها منجر شد(Ghobadian et al. 2012). سرریزهای تاج‌دایره‌ای و استوانه‌ای در مقایسه با دیگر سازه‌های اندازه‌گیری و کنترل جریان، از مزایای قابل توجهی برخوردار هستند. این مزایا شامل ایجاد الگوی پایدار جریان ریزشی، عبور آسان اجسام شناور بدون اختلال در جریان، طراحی ساده و قابلیت انعطاف، ضریب جریان بالا و پایداری جریان است. از دیگر ویژگی‌های برجسته این سرریزها می‌توان به قابلیت طراحی در مقیاس‌های بزرگ و کاهش هزینه‌های ساخت و بهره‌برداری اشاره کرد. به همین دلیل، سرریزهای تاج‌دایره‌ای به‌طور گسترده‌ای در مهندسی هیدرولیک به عنوان ابزارهایی برای اندازه‌گیری جریان و سازه‌های تخلیه آب در پروژه‌های مختلف به کار می‌روند. این سرریزها به‌ویژه در کنترل سطح آب در کانال‌ها، مخازن و همچنین در سیستم‌های آبیاری و مدیریت منابع آب بسیار موثر هستندBeigi Poor & Bina, 2005)). به طور خلاصه، سرریز تاج‌دایره‌ای شامل تاجی است که به شکل یک قطاع از دایره با شعاع طراحی شده و به صورت عمود بر جریان قرار گرفته است، به‌طوری‌که دیواره‌های بالادست و پایین‌دست آن با افق زاویه‌های α و β تشکیل می‌دهند. رابطه ارتفاع - دبی در واحد عرض جریان (q) برای این سرریزهای لبه‌باریک به صورت زیر نوشته می‌شود. در اکثر تحقیقات مرتبط با سرریزها، ضریب دبی (Cd) از رابطه (1) محاسبه می‌شود.

در رابطه (1)، q دبی در واحد عرض در بالادست سرریز، g شتاب ثقل، بار آبی کل نسبت به تاج سرریز و ضریب دبی سرریز است. شکل (1)، جریان عبوری از روی سرریز را نشان می‌دهد که در آن قطاعی از دایره به شعاع دیواره‌های بالادست و پایین‌دست آن است که با افق زوایای و می‌سازد و به صورت عمود بر جهت جریان نصب می‌شود. بار آبی کل نسبت به تاج سرریز و D بازشدگی دریچه است.

شکل 1- مقطع یک سرریز تاج‌دایره‌ای

Fig 1. Cross-section of a circular-crested weir

یکی از روش‌های مهم و رایج برای تحلیل رفتار سیال در سرریزها، استفاده از مدل‌های فیزیکی است. این مدل‌ها می‌توانند در صورتی که شرایط مدل‌سازی به‌دقت شبیه‌سازی شوند، نتایج دقیقی ارائه دهند که در بسیاری از موارد قابل استناد هستند. با این حال، یکی از چالش‌های اصلی استفاده از مدل‌های فیزیکی در پروژه‌های مهندسی، زمان طولانی و هزینه‌های بالا است. به‌عنوان‌ مثال، برای بررسی پدیده‌هایی پیچیده مانند کاویتاسیون، ممکن است مدت زمان زیادی نیاز باشد و حتی در برخی از مسائل، مدل‌سازی شرایط محیطی با دماهای بسیار بالا، مانند 4000 درجه سانتی‌گراد، ممکن است دشوار یا حتی غیرممکن باشد. این موضوع باعث می‌شود که استفاده از روش‌های عددی به‌عنوان یک راه‌حل کارآمد، امکان‌پذیر شود. روش‌های عددی، به‌ویژه روش‌های دینامیک سیالات محاسباتی(CFD)، این امکان را فراهم می‌آورند که در مدت زمان کوتاه‌تری به تحلیل سیستم‌های پیچیده پرداخته شود. این روش‌ها به محققان اجازه می‌دهند که گزینه‌های مختلف را بررسی کرده و آن‌ها را با هزینه‌های بسیار کمتر نسبت به آزمایش‌های فیزیکی ارزیابی کنند. علاوه‌بر سرعت بیشتر، روش‌های عددی می‌توانند اطلاعات دقیق‌تری از جمله تغییرات سرعت، فشار و دما را در سراسر محدوده موردنظر ارائه دهند. این در حالی است که جمع‌آوری چنین داده‌هایی از طریق آزمایش‌های فیزیکی معمولاً پیچیده، زمان‌بر و گاهی غیرعملی است.

یکی از مزایای برجسته استفاده از روش‌های CFD نسبت به روش‌های تجربی عبارت است از:

کاهش چشم‌گیر در زمان و هزینه طراحی‌های جدید: با استفاده از شبیه‌سازی‌های عددی، می‌توان طراحی‌ها و آزمایش‌های مختلف را بدون نیاز به ساخت مدل‌های فیزیکی و انجام آزمایش‌های پیچیده و پرهزینه انجام داد.

امکان بررسی سیستم‌هایی که آزمایش بر روی آن‌ها دشوار یا غیرممکن است: برای مثال، تحلیل سیستم‌های بسیار بزرگ یا سیستم‌هایی که در شرایط خاصی فعالیت می‌کنند، با استفاده از شبیه‌سازی‌های عددی امکان‌پذیر است.

قابلیت مطالعه سیستم‌ها در شرایط غیرعادی یا تصادفی: به کمک روش‌های CFD می‌توان به شبیه‌سازی شرایط غیرمنتظره و یا حالت‌های ناپایدار که در آزمایش‌های فیزیکی به‌راحتی قابل بررسی نیستند، پرداخته و رفتار آن‌ها را تحلیل کرد.

جزئیات بسیار بالای نتایج: روش‌های عددی به‌ویژه به محققان این امکان را می‌دهند که در سطح بسیار دقیقی از نتایج، از جمله پروفایل‌های سرعت، فشار و دما، تجزیه‌وتحلیل انجام دهند.

تحقیقات زیادی تاکنون روی سرریزهای تاج‌دایره‌ای انجام شده است. (Ramamurthy & Vo, 1993) به بررسی تأثیر زاویه وجه بالادست و پایین‌دست و شعاع تاج سرریز بر ضریب دبی پرداختند. (Negm et al., 2002) ضریب دبی و ویژگی‌های هیدرولیکی سرریز لبه‌تیز دریچه‌دار با فشردگی جانبی را بررسی کردند. (Bhajantri et al., 2006) مدل دوبعدی از سرریز تاج‌دایره‌ای را مورد مطالعه قرار دادند و توزیع سرعت، فشار و ضریب دبی را با نتایج مدل‌های فیزیکی مقایسه کردند که تطابق خوبی داشت. (Samani, & Mazaheri, 2009) جریان مستغرق و نیمه‌مستغرق را برای سرریز لبه‌تیز دریچه‌دار بدون انقباض جانبی بررسی کردند و نتایج دقیق مدل‌های آزمایشگاهی را برای تعیین عمق آب و رابطه دبی تاج و دریچه مشاهده کردند. (Schmocker et al., 2011) تأثیر زاویه وجه بالا و پایین‌دست، جریان آزاد و مستغرق، و شعاع تاج سرریز بر ضریب دبی و ویژگی‌های هیدرولیکی را با استفاده از 20 مدل آزمایشگاهی بررسی کردند. (Severi et al., 2014) نیز تغییرات ضریب دبی سرریز استوانه‌ای دریچه‌دار در شرایط مختلف جریان آزاد، ترکیبی، و زیرگذری را مطالعه کردند.

کالورت یا آبرو زیرزمینی سازه متداولی است که برای عبور جریان از زیر جاده، بزرگراه، خاک‌ریز، مناطق مسکونی و سرریز استفاده می‌شود. کالورت‌ها از لحاظ شکل مقطع به پنج دسته مستطیلی، مربعی، دایره‌ای، بیضی و نیم‌بیضی تقسیم‌بندی می‌شوند. انتخاب شکل کالورت به عواملی مانند نوع و محل کاربرد و محدودیت اجرا بستگی دارد. متداول‌ترین نوع کالورت مستطیلی است. طراحی کالورت متاًثر از عوامل مختلفی از جمله عملکرد هیدرولیکی کالورت، استحکام بنا، هزینه‌های کلی سازه و نگهداری و سایر موارد است؛ به همین منظور همیشه طراحی یک کالورت ایده‌آل مد نظر است. هدف این تحقیق بررسی شرایط هیدرولیکی جریان با استفاده از مدل‌سازی عددی است.

2- مواد و روش‌ها

2-1- شبیه‌سازي جريان به روش CFD

نرم‌افزار مورد استفاده از یک الگوریتم کلی برای شبیه‌سازی پیروی می‌کند که شامل مراحل زیر است:

تعیین هندسه‌ی جریان: در این مرحله، سیال و مسیر حرکت آن مدل‌سازی می‌شود. این مدل‌سازی می‌تواند در محیط‌هایی مانند ANSYS Design Modeler یا گمبیت انجام شود.

تولید شبکه یا المان: برای حل عددی معادلات سیال از تکنیک‌های ریاضی پیشرفته استفاده می‌شود. نرم‌افزار فلوئنت از روش حجم محدود یا کنترل حجم بهره می‌برد. تولید شبکه به دقت مسأله و نوع هندسه جریان وابسته است و نقش کلیدی در حل دقیق مسأله ایفا می‌کند. این مرحله می‌تواند در نرم‌افزارهایی مانند ANSYS Icem CFD یا گمبیت انجام شود. خروجی این نرم‌افزارها با فرمت msh به عنوان ورودی در فلوئنت استفاده می‌شود.

اعمال شرایط مرزی: در این مرحله مدل جریان (آرام یا آشفته)، ویژگی‌های فیزیکی سیال، فازهای جریان، شرایط مرزی روی سطوح و تنظیمات مربوط به مش دینامیکی مشخص می‌شود.

حل مسأله: تنظیمات مربوط به فرآیند حل در این مرحله انجام می‌شود.

پردازش نتایج: در این مرحله، تمام نتایج محاسباتی و اهداف موردنظر قابل بررسی هستند. این شامل مشاهده‌ی کانتورهای مختلف، خطوط جریان، بردارهای نرمال، استخراج نمودارها و ساخت انیمیشن از مدل جریان است. پردازش نتایج می‌تواند در نرم‌افزار فلوئنت یا محیط‌های ANSYS Workbench و CFD Post انجام شود.

در تحقیق حاضر، برای تعیین سطح آزاد جریان و محاسبه عمق آب بالادست، از مدل به دلیل سرعت همگرایی و دقت بالاتر استفاده شد. همچنین، برای برداشت پروفیل‌های سرعت، از مدل Reynolds Stress بهره گرفته شد. مقاطع پروفیل سرعت در طول کانال و سرریز با استفاده از قابلیت Iso Surface نرم‌افزار تعیین شدند و فاصله‌ی این مقاطع با نزدیک شدن به سرریز، به دلیل افزایش تغییرات، کاهش یافت. اندازه‌گیری دبی جریان عبوری از روی تاج و کالورت سازه با استفاده از پروفیل سرعت و انتگرال‌گیری انجام شد(Izadinia & Heidari., 2023).

مدل یکی از پرکاربردترین و شناخته‌شده‌ترین مدل‌های آشفتگی در دینامیک سیالات محاسباتی (CFD) است که برای شبیه‌سازی جریان‌های آشفته به کار می‌رود. این مدل بر اساس حل دو معادله برای توصیف رفتار آشفتگی در جریان طراحی شده است:

K: نشان‌دهنده انرژی جنبشی آشفتگی است که معیاری برای اندازه‌گیری شدت آشفتگی در جریان است.

:(ε)نرخ اتلاف انرژی جنبشی آشفتگی است که نرخ تبدیل انرژی جنبشی آشفتگی به انرژی داخلی را نشان می‌دهد(Izadinia et al., 2022).

2-1-1- تأثیر مدل ()

دقت بیشتر در شبیه‌سازی جریان‌های آشفته: این مدل با حل معادلاتی برای انرژی آشفتگی و نرخ اتلاف آن، توزیع آشفتگی را در جریان با دقت بیشتری محاسبه می‌کند. در جریان‌هایی با تغییرات شدید سرعت، مانند جریان روی سرریزها یا پشت موانع، مدل به دلیل دقت بیشترش می‌تواند نتایج دقیقی ارائه دهد.

کاربرد در جریان‌های آزاد و نیمه‌مستغرق: این مدل برای شبیه‌سازی جریان‌های آزاد، مانند جریان‌های عبوری از سرریزها، که شامل ناحیه‌های مختلف از جمله ناحیه‌های جداشده و بازگردشی هستند، مناسب است. استفاده از مدل به محققان امکان می‌دهد که رفتار جریان در نزدیکی سازه‌ها مانند سرریزها را با دقت شبیه‌سازی کنند.

همگرایی سریع‌تر: مدل نسبت به برخی دیگر از مدل‌های آشفتگی، مانند مدل‌های پیچیده‌تر Reynolds Stress، از همگرایی سریع‌تری برخوردار است. این ویژگی به معنای کاهش زمان محاسبات است، به‌ویژه در مسائلی که حجم محاسبات زیاد است.

سادگی و پایداری: مدل به دلیل ساختار ساده‌تر خود نسبت به سایر مدل‌های آشفتگی، پایداری بالاتری دارد و به‌راحتی می‌تواند در شرایط مختلف جریان به کار گرفته شود. این ویژگی باعث می‌شود که این مدل در شبیه‌سازی‌های مهندسی سیالات به طور گسترده استفاده شود.

محدودیت‌ها: مدل در شبیه‌سازی جریان‌های نزدیک به دیواره (جریان‌های لایه مرزی) یا جریان‌های دارای تغییرات شدید در مقیاس کوچک، مانند جریان‌های با گردابه‌های پیچیده یا جریان‌هایی که به شدت وابسته به شرایط مرزی هستند، دقت کمتری دارد. برای این‌گونه مسائل، مدل‌های پیچیده‌تری مانند k-ℇ یاReynolds Stress ترجیح داده می‌شوند.(Izadinia et al., 2023)

2-2- مدل‌های مورد بررسی

(Ramamurthy & Vo, 1993) و (Schmocker et al., 2011) در پژوهش‌های خود نشان دادند که زاویه وجه بالادست سرریز تاج‌دایره‌ای تاُثیر چشم‌گیری در تغییرات ضریب دبی سرریز ندارد از این رو در تحقیق حاضر شعاع تاج سرریز و ارتفاع دریچه به‌عنوان متغیر در نظر گرفته شده‌اند. گروه A با سه شعاع تاج مختلف 5، 10 و 15 سانتی‌متر و 6 دبی 180، 220، 260، 300، 340 و 380 لیتر بر ثانیه و گروه B با 4 ارتفاع دریچه مختلف 3، 5، 5/7 و 10 سانتی‌متر و 6 دبی در واحد عرض مختلف مورد ارزیابی قرار گرفته است. پارامتر ثابت در این گروه شعاع تاج 10 سانتی‌متر و زاویه وجه پایین‌دست 45 درجه است، مورد ارزیابی قرار گرفته‌اند. پارامتر ثابت در مدلهای A ارتفاع کالورت 5 سانتی‌متر و زاویه وجه پایین‌دست 45 درجه است. مشخصات مدلها در جدول 1 آورده شده است.

2-3- طراحی هندسه سازه در نرم افزار گمبیت

برای طراحی هندسه سرریز-دریچه و کانال در نرم‌افزار گمبیت ابتدا از منوی Operation در سمت راست نرم‌افزار گزینه Geometry را انتخاب کرده و سپس برای طراحی دو بعدی از سه گزینه Vertex، Edge و Face به ترتیب برای ایجاد نقطه، خط و سطح استفاده شده است. در زیر مجموعه گزینه Face که برای طراحی سطح استفاده شده گزینه‌های تبدیل خطوط به سطح، ایجاد اشکال هندسی و اجتماع و اشتراک‌گیری از سطوح نیز در طراحی دو بعدی مورد استفاده قرار گرفته است. در شکل (2) کلیه گزینه های مورد استفاده نشان داده شده است.

شکل2- جزئیات پنجره Geometry برای طراحی هندسی سرریز-دریچه و حجم کنترل

(1)

Fig 2. Details of the Geometry window for designing the geometry of the weir-gate and control volume

2-3-1- مش‌بندی هندسه سازه در نرم افزار گمبیت

برای مش‌بندی هندسه سرریز-دریچه و کانال در نرم‌افزار گمبیت ابتدا از منوی Operation در سمت راست نرم‌افزار گزینه Mesh را انتخاب کرده و سپس از گزینه‌ Face برای مش‌بندی سطح استفاده شده است. با انتخاب گزینه Mesh Face در زیر مجموعه گزینه Face شکل و نوع مش‌بندی در قسمت Element و Type به ترتیب Quad و Pava انتخاب شده و اندازه‌ی چشمه‌ها در قسمت Spacing تعیین شده است. در شکل (3) کلیه گزینه‌های مورد استفاده نشان داده شده است.

2-3-2- تعیین شرایط مرزی در نرم‌افزار گمبیت

برای تعیین شرایط مرزی حجم کنترل در نرم افزار گمبیت ابتدا از منوی Operation در سمت راست نرم‌افزار گزینه Zones را انتخاب کرده و سپس از گزینه‌ Specify Boundary Types برای تعیین شرایط مرزی استفاده شده است. در مرحله بعد ابتدا لبه‌ی مورد نظر را انتخاب کرده و سپس در قسمت Type نوع شرط مرزی برای لبه‌ی مورد نظر تعیین شده است. در شکل (4) کلیه گزینه‌های مورد استفاده نشان داده شده است.

شکل3- جزئیات پنجره Mesh برای مش‌بندی سرریز- دریچه و حجم کنترل

شکل4- جزئیات پنجره Zones برای تعیین شرایط مرزی حجم کنترل

Fig 3. Details of the Mesh window for meshing the weir-gate and control volume

Fig 4. Details of the Zones window for defining boundary conditions of the control volume

2-4- صحت‌سنجی مدل عددی

برای صحت‌سنجی این پژوهش از مدل‌سازی عددی 5 مدل آزمایشگاهی مورد بررسی در مقالات بین المللی استفاده شده است. مدل‌های آزمایشگاهی مورد بررسی در این مقالات در نرم‌افزار شبیه‌سازی شده است. برای صحت‌سنجی دقت محاسبه ضریب دبی سرریز تاج‌دایره‌ای از دو مدل آزمایشگاهی مورد بررسی توسط) (Ramamurthy & Vo, 1993استفاده شد. دو سرریز تاج‌دایره‌ای به شعاع‌های 2/15 و 54/2 سانتی‌متر به صورت عددی در نرم‌افزار فلوئنت مدل‌سازی شد. مقایسه نتایج مدل عددی با مدل آزمایشگاهی مطابقت خوبی را نشان می‌دهد.

برای صحت‌سنجی توانایی مدل‌سازی سرریز و دریچه توسط نرم افزار فلوئنت از مدل آزمایشگاهی مورد بررسی توسط (Samani & Mazaheri, 2009) استفاده شد. در این بررسی سرریز لبه‌تیز دریچه‌دار بدون فشردگی جانبی در نرم افزار فلوئنت مدل‌سازی شد. مقایسه عمق آب بالادست محاسبه شده در مدل عددی با مدل آزمایشگاهی مطابقت خوبی را نشان می‌دهد. عمق جریان در بالادست و پایین‌دست سرریز به ترتیب در فاصله‌های دو برابر و هشت برابر نسبت به مرکز سرریز اندازه‌گیری شدند (Sing and Kumar, 2022).

برای صحت‌سنجی دقت محاسبه پروفیل سرعت و به دنبال آن پروفیل پارامتر‌های آشفتگی و فشار از دو مدل آزمایشگاهی مورد بررسی توسط (Bagheri & Heydarpour, 2010) استفاده شد. در این بررسی سرریز تاج‌دایره‌ای به شعاع 065/8 سانتی‌متر با سه مدل آشفتگی متفاوت در نرم افزار فلوئنت مدل‌سازی شد. در این مدل‌سازی پروفیل سرعت تاج سرریز برای دو نسبت متفاوت 85/0 و 2/1 استخراج شد. مقایسه نتایج مدل عددی با مدل آزمایشگاهی مطابقت خوبی را نشان می‌دهد. پارامترهای شکل (5)، فاصله تاج سرریز تا کالورت،P ارتفاع سرریز، D بازشدگی دریچه سرریز و R شعاع تاج دایره‌ای است.

شکل5- مشخصات مدل‌های مورد ارزیابی

Fig 5. Specifications of the evaluated models

جدول1- مشخصات هندسی مدل‌ها

Table 1. Geometric specifications of the models

مدل	R(cm)	D(cm)	β	Q(lit/s)
A1	5	5	45	180،220،260،300،340،380
‌A2	10	5	45	180،220،260،300،340،380
A3	15	5	45	180،220،260،300،340،380
B1	10	3	45	180،220،260،300،340،380
‌B2	10	5	45	180،220،260،300،340،380
‌B3	10	5/7	45	180،220،260،300،340،380
B4	10	10	45	180،220،260،300،340،380

3- نتایج و بحث

3-1- بررسی پروفیل سطح آب

پارامتر متغیر در سرریزها در گروه A شعاع تاج سرریز است. در شکل (6) پروفیل سطح آب برای دبی 260 لیتر بر ثانیه نشان داده شده است. همان‌گونه که مشخص است، سطح آزاد جریان در بالادست با افزایش شعاع تاج افزایش یافته و پروفیل سطح آب در پایین‌دست از سرریز فاصله گرفته است. در پایین‌دست (در فاصله 8/0 متر تا انتها) سطح آب برای کلیه شعاع‌های سرریز برابر است.

شکل6- تغییرات پروفیل سطح آب در کالورت - سرریز تاج دایرهای بر اساس شعاع تاج سرریز

Fig 6. Variations in the water surface profile in the culvert-circular-crested weir based on the radius of the weir crest

پارامتر متغیر در سرریزها در گروه B مقدار بازشدگی دریچه است. در شکل (7) پروفیل سطح آب برای دبی 260 لیتر بر ثانیه نشان داده شده است. سطح آب در بالادست با افزایش بازشدگی دریچه مقدار ناچیزی افزایش یافته است. به دلیل افزایش بازشدگی دریچه و به دنبال آن افزایش دبی دریچه و کاهش دبی سرریز پروفیل سطح آب در پایین‌دست از سرریز فاصله گرفته است. در پایین‌دست (در فاصله 8/0 متر تا انتها) سطح آب برای کلیه مقادیر بازشدگی دریچه برابر است.

شکل7- تغییرات پروفیل سطح آب در کالورت- سرریز تاج‌دایرهای براساس بازشدگی دریچه

Fig 7. Variations in the water surface profile in the culvert-circular-crested weir based on the gate opening

3-2- بررسی پروفیل‌های سرعت در طول کانال و سرریز

محور افقی مقدار سرعت بی‌بعد شده با میانگین سرعت بالادست () و محور قائم ارتفاع جریان (y) نسبت به ارتفاع بیشینه پروفیل () مورد نظر است. محور افقی بالای هر شکل فاصله هر پروفیل از ابتدای سرریز تقسیم بر طول سرریز (L) را نشان می‌دهد. همان‌گونه که اشاره شده است، پارامتر متغیر در کالورت - سرریز تاج‌دایرهای گروهA شعاع تاج سرریز است. شکل (8- الف) توزیع سرعت را در بالادست سرریز نشان می‌دهد. با توجه به شکل جریان دارای مقدار سرعت بیشینه در نزدیک سطح آزاد است و با حرکت به سمت‌ پایین‌دست جریان یکنواخت و توسعه‌یافته می‌شود. با توجه به اینکه در سرریز با شعاع تاج 15 سانتی‌متر سرعت جریان در نزدیکی بستر کاهش و با نزدیک شدن به سطح آب نسبت به دو شعاع دیگر 5 و 10 سانتی‌متر، 17 درصد اختلاف داشته می‌توان بیان نمود افزایش شعاع تاج سرریز باعث اثر گذاشتن در جریان بالادست سرریز شده است در حالیکه در دو شعاع دیگر پروفیل‌های سرعت تقریباً یکسان هستند. با نزدیک شدن به سرریز جریان از حالت یکنواخت خارج شده و دارای مقدار بیشینه‌ای در ورودی کالورت است. سرعت ورودی آب به کالورت با افزایش شعاع سرریز کمتر از 4 درصد افزایش یافته است و با نزدیک شدن به سطح آزاد این روند معکوس شده و بیشترین سرعت مربوط به پروفیل با کمترین شعاع تاج سرریز است. شکل (8-ب) تغییرات سرعت را در طول کالورت - سرریز نشان می‌دهد. همان‌گونه که مشخص است بیشینه سرعت جریان در دریچه نزدیک کف کانال و در ابتدای کالورت است. بیشینه سرعت جریان کالورت با کاهش شعاع تاج کمتر از 3 درصد کاهش می‌یابد. تغییرات سرعت روی سرریز نشان می‌دهد بیشینه سرعت جریان روی تاج سرریز است که با کاهش شعاع 13 درصد افزایش می‌یابد که نشان دهنده تأثیر کاهش شعاع در افزایش سرعت جریان تاج سرریز است. بیشینه سرعت جریان قبل از تاج سرریز نزدیک سطح آزاد است و با نزدیک شدن به تاج سرریز و پایین‌دست تاج این روند معکوس شده و بیشینه سرعت جریان در نزدیک سطح سرریز است. در پایین‌دست سرریز با کاهش شعاع سرعت کاهش یافته است. شکل (8-ج) تغییرات سرعت را در پایین‌دست سازه نشان می‌دهد. مشاهده می‌شود که با فاصله گرفتن از سازه جریان یکنواخت شده و بیشینه سرعت جریان در نزدیک سطح آزاد است. درصد اختلاف سرعت جریان پایین‌دست سازه بین شعاع 5 و 10 سانتی‌متر، 7 درصد و بین شعاع 10 و 15 سانتی‌متر کمتر از 3 درصد است که نشان دهنده تأثیر کم شعاع در سرعت پایین‌دست سازه است.

پارامتر متغیر در کالورت-سرریز تاج دایره‌ای مقدار بازشدگی دریچه است. شکل (9- الف) توزیع سرعت را در بالادست سرریز نشان می‌دهد. با توجه به شکل جریان دارای مقدار سرعت بیشینه در نزدیک سطح آزاد است و با حرکت به سمت‌ پایین‌دست جریان یکنواخت و توسعه‌یافته می‌شود. سرعت جریان با کاهش مقدار بازشدگی دریچه 10 درصد افزایش می‌یابد که علت آن کاهش سطح مقطع جریان با داشتن یک دبی ثابت است. با نزدیک شدن به سرریز جریان از حالت یکنواخت خارج شده و دارای مقدار بیشینه‌ای در ورودی دریچه و نزدیک کف کانال است. مقدار بیشینه سرعت ورودی آب به دریچه با افزایش بازشدگی از کف کانال فاصله گرفته و 11 درصد افزایش یافته است. بیشترین مقدار بیشینه سرعت مربوط به بازشدگی 10 سانتی‌متر و در فاصله از کف کانال اتفاق افتاده است. شکل (9-ب) تغییرات سرعت را در طول سرریز-دریچه نشان می‌دهد. همان‌گونه که مشخص است بیشینه سرعت جریان در دریچه نزدیک کف کانال و در ابتدای دریچه است. با افزایش بازشدگی دریچه سرعت جریان در دریچه در حدود 11 درصد افزایش یافته است. مقدار سرعت تا مقطع بر روی سرریز با افزایش بازشدگی دریچه کاهش یافته که نشان دهنده کاهش دبی عبوری روی سرریز است و بعد از آن با توجه به کاهش سطح مقطع جریان بیشینه سرعت تقریباً ثابت شده است. با افزایش بازشدگی دریچه سرعت جریان روی تاج سرریز 5/11 درصد کاهش یافته است. شکل (9-ج) تغییرات سرعت را در پایین‌دست سرریز نشان می‌دهد. مشاهده می‌شود که با فاصله گرفتن از سرریز جریان یکنواخت شده و بیشینه سرعت جریان در نزدیک سطح آزاد است. سرعت جریان در پایین‌دست سرریز با افزایش بازشدگی دریچه 3 درصد افزایش یافته که نشان دهنده تأثیر ناچیز دریچه در تغییرات سرعت پایین‌دست سازه است.

شکل 8- تغییرات پروفیل‌های سرعت بر اساس شعاع تاج سرریز الف) در بالادست، ب) در طول سرریز-دریچه و ج) در پایین‌دست کالورت- سرریز تاج‌دایرهای

Fig 8. Variations in velocity profiles based on the radius of the weir crest: a) upstream, b) along the weir-gate, and c) downstream of the culvert-circular-crested weir

شکل 9- تغییرات پروفیل‌های سرعت بر اساس بازشدگی دریچه الف) در بالادست، ب) در طول سرریز-دریچه و ج) در پایین‌دست کالورت- سرریز تاج‌دایره‌ای

Fig 9. Variations in velocity profiles based on the gate opening: a) upstream, b) along the weir-gate, and c) downstream of the culvert-circular-crested weir

4- نتیجه‌گیری

نتایج بدست آمده از این تحقیق نشان می‌دهد که ضریب دبی در سرریزهای دریچه‌دار بیشتر از سرریزهای بدون دریچه است. همچنین با کاهش شعاع تاج سرریز، ضریب دبی افزایش می‌یابد. در سرریزهای گروه A، کمترین ضریب دبی مربوط به شعاع تاج 15 سانتی‌متر است که ضریب دبی در محدوده 2/1 تا 34/1 تغییر می‌کند، در حالی که بیشترین ضریب دبی مربوط به شعاع 5 سانتی‌متر است که ضریب دبی در محدوده 42/1 تا 52/1 تغییر دارد. علاوه بر این، با افزایش بازشدگی دریچه در نسبت H1/R، مقدار ضریب دبی افزایش یافته و اثر دریچه بر ضریب دبی تقویت می‌شود. در سرریزهای گروه B، بیشترین ضریب دبی مربوط به بازشدگی دریچه 10 سانتی‌متر است که ضریب دبی در محدوده 3/1 تا 46/1 تغییر می‌کند و کمترین ضریب دبی مربوط به بازشدگی 3 سانتی‌متر است که ضریب دبی در محدوده 27/1 تا 4/1 تغییر دارد. این نتایج نشان می‌دهند که تغییرات در شعاع تاج و بازشدگی دریچه تاثیر قابل توجهی بر ضریب دبی سرریزها دارند و می‌توانند در بهینه‌سازی طراحی‌های هیدرولیکی موثر باشند.

5- تضاد منافع نویسندگان

نویسندگان این مقاله اعلام می‌دارند که هیچ تضاد منافعی در رابطه با نویسندگی و یا انتشار این مقاله ندارند.

6- منابع

Bagheri, S., & Heydarpour, M. (2010). Application of Euler equations in determining the discharge coefficient of circular crest spillways. Iranian Journal of Water and Soil Research, 40(1)67-73. (In Persian)

Beigi Poor, G., & Bina, M. (2005). Hydraulic flow in circular and cylindrical crested weirs. 5th Iranian Hydraulic Conference, Shahid Bahonar University of Kerman, Kerman, Iran, 1-8.

Bhajantri, M.R., Eldho, T.I., & Deolalikar, P.B. (2006). Hydrodynamic modelling of flow over a spillway using a two-dimensional finite volume-based numerical model. Sadhana, 31(6), 743-754. https://doi.org/10.1007/BF02716893

Ghobadian, R., Farmani Fard, M., & Karami, A. (2012). Effect of crest roughness and body on discharge coefficient and hydraulic conditions of flow in circular crested weirs using a physical model. Journal of Soil and Water Sciences, 22(1), 1-8.

Izadinia, E., & Heidari, D. (2023). Flow simulation in a combined culvert-circular crown spillway structure, 22nd Iranian Hydraulic Conference, Maragha. https://civilica.com/doc/1963858

Izadinia, E., & Heidari, D. (2022). Flow simulation in a combined structure of culvert and circular crown spillway, 21st National Conference of Iranian Hydraulics, Ahvaz. https://civilica.com/doc/1694896

Negm, A.M., Al-Brahim, A.M., & Alhamid, A.A. (2002). Combined-free flow over weirs and below gates. Journal of Hydraulic Research, 40(3), 359-365. https://doi.org/10.1080/00221680209499950

Ramamurthy, A.S., & Vo, N.D. (1993). Characteristics of circular crested weir. Journal of Hydraulic Engineering, 119(9), 1055–1062.

Samani, J., & Mazaheri, M. (2009). Combined flow over weir and under gate. Journal of Hydraulic Engineering, 135(3), 224–227. https://doi.org/10.1061/(ASCE)0733-9429(2009)135:3(224)

Schmocker, L., Halldórsdóttir, B., & Hager, W. (2011). Effect of weir face angles on circular-crested weir flow. Journal of Hydraulic Engineering, 137(6), 637–643. https://doi.org/10.1061/(ASCE)HY.1943-7900.0000346

Severi, A., Masoudian, M., Kordi, E., & Roettcher, K. (2014). Discharge coefficient of combined-free over-under flow on a cylindrical weir-gate. ISH Journal of Hydraulic Engineering. 42-52. https://doi.org/10.1080/09715010.2014.939503

شارک

عنوان URL للمقالة

تحلیل عددی هیدرولیک جریان در سیستم کالورت-سرریز تاج‌دایره‌ای همراه با دریچه

سند

الروابط

المراكز ذات الصلة

دعامة

الصفحات الرسمية