بررسی روش های مختلف گسسته سازی بر جواب معادله پخش یک بعدی ناپایا

بررسی روش های مختلف گسسته سازی بر جواب معادله پخش یک بعدی ناپایا

بررسی روش های مختلف گسسته سازی بر جواب معادله پخش یک بعدی ناپایا در حالت نفوذ لایه مرزی با تحریک اولیه در مرکز

در این پژوهش به بررسی روش­های مختلف گسسته سازی بر جواب معادله پخش یک­بعدی ناپایا  در حالت نفوذ لایه مرزی با تحریک اولیه در مرکز پرداختیم. برای این منظور با گسسته سازی معادلات فوق توسط روش­های تفاضل محدود و حجم محدود و استفاده از کد نوشته شده در MATLAB به بررسی روش­های مختلف حل عددی پرداخته­ایم.  با افزایش پارمتر r (عدد کورانت) که متناسب نسبت گام زمانی بر توان دوم گام مکانی است، و در نتیجه افزایش گام زمانی یا افزایش تعداد گره­ها حل از حل واقعی فاصله می­گیرد. با افزایش عدد کورانت حل در روش صریح به سرعت ناپایدار می­شود. روش دوفورت فرانکل با اینکه در بررسی پایداری نیومن همواره پایدار است ولی با افزایش عدد کورانت حل آن از حل واقعی فاصله می­گیرد و حل نوسانی می­شود. دو روش صریح و کرانک نیکلسون نیز همواره پایدارند. متدهای صریح مرتبه اول و دوفورت فرانکل مقدار تابع را کمتر از مقدار واقعی در هر گام زمانی نشان می­دهند. روش صریح نیز مقدار را برای r های بزرگ بیش از مقدار واقعی تخمین می­زند.

 

فهرست

۱    مقدمه. ۴

۱‌.۱‌   لایه مرزی.. ۴

۲    تئوری.. ۵

۳   تعریف مسائل.. ۸

۴   نتایج حل نفود لایه مرزی.. ۸

۴‌.۱‌   رشد لایه مرزی در کورانت ثابت… ۱۲

۵   بحث و نتیجه گیری.. ۱۵

پیوست‌ها ۱۷

 

 

فهرست شکل ها

شکل ‏۴‌.‌‌۱ نمودار تغییرات سرعت بر حسب مکان t=0.5 s , r=0.1. 8

شکل ‏۴‌.‌‌۲ نمودار تغییرات سرعت بر حسب مکان  t=0. 5 s , r=0.5 0. 9

شکل ‏۴‌.‌‌۳ نمودار تغییرات سرعت بر حسب مکان t=0.5 s , r=0.6. 9

شکل ‏۴‌.‌‌۴ نمودار تغییرات سرعت بر حسب مکان t=0.5 s ,r= 1. 10

شکل ‏۴‌.‌‌۵ نمودار تغییرات سرعت بر حسب مکان t=0.5 s , r=2.5. 10

شکل ‏۴‌.‌‌۶ نمودار تغییرات سرعت بر حسب مکان  t=0.5 s , r=4. 11

شکل ‏۴‌.‌‌۷ نمودار تغییرات سرعت بر حسب مکان t=0.01 s , r=0.3. 12

شکل ‏۴‌.‌‌۸ نمودار تغییرات سرعت بر حسب مکان t=0.05 s , r=0.3. 12

شکل ‏۴‌.‌‌۹ نمودار تغییرات سرعت بر حسب مکان t=0.1 s , r=0.1. 13

شکل ‏۴‌.‌‌۱۰ نمودار تغییرات سرعت بر حسب مکان t=0.2 s , r=0.1. 13

شکل ‏۴‌.‌‌۱۱ نمودار تغییرات سرعت بر حسب مکان  t=0.4 s , r=0.1. 14

شکل ‏۴‌.‌‌۱۲ نمودار تغییرات سرعت بر حسب مکان  t=2 s , r=0.1. 14

شکل ‏۴‌.‌‌۱۳  نمای سه بعدی پخش لایه مرزی x,t,u   برای t=2 s , r=0.1 با روش کرنک نیکلسون.. ۱۵

تمام کدهای لازم این پروژه به پیوست می باشد.

خرید کامل پروژه

قیمت :۲۰۰۰۰تومان

 

درباره نویسنده

88مطلب نوشته است .

نوشتن دیدگاه

شما میتوانید از تصاویر مخصوص خود در قسمت نظرات استفاده نمایید برای اینکار از وب سایت آواتارکمک بگیرید .

تمام حقوق این سایت برای © 2019 پرتال جامع مهندسی مکانیک. محفوظ است.
قدرت گرفته از وردپرس فارسی
با توجه به مشکلاتی که در سایتهای بانکی به وجود آمده است لطفا در صورت مشکل در پرداخت از فیلتر شکن یا VPN استفاده نمایید.

پرتال جامع مهندسی مکانیک