مدلسازی CFD حالت پایدار و تجزیه و تحلیل تجربی ریزاقلیم محلی در دبی - ساختمان سبز ایرانیان

  • Home
  • مدلسازی CFD حالت پایدار و تجزیه و تحلیل تجربی ریزاقلیم محلی در دبی

© SF Fatima و HN Chaudhry، منتشر شده توسط EDP Sciences، 2017

مجوز Creative Commonsاین یک مقاله با دسترسی آزاد است که تحت شرایط Creative Commons Attribution License ( http://creativecommons.org/licenses/by/4.0 ) توزیع شده است، که استفاده، توزیع و بازتولید نامحدود در هر رسانه را مجاز می‌کند، مشروط بر اینکه کار اصلی باشد. به درستی ذکر شده است.

1. معرفی

در چند سال گذشته، جهان به سرعت در حال توسعه به سمت شهرنشینی بوده است که بر محیط زیست تأثیر منفی می گذارد و منجر به افزایش قابل توجهی در نیازهای سرمایش می شود. یکی از پیامدهای این پدیده اثر جزیره گرمایی است. افزایش دمای سطحی سازه های شهری منجر به تغییر ریزاقلیم شهری در نتیجه تغییرات در تبادل حرارتی همرفتی از این سطوح می شود.1 ]. این منجر به افزایش دمای هوا می‌شود، جایی که کلان‌شهرها با ساختارهای سخت ساخته شده دارای میانگین دمای هوای بالاتری نسبت به مناطق روستایی اطراف آن هستند.

باشه [2 ] پدیده جزیره گرمایی را به عنوان گرمای به دام افتاده در مناطق شهری متراکم که منجر به افزایش دمای هوای محیط در مرکز شهر در مقایسه با مناطق دورافتاده اطراف می‌شود، به تصویر کشید. تام و همکاران [3 ] از آن به عنوان یک پدیده انسان زا به دلیل تغییرات سطوح شهری توسط سکونت انسان یاد کرد که منجر به افزایش دمای هوا در یک مرکز شهری متراکم می شود. افزایش دمای ناشی از این پدیده منجر به تقاضای انرژی بیشتر برای سرمایش، افزایش غلظت آلودگی هوا و گازهای گلخانه ای و کاهش آسایش حرارتی می شود.

شرایط جوی مانند سرعت و جهت باد یکی از عوامل اصلی تاثیرگذار بر شدت جزایر گرمایی در محیط ساخته شده است.4 ]. سرعت باد بالا امکان کاهش پدیده را فراهم می‌کند در حالی که شرایط راکد آن را افزایش می‌دهد، جایی که الگوی جریان باد و سرعت تحت تأثیر پیکربندی شهری است که به نوبه خود بر شدت تبادل حرارت تأثیر می‌گذارد.5 ]. آلگرینی و همکاران [6 ] افزایش انتقال حرارت همرفتی با نمای ساختمان و حذف گرما به دلیل افزایش جریان باد در دره خیابان مشاهده شد. رنگ قهوه ای [7 ] نشان داد که افزایش دما با توجه به سرعت باد متفاوت است، جایی که شدت آن در شرایط باد آرام قابل توجه است و در شرایط باد شدید کم است. فیشر و همکاران [8 ] هوای آرام را زمانی که سرعت باد کمتر از 3 متر بر ثانیه است طبقه بندی کرد. از آنجایی که سرعت باد بیش از 3 متر بر ثانیه است، هوای گرم را از شهر خارج می کند و شدت پدیده جزیره گرمایی را کاهش می دهد.

این مطالعه یک ارزیابی دمای حالت پایدار را در زمینه الگوی جریان باد، سرعت میکرو اقلیم، هندسه ساختمان و دمای سطح شهری با استفاده از پردیس دانشگاه هریوت وات دبی (HWUDC) به عنوان مطالعه موردی انجام می‌دهد. پارامترهای حرارتی و آیرودینامیکی مؤثر بر جزایر گرمایی با استفاده از دینامیک سیالات محاسباتی شبیه‌سازی عددی شدند که نتایج آن از طریق اندازه‌گیری‌های میدان فیزیکی تأیید شد. مکان‌های کانونی خاص در اطراف زیرساخت ساختمان در پاسخ به آب و هوای محلی خارجی برای تعیین کمیت بارهای خنک‌کننده به دلیل افزایش حرارت رسانای محسوس از طریق خارجی تعیین شدند.

2 کارهای مرتبط قبلی

مطالات گذشته [9 –12 ] تاثیر دمای سطح حالت پایدار و جزایر گرمایی را به عنوان تابعی از هندسه/جهت شهری، ماهیت سطوح شهری، شرایط هواشناسی، حرکت هوا و سرعت هوا بررسی کرده‌اند. شدت اثر جزیره گرمایی را می توان با اصلاح این پارامترها به صورت استراتژیک کاهش داد. روش‌های ارزیابی پیاده‌سازی‌شده در ادبیات، مبنایی اساسی در رویکرد ارزیابی مشخصات دما و بارهای خنک‌کننده در محیط ساخته‌شده واقع در آب‌وهوای گرم و خشک فراهم می‌کند. مطالعه انجام شده توسط راضی [13 ] نشان داده است که محیط های حرارتی تحت تأثیر محیط ساخته شده قرار می گیرند و افزایش دما می تواند به خوبی بین 2 تا 4 درجه سانتیگراد باشد. استفاده از ابزارهای عددی مانند CFD در پیش‌بینی توزیع جریان و حرارت پایدار و دینامیکی در اطراف ساختمان‌ها نیز به طور گسترده در ادبیات قبلی مورد بررسی قرار گرفته است.14 –16 ] همانطور که در بخش های بعدی توضیح داده شده است.

2.1 تأثیر، علل و استراتژی های کاهش جزایر گرمایی حالت پایدار

ابید و همکاران [ 9 ] تأثیر تابش خورشیدی و پراکندگی آلودگی ریزاقلیم شهری در شهر مصدر، ابوظبی را در سه زمان از روز بررسی کرد. شبیه‌سازی‌های CFD در ANSYS FLUENT با استفاده از مدل k – e با معادلات ناویر-استوکس ثابت انجام شد . نتایج شبیه‌سازی دمای بالاتر هوا را در پنجره‌های رو به خورشید نشان می‌دهد. مشخص شد که افزایش دمای هوای محیط به دلیل دمای بالاتر سطح زمین و سقف است که از طریق تابش خورشیدی گرم می شود. کار به این نتیجه رسید که نسبت ابعاد بالاتر از گرم شدن تابش مستقیم خورشیدی سطوح شهری و کاهش اثرات منفی بر اقلیم خرد شهری جلوگیری می کند.

چو و همکاران [10 ] برای بررسی پدیده جزیره گرمایی، محیط حرارتی را در مقیاس وسیع در یک منطقه اطراف (5 کیلومتر مربع و 10 کیلومتر مربع ) از خلیج ساحلی توکیو تجزیه و تحلیل کرد. شبیه‌سازی‌های CFD با استفاده از مدل آشفتگی k – ε انجام شد  و در برابر اندازه‌گیری‌های مشاهده‌ای اعتبارسنجی شد. نتایج نشان داد که دمای هوا در فاصله 2 متری از سطح زمین به میزان 1 تا 2 درجه سانتی گراد در بالای مناطق و رودخانه های بزرگ رویشی کاهش یافته است. این مطالعه نشان داد که تکانه رو به بالا (محرک شناوری) سرعت هوا در منطقه ای با دمای هوا بالاتر رخ می دهد. کار تحقیقاتی به این نتیجه رسید که وجود جفت حرکت به سمت بالا و پایین سرعت هوا بر تغییر دمای هوای محیط تأثیر می‌گذارد.

طالب و ابوحجله [11 ] تأثیرات پیکربندی های مختلف شهری در دبی را بر تغییرات دما مورد مطالعه قرار داد. 27 شبیه‌سازی با ترکیب‌های مختلف چیدمان شهری، سرعت باد (0.1 متر بر ثانیه برای سناریوی بدون باد، 3.6 متر بر ثانیه به عنوان سرعت اولیه باد و 7 متر بر ثانیه برای تابستان تنها برای درک تأثیر افزایش سرعت باد بر ثانیه انجام شد. تغییرات دما) و زمان سال (ژوئن، دسامبر و سپتامبر) با ورودی دمای اولیه ثابت 32 درجه سانتیگراد. نتایج نشان داد که پیکربندی Bastakiya ( شکل 1) در تابستان دمای کمتری داشت. در مقایسه نتایج با حضور و عدم حضور باد، کاهش نوسانات دما در حضور باد وجود داشت که منجر به کاهش وجود نقاط داغ شد. این مطالعه به این نتیجه رسید که پیکربندی بستکیا از نظر موقعیت مکانی و اتحاد بادهای غالب به موازات جاده، بهترین عملکرد حرارتی تحت هر شرایطی است.

پریادرشینی و همکاران [ 12 ] عوامل اساسی ایجاد پدیده جزیره گرمایی در سنگاپور را برجسته کرد. شبیه‌سازی‌های CFD از  معادله مدل آشفتگی k – ε تحت تکنیک حجم محدود استفاده کردند. شرایط مرزی از داده های هواشناسی هواشناسی و اندازه گیری های میدانی به دست آمد. آزمایش از طریق اندازه گیری میدان فیزیکی تایید شد. مشخص شد که دما با افزایش سرعت هوای متناظر تا 35 درصد به دلیل وجود برج‌های بلند با نسبت H/W بهینه، 0.7 درجه سانتی‌گراد کاهش یافت. این مطالعه به این نتیجه رسید که با مکان‌یابی استراتژیک برج‌های مرتفع با نسبت H/W بهینه، جریان هوا را می‌توان در دره‌های خیابانی افزایش داد و در نتیجه شدت جزیره گرمایی را کاهش داد.

بند انگشتی عکس. 1Bastakiya (چپ)، متعامد (مرکز)، حجم عمودی (راست) مدل سازی شده در ENVI-met [ 11 ].

2.2 کاربرد CFD در ارزیابی پروفیل های حرارتی و جریان

آلگرینی و همکاران [ 6 ] تأثیر ریزاقلیم شهری را بر تقاضای انرژی ساختمان در یک محله شهری بررسی کرد. حوزه محاسباتی متشکل از 14 ساختمان است که در نسبت‌های ابعادی مختلف چیده شده‌اند که با استفاده از معادلات سه بعدی رینولدز میانگین ناویر-استوکس (RANS) و  k – ε شبیه‌سازی شده‌اند.مدل آشفتگی در کد عددی OpenFOAM. در مجموع چهار شرایط آب و هوایی مختلف با سرعت باد از 0.3 متر بر ثانیه تا 4.7 متر بر ثانیه مدل‌سازی شد. یافته‌ها نشان داد که تفاوت در دمای هوا در مورد جریان نزدیک و هوای داخل یک دره شهری یا شدت جزیره گرمایی محلی می‌تواند تا 2.5 درجه سانتی‌گراد باشد که باعث افزایش قابل توجهی در تقاضای خنک‌کننده می‌شود. این مطالعه نتیجه گرفت که رویکرد CFD ابزار مفیدی برای پیش‌بینی شدت جزیره گرمایی برای مناطق شهری است.

توپارلار و همکاران [ 14 ] تحلیلی از ریزاقلیم شهری در شهر روتردام، برگپولدر زوید با انجام شبیه‌سازی‌های CFD با استفاده از معادله RANS ناپایدار سه بعدی تحت  مدل آشفتگی k – ε بر روی یک شبکه با وضوح بالا ارائه کرد. نتایج نشان می‌دهد که در شرایط باد آرام، دمای سطح عناصر شهری نسبتاً بالا بوده است. خیابان های عمود بر سرعت باد اثر خنک کننده تهویه را کاهش داده است. CFD دمای سطح شهری مشابهی را در مقایسه با تصاویر حرارتی ماهواره‌ای با انحراف 7.9 درصد تولید می‌کند. این مطالعه به این نتیجه رسید که CFD می تواند به عنوان یک ابزار عددی برای ارزیابی میکرو اقلیم شهری و شناسایی مناطق مشکل ساز استفاده شود.

آهنگ و لیو [15 ] توزیع رطوبت حرارتی و سرعت را نسبت به اقلیم ریز شهری در حضور تابش خورشیدی و انتقال رطوبت ارزیابی کرد. شرایط مرزی از داده های به دست آمده از طریق مرکز هواشناسی تعیین شد. معادله اغتشاش استاندارد  k – ε به عنوان مدل آشفتگی استفاده شد. نتایج نشان می دهد که پیکربندی ساختمان به طور قابل توجهی دمای هوا و رطوبت را تغییر داده است. با افزایش ارتفاع ساختمان، اثر خنک کنندگی ویژگی آب به دلیل کاهش قرار گرفتن در معرض خورشید مسدود شده توسط ساختمان های بلند کاهش یافت. این تحقیق نشان داد که شبیه‌سازی‌های CFD می‌تواند به عنوان یک ابزار موثر برای پیش‌بینی محیط حرارتی یک میکرو اقلیم شهری مورد استفاده قرار گیرد.

ستایجه و همکاران [28 ] از ANSYS FLUENT در یک مطالعه موردی در شهر مدینه، عربستان سعودی برای تجزیه و تحلیل آسایش حرارتی در خیابان‌های شهری استفاده کرد. در شبیه سازی از معادله آشفتگی RANS با  مدل k – ε و تابع دیوار استاندارد استفاده شد. نتایج شبیه‌سازی با اندازه‌گیری‌های میدانی مقایسه شد، که در آن سرعت هوای به‌دست‌آمده از شبیه‌سازی 6.5 متر بر ثانیه و 7 متر بر ثانیه از اندازه‌گیری‌های میدانی بود که نشان‌دهنده خطای کمتر از 20 درصد بود. این مطالعه بیان کرد که با توجه به Wilkinson و همکاران. [17 ]، خطای قابل قبول برای نتایج CFD تا 20 درصد است، بنابراین شبیه سازی انجام شده اعتبار می یابد. نتایج شتابی را در سرعت باد 1.3 متر بر ثانیه با ساختمان‌های بلندتر نشان می‌دهند که حرکت باد را در امتداد خیابان‌ها افزایش می‌دهند و باعث کاهش دمای هوا به میزان 2 درجه سانتی‌گراد می‌شوند. این مطالعه به این نتیجه رسید که هندسه و جهت گیری شهری بر سرعت و دمای هوا تأثیر می گذارد که بر آسایش حرارتی در امتداد خیابان ها تأثیر می گذارد. این مطالعه همچنین پیشنهاد کرد که CFD می تواند برای ارزیابی مسائل پیچیده ریزاقلیمی مانند آسایش حرارتی، انتشار آلودگی، کیفیت الگوهای جریان هوا و باد همراه با اثرات پوشش گیاهی و رطوبت نسبی مورد استفاده قرار گیرد.

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

3 حوزه فیزیکی

پردیس دانشگاه هریوت وات دبی به عنوان مطالعه موردی برای تحلیل انتقادی وقوع پارامترهای حرارتی و جریان انتخاب شد. پردیس شامل ساختمان های فاز 1 و فاز 2 با مساحت کلی 28000 متر مربع است که در شهر آکادمیک بین المللی دبی واقع شده است. ساختمان دارای هندسه متقارن با دیوارهای بیرونی برافراشته و فرورفته به شکل یکنواخت است. فاز 1 ساختمان 51 متر عرض دارد و در طول 132 متر گسترش یافته و ارتفاع آن به 20 متر می رسد. فاز 2 ساختمان شامل سه بلوک متشکل از پنج طبقه با طول، عرض و ارتفاع به ترتیب 32 متر، 12 متر، 16 متر است. یک بلوک برای سالن با ابعاد 25 متر × 35 متر × 5 متر. این دو فاز در فاصله 15 متری از بلوک سالن مرکزی قرار دارند. شکل 2HWUDC را به صورت سه بعدی در مقیاس 1:100 نشان می دهد.

نسبت پنجره به دیوار (WWR) در هندسه 37٪ بود. مناطق لعاب با WWR بالا عمدتاً برای اهداف ساده‌سازی باعث افزایش قابل توجه دما می‌شوند. به منظور تنظیم شرایط مرزی ورودی، یک محفظه 212 × 152 × 40 متر در اطراف هندسه ایجاد شد به طوری که 2.5 برابر طول حوزه فیزیکی در جهت x، y و z بود تا از اثر برگشتی جلوگیری شود.

این مطالعه فعالیت تبادل حرارتی را در میان دو فاز ساختمانی (فضای گردش انتقالی) که سطح نسبتاً بالایی از حرکت عابر پیاده دارد، ارزیابی کرد. دمای بالای ناشی از این پدیده منجر به افزایش ناراحتی حرارتی می‌شود، بنابراین شناسایی فعالیت در این منطقه بسیار مهم تلقی می‌شود. ساختمان برای اهداف تصویری به سه منطقه تقسیم شده است که در شکل 3 نشان داده شده است .

بند انگشتی شکل 2(الف) پردیس دانشگاه هریوت وات دبی، (ب) نسخه مدل شده.

بند انگشتی شکل 3(الف) تقسیم مناطق و فازهای طرح HWUDC مدل شده، (ب) سه منطقه ایجاد شده از جهت شمال.

3.1 راه اندازی آزمایشی

اندازه‌گیری‌های میدانی در اطراف HWUDC مورد نیاز برای شرایط مرزی ورودی و اعتبارسنجی تجربی برای نشان دادن نمایشی نزدیک‌تر از سناریوی مورد واقعی انجام شد. اندازه‌گیری‌ها در 2 فوریه 2015 در ساعت 16:00 ثبت شد. از آنجایی که باد از جهت شمال غربی (NW) غالب بود، پارامترهای سرعت و دما در سمت شمال غربی محفظه تقارن به عنوان شرایط مرزی ورودی تنظیم شدند. به منظور ارزیابی اثر دمای بالای نمای پردیس، دمای سطح شیشه و دیوار خارجی نیز به عنوان شرایط مرزی وارد شده است. تنش سنج حرارتی Extech (HT30) با دقت 1±٪ با وضوح 0.1 درجه سانتیگراد برای ثبت دمای حباب خشک استفاده شد. دماسنج مادون قرمز Meterman IR 608 با دقت 2% ± با وضوح 0.2 درجه سانتی گراد و بادسنج سیم داغ (C.

شکل 4a شرایط مرزی را نشان می دهد که بر روی شیشه ساختمان ها مشخص شده است. پنجره های رو به جنوب با هم گروه بندی شدند و در دمای سطح 42 درجه سانتی گراد به دست آمده از اندازه گیری های میدانی تنظیم شدند. شکل 4b تصویر واقعی موقعیت ورودی را نشان می دهد.

به منظور اعتبارسنجی نتایج تولید شده از مدل CFD، 10 نقطه اندازه گیری گسسته در امتداد منطقه گردش ایجاد شد که دو فاز را در سه منطقه به هم متصل می کند. شکل 5 نقاطی را نشان می دهد که پارامترهای مورد بررسی در آنها اندازه گیری شده اند. دمای حباب خشک و سرعت هوا ثبت‌شده در نقاط تعیین‌شده به‌عنوان پارامترهای اندازه‌گیری برای تأیید یافته‌های CFD استفاده شد.

بند انگشتی شکل 4(الف) شرایط مرزی لعاب، (ب) دمای سطح، دمای محیط و سرعت اندازه‌گیری شده در سمت جنوب برای شرایط مرزی ورودی.

بند انگشتی شکل 5نقاط نشان دهنده اندازه گیری میدان فیزیکی انجام شده در HWUDC (ابعاد بر حسب متر).

3.2 کاهش داده ها

قرائت‌های به‌دست‌آمده از اندازه‌گیری‌های میدان فیزیکی برای محاسبه افزایش حرارت محسوس ناشی از انتقال حرارت رسانا از طریق شیشه به منظور تجزیه و تحلیل تأثیر بر بار خنک‌کننده ساختمان اعمال شد. ضریب انتقال حرارت برای شیشه های دوجداره ساختمان 2.26 W/m 2  K [.18 ]. مساحت سطح لعاب 87.31 متر مربع برآورد شد . تفاوت دما را می توان از طریق داده های اندازه گیری شده یا نتایج محاسباتی به دست آورد. این متغیرها برای بدست آوردن درجه انتقال حرارت همرفتی ناشی از لعاب در رابطه (1) ادغام شدند .(1)که در آن لعاب Q نشان دهنده بار محسوس ناشی از تبادل حرارتی همرفتی از طریق لعاب (W)، U نشان دهنده عبور حرارتی لعاب (W/m 2  K)، A نشان دهنده سطح شیشه خارجی (m2 ) و ΔT نشان دهنده ضریب انتقال حرارتی لعاب است. تفاوت دما بین محیط بیرون و داخل در درجه سانتیگراد.

4 حوزه محاسباتی

4.1 روش های حل

با توجه به Chaudhry [ 16 ]، ناویر استوکس سه بعدی رینولدز متوسط ​​(RANS) با استفاده از  مدل آشفتگی k – ε یکی از مناسب ترین روش های عددی برای شبیه سازی محیط ساخته شده است. این مدل دو معادله حاکم، انرژی جنبشی و اتلاف آشفته را حل می کند. ابید و همکاران [ 9 ] از RANS همراه با مدل انرژی جنبشی آشفته برای شبیه‌سازی اثر تابش خورشیدی در ریزاقلیم شهری استفاده کرد. پریادرشینی و همکاران [ 12 ] از  مدل آشفتگی k – ε برای شبیه سازی پدیده جزیره گرمایی در سنگاپور استفاده کرد. علاوه بر این، آلگرینی و همکاران. [ 6 ]، سونگ و لیو [ 15 ]، توپارلار و همکاران. [14 ]، و غیره، همان را برای تحلیل و شبیه سازی محیط های حرارتی شهری همانطور که در بخش 2.2 بررسی شد، اعمال کرده اند . بنابراین مطالعه پژوهشی کنونی از معادله رینولدز-متوسط ​​ناویر-استوکس (RANS) حالت پایدار سه بعدی با استفاده از مدل  آشفتگی k – ε برای حل دامنه محاسباتی گسسته استفاده کرد.

بنابراین، معادلات RANS به همراه معادلات تکانه و تداوم با استفاده از کد تجاری CFD برای شبیه‌سازی میدان سرعت و فشار حل شد. این مدل از روش کنترل حجم و روش نیمه ضمنی برای معادلات پیوندی با فشار (SIMPLE) الگوریتم جفت شدن سرعت-فشار با گسسته سازی مرتبه دوم رو به باد استفاده می کند. انرژی جنبشی آشفتگی، k و سرعت اتلاف آن، e ، از معادلات انتقال زیر فرمول‌بندی شده در معادلات (2) و (3) به دست می‌آیند .(2)(3)که در آن G k تولید انرژی جنبشی آشفتگی ناشی از گرادیان‌های میانگین سرعت را نشان می‌دهد، G b نشان‌دهنده تولید انرژی جنبشی آشفتگی ناشی از شناوری است. Y M نشان دهنده سهم اتساع نوسانی در آشفتگی تراکم پذیر در نرخ اتلاف کلی است. C 1 e ، C 2 e و C 3 e ثابت هستند. σ k و σ e اعداد پراندتل آشفته برای k و e هستند .

4.2 تولید مش

مطالعات قبلی توسط Shengwei [19 ]، دندگاوال [20 ]، هودال و همکاران. [21 ] نشان داده اند که یک شبکه چهار وجهی قادر است جواب های عددی را با دقت قابل قبولی در هزینه نسبتا کم و زمان محاسباتی محاسبه کند. بنابراین یک شبکه چهار وجهی بدون ساختار به دلیل پیکربندی هندسی پیچیده آن که دارای مجموعه ای از دیوارهای بیرونی برافروخته و فرورفته است، اعمال شد. مرکز ربط شبکه بدون ساختار به منظور حل کردن میدان‌های جریان و شارهای دما در اطراف گوشه‌های دیوارهای فرورفته و فلاش شده، روی «ریز» با اندازه‌های هموارسازی «بالا» تنظیم شد. مدل مش شده شامل 113 040 گره و 580 464 عنصر برای محاسبه کدهای عددی حوزه فیزیکی است. شکل 6 مدل مشبک شده را با شرایط ورودی نشان می دهد.

بند انگشتی شکل 6نمایش دامنه جریان.

4.3 شرایط مرزی

شرایط مرزی اعمال شده ( صفحه 1 ) شامل سرعت ورودی 0.4 متر بر ثانیه است که از جهت شمال غربی نزدیک می شود. هندسه به عنوان یک منطقه جامد مدل شد در حالی که محفظه به عنوان یک منطقه سیال برای تجزیه و تحلیل مدل شد. شرایط مرزی در طول بررسی عددی برای همه مدل‌های تحلیل‌شده یکسان نگه داشته شد.

میز 1شرایط مرزی.

5 نتایج و تجزیه و تحلیل

5.1 اعتبار سنجی مدل

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

همانطور که از شکل 7 مشاهده شد، درصد خطا برای متغیر دما از 2٪ تا 17٪ با خطای متوسط ​​11٪ متغیر بود، در حالی که پارامتر سرعت دارای یک تغییر میانگین در محدوده 19٪ بود. نوارهای خطا برای ضبط‌های CFD با 5% رسم می‌شوند. به گفته ویلکینسون و همکاران. [ 17 ]، یک خطای قابل قبول برای نتایج CFD تا 20٪ است. بنابراین می توان ثابت کرد که مدل عددی اعتبار سنجی شده است.

بند انگشتی شکل 7مقایسه بین CFD و نتایج تجربی برای (الف) دما و (ب) سرعت.

5.2 پروفیل های جریان و حرارتی

شکل 8 توزیع دما را در جهت شرایط ورودی به دمای محیط 27 درجه سانتیگراد نشان می دهد. دماهای بالاتر در اطراف مناطق مجاور نمای لعابدار مشاهده شد. بخشی از منطقه 1 که در بین هندسه شسته شده و فرورفته ساختمان فاز 1 قرار دارد، همانطور که نشان داده شده است، دمای سطح نسبتاً گرمتر (هات اسپات) را نشان می دهد که منجر به شدت تبادل حرارتی احتمالی بالاتر در آن ناحیه می شود. دمای سطح پایین‌تر در منطقه 2 مشاهده شد، اگرچه شارهای حرارتی همرفتی بالاتر که معمولاً با هر نمای شیشه‌ای ساختمان رخ می‌دهد منجر به افزایش دما می‌شود. بالاترین دمای سطح کلی که در سطح نسبتاً بزرگتری توزیع شده است در منطقه 3 است.

شکل 9سطوح سرعت میدان جریان را نشان می دهد که در اطراف 3 منطقه پردیس به هم می زند. با سرعت ورودی 0.4 متر بر ثانیه، سطوح کانتور در منطقه 1 نشان می دهد که این منطقه بالاترین سرعت هوا را در آن منطقه تا 0.64 متر بر ثانیه تجربه کرده است. به دنبال تماس بین ساختمان و جریان باد نزدیک، سرعت مجاور ساختمان ها به میزان قابل توجهی کاهش یافت. پروفیل جریان در اطراف هندسه فرورفته کاهش قابل توجهی در سرعت باد نشان می دهد که تقریباً می تواند به عنوان شرایط باد راکد در نظر گرفته شود. با این حال، راهروی بیرونی بین فاز 1 و فاز 2 ساختمان، جریان باد را در مقایسه با جریان بین هندسه ساختمان تسریع کرده است. نزدیک شدن باد از منطقه 1 به راهروی باریک منطقه 2 وارد شد.

تغییرات دما بین پیکربندی ساختاری ساختمان ها و میکرو اقلیم اطراف در شکل 10 نشان داده شده است . افزایش دما 28 درصد در دیوارهای خارجی و لعاب مشاهده شد که با افزایش دمای سطح نشان داده شد. افزایش دمای 7٪ در مناطق مجاور به هر دو فاز ساختمان مشاهده شد که تأثیر محلی جزایر گرمایی را برجسته می کند.

به منظور ارزیابی کمی یافته ها، خطی ایجاد شد که از طریق سه منطقه در ارتفاع 2 متری از سطح زمین متقاطع می شود تا شرایط تجربه شده توسط عابران پیاده را شبیه سازی کند. شکل 11a شکل گیری دما را در سراسر هندسه ساختمان نشان می دهد. همانطور که مشاهده شد، دما در منطقه 3 با مقدار 28.62 درجه سانتیگراد به بالاترین حد خود رسید. این روند بر خلاف سرعت هوا در شکل 11b بود که نشان دهنده کاهش به شرایط تقریباً راکد بود، بنابراین تناسب معکوس بین دو پارامتر را تأیید کرد.

دو خط عمود بر هم (خطوط 3 و 4) در ناحیه 1 و 3 بر روی نمودار دمایی که دو فاز ساختمان را به هم متصل می کند رسم شد. شکل 12a نشان می دهد که خط 4 در منطقه 3 دارای بالاترین دمای حداکثر 33.8 درجه سانتی گراد در مقایسه با خط 3 در منطقه 1 است. روند کلی نمودار همچنین نشان می دهد که منطقه 3 دارای شدت دمای بالاتری نسبت به منطقه 1 است. بعلاوه دما در ناحیه نزدیک به نمای ساختمان افزایش می یابد و تقریباً در ناحیه مرکزی خطوط در هر دو ناحیه ثابت می ماند. کمترین دما در این مناطق مرکزی خط 3 و 4 به ترتیب 26.8 درجه سانتیگراد و 27.1 درجه سانتیگراد ثبت شد که منطقه 1 دارای کمترین دما بود.

شکل 12b سرعت هوای کمتری را در نزدیکی ساختمان فاز 1 نسبت به فاز 2 (منطقه 1) نشان می دهد. حداکثر سرعت هوا 0.17 متر بر ثانیه در خط 4 است که در زیر منطقه 3 نزدیک به ساختمان فاز 2 قرار دارد. با این حال، سرعت کلی هوا در خط 3 در منطقه 1 بیشتر بود. کمترین حداقل سرعت در سمت فاز 1 منطقه 3 با 0.005 متر بر ثانیه مشاهده شد که تقریباً می تواند به عنوان شرایط باد راکد در نظر گرفته شود.

بند انگشتی شکل 8دمای سطح در صفحه XY در 2 متر بالاتر از سطح زمین.

بند انگشتی شکل 9توزیع سرعت در صفحه XY در ارتفاع 2 متری از سطح زمین.

بند انگشتی شکل 10بازنمایی تشکیلات دمایی محلی بر روی هندسه ساختمان و محیط اطراف.

بند انگشتی شکل 11تغییرات آنی در محوری (الف) دما، (ب) پارامترهای سرعت.

بند انگشتی شکل 12تغییرات آنی در شعاعی (الف) دما، (ب) پارامترهای سرعت.

5.3 بارهای خنک کننده محسوس

نقاط کاوشگر در تمام مناطق مجاور پنجره ها در راهروی بیرونی در رابط نتایج ANSYS ترسیم شد و با استفاده از ماشین حساب عملکرد میانگین در نرم افزار محاسبه شد. افزایش گرمای رسانا ناشی از لعاب با حداکثر دمای شناسایی شده (28.22 درجه سانتیگراد) برای یافتن سهم احتمالی افزایش حرارت محسوس بر بار خنک کننده ساختمان به عنوان بدترین سناریو محاسبه شد. شکل 13 نمایش گرافیکی افزایش بارهای خنک کننده قبل و بعد از اثر جزیره گرمایی حالت پایدار را نشان می دهد. افزایش 1.22 درجه سانتیگراد در دمای هوای محیط در مجاورت لعاب باعث شد که افزایش گرمای رسانا از طریق لعاب 19.61٪ افزایش یابد. این در بخش 6.3 بیشتر مورد بحث قرار گرفته است .

بند انگشتی شکل 13بار خنک کننده معقول ناشی از لعاب قبل و بعد از اثر جزیره گرمایی.

6 بحث

6.1 وجود نقطه داغ و رابطه بین دما و سرعت هوا

رایج ترین روندی که می توان از نتایج آن تحلیل کرد، افزایش دمای محلی در مناطق با سرعت باد کم بود. با توجه به یوان و نگ [22 ] سرعت باد کمتر از 0.3 متر بر ثانیه به عنوان تهویه طبیعی در سطح عابر پیاده راکد در دره خیابان در نظر گرفته می‌شود که به معنای حذف ناکافی گرما است که در طول زمان افزایش می‌یابد. گرمای دفع شده از سطوح شهری (پنجره ها، زمین و غیره) به دلیل تبادلات حرارتی همرفتی کمتر به دلیل حرکت محدود هوا، به افزایش دمای محیط می افزاید. از این رو شدت دمای بالاتر در این مناطق منجر به نقاط داغ تجربه شد. به دلیل وجود این هات اسپات ها، دمای هوای مجاور با دمای بالاتر کانون ها در تعامل است. قانون صفر ترمودینامیک بیان می کند که دمای بالاتر از یک سیال یا جسم به جسم یا سیالی با دمای پایین تر برای رسیدن به تعادل حرارتی منتقل می شود.23 ]، بنابراین دمای هوای محیط نسبتاً پایین تر (27 درجه سانتیگراد) از طریق تبادل گرما با نقاط داغ (29.7 درجه سانتیگراد) افزایش می یابد و در نتیجه در سراسر حوزه فیزیکی پخش می شود. رابطه بین دما و سرعت باد مشاهده شده نتیجه می گیرد که دما با سرعت هوا نسبت معکوس دارد.

همانطور که در شکل 6 مشاهده می شود، منطقه 2 افزایش دما 0.74 درجه سانتی گراد (27.7 درجه سانتی گراد) را تجربه کرد . به گفته ووگت [24 ]، جزایر گرمایی زمانی ایجاد می شوند که افزایش دما بیش از 1 درجه سانتی گراد از دمای واقعی محیط باشد. شکل 7 سرعت نسبتاً بالاتر هوا (0.32-0.4 متر بر ثانیه) را نشان می دهد که به تخلیه گرما کمک می کند. اگرچه وجود لعاب باعث افزایش دما به دلیل افزایش تبادل حرارت همرفتی از دمای نسبتاً بالاتر سطح لعاب در منطقه 2 می شود، راهروی باریک و ساختمان مرکزی فاز 2 پلکانی امکان افزایش اثر جریان را فراهم می کند که به افزایش سرعت هوا کمک می کند و امکان حذف گرما را فراهم می کند. . برای حمایت از این پدیده، شواهدی از ادبیات نشان می دهد که افزایش سرعت باد باعث افزایش انتقال حرارت همرفتی با نمای ساختمان و حذف گرما می شود [ 6]]. در نتیجه، منطقه 2 و محل ورودی دمای سطح پایین تری را نشان می دهد که نشان دهنده انتقال حرارت ناچیز است.

6.2 رابطه بین هندسه ساختمان و سرعت/الگوی جریان باد

همانطور که از یافته‌ها مشاهده شد، باد فشار بیشتری را بر سمت بادگیر ساختمان (منطقه 1) وارد کرد و باعث شد جریان باد به سمت بادگیر (منطقه 3) حرکت کند. سرعت اولیه باد غالب 0.4 متر بر ثانیه اندازه گیری شد، با این حال، پس از ضربه، باد به صورت جانبی در اطراف ساختمان قیچی کرد و باعث کاهش سرعت در اطراف ساختمان به دلیل کشش اصطکاکی ناشی از ناهمواری سطح ساختمان شد. به دلیل شتاب ناشی از انحراف باد هنگام برخورد با ساختمان، سرعت باد در فضای خالی بیرونی افزایش می یابد. بنابراین، حداکثر سرعت هوا در منطقه یک منطقه به 0.64 متر بر ثانیه رسید.

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

6.3 پیامدهای افزایش دما بر تقاضای سرمایش

مشخصات حرارتی در شکل های 6 و 8 دمای بالاتری را در مجاورت لعاب نشان می دهد. حداکثر دما در منطقه 2 در حدود 35.08 درجه سانتیگراد است که 23٪ بالاتر از دمای محیط است. افزایش دمای لعاب بر دمای هوای حاصل در مجاورت این نماهای لعابدار تأثیر می گذارد. این افزایش دما باعث شد که بار خنک کننده معقول از طریق لعاب 19.61٪ افزایش یابد. به دنبال اثر جزیره گرمایی، دمای محیط 4.32 درصد افزایش یافت. از آنجایی که بار خنک کننده محسوس به دمای خارجی بستگی دارد، هرگونه تغییر در دمای خارجی منجر به تغییر در بارهای خنک کننده محسوس از طریق لعاب می شود. بنابراین مشاهده شد که بارهای خنک کننده از 986 وات به 1227 وات افزایش یافته است.

از یک مطالعه انجام شده در یونان می توان دریافت که با فراتر رفتن میانگین جزیره گرمایی حالت پایدار از 10 درجه سانتیگراد، بار خنک کننده ساختمان ها دو برابر می شود.25 ].

بر اساس EPA [26 ] به ازای هر 0.6 درجه سانتیگراد افزایش دمای هوا، مصرف برق به دلیل افزایش تقاضای خنک کننده بین 1.5٪ تا 2٪ افزایش می یابد. علاوه بر این، تغییرات در دمای هوای محیط خارجی منجر به عواقب قابل توجهی بر عملکرد حرارتی ساختمان می شود [ 13 ]. الصلال [27 ] تخمین زده است که به دلیل شرایط آب و هوایی سخت امارات متحده عربی، 60 درصد از انرژی مصرفی ساختمان ناشی از تقاضای بالای سرمایش است. بنابراین حتی یک درجه افزایش جزئی در دمای هوا می تواند افزایش فاجعه باری در کل تقاضای انرژی یک ساختمان داشته باشد. از این رو HI یکی از پدیده هایی است که بر اقلیم میکرو محلی تأثیر منفی می گذارد و منجر به افزایش تقاضای خنک کننده فضا می شود که به نوبه خود باعث افزایش انتشار کربن به دلیل افزایش اوج بار الکتریکی ناشی از تهویه مطبوع می شود.

7 نتیجه گیری

این مطالعه جریان باد و توزیع دما را برای ارزیابی دمای حالت پایدار و مشخصات جریان در اطراف ساختمان‌ها در امارات تجزیه و تحلیل کرد. بررسی عددی با استفاده از CFD انجام شد که یافته‌های آن با استفاده از آزمایش میدان دور تأیید شد. یافته‌ها حاکی از آن است که آن مکان‌های خاص در منطقه 1 و 3، 9.1 درصد افزایش دمای محیط را تجربه کرده‌اند که باعث افزایش دمای 2.7  درجه سانتی‌گراد شده و منجر به تشکیل کانون‌ها در این مناطق شده است. یافته‌ها نشان داد که الگوی جریان باد و سرعت یکی از عوامل اصلی شکل‌گیری دماهای بالاتر همراه با هندسه ساختمان است که سرعت باد و تبادل حرارت همرفتی بالا از نمای شیشه‌ای با دمای سطح نسبتاً بالا را تغییر می‌دهد.

مناطق با دمای هوای بالا سرعت نسبتاً کمتری در مناطق مربوطه داشتند. بنابراین نتیجه گیری شد که دمای هوا با سرعت باد در اطراف محوطه دانشگاه نسبت معکوس دارد. بالاترین دمای حداکثر 28.22 درجه سانتیگراد در مجاورت پنجره نزدیک منطقه 1 به دلیل افزایش 4.32 درصدی در دمای سطح مجاور مشاهده شد. افزایش دما باعث افزایش 19.61 درصدی بار محسوس شد که منجر به تأثیر نامطلوب بر تقاضای برق شد. بر اساس تجزیه و تحلیل گسترده انجام شده در این کار، این مطالعه پتانسیلی را برای ارزیابی و کمی کردن عواقب جزیره گرمایی حالت پایدار در زمینه الزامات خنک کننده معقول ایجاد کرده است.

