کد متلب حل معادله انتقال حرارت دو بعدی ناپایدار(دو بعد مکان و یک بعد زمان) به روش ADI با شرایط مرزی دما ثابت با توضیحات خط به خط در ورد
این فایل حاوی کد متلب حل معادله پارابولیک نا پایدار دو بعدی انتقال حرارت به روش ADI می باشد
مشخصات فایل
| تعداد صفحات |
3 |
| حجم |
0 کیلوبایت |
| فرمت فایل اصلی |
rar |
| دسته بندی |
مهندسی شیمی |
توضیحات کامل
به نام خدا
این فایل حاوی کد متلب حل معادله پارابولیک نا پایدار دو بعدی انتقال حرارت به روش ADI با شرایط مرزی دما ثابت می باشد که در زیر آمده است:
T_{t} = alpha*( T_{xx} + T_{zz} )
Boundary and initial condition is constant temperature
این برنامه ابتدا ابعاد و تعداد مش ها در سه جهت را از کاربر دریافت کرده و با تشکیل ماتریس ضرایب A و ماتریس ثابتها b آنها را به روش ماتریس معکوس حل و نتایج هر گام زمان را ضمن نمایش در command window با دستور mesh رسم می کند.
در انتها نیز حل عددی را با حل تحلیلی مقایسه و خطا را با دستور contour رسم می کند.
توضیحات برنامه در دو بخش توضیحات کلی و جزئی در 4 صفحه ورد ارائه شده است.
توضیحات کلی شامل نحوه ی مش بندی، نام گذاری نود ها، درایو کردن معادلات، تشکلی ماتریس ضرایب و ثابت ها A و b ، معرفی پارامتر ها می باشد.
توضیحات جزئی نیز شامل توضیحات خط به خط برنامه است که در زیر آورده شده است:
ابتدا برنامه در خطوط 16 تا 19 مقادیر مربوط به ابعاد و تعداد گامهای زمان و زمان انتهایی را از کاربر دریافت می کند. سپس با این مقادیر تعداد گره ها در جهت x و z را محاسبه می کند و دو رشته با فواصل محاسبه شده با نامهای x و z می سازد.
در بلوک for مربوط به خطوط 31 تا 36 مقادیر شرایط مرزی را از تابع آن (Bound_Cond) استخراج کرده و در ماتریس T1 قرار می دهد.
در خط 28 هم شرط اولیه را از کاربر دریافت و در ماتریس T1 قرار می دهد.
دقت شود که شرایط مرزی بعد از دریافت شرط اولیه اعمال می شود (در خط 34) زیرا؛ لبه های صفحه در لحظه شروع انتقال حرارت هم در دمای ثابت قرار دارند.
در بلوک مربوط به خطوط 42 تا 135 که بلوک اصلی برنامه است ماتریسها را تشکیل و در خطوط76 و 111 آنها را به روش ماتریس معکوس حل می کند.
در این بلوک مسئله به دو قسمت تقسیم شده و برای هر بعد جداگانه ماتریس ها تشکیل و حل می شوند.
خطوط 47 تا 81 برای بعد x و خطوط 85 تا 119 برای بعد z می باشد.
در خطوط 123 تا 133 ضمن نشان دادن ماتریس توزیع دما در گام جاری زمان این ماتریس توسط دستور mesh در نمودار مربوطه رسم می شود.
در خطوط 136 تا 140 هم حل تحلیلی به دست می آید و اختلاف آن با حل عددی که در خط 141 به دست می آید؛ در فیگوری جدا با دستور contour رسم می شود.
در انتها در command window ماتریس های مربوط به توزیع دما در هر گام زمان نمایش داده می شود.
توضیحات بیشتر و دانلود
صدور پیش فاکتور، پرداخت آنلاین و دانلود
