بخشی از مقاله
خلاصه
پرتابکننده جامی شکل، از جمله سازههای استهلاک انرژی استکه معمولاً در سدهای بلند و در مواردی که سرعت جریان بیش از 15 تا 20 متر بر ثانیه باشد مورد استفاده قرار میگیرد. از آن جا که جریان در پرتابکننده جامی در عمل، به شدت آشفته، غیر دائمی و سه بعدی میباشد، بنابراین در این پژوهش، با در نظر گرفتن هندسهای متفاوتِ جام - باکت - ، شامل مقادیر مختلف شعاع انحناء، زاویه انحراف جریان و ارتفاع جام و نیز اعمال شرایط مختلف جریان نظیر مقادیر مختلف عمق، سرعت و عدد فرودِ ورودی، شبیهسازی عددی جریان آشفته با استفاده از مدل FLUENT صورت گرفته است.
این مدل عددی برای حل جریان آشفته از انواع مدل های مختلف مدلهای آشفتگی استفاده مینماید. در این تحقیق با شبیهسازی جریان آشفتهی 3 بعدی و 2 فازی - آب و هوا - به روش VOF ، در جام و محاسبه حداکثر هد فشار دینامیکی، مشاهده گردید که مدل آشفتگی - k استاندارد بیشترین انطباق را با نتایج حاصل از مدل آزمایشگاهی نشان می دهد. بر این اساس، محل وقوع حداکثر فشار دینامیکی، به هندسه پرتابکننده جامی بستگی دارد و مقدار آن نیز وابسته به عدد فرود جریان ورودی، عمق جریان، زاویه انحراف جریان در جام و همچنین شعاع جام بوده و تقریباً در قسمت میانیِ گودیِ جام رخ میدهد.
.1 مقدمه
مستهلک کنندههای انرژی، سازههایی هستند که در انتهای سرریزها و تخلیهکنندههای تحتانی سدها به منظور استهلاک انرژی مازاد جریان ساخته میشوند. پرتابکنندهی جامی در سدهای بلند و در مواردی که سرعت جریان بیش از 15 تا 20 متر بر ثانیه باشد مورد استفاده قرار میگیرد. عملکرد هیدرولیکی و نیز ارزانتر بودن این سازه باعث شده است که در صورت مناسب بودن شرایط زمین شناسی بستر رودخانه در پایاب، در بیشتر سدها از پرتابکننده جامی به عنوان سازه مستهلک کننده انرژی استفاده شود.[2]
به طور کلی حوضچههای پرتابی جریان، شامل پنج قسمت به شرح زیر میباشد: -1 تنداب ورودی، -2 جام یا قسمت پرتاب کننده جریان - دفلکتور - ، -3 پخش و پراکنده شدن جت جریان در هوا، -4 برخورد جت جریان با کف بستر و ایجاد فرسایش و آبشستگی، -5 محدوده پایاب. در شکل 1، نمونهای از طرح شماتیک یک پرتابکنندهی جامی - فلیپ باکت - نشان داده شده است.
جریان در پرتابه جامی شکل در طبیعت، به صورت سه بعدی، غیر ماندگار و به شدت آشفته میباشد. مطالعه عددی بر روی جریان آشفته در سرریز سدها از گذشته تا کنون مورد توجه بسیاری از محققین بوده است. به عنوان نمونه، کراگر و راچمن در سال 2006 با توسعه معادلات آبهای کم عمق، مدل سه بعدی برای شبیهسازی جریانهای فوق بحرانی ارائه نمودند. لارس و همکاران در سال 2006 به بررسی منحنی مسیر پرتابه جت جریان و نیز مقادیر فشار در کف پرتابکننده جامی در پایانه تنداب افقی پرداختند.[3]
آنها تحلیل خود را در دو حالت، دو بعدی و سه بعدی، و با استفاده از مدل عددی اجزای محدود جزئی - PFEM - انجام داده و تأثیر افزایش عدد فرود و دبی جریان را بر روی طول پرتابه و فشار کف جام بررسی نمودند. صباغ یزدی و رستمی - 1386 - در پژوهشی به بررسی عمق، سرعت و توزیع فشار بر روی بسترهای دارای انحنای قائم، توسط مدل FLOW-3D پرداختند.
پیروز و همکاران - 1388 - کاربرد مدلهای آشفتگی را در شبیهسازی جریان هوادهی شده در تخلیهکننده تحتانی سد سفیدرود، با استفاده از مدل عددی FLUENT بررسی نمودند.[4] در تحقیق حاضر، با بیان معادلات حاکم بر انواع مدلهای آشفتگی موجود در مدل FLUENT ، جریان آشفتهی سه بعدی در سازه پرتابه جامی شکل، شبیهسازی گردیده و پارامتر حداکثر هد فشار دینامیکی در جام، مبنای صحتسنجی مدل عددی قرار گرفته است.
.2 مدل آزمایشگاهی
به منظور صحت سنجی مدل عددی و یافتن بهترین مدل آشفتگی برای مقدار حداکثر هد فشار در پرتابکننده جامی، از آزمایشات هلر و همکاران - 2005 - استفاده شده است. این آزمایشات در یک کانال مستطیلی به عرض 0/5 متر، ارتفاع 0/7 متر و طول 7 متر انجام شده است. کف کانال و دیواره سمت راست از جنس PVCو دیواره چپ کانال، از جنس شیشه است.
جریان تقربِ فوق بحرانی تا مقدار 0/186 مترمکعب بر ثانیه اندازهگیری شده است. در این آزمایشات، زوایای مختلف انحراف جریان، اعداد فرود مختلف و اعماق مختلف جریان مورد بررسی قرار گرفته است. شکل 2، طرح مدل آزمایشگاهی مذکور را نشان میدهد.[7] همچنین جدول 1، نشاندهنده شرایط هندسی آزمایشات صورت گرفته در جام پرتابی میباشد.
.3 مدل عددی
در این تحقیق به منظور شبیهسازی جریان آشفته از مدل عددی FLUENT 6.3.26 استفاده گردیده است. FLUENT یک برنامه CFD جهانی است که با زبان برنامهنویسی C نوشته شده است و برای محدوده زیادی از کابردهای سیال مورد استفاده قرار میگیرد. با توجه به ابعاد مدل آزمایشگاهی، شبکه بندی مدل عددی در نرم افزار Gambit با ابعاد سلولی 3 سانتیمتر یک بررسی بر روی انواع مشبندی ممکن در Gambit انجام گردید و بر این اساس، در قسمت کانال تنداب از مش مربعی با الگوی Pave و در نواحی نزدیک به انحنای جام و دفلکتور پرتابی و ناحیه خروجی جریان و همچنین در قسمت مرز هوای مسئله، از مش چهار وجهی - Tetrahedron - استفاده گردید.
بنابراین، تعداد مش در شبکهی منتخب مسئله، 500564 عدد میباشد. به منظور انجام صحتسنجی مدل عددی با مدل آزمایشگاهی - برای پارامتر حداکثر هد فشاری - ، مدلسازی جریان در 4 حالت - آزمایش شماره 1 تا - 4 که در جدول 2، ذکر شده اجرا گردید. در هر یک از این حالتها، تمامی مدلهای آشفتگی موجود در نرمافزار FLUENT تا رسیدن به همگرایی کامل اجراء و نتایج نهایی در جدول ذیل بیان گردیده است. همچنین در شکلهای5 تا 9، بترتیب بردارهای سرعت جریان و نحوه توزیع سرعت و فشار جریان در پرتابه جامی و نیز کانتور شدت آشفتگی در جریان ورودی به پرتابکننده جامی و کانتور نسبت حجمی فازهای سیال در جریان 2 فازی آب و هوا، نشان داده شده است.[8]