منابع

  1. EPA، اثر جزیره گرمایی [آنلاین] (2014)، موجود: http://www.epa.gov/heatisland/ (دسترسی: 09/15/2014) [Google Scholar]
  2. TR Oke، مبنای انرژی جزیره گرمایی شهری، JR Meteorol. Soc. 108 ، 1 (1982) [Google Scholar]
  3. توسط تام، وا گوگ، تی محسن، تأثیر شهرنشینی و گرمای شهری و تأثیر آن بر تغییرات دمایی روزانه، اقلیم شهری. 12 ، 1 (2014) [CrossRef] [Google Scholar]
  4. P. Rajagopalan، KC Lim، E. Jamei، جزیره گرمایی شهری و ویژگی‌های جریان باد یک شهر گرمسیری، انرژی خورشیدی 107 ، 159 (2014) [CrossRef] [Google Scholar]
  5. گروه تحقیقات آب و هوا، پروژه تحقیقاتی جزیره گرمایی شهری [آنلاین] (2012)، موجود: http://www.uta.edu/faculty/awinguth/uhi/dfw_uhi.html (دسترسی: 10/04/2014) [Google Scholar ]
  6. J. Allegrini، V. Kämpf، J. Dorer، J. Carmeliet، مدل‌سازی ریزاقلیم شهری و تأثیر آن بر تقاضای انرژی ساختمان در یک محله شهری، در CISBAT 2013 (2013) [Google Scholar]
  7. M. Brown, Urban parameterizations for mesoscale meteorological models [آنلاین] (2014)، موجود: http://www.researchgate.net/publication/258703466_Urban_parameterizations_for_mesoscale_meteorological_models (دسترسی: 2014/10) [Google/10]
  8. B. Fisher, J. Kukkonen, M. Piringer, MW Rotach, M. Schatzmann, Meteorological اعمال شده در مسائل آلودگی هوای شهری: مفاهیمی از COST 715, Atmos. شیمی. فیزیک بحث و گفتگو. 5 ، 7903 (2005) [CrossRef] [Google Scholar]
  9. N. Abid، I. Kolo، I. Janajreh، مطالعه دینامیک سیالات محاسباتی بر روی باد، اثر خورشیدی و پراکندگی آلاینده در شهر مصدر، در نوزدهمین کنفرانس مکانیک سیالات استرالیا ملبورن، 2014 (2014) [Google Scholar]
  10. K. Cho, T. Kono, Y. Ashie, شبیه‌سازی CFD در مقیاس بزرگ پدیده جزیره گرمایی و اقدامات متقابل در توکیو، در گزارش سالانه مرکز شبیه‌ساز زمین، آوریل 2007 تا مارس 2008 (2008) [Google Scholar]
  11. D. Taleb، B. Abu-Hijleh، جزایر گرمایی شهری: تأثیر بالقوه پیکربندی های شهری ارگانیک و ساختار یافته بر تغییرات دما در دبی، امارات متحده عربی، Renew. Energy 50 , 747 (2012) [Google Scholar]
  12. R. Priyadarshini، WN Hien، CKW David، مدل‌سازی ریزاقلیمی محیط حرارتی شهری سنگاپور برای کاهش جزیره گرمایی شهری، انرژی خورشیدی 82 ، 727 (2008) [Google Scholar]
  13. ح. راضی، در مورد تأثیر گرمایش جهانی و محیط ساخت امارات متحده عربی، گرمایش جهانی . استوارت آرتور هریس [آنلاین] (2010)، موجود: http://www.intechopen.com/books/global-warming/on-the-effect-of-global-warming-and-the-uae-builtenvironment (دسترسی در : 2015/29/03) [Google Scholar]
  14. Y. Toparlar، B. Blocken، P. Vos، GJF van Heijst، WD Janssen، T. van Hooff، H. Montazeri، HJP Timmermans، شبیه‌سازی CFD و اعتبارسنجی میکرو اقلیم شهری: مطالعه موردی برای Bergpolder Zuid، روتردام، ساختمان. محیط زیست 83 ، 1 (2014) [Google Scholar]
  15. X. Song, J. Liu, CFD simulation of micro-climate in waterfront (2013), فرانسه [آنلاین]، موجود: http://www.ibpsa.org/proceedings/BS2013/p_1066.pdf (دسترسی: 2014/ 15/12) [Google Scholar]
  16. HN Chaudhry، مطالعه ای در مورد بهینه سازی پارامترهای هندسی لوله حرارتی برای خنک سازی غیرفعال پایدار در محیط ساخته شده، برنامه کاربردی. حرارت مهندس 93 , 486 (2016) [Google Scholar]
  17. اس. ویلکینسون، اس. هانا، ال. هسلگرن، وی. مولر، محاسبات و عملکرد آیرودینامیک القایی، در مجموعه مقالات سی و یکمین کنفرانس eCAADe 2013 (2013)، جلد. 2، صفحات 29–48 [آنلاین] [Google Scholar]
  18. راهنمای CIBSE A، طراحی محیطی (موسسه مجاز مهندسین خدمات ساختمان، لندن، 2006)، ویرایش هفتم، ص. 15 [Google Scholar]
  19. M. Shengwei، 5 سوء تفاهم برتر در مورد مش (خوب) [آنلاین] (2011)، موجود: http://caewatch.com/top-5-misunderstandings-on-good-mesh/ (دسترسی: 02/05/2015) ) [Google Scholar]
  20. P. Dandgawal، برآورد شرایط محیطی طبیعی با CFD [آنلاین] (2013)، موجود: http://www.cctech.co.in/live/2013/06/19/estimate-the-natural-environment-conditions- with-cfd/ (دسترسی در: 02/05/2015) [Google Scholar]
  21. S. Houdal, N. Zemmouri, A. Hasseine, R. Atthani, R. Belarbi, A CFD model for simulating flow شهری در شبکه خیابانی پیچیده مورفولوژیکی, TOJSAT: Online J. Sci. تکنولوژی 2 ، 1 (2012) [Google Scholar]
  22. C. Yuan، E. Ng، تخلخل ساختمان برای تهویه بهتر شهری در شهرهای با تراکم بالا: یک مطالعه پارامتری محاسباتی، ساختمان. محیط زیست 50 ، 176 (2011) [CrossRef] [Google Scholar]
  23. Y. Çengel، M. Boles، ترمودینامیک: یک رویکرد مهندسی (McGraw-Hill، ایالات متحده، 2007) [Google Scholar]
  24. JA Voogt, Urban Heat Islands: hotter citys [آنلاین] (2004)، موجود: http://www.actionbioscience.org/environment/voogt.html (دسترسی: 2014/04/10) [Google Scholar]
  25. M. Santamouris, N. Papanikolaou, I. Livada, I. Koronakis, C. Georgakis, A. Argiriou, DN Assimakopoulos, On تاثیر اقلیم شهری بر مصرف انرژی ساختمانها, Sol. Energy 70 , 201 (2001) [Google Scholar]
  26. EPA، اثر جزیره گرمایی [آنلاین] (2013)، موجود: http://www.epa.gov/heatislands/impacts/index.htm (دسترسی: 03/29/2015) [Google Scholar]
  27. KA Al-Sallal، ساختمان های برج در دبی – آیا آنها پایدار هستند؟ در کنفرانس CTBUH، سئول، 2004 [آنلاین] (2004)، موجود : http://ctbuh.org/Portals/0/Repository/Al-Sallal_2004_TowerBuildingsinDubai.b82e3a40-df34-41f1.8759-4804 /01/04) [Google Scholar]
  28. K. Setaih، N. Hamza، MA Mohammed، S. Dudek، T. Townshend، مدل سازی CFD به عنوان ابزاری برای ارزیابی شرایط آسایش حرارتی در فضای باز در محیط های شهری در آب و هوای خشک گرم، J. Inf. تکنولوژی ساخت و ساز 19 , 248 (2014) [Google Scholar]

Leave A Comment