УСЛОВИЕ КУРАНТА
ДЛЯ ЯВНОЙ СХЕМЫ ЭЙЛЕРА
(ЛИНЕЙНОЕ УРАВНЕНИЕ ДИФФУЗИИ)
-
коэффициент диффузии:
-
-
источники тепла отсутствуют:
-
-
поддержание постоянной температуры на краях стержня:
-
-
начальное распределение температуры (см. график):
шаг по времени:
число точек по оси X:
соответственно, пространственный шаг:
дискретизация начального условия на сетке
Разностное решение для этих шагов можно увидеть здесь.
Следующее выражение для шагов сетки называют соотношением Куранта.
Можно показать, что явная разностная схема Эйлера для задачи диффузии является устойчивой, если соотношение Куранта меньше 1.
Если же оно превышает 1, то разностная схема неустойчива, и ее ответ неверен.
Рассмотрим увеличенный шаг по времени, который даст соотношение Куранта большее 1:
А вот, как выглядит решение по этой схеме на первых же шагах по времени:
Из физических соображений очевидно, что это решение неверно. О неустойчивости схемы говорит характерная "разболтка" сеточного решения.
Советую читателю поэкспериментировать, посмотрев, что будет происходить на следующих шагах по времени, а также взяв другие значения пространственного и временного шага и коэффициента диффузии.