سعيد رضا صباغ يزدی | نرم افزارهاي علمي | صفحه اصلي | |||||
نرم افزار
تحلیل مسائل علمی و
صنعتی نصير
|
|||||||
کاربرد
نرم افزار احجام محدود NASIR تحلیل مسائل علمی و صنعتی
در مدلسازی جریان میانگین عمقی با سطح آزاد
(در کانالها و رودخانه ها)
سعيدرضا صباغ يزدي
دانشیار دانشکده مهندسی عمران
دانشگاه صنعتي خواجه نصير الدين طوسي
SYazdi@kntu.ac.ir
بهار 1386
معرفی نرم افزار حجم محدود تحلیل مسائل علمی و صنعتی NASIR
1. تحلیل مسائل علمی و صنعتی در مدلسازی جریان میانگین عمقی با سطح آزاد (در کانالها و رودخانه ها)
2. تحلیل مسائل علمی و صنعتی در مدلسازی جریان میانگین عمقی با سطح آزاد (فرسایش، انتقال و ته نشست رسوب)
3. تحلیل مسائل علمی و صنعتی در مدلسازی جریان غیردائمی در لوله ها
4. تحلیل مسائل علمی و صنعتی در مدلسازی تولید و انتقال حرارت در بتن
5. تحلیل مسائل علمی و صنعتی در مدلسازی جریان در محیط متخلخل
7. تحلیل مسائل علمی و صنعتی در مدلسازی جریان دو بعدی باد و انتقال آلایندها
8. تحلیل مسائل علمی و صنعتی در مدلسازی جریان سه بعدی باد
مقدمه
در حال حاضر این ماجولها برای حل عددی مسائلی همچون موضوعات زیر مناسب میباشند:
1. جریان با سطح آزاد آشفته در مجاری با پیچیدگیهای هندسی
2. جریان غیرماندگار آشفته دارای سطح آزاد در اطراف پایه ها با شکل دلخواه
در نوشتار حاضر نمونه هائي از نتايج آزمونهای صحت (بكمك مقايسه با آزمونهاي تحليلي و يا اندازهگيريهاي آزمايشگاهي) و کاربرد ماجولهای نرم افزار مذکور در مسائل مرتبط با جریان میانگین عمقی با سطح آزاد (در کانالها و رودخانه ها) ارائه میگردد. لازم بذکر است که جزئیات فنی هر کدام از موارد در مقالات منتشر شده توسط نگارنده و همکاران مورد تشریح قرار گرفته است. در این نوشتار برای هر یک از موارد، پس از تشریح مسئله مدلسازی شده، به اشاره مختصر به نمونه ای از نتایج بسنده شده است. لذا جهت حفظ اختصار از ذکر توضیحات مربوط به مسئله و مدلسازی آن شامل فنون محاسباتی خاص مورد مذکور، مراحل صحت سنجی کار مدلسازی و نتایج تفصیلی بدست آمده پرهیز شده است. علاقمندان میتوانند جزئیات کامل فنون و روند صحت سنجی و شبیه سازی هر آزمون یا کاربرد را در مقاله و یا مقالات مربوط به آن مسئله دنبال نمایند. مشخصات مقاله و یا مقالات مربوط به هر مسئله تحلیل شده در پاورقی صفحه مربوط به همان مورد آورده شده است.
جريان با سطح آزاد زيربحراني بياستهلاك در كانال با خم ه180
قبل از مدلسازی یک جریان واقعی ابتدا سعی میشود از حل بدون خطای تک تک عبارتهای معادلات حاکم بر مسئله اطمینان حاصل شود. لذا در اینجا برای اطمینان از حل دقیق عبارتهای انتقالی در غیاب عبارتهای استهلاکی (میرا کننده خطاهای عددی) در معادلات جریان با سطح آزاد، جريان زيربحراني در يك كانال افقي بدون اصطكاك داراي يك خم 180 درجه (که برگشت کامل جریان را در بر دارد) مداسازی شده است. لازم به ذكر است كه قبلاً توسط (De Vreind, 1977) آزمايشاتي بر روي كانال فوق انجام شده بود. عرض كانال مورد استفاده 7/1 متر و شعاع مياني خم نيز برابر با 25/4 متر است و براي اينكه جريان پيش از وارد شدن به خم و پس از خروج از آن فرصت كافي براي شكلگيري داشته باشد دو كانال مستقيم بطول 6 متر در ابتدا و انتهاي خم تعبيه شده است. مقدار دبي جريان در بالادست كانال برابر با 19/0 متر مكعب بر ثانيه و مقدار عمق نيز در پاييندست برابر با 18/0 متر تنظيم شده است. در اين خم نسبت شعاع ميانگين به عرض كانال برابر با 5/2 ميباشد(Jia & Wang, 1999).
براي شبيهسازي عددي جريان در اين كانال از يك شبكه بدونساختار مثلثي استفاده شد كه اندازه اين شبكه در محل خم تا حدي ريزتر شده است. اين شبكه از 1574 گره، 2742 مثلث و 4315 يال تشكيل شده است. نتايج بدست آمده از مدل در اين مرحله در شكل زير ارائه شده است.
الف) نواحي هم عمق محاسبه شده، ب) بردارهاي سرعت محاسبه شده
(در نواحي هم رنگ مقادير عمق و سرعت يكسان است).
بطوريكه در شكل فوق مشاهده ميشود، در محدوده خم موجود در كانال، در نزديكي ديواره داخلي عمق كاهش و سرعت افزايش يافته است. همچنين در نزديكي ديواره خارجي عمق افزايش و مقدار سرعت كاهش يافته است. علاوه بر اين جريان پديد آمده در خم از تقارن بسيار كاملي برخوردار است و منحنيهاي همعمق نيز در اين محدوده با يكديگر موازيند.
در قسمتهاي مستقيم كانال مقادير عمق و سرعت ثابت شدهاند، بطوريكه در اين نواحي سطح آب تقريباً ثابت شده است. نكته قابل توجه ديگر اين است كه مقدار پارامترهاي عمق و سرعت در قسمتهاي ابتدايي و انتهايي كانال يكسان هستند و به همين علت ميتوان گفت كه هيچگونه استهلاك ناخواسته عددي وارد محاسبات نشده است و انطباق خوب نتايج مدل عددي با حل تحليلي نشان از دقت مناسب مدل تدوين شده دارد[8] .
جريان باسطح آزاد زيربحراني بدوناستهلاك در مجراي پيچان
در اين مرحله شبيهسازي عددي جريان در يك كانال مارپيچ بدون شيب كه قبلا توسط (Chang,1971) برروي آن آزمايشاتي انجام شده بود، انجام گرديد. كانال مورد استفاده در اين مرحله از دو خم متناوب 90 درجه در دو جهت متفاوت تشكيل شده است كه اين دو خم توسط يك كانال مستقيم به يكديگر متصل شدهاند. قبل از خم اول و پس از خم دوم نيز كانالهاي مستقيمي تعبيه شده است تا جريان پيش از وارد شدن به خمها كاملاً شكل گرفته باشد. در اين كانال با زبري كم عمق اوليه برابر با 115/0 متر و سرعت اوليه برابر با 366/0 متر بر ثانيه گزارش شده است و لذا عدد فرود برابر با 35/0 بدست ميآيد. در خمهاي اين كانال عرض كانال برابر با33/2 متر و شعاع ميانگين آنها نيز برابر با 53/8 متر است. بنابراين نسبت شعاع ميانگين به عرض كانال در اين خمها برابر با 66/3 خواهد بود و ميتوان گفت كه جريانهاي ثانويه در آن ناچيز بوده و جريان غالب در خمها دو بعدي است(Ye & Mc.Corquodale, 1997) .
شبكه مورد استفاده در اين مرحله، يك شبكه بدون ساختار مثلثي است كه به روش دلاوني توليد گرديده است و اندازه شبكه در شعاع داخلي خمها تا حدي ريز شده است. اين شبكه از 1748 گره، 3069 مثلث و 4816 يال تشكيل شده است. در اين آزمون تمامي استهلاكهاي فيزيكي ناچيز بودهاند و جرياني غيرلزج و بدون اصطكاك شبيهسازي شده است. نتايج بدست آمده از مدل در اين مرحله در اشكال زير ارائه شده است. نتايج محاسباتي جريان متقارني را نشان ميدهد. همانطور كه انتظار ميرود در نزديكي ديوارهاي داخلي خمها مقدار سرعت افزايش و عمق كاهش يافته است و در نزديكي ديوارهاي خارجي خمها مقدار سرعت كاهش و عمق افزايش يافته است. در قسمتهاي مستقيم كانال سطح آب افقي و سرعت در عرض كانال يكسان است. ثابت بودن مقادير عمق و سرعت در محور مياني كانال از عدم وجود خطا در حل جريان بياستهلاك حكايت دارد [9].
نواحي همعمق محاسبه شده (در نواحي همرنگ مقادير عمق يكسان است)
بردارهاي سرعت محاسبه شده (در نواحي همرنگ مقادير سرعت يكسان است)
شبيهسازي جريان چرخشي در بازشدگي بيزبري
در اين كار نتايج مدل عددي براي شبيهسازي جريانهاي بي استهلاک بستر داراي چرخه شده است. چنين جرياناتي در کانالهاي با زبري بستر ناچيز و يا داراي عمق زياد بوقوع ميپيوندند. حل صريح جريانهاي فاقد استهلاک فيزيکي با مشکل رشد نوسانات ناخواسته عددي و ناپايداري روند حل مواجه است. در مدل تدوين شده براي غلبه بر نوسانات عددي پديد آمده در حل صريح معادلات در غياب عبارتهاي ميراگر نوسانات (عبارتهاي بيانگر تاثير اصطكاك، لزجت و آشفتگي) از روش لزجت مصنوعي استفاده شده است. روش استفاده شده براي اعمال لزجت مصنوعي مناسب شبكههاي بيساختار بوده و نوسانات عددي ناخواسته را بگونهاي از بين ميبرد كه صدمهاي به دقت نتايج حل وارد نگردد. با استفاده از اين مدل جريان چرخشي در يك كانال با بازشدگي ناگهاني تحت تاثير آشفتگي شبيهسازي شده است و نتايج حاصل از مدل با اندازه گيريهاي آزمايشگاهي مقايسه شدهاند[12].
شبكه بي ساختار مثلثي كانال با بازشدگي ناگهاني
خطوط جريان محاسبه شده نشانگر جريان با چرخه هاي پايدار در پشت بازشدگي کانال
مقايسه سرعت آب در امتداد عرض كانال بدست آمده از مدل عددي و اندازه گيريهاي آزمايشگاهي در سه مقطع از کانال
جريان چرخشي در بازشدگي كانال بازبري بستر
چرخابه هاي جريان کم عمق که در کانالهاي داراي باز شدگي ناگهاني عمدتا از زبري بستر تاثير ميگيرند. در مدلسازي جرياناتي که تاثير استهلاک در آنها غالب ميباشد؛ روش مدلسازي آشفتگي اثر قابل ملاحظه اي در کيفيت نتايج حل عددي دارد. در جريانهاي داراي عمق کم براي مدلسازي رياضي جريان دو بعددي ميتوان از معادلات ميانگين عمقي به همراه مدل آشفتگي مناسب استفاده نمود. معادلات مورد استفاده با بكارگيري احجام محدود قابل گسسته سازي مي باشند. مدلهاي چسبندگي گردابه اي ثابت و شبيه سازي گردابه بزرگ براي مدلسازي آشفتگي جريان مورد آزمون قرار گرفت. جهت صحت يابي، نتايج مدل عددي با اندازه گيريهاي آزمايشگاهي جريان در يك كانال با گشودگي ناگهان مقايسه شده است[13].
خطوط جريان محاسبه شده روي شبكه بيساختار مثلثي توسط مدل عددي
خطوط جريان محاسبه شد براي چند آزمون که اندازه گيريهاي آزمايشگاهي براي آن گزارش شده است
مقايسه نتايج مدل عددي با اندازه گيريهاي آزمايشگاهي نسبت عرض بازشدگي به طول چرخه براي چند آزمون مختلف
شبيهسازي جريان چرخشي درکانال شیب دار با بازشدگي ناگهانی
در اين آزمون نتايج مدل عددي براي شبيهسازي جريانهاي آشفته چرخابی در کانال با بستر با شیب 1000/1 ارزیابی شده است. اين جريان در کانالي با زبري بستر ناچيز و يا داراي عمق زياد بوقوع ميپيوندد. حل صريح جريانهاي فاقد استهلاک فيزيکي با مشکل رشد نوسانات ناخواسته عددي و ناپايداري روند حل مواجه است. در مدل تدوين شده براي غلبه بر نوسانات عددي پديد آمده در حل صريح معادلات در غياب عبارتهاي ميراگر نوسانات (عبارتهاي بيانگر تاثير اصطكاك، لزجت و آشفتگي) از روش لزجت مصنوعي استفاده شده است. روش استفاده شده براي اعمال لزجت مصنوعي مناسب شبكههاي بيساختار بوده و نوسانات عددي ناخواسته را بگونهاي از بين ميبرد كه حتی در جریانات آشفته صدمهاي به دقت نتايج حل وارد نگردد. با استفاده از اين مدل جريان چرخشي در يك كانال با بازشدگي ناگهاني تحت تاثير آشفتگي شبيهسازي شده است و نتايج با اندازه گيريهاي آزمايشگاهي مقايسه شدهاند[14].
مشخصات هندسی آزمون جریان عمیق در کانال آزمایشگاهی با بازشدگي ناگهاني
خطوط جريان محاسبه شده نشانگر جريان با چرخه هاي پايدار در پشت بازشدگي کانال
مقايسه سرعت آب در امتداد عرض كانال بدست آمده از مدل عددي و اندازه گيريهاي آزمايشگاهي در سه مقطع از کانال
(بالا : و پایین: )
جريان چرخشي در پشت آبشكن عمود بر جريان
يكي از مهمترين خصوصيت جريان تحت تاثير سازههايي نظير آبشكنها، توسعة يك ناحية چرخشي در پايين دست ميباشد. جريان كمعمق چرخشي عمدتا در آبراهههاي طبيعي به صورت انفصال جريان از خط ساحلي مشاهده ميشود. جريان آب وقتي به نوك آبشكن ميرسد به طرف ساحل مقابل شتاب ميگيرد. سپس در فاصلهاي در پاييندست آبشكن در ناحية اتصال مجدد جريان مجددا به ساحل اوليه ميرسد. در نتيجه در پاييندست آبشكن يك منطقة چرخشي شكل ميگيرد كه طول آن تا چندين برابر طول آبشكن امتداد مييابد. طول اين ناحيه چرخشي با تغيير شرايط هيدروليكي جريان (عمق، سرعت و زبري) و تغيير هندسة كانال (عرض كانال و طول آبشكن) تغيير پيدا ميكند. در کار حاضر اثر زبري بستر و لزجت آشفتگي در قسمت نيروهاي كالبدي معادلات حرکت لحاظ ميشود. ضريب مقاومت بستر توسط رابطه مانينگ و رابطة اصلاحي ضريب زبري محاسبه شده و تاثير آن بر طول گردابه پاييندست آبشكن مورد بررسي قرار گرفته است. براي صحتسنجي مدل عددي از نتايج آزمايشگاهي مقالة آقايان تينگ سانچالي و ماهسواران 1990 استفاده شده است . مدل براي شرايط مختلف هيدروليکي و براي آبشکنها با طولهاي مختلف اجرا شده است و طول گردابه ايجاد شده در پاييندست اندازهگيري شده است[15].
شبكه بي ساختار مثلثي كانال با آبشکن عمود بر جريان
خطوط جريان محاسبه شده نشانگر جريان با چرخه هاي پايدار در پشت آبشکن عمود بر جريان
خطوط هم عدد فرود محاسبه شده جريان با چرخه هاي پايدار در پشت آبشکن
مقايسه طول گردابه پشت آبشكن بدست آمده از مدل عددي و اندازه گيريهاي آزمايشگاهي
جريان چرخشي در پشت آبشكن مايل
در اينجا چرخه هاي جريان در پايين دست آب شکنها و بندرگاهها بوسيله نرم افزار حل احجام محدود معادلات آبهاي کم عمق مدلسازي شده است. در جريانها کم عمق در برابر طول و عرض ناچيز بوده و تغييرات سطح آب نيز تدريجي است و ميتوان فشار را در عمق هيدرواستاتيک فرض کرد. اثر زبري بستر و لزجت آشفنگي در قسمت نيروهاي كالبدي معادلات حرکت لحاظ ميشود. ضريب مقاومت بستر توسط رابطه مانينگ و رابطة اصلاحي ضريب زبري محاسبه شده و تاثير آن بر طول گردابه پاييندست آبشكن مايل مورد بررسي قرار گرفته است. براي صحتسنجي مدل عددي از نتايج آزمايشگاهي مدلسازي فيزيكي انجام شده در مركز تحقيقات آب وزارت نيرو استفاده شده است. مدل براي شرايط مختلف هيدروليکي و براي آبشکنها با طولهاي مختلف اجرا شده است و طول گردابه ايجاد شده در پاييندست آبشکن براي واسنجي نتايج استفاده شده است[16].
شبكه بي ساختار مثلثي كانال با آبشکن مايل
خطوط جريان محاسبه شده نشانگر جريان با چرخه هاي پايدار در پشت آبشکن مايل
خطوط هم عدد فرود محاسبه شده جريان با چرخه هاي پايدار در پشت آبشکن مايل
مقايسه طول گردابه پشت آبشكن بدست آمده از مدل عددي و اندازه گيريهاي آزمايشگاهي
جريان چرخشي در پشت آبشكن های مختلف
در اين قسمت نتایج شبیه سازی چرخه هاي جريان در پايين دست آب شکنها بوسيله حل احجام محدود جهت مند معادلات آبهاي کم عمق ارائه شده است. اینکار با استفاده از شیوه مرکز سلول روش تحلیلگر ریمان (گودونوف) برای تضمین پایداری و دقت محاسبات انجام شده است. در جريانها کم عمق در برابر طول و عرض ناچيز بوده و تغييرات سطح آب نيز تدريجي است و ميتوان فشار را در عمق هيدرواستاتيک فرض کرد. مقاومت بستر با استفاده از رابطه مانينگ محاسبه شده است. براي صحتسنجي مدل عددي از اندازه گیریهای آزمايشگاهي گزارش شده در مراجع استفاده شده است. مدل براي براي آبشکنها با طولهاي مختلف اجرا شده است و طول و زوایای مختلف باجرا درآمد و طول گردابه پاييندست و سرعت در نوک آبشکن براي واسنجي استفاده شد[17].
شبكه بي ساختار مثلثي كانال با آبشکن در جهات مختلف (راست عمود بر جریان، وسط:مایل به پایین دست، چپ: مایل به بالادست)
خطوط جريان محاسبه شده در كانال با آبشکن در جهات مختلف (راست عمود بر جریان، وسط:مایل به پایین دست، چپ: مایل به بالادست)
نقشه مناطق هم عدد فرود محاسبه شده در كانال با آبشکن در جهات مختلف (راست عمود بر جریان، وسط:مایل به پایین دست، چپ: مایل به بالادست)
زاويه نصب موج شکن |
طول موج شکن(cm) |
مدل |
Vapp(m/s) |
Vtip(m/s) |
L(m) |
L/l |
خطايL/l (درصد) |
خطايVtip (درصد) |
90 |
20 |
فيزيکي |
25/0 |
31/0 |
35/2 |
75/11 |
- |
- |
عددي* |
25/0 |
29/0 |
40/2 |
12 |
1/2 |
5/6 |
||
عددي** |
25/0 |
285/0 |
75/1 |
75/8 |
5/25 |
1/8 |
||
90 |
30 |
فيزيکي |
25/0 |
35/0 |
60/3 |
12 |
- |
- |
عددي* |
25/0 |
34/0 |
75/3 |
5/12 |
2/4 |
9/2 |
||
عددي** |
25/0 |
33/0 |
52/3 |
73/11 |
5/7 |
2/3 |
||
45 |
30 |
فيزيکي |
4/0 |
49/0 |
0/3 |
10 |
- |
- |
عددي* |
4/0 |
47/0 |
82/2 |
4/9 |
0/6 |
1/4 |
||
عددي** |
4/0 |
45/0 |
79/2 |
3/9 |
0/7 |
2/8 |
||
135 |
30 |
فيزيکي |
4/0 |
51/0 |
0/3 |
10 |
- |
- |
عددي* |
4/0 |
48/0 |
94/2 |
8/9 |
0/2 |
1/5 |
||
عددي** |
4/0 |
47/0 |
82/2 |
4/9 |
0/6 |
8/7 |
مقايسه طول گردابه پشت آبشكن بدست آمده از مدل عددي و اندازه گيريهاي آزمايشگاهي
(* حل عددي بدون اعمال عبارتهای آشفتگی، ** حل عددي با اعمال عبارتهای آشفتگی)
جريان چرخشي در محل انشعاب خروجي از كانال اصلي
چنين حالتي از جريان زماني حاصل ميشود که يک کانال از يک کانال اصلي منشعب ميشود. طول ناحيه چرخشي جريان با تغيير شرايط هيدروليکي جريان و تغيير هندسه کانال تغيير پيدا ميکند. در اينجا از معادلات آبهاي کم عمق که مناسب جريانهايي است که مقدار عمق در برابر عرض کانال ناچيز بوده و تغييرات تراز سطح آب نيز تدريجي ميباشد. براي صحت سنجي نتايج مدل عددي از نتايج آزمايشگاهي محققان پيشين استفاده ميشود. مدلسازي آشفتگي جريان با استفاده از مدل آشفتگي لزجت گردابه اي ثابت و زیرمقیاس شبکه تجربه می شود و بر طبق مشاهدات، با مقايسه نتايج حاصل از مدل عددي با اندازه گيريهاي آزمايشگاهي محققان پيشين، نتايج حاصل از مدل عددي با مدل آشفتگي اسماگورينسکي در مقايسه با نتايج حاصل از مدل عددي با مدل آشفتگي لزجت گردابه اي ثابت از صحت و درستي مطلوبي برخوردار است[18].
|
شکل شماتيک انشعاب خروجي از كانال اصلي و پارامترهاي جريان مورد استفاده در ارزيابي
خطوط جريان (راست) و خطوط هم عدد فرود (چپ) محاسبه شده نشانگر جريان با چرخه هاي پايدار در انشعاب خروجي از كانال اصلي
مقايسه پارامترهاي گردابه بدست آمده از مدل عددي و اندازه گيريهاي آزمايشگاهي
جريان چرخشي در جریان بعد از ورودي فرعي به كانال اصلی
در اين کار از نرم افزار تحلیلگر احجام محدود NASIR برای بررسی جريانهای آشفته زير بحراني کم عمق دو بعدی در كانالهاي با سطح آزاد دارای ورودی فرعی استفاده شده است. نرم افزار مذکور برای بررسی تاثیر تغییرات عدد فرود جریان خروجی از کانال اصلی بر کاهش عرض موثر جریان پایین دست محل تقاطع با جریان از یک ورودی فرعی بکار گرفته شده است. هردو كانال اصلی و فرعی با مقطع مستطيلي و بدون شيب هستند، نتايج بدست آمده با توجه به متغيرهاي جريان شامل تغييرات نسبت دبي دو شاخه اتصال و عدد فرود پايين دست، تأثير نسبتهاي عرضي مختلف دو شاخه اتصال بر هيدروليك و ويژگيهاي جريان مورد بررسي قرار گرفته است. نتايج مدلسازي عددي مسئله بر روی شبكه بندي مثلثي بدون ساختار ارائه و دقت نتايج حاصل از مدل رايانه اي به وسيله جداول و نمودارهايي براي مقايسه پارامترهاي عمق و سرعت با اندازه گيريهای آزمايشگاهی در دسترس آزمونهای مذکور مقايسه شده اند. نتايج تحقيق حاضر نشان می دهد که مي توان بجای مدلهای فيزيکی با اطمينان خوبی مدل رايانه ای مورد استفاده را برای کنترل طراحی پروژه هاي مهندسي اتصال کانالها بکاربرد. باکاربرد اين ابزار پيشگويی و مطالعه هيدروليکی مسائل واقعی مي تواند با كمترين هزينه انجام گردد و نتايج رضايت بخشي حاصل شود[19].
ناحيه اتصال دو شاخه كانال اصلي و فرعي
شبکه مورد استفاده (راست) و خطوط جریان (چپ) محاسبه شده نشانگر جريان با چرخه هاي پايدار در پایین دست انشعاب خروجي از كانال اصلي
تغييرات با مقادير در آزمون هاي آزمايشگاهيو عددي به ازاي اعداد فرود مختلف براي حالت
جريان سطح آزاد با رژیمهای مختلف در کانال با برآمدگی در کف
در این بخش نتایج حل احجام محدود جریان در یک کانال 25×1m با زبری ناچیز با برآمدگی در کف برای رژیمهای مختلف ارائه میشود. در آزمون اول جریان زیربحرانی عبوری از روی برآمدگی با منظور نمودن و عمق در پایین دست مدلسازی میشود. در آزمون دوم جریان انتقالی از زیربحرانی (قبل از برآمدگی) به فوق بحرانی (بعد از برآمدگی) با لحاظ کردن و عمق در پایین دست مدلسازی میشود. در آزمون سوم جریان زیر بحرانی با عبور از رژیم بحرانی در عبور از روی برآمدگی با منظور نمودن و عمق در پایین دست کانال مدلسازی میشود[24] .
معادله و شکل برآمدگی در کف در مقطع قائم
|
|
مقايسه تراز سطح آب بدست آمده از مدل عددي در محور مرکزي کانال با برآمدگی در کف با حل تحلیلی
(راست: جریان زیربحرانی، وسط: جریان انتقالی از زیربحرانی به فوق بحرانی، چپ: جریان زیر بحرانی با عبور از رژیم بحرانی)
جريان سطح آزاد در کانال با تنگ شدگی در محل برآمدگی کف
در اینجا جریان در یک کانال 25×1m با زبری ناچیز با تنگ شدگی عرضی در محل برآمدگی کف شبیه سازی میشود. در این آزمون جریان انتقالی از زیربحرانی (قبل از تنگ شدگی در محل برآمدگی) به فوق بحرانی (بعد از تنگ شدگی در محل برآمدگی) با منظور کردن و عمق در پایین دست مدلسازی میشود. نتایج حل عددی با استفاده از روش احجام محدود با حل تحلیلی و نتایج محاسباتی محققین قبلی مقایسه شده است[25] .
معادله عرض تنگ شدگی و برآمدگی در کف در مقطع قائم
مقايسه تراز سطح آب بدست آمده از مدل عددي در محور مرکزي کانال با تنگ شدگی در محل برآمدگی در کف با حل تحلیلی
نمای سه بعدی از سطح آب محاسبه شده روی بستر کانال با تنگ شدگی در محل برآمدگی در کف با حل تحلیلی
جريان باسطح آزاد غير چرخشي در اطراف پايه پلها
با استفاده از اين مدل ميتوان جريان اطراف پايه پلها را شبيهسازي نمود. اگر هندسه پايه پل منطبق با خطوط جريان طراحي شود، ميتوان جريان عاري از چرخابه ها را انتظار داشت. اما جريان در پايين دست پايه هاي داراي گوشه هاي تيز قطعا داراي چرخابهاي ناپايدار خواهد بود. در اين مدلسازي، محيط حل توسط يك شبكه مثلثي بيساختار گسسته شده است تا امكان مدلسازي مرزهاي پيچيده هندسي بسادگي فراهم گردد. در کار حاضر براي ارزيابي دقت مدل، ابتدا جريان بي استهلاک (با زبري ناچيز بستر) شبيه سازي شده است. براي غلبه بر نوسانات عددي پديد آمده در حل صريح معادلات در غياب عبارتهاي ميراگر نوسانات (عبارتهاي بيانگر تاثير اصطكاك، لزجت و آشفتگي) از روش لزجت مصنوعي استفاده شده است. روش استفاده شده براي اعمال لزجت مصنوعي مناسب شبكههاي بيساختار بوده و نوسانات عددي ناخواسته را بگونهاي از بين ميبرد كه صدمهاي به دقت نتايج حل وارد نگردد.
خطوط جريان روي شبكه بي ساختار مثلثي، جريان پايدار پيرو مرز (راست: پايه با انتهاي گرد، چپ: پايه با انتهاي دو قوسي)
تقارن نتايج محاسباتي در بالادست و پايين دست پايه پل منطبق با خطوط جريان عدم ورود لزجت غير واقعي را در مدل نشان ميدهد. در مرحله بعد به مدلسازي چرخابه ها در پايين دست پايه هاي تيز گوشه پرداخته ميشود [10].
خطوط جريان روي شبكه بي ساختار مثلثي، جريان با چرخه هاي ناپايدار (راست: پايه با انتهاي تيز، چپ: پايه با انتهاي صاف)
خطوط جريان روي شبكه بي ساختار مثلثي، جريان با چرخه هاي ناپايدار نوساني در پشت چند پايه با انتهاي تخت
جريان با سطح آزاد و آبشستگي در اطراف پايه هاي پل
با استفاده از اين مدل مي توان مدلسازي جريان اطراف انواع پايه ها را انجام و بدنبال آن نحوه شكلگيري آبشستگي در هر لحظه و عمق آبشستگي بیشینه را محاسبه نمود. اين مدل بوسيله حل معادلات جريان در اطراف پايه ها، مقادير سرعت در دو جهت افقی و عمق جريان را در هر گام زماني محاسبه مي کند. همزمان معادلات مدل آشفتگي استاندارد برای آبهای کم عمق نیز حل مي شوند. سپس با کمک انرژي آشفتگی مقدار سرعت برشي و تنش برشي برآورد مي شود. با استفاده از رابطه بين تنش برشي و نرخ آبشستگي، عمق آبشستگی در هر نقطه تعيين مي شود. اين روش براي نوع پايه قابل استفاده مي باشد. جهت صحت سنجي، از اندازه گیریهای آزمایشگاهی موجود برای آبشستگی در پایه های روی بستر دانه ای و رسی استفاده شده است[11].
نتايح مدل عددي و اندازه گیریهای آزمایشگاهی آبشستگی، بالا: در بستر دانه ای، پایین: در بستر رسی
[1] Numerical Analyzer for Scientific and Industrial Requirements (NASIR)
[2] صباغيزدي، سعيدرضا ،1381، تدوين مدل نيمه سهبعدي جريان درخليج فارس و پراكنش آلودگي در آن، گزارش طرح تحقيقاتي درون دانشگاهي شماره 1668، معاونت پژوهشي دانشگاه صنعتي خواجه نصيرالدين طوسي.
[3] Chorin, A. , 1967, A Numerical Method for Solving Incompressible Viscous Flow Problems, Journal of Computational Physics, Vol. 2, 12-26.
[4] Subgrid-scale
[5] صباغيزدي، سعيدرضا ، 01384 بررسی تاثير آشفتگی بر مدلسازی عددي ميدان فشار ناشی از جريان باد در اطراف ساختمانهای بزرگ ، گزارش طرح تحقيقاتي درون دانشکاهي شماره 1390، معاونت پژوهشي دانشگاه صنعتي خواجه نصيرالدين طوسي.
[6] صباغيزدي، سعيدرضا ،1383، تدوين مدل عددي سهبعدي توليد و انتقال حرارت هيدراتاسيون سيمان در مونوليت سد بتني، گزارش طرح تحقيقاتي درون دانشکاهي شماره 1846، معاونت پژوهشي دانشگاه صنعتي خواجه نصيرالدين طوسي.
[7] Sabbagh-Yazdi, Saeed-Reza, 1997, Simulation of the Incompressible Flow Using Artificial Compressibility Method, Ph.D. Thesis, University of Wales, Swansea.
[8] Sabbagh-Yazdi, R. S. & Mohammadzadeh-Qomi, M. , 2004, Finite Volume Solution of Two-Dimensional Convection Dominated Sub-Critical Flow Using Unstructured Triangular Meshes, International Journal of Civil Engineering, Vol. 2 , No. 3, pp 78-91.
[9] Sabbagh Yazdi, Saeed Reza & Mohamadzadeh Qomi, Mahyar, 2003, Using 2D Unstructured Mesh for Numerical Simulation of Free Surface Flow in Meandering Canal, International Symposium on Shallow Flows, Delft University of Technology, The Netherlands.
[10] سعيدرضا صباغ يزدي، 1380، مدلسازي جريان دو بعدي در اطراف پايه پل با شكل دلخواه بوسيله معادلات ميانگين عمقي ، سومين كنفرانس هيدروليك ايران، دانشگاه تهران.
[11] رضا دهقان نيّري، سعيدرضا صباغيزدي، " مدلسازي دوبعدي جريان و آبشستگي بستر رسي اطراف پايههاي دايره اي با کاربرد احجام محدود همپوشان بي ساختار" پذیرفته شده برای چاپ در مجله امیرکبیر (دانشگاه صنعتی امیرکبیر).
[12] Sabbagh Yazdi, Saeed Reza & Kermani, A., 2005, Hydrodynamic Modeling of Tidal Currents
in Persian Gulf Using Finite Volume Method, XXXI IAHR Congress, Seoul, Korea.
[13] سعيدرضا صباغ يزدي، 1379، حل عددي جريان آشفته دويعدي اطراف پايههاي پل با شكل دلخواه ، چهارمين كنفرانس بينالمللي سواحل، بنادر و سازههاي دريايي ICOMPAS2001، بندر عباس.
[14] Sabbagh Yazdi, Saeed Reza & SeediFar, B.., 2007, Analysis of Two-Dimensional Turbulence Flow in Channel with Sudden Expantion by Using Cell Centre Finite Volume Method, XXXII IAHR Congress, Veins Italy.
[15] سعيدرضا صباغ يزدي، و ذوقي، رضا ، 1385، شبيهسازي عددي جريان کمعمق گردابي در پاييندست آبشکن قائم با استفاده از ضريب زبري اصلاحشده مانينگ، هفتمين سمینار بین المللی مهندسی رودخانه، اهواز، ايران.
[16] صباغ يزدي ، سعيدرضا و ذوقي، رضا ، 1384، شبيه سازي عددي پديده جريان کم عمق چرخشي در پايين دست يک آبشکن، دومين كنگره ملي مهندسيعمران، دانشگاه علم و صنعت ايران.
[17] فرسیو، عطا و صباغ يزدي ، سعيدرضا ، 1385، بررسي صحت نتايج مدل حجم محدود جهتمند با الگوي مرکزيت سلول براي شبيه سازي جريان گردابه اي دو بعدي مجاور موج شکن منفرد ، هفتمين كنگره بین المللي مهندسيعمران، دانشگاه تربیت مدرس ایران.
[18] صباغ يزدي ، سعيدرضا و تقي زاده شياده، يوسف ، 1384، شبيه سازي عددي جريان گردابه اي با سطح آزاد در انشعاب قائم از کانال، دومين كنگره ملي مهندسيعمران، دانشگاه علم و صنعت ايران.
[19] صباغ يزدي ، سعيدرضا و بیرجندی، وحید 1385، شبيه سازي گردابه هاي كم عمق با استفاده از تغييرات عدد فرود خروجي در كانال با ورودي فرعي بوسيله مدل احجام محدود بي ساختار، هفتمين سمینار بین المللی مهندسی رودخانه، اهواز، ايران.
[20] Sabbagh-Yazdi, R. S, the paper is under preparation
[21] صباغ يزدي، سعيدرضا و محمدزاده، مهيار، 1383، شبيهسازي جريان دوبعدي در تندآب با ديوارهاي غيرموازي با روش احجام محدود بر روي شبكه بيسازمان ريز شده تطبيقي، نهمين کنفرانس ديناميک شاره های ايران، دانشکده مکانيک دانشگاه شيراز.
[22] Sabbagh-Yazdi, R. S. & AmiNejad, B. 2005, Simulation of High Velocity Flows in Chutes with Non-Parallel Sidewalls using Upwind and Central Over-Lapping Finite Volumes, 73rd Annual Meeting of ICOLD, Tehran, IRAN, Paper No.: 072-OT.
[23] Sabbagh-Yazdi,
R. S, Comparison of Central and Upwind Flux Averaging in Overlapping Finite
Volume Methods for Simulation of Super-Critical Flow with Shock Waves, The 9th WSEAS International Conference on
Applied Mathematics (MATH 06) , Istanbul, Turkey May 2006.
[24] Sabbagh-Yazdi, R. S. & SaeediFar, B. (2007), Accurate and Efficient Numerical Solution for Trans-Critical Steady Flow in a Channel with Variable Geometry, 2nd IASME/WSEAS International Conference on Continuum Mechanics, Portoroz (Portorose), Slovenia
[25] Sabbagh-Yazdi, R. S. & SaeediFar, B. (2007), Accurate and Efficient Numerical Solution for Trans-Critical Steady Flow in a Channel with Variable Geometry, 2nd IASME/WSEAS International Conference on Continuum Mechanics, Portoroz (Portorose), Slovenia
[26] Sabbagh-Yazdi, R. S. & Mohammadzadeh-Qomi, M. (2004), Finite Volume Solution of Mixed Sub & Super Critical 2D Free Surface Flow Using Unstructured Triangular Meshes, 9th International Conference on Hydro-Informatics (IAHR), Singapore.
[27] . Sabbagh-Yazdi, Saeed-Reza, (the paper is to be prepared)
[28] فرسیو، عطا ، 1384، تدوین مدل عددی به منظور بررسی تاثیر شکل دیوار هدایت جریان بالادست سرریز، بر جریان عبوری از سرریز، پایان نامه کارشناسی ارشد سازه های هیدرولیکی، دانشکده مهندسی عمران، دانشگاه صنعتی خواجه نصیرالدین طوسی
[29] صباغ يزدي، سعيدرضا، شار انتقالي جهتمند در مرزهاي احجام محدود بي ساختار با گام محاسباتي چند مرحله اي براي مدلسازي جريان در سرريز تندآب، ارائه شده براي چاپ در مجله انجمن هيدروليك (تحت بررسي).
[30] Sabbagh-Yazdi, R. S. & Zounemat Kermani, M. (2007), Velocity Profile over Spillway by Finite Volume Solution of Slopping Depth Averaged Flow, 2nd IASME/WSEAS International Conference on Continuum Mechanics, Portoroz (Portorose), Slovenia.
[31] Sabbagh-Yazdi, R. S. & AminNejad, B., 2005, Finite Volume Solution of Mixed Sub & Super Critical 2D Free Surface Flow Using Unstructured Triangular Meshes, XXXI IAHR Congress, Seoul, Korea.
[32]
Sabbagh-Yazdi, R. S. (2007), Spillway Flow Modeling by Finite Volume Solution
of Slopping Depth Averaged Equations on Triangular Mesh; Application to
KAROUN-4, 10th WSEAS International Conference on Applied Mathematics, Dallas
(Texas), USA.
[33] صباغ يزدي، سعيدرضا ، 1386، شبيهسازي جريان فوقبحراني برروي سرريز تندآب با منظور نمودن پروفيل عمودي سرعت و ورود هوا از سطح آب، طرح تحقیقاتی شماره Dam2 81170 معاونت پژوهشی وزارت نیرو.
[34] Sabbagh-Yazdi, R. S. & Safaeih, R., B. 2007, Post-Processing of Air Entrainment on NASIR Flow Solver Results for Skimming Flow over Stepped Chutes, 9th WSEAS International Conference onAutomatic Control, Modeling and Simulation, Istanbul, Turkey.
[35] صباغ يزدي، سعيدرضا، صفائیه ، رزیتا، مقاله در دست تهیه است.
[36] فرسیو، عطا، فرسیو، عطا ، 1384، تدوین مدل عددی به منظور بررسی تاثیر شکل دیوار هدایت جریان بالادست سرریز، بر جریان عبوری از سرریز، پایان نامه کارشناسی ارشد سازه های هیدرولیکی، دانشکده مهندسی عمران، دانشگاه صنعتی خواجه نصیرالدین طوسی
[37] Sabbagh-Yazdi, Saeed-Reza and SaeediFar, Behzad, , 2007, A Cell Center Finite Volume Model for 2D Numerical Simulation of Deposition-Erosion and Transport of Suspended Sediment in Free Surface Flows, 9th WSEAS Int.Conf. on Automatic Control, Modeling & Simulation , Istanbul, Turkey
[38] Sabbagh-Yazdi, Saeed-Reza and SaeediFar, Behzad, , 2007, A Cell Center Finite Volume Model for 2D Numerical Simulation of Deposition-Erosion and Transport of Suspended Sediment in Free Surface Flows, 9th WSEAS Int.Conf. on Automatic Control, Modeling & Simulation , Istanbul, Turkey
[39] Sabbagh-Yazdi, Saeed-Reza and SaeediFar, Behzad, , 2007, Simulation of Suspended Sediment Transport in Dam Reservoir Using a Cell Center Finite Volume Solver, 11th WSEAS Int. Conf. on COMPUTERS , Vancouver, Canada
[40] صباغ يزدي، سعيدرضا، و عربي، حامد، اعمال اثر تروخشک شدن سواحل در مدلسازي احجام محدود جريانات دو بعدي ناشي از باد در درياچه ها، مجله علوم و تكنولوژي محيط زيست شماره 26، صفحه 37-46.
[41] صباغ يزدي، سعيدرضا، و عربي، حامد، اعمال اثر تر و خشک شدن سواحل در مدلسازي احجام محدود جريانات دو بعدي ناشي از باد در درياچه ها، مجله علوم و تكنولوژي محيط زيست شماره 26، صفحه 37-46.
[42] صباغ يزدي، سعيدرضا و صادقي گويا، محمد ، 1383، مدلسازي عددي دقيق جريان جزرومدي انتقالي دو بعدي در خورها ، اولين كنگره ملي مهندسيعمران، دانشگاه صنعتي شريف.
[43] صباغ يزدي، سعيدرضا و کرماني، علي، 1383، كاربرد لزجت مصنوعي مناسب شبكه بيساختار در حل عددي جريان دوبعدي غيرماندگار روي بستر با شيب متغير و اصطكاك ناچيز، اولين كنگره ملي مهندسيعمران، دانشگاه صنعتي شريف.
[44] صباغ يزدي، سعيدرضا، و عربي، حامد، بهمن ماه 1385 آزمون مدل احجام محدود ميانگين عمقی جهت شبيه سازی امواج عمود بر نوار ساحلی، هفتمين سمينار بين المللي مهندسي رودخانه اهواز. اهواز . دانشگاه شهيد چمران.
[45] صباغ یزدی، سعیدرضا و صادقی گویا، ارسلان، 1383، مدلسازي عددي دقيق جريان جزرومدي انتقالي دو بعدي در خورها، اولین کنگره ملی مهندسی عمران، دانشگاه صنعتی شریف.
[46] Sabbagh-Yazdi, R. S. & Zounemat Kermani, 2007, Computer Simulation of Tidal Currents in a Marine Canal, Case Study of QESHM Canal in HOURMOZ Strait, Iran, , XXXII IAHR Congress, Veins Italy.
[47] Sabbagh-Yazdi, R. S., Zounemat Kermani, M. & Kermani, A., Dec. 2006, Solution of Depth Averaged Tidal Currents in Persian Gulf on Unstructured Overlapping Finite Volumes. International Journal of Numerical Methods in Fluids.
[48] Sabbagh-Yazdi, S.R., 2004, Numerical Modeling of Persian Gulf Salinity Variations Due to Tidal Effects, Journal of Environmental Science & Technology. Vol. 1, No. 1, 41-48.
[49] Sabbagh-Yazdi, S.R., 2005, Coupled Solution of Oil Slick and Depth Averaged Tidal Currents on Three-Dimensional Geometry of the Persian Gulf, Journal of Environmental Science & Technology, Vol. 2 , No. 4, pp. 309-317.
[50] S.R Sabbagh-Yazdi, “Investigation of Coriolis Forces on the Caspian Sea Currents Using a Two Dimensional Numerical Model “, 4th International of Coasts, Ports and Marin Structures (ICOMPAS2002) ,Ramsar, Iran.
[51] صباغ يزدي، سعيدرضا، پاييز 1379، مدل عددي دوبعدي مقدماتي تحليل پراکنش آلودگي در گستره درياي خزر با در نظر گرفتن اثر باد، مجله علوم و تكنولوژي محيط زيست شماره 6و5، ص 35-52.
[52] Hamidi M and Sabbagh-Yazdi S. R., (2006) " Numerical Modeling of Seawater Intrusion in Coastal Aquifer Using Finite Volume Unstructured Mesh Method ", WSEAS Transactions on Mathematics, Vol. 6, No. 5, pp. 648-655.
[53] Hamidi M & Sabbagh-Yazdi S. R, (2006), " Numerical Study on Effect of Subsurface Dam on Controlling Seawater Intrusion in Confined Coastal Aquifer ", 7th International Civil Engineering Conference TM University -Iran, on CD.
[54] صباغ یزدی، سعید رضا و بیات، بابک ، 1385، تحلیل تراوش دوبعدی در پي و محاسبه زیرفشار در سدهاي بتني وزني، هفتمین کنگره بین المللی مهندسی عمران، دانشگاه تربیت مدرس.
[56] عباسی، علی و صباغ يزدي،
سعيدرضا، 1385، تحليل پديده ضربه قوچ در تونلهای آب بر نيروگاهها با روش حجم محدود
جهتمند، هفتمین کنگره بین المللی مهندسی عمران، دانشگاه تربیت
مدرس.
[57] صباغ يزدي، سعيدرضا، 1383 ، حل عددي دقيق جريان دوبعدي تراكمناپذير غيرلزج دائمي با روش احجام محدود بر روي شبكه مثلثي، مجله بينالمللي علوم مهندسي (دانشگاه علم و صنعت ايران)، جلد 15، شماره 2، ص 31-43.
[58] صباغ يزدي، سعيدرضا، 1381 ، استفاده از استهلاك مصنوعي براي شبكه اجزاء بدونساختار در حل عددي معادلات دائمي غير لزج تراكمناپذير ، چهارمينهمايشانجمنهوافضاي ايران، دانشگاه صنعتي اميرکبير.
[59] صباغ يزدي، سعيدرضا، 1383 ، حل عددي دقيق جريان دوبعدي تراكمناپذير غيرلزج دائمي با روش احجام محدود بر روي شبكه مثلثي، مجله بينالمللي علوم مهندسي (دانشگاه علم و صنعت ايران)، جلد 15، شماره 2، ص 31-43.
[60] Sabbagh Yazdi, Saeed Reza & Hadian, Abbas., 2003, Simulation of Air Pollution Transport Around Tall Buildings Using a Two-Dimensional Numerical Flow Solver, 4th National Aero Space Conference, Isfehan University of Technology, Isfehan, Iran.
[61] Sabbagh-Yazdi S. R. and Hadian A., (2002) "Two-Dimensional Numerical Simulation of Wind Effect on Pollutant Transport around Urban Buildings", ", Journal of Environmental Science and Technology Vol.13, pp. 33-46
[62] Sabbagh-Yazdi S. R. and Hadian A., 2003, Numerical Simulation of 2D Air Pollution Transport over Walls, 4th Conference of Iranian Aerospace Society, Amir-Kabir University of Technology, Tehran-Iran, pp. 208-212
[63] Sabbagh-Yazdi S. R, Torbati M, MeysamiAzad F, and Haghighi B, 2007, Computer Simulation of Changes in the Wind Pressure Due to Cooling Towers-Buildings Interference, WSEAS Transactions on Mathematics, Vol. 6, No.1, pp.205-214 .
[64] Sabbagh Yazdi, Saeed Reza, 1997, Simulation of the Incompressible Flow Using Artificial Compressibility Method, Ph.D. Thesis, University of Wales, Swansea.
[65] Sabbagh Yazdi, Saeed Reza, 1997, Simulation of the Incompressible Flow Using Artificial Compressibility Method, Ph.D. Thesis, University of Wales, Swansea.
[66] Sabbagh Yazdi, Saeed Reza, 1997, Simulation of the Incompressible Flow Using Artificial Compressibility Method, Ph.D. Thesis, University of Wales, Swansea.
[67] Sabbagh-Yazdi S.R, Haghighi B. and Torbati M., 2006, Application of NASIR[67] Flow-Solver for Finite Volume Solution of Wind Effects on a Set of Cooling Towers, 10th WSEAS International Conference on Applied Mathematics (MATH 06), Dallas(Texas)-USA.
.
[68] Sabbagh-Yazdi S. R, Torbati M, MeysamiAzad F, and Haghighi B, 2007, Computer Simulation of Changes in the Wind Pressure Due to Cooling Towers-Buildings Interference, WSEAS Transactions on Mathematics, Vol. 6, No.1, pp.205-214 .
[69] Sabbagh Yazdi, Saeed Reza & Hadian, Abbas., 2004, Pressure-Velocity Coupled Finite Volume Solution of Steady Incompressible Flow Using Artificial Compressibility Technique, International Journal of Engineering, Transaction A: Basics, Vol. 17, No. 2, pp 109-118.
[70] صباغ يزدي، سعيدرضا، 1379، شيوهاي بر مبناي ضلع اجزاء شبكه بدون ساختار براي حل عددي جريان غير چسبنده تراكم ناپذير، اولين كنفرانس بينالمللي هوا فضاي ايران ( (Aero2001، دانشگاه صنعتي شريف.
[71] صباغ يزدي، سعيدرضا و صالحيگورابسري، وحيد، 1382، كاربرد مدل تحليل عددي توليد و توزيع دما در سازههاي بتني، ششمين كنفرانس مهندسي عمران، دانشگاه صنعتي اصفهان.
[72] صباغ يزدي، سعيدرضا و دولوو، كاميار ، 1383، بررسي شيوههاي مختلف احجام محدود براي شبكههاي بيساختار در حل عددي توليد و انتتقال حرارت، اولين كنگره ملي مهندسيعمران، دانشگاه صنعتي شريف.
[73] Sabbagh Yazdi, S.R. & Bagheri, A.R., 2004, Computer Simulation of Cement Heat Generation and Temperature Profiles in Mass Concrete Structures, International Journal of Engineering Science, Vol. 15, No. 2, pp. 65-71.
[74] صباغ يزدي، سعيدرضا ، 1384، تدوين مدل عددي سهبعدي توليد و انتقال حرارت هيدراتاسيون سيمان در مونوليت سد بتني، گزارش طرح تحقيقاتي درون دانشکاهي، معاونت پژوهشي دانشگاه صنعتي خواجه نصيرالدين طوسي.
[75] صباغ يزدي، سعيدرضا و محمودي خالد آبادي، محمد ، 1384، ارزيابي مدل حل عددي سه بعدي پديده توليد و انتقال حرارت در بتن، دومين كنگره ملي مهندسيعمران، دانشگاه علم و صنعت ايران.
[76] صباغ يزدي، سعيدرضا و باقري، عليرضا، 1379، توسعه و كاربرد نرمافزار تحليل توزيع حرارت ناشي از هيدراسيون سيمان در ساخت لايهاي سدهاي RCC، چهارمين كنفرانس سد سازي (IRCOLD2001 ).
[77] Sabbagh-Yazdi S. R. and Bagheri, A., 2001, Thermal Numerical Simulation of the Laminar Construction of RCC Dams, 6th Inter. Conf. on Computational Modeling of Moving Boundary Problems, Lemnos, Greece.
[78] صباغ يزدي، سعيدرضا ، 1384، تدوين مدل عددي سهبعدي توليد و انتقال حرارت هيدراتاسيون سيمان در مونوليت سد بتني، گزارش طرح تحقيقاتي درون دانشکاهي، معاونت پژوهشي دانشگاه صنعتي خواجه نصيرالدين طوسي.
[79] Sabbagh-Yazdi S. R. and Davallo k., 2005, Two Dimensional Numerical Simulation of Temperature Field During the Construction of the Gravity Dams with Cooling Pipe, to be Appeared in Proceeding of 7th Inter. Conf. of Civil Engineering, Tarbiat Modares University, Tehran, Iran.
[80] Numerical Analyzer for Scientific and Industrial Requirements (NASIR)
معرفی نرم افزار حجم محدود تحلیل مسائل علمی و صنعتی NASIR
1. تحلیل مسائل علمی و صنعتی در مدلسازی جریان میانگین عمقی با سطح آزاد (در کانالها و رودخانه ها)
2. تحلیل مسائل علمی و صنعتی در مدلسازی جریان میانگین عمقی با سطح آزاد (فرسایش، انتقال و ته نشست رسوب)
3. تحلیل مسائل علمی و صنعتی در مدلسازی جریان غیردائمی در لوله ها
4. تحلیل مسائل علمی و صنعتی در مدلسازی تولید و انتقال حرارت در بتن
5. تحلیل مسائل علمی و صنعتی در مدلسازی جریان در محیط متخلخل
7. تحلیل مسائل علمی و صنعتی در مدلسازی جریان دو بعدی باد و انتقال آلایندها