مقدمه
معادلات سهموی (Parabolic Equations) به طور گسترده‌ای در مدل‌سازی پدیده‌های فیزیکی و مهندسی استفاده می‌شوند. یکی از روش‌های مؤثر برای حل این معادلات، روش ( ADI (Alternating Direction Implicit) است. این روش به ویژه در مسائل چند بعدی و زمانی که پایداری و دقت اهمیت دارد، کاربردی است.
فرم کلی معادله سهموی
معادله سهموی به شکل زیر است:
که در آن:
• u: تابع ناشناخته (مثلاً دما یا غلظت)
• t : زمان
• x: مکان
• α: ثابت انتشار
• ( f(x, t: تابع منبع یا منبع حرارتی
روش ADI
روش ADIیک تکنیک عددی است که به حل معادلات دیفرانسیلجزئی سهموی کمک می‌کند. این روش به دو مرحله تقسیم می‌شود:
1. مرحله اول : حل معادله در جهت اول (به عنوان مثال، محور x)
2. مرحله دوم : حل معادله در جهت دوم (به عنوان مثال، محور y)
این رویکرد به ما اجازه می‌دهد تا از مزایای روش‌های ضمنی استفاده کنیم بدون اینکه نیاز به حل یک سیستم بزرگ از معادلات خطی داشته باشیم.
مراحل اجرای روش ADI
1. گسسته‌سازی دامنه : دامنه زمانی و مکانی را به شبکه‌ای گسسته تقسیم کنید.
2. تبدیل معادله به سیستم خطی : با استفاده از روش ضمنی، معادله را به یک سیستم خطی تبدیل کنید.
3. حل سیستم خطی : از روش‌های عددی برای حل سیستم خطی استفاده کنید.
4. به‌روزرسانی مقادیر : مقادیر جدید را در شبکه گسسته به‌روزرسانی کنید و مراحل را تکرار کنید.
کد نویسی با پایتون
کد پایتون برای حل معادله سهموی با روش ADI
در اینجا یک کد ساده برای حل معادله انتقال حرارت در یک میله یک بعدی با استفاده از روش ADI آورده شده است:
