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