معادلات سهموی (parabolic equations) به نوعی از معادلات دیفرانسیل جزئی اطلاق می‌شود که به‌خصوص در مدل‌سازی پدیده‌های انتقال حرارت و دیگر فرآیندهای دینامیکی کاربرد دارند. یکی از مهم‌ترین ویژگی‌های این معادلات، وجود شرایط مرزی است که می‌تواند بر حل آن‌ها تأثیر بگذارد. در اینجا، ما به بررسی روش‌های عددی ضمنی برای حل معادلات سهموی با استفاده از شرط مرزی نیومن (Neumann boundary condition) می‌پردازیم.
شرط مرزی نیومن
شرط مرزی نیومن به معنای تعیین مقدار مشتق تابع در مرزهای دامنه است. به عنوان مثال، اگر u تابع مورد نظر باشد، شرط مرزی نیومن به صورت زیر بیان می‌شود:
که در آن n جهت نرمال به مرز و (g(x, t تابعی است که می‌تواند به زمان و فضا وابسته باشد.
روش Crank-Nicolson
روش Crank-Nicolson یک روش عددی ضمنی است که برای حل معادلات سهموی بسیار مناسب است. این روش با استفاده از میانگین مقادیر در زمان‌های n و n+1 کار می‌کند و معمولاً دقت بالایی دارد.
پیاده‌سازی در فرترن
کد فرترن
توضیحات کد فرترن
- • مشابه متلب، پارامترها تعریف شده و آرایه‌ها تخصیص داده می‌شوند.
- • شرایط اولیه تعیین می‌شود.
- • در حلقه اصلی، مقادیر دما با استفاده از روش Crank-Nicolson محاسبه می‌شود.
- • شرایط مرزی نیومن در انتهای دامنه اعمال می‌شود.
- • نتایج باید با استفاده از کتابخانه‌های گرافیکی نمایش داده شود.
تصویر
