مقدمه
معادلات سهموی (Parabolic Equations) در بسیاری از زمینه‌های علمی‌و مهندسی، به ویژه در مدل‌سازی انتقال حرارت و diffusion، کاربرد دارند. یکی از روش‌های عددی موثر برای حل این معادلات، روش Crank-Nicolson است که به عنوان یک روش ضمنی شناخته می‌شود. این روش ترکیبی از روش‌های پیشرو ( Explicit) و پسرو ( Implicit) است و به دلیل پایداری بالای آن، به ویژه برای مسائل با زمان طولانی، بسیار محبوب است
شرط مرزی نیومن
شرط مرزی نیومن به معنای تعیین مقدار مشتق تابع در مرزهای دامنه است. به عنوان مثال، برای تابع ( u(x, t، شرط مرزی نیومن به صورت زیر بیان می‌شود:
که در آن nجهت نرمال به مرز و ( g(x, tتابعی است که می‌تواند به زمان و فضا وابسته باشد.
روش Crank-Nicolson
روش Crank-Nicolson یک روش عددی ضمنی است که برای حل معادلات سهمویاستفاده می‌شود. این روش با استفاده از میانگین مقادیر در دو زمان nو n+1کار می‌کند. معادله عمومی‌برای این روش به صورت زیر است:
uᵢⁿ⁺¹ - uᵢⁿ / Δ t = 1 / 2 (( ∂² u / ∂ x² |ᵢⁿ + ∂² u / ∂ x² |ᵢⁿ⁺¹ ))
پیاده‌سازی در متلب
کد متلب
توضیحات کد متلب
- • پارامترها : طول دامنه، زمان نهایی، تعداد نقاط فضایی و زمانی و ضریب انتشار تعریف شده‌اند.
- • ماتریس‌های A : برای حل ضمنی معادله تشکیل می‌شوند.
- • شرایط مرزی نیومن: با تنظیم مقادیر مناسب در ماتریس Aاعمال می‌شود.
- • حل معادله: با استفاده از یک حلقه، مقادیر دما در زمان‌های مختلف محاسبه می‌شود.
- • نمایش نتایج: نتایج به صورت سه‌بعدی نمایش داده می‌شود.
نصویر
