Схема ФТКС
В численном анализе метод FTCS (прямое время-центрированное пространство) представляет собой метод конечных разностей, используемый для численного решения уравнения теплопроводности и аналогичных параболических уравнений в частных производных . [ 1 ] Это метод первого порядка по времени, явный по времени и условно устойчивый при применении к уравнению теплопроводности. Когда он используется в качестве метода для уравнений адвекции или, в более общем плане, гиперболических уравнений в частных производных , он нестабилен, если не учитывать искусственную вязкость. Аббревиатуру FTCS впервые использовал Патрик Роуч. [ 2 ] [ 3 ]
Метод
[ редактировать ]Метод FTCS основан на прямом методе Эйлера во времени (отсюда «прямое время») и центральной разности в пространстве (отсюда и «центрированное пространство»), что обеспечивает сходимость первого порядка во времени и сходимость второго порядка в пространстве. Например, в одном измерении, если уравнение в частных производных имеет вид
затем, позволяя прямой метод Эйлера определяется следующим образом:
Функция должны быть пространственно дискретизированы с помощью центральной разностной схемы. Это явный метод , который означает, что могут быть вычислены явно (нет необходимости решать систему алгебраических уравнений), если значения на предыдущем временном уровне известны. Метод FTCS не требует больших вычислительных затрат, поскольку метод является явным.
Иллюстрация: одномерное уравнение теплопроводности
[ редактировать ]Метод FTCS часто применяется для решения диффузионных задач. Например, для одномерного уравнения теплопроводности :
Схема FTCS определяется следующим образом:
или, позволяя :
Стабильность
[ редактировать ]Как получено с помощью анализа устойчивости фон Неймана , метод FTCS для одномерного уравнения теплопроводности является численно устойчивым тогда и только тогда, когда выполняется следующее условие:
То есть выбор и должен удовлетворять вышеуказанному условию, чтобы схема FTCS была стабильной. В двумерном измерении условие становится
Если мы выберем , то условия устойчивости становятся , , и для одно-, двух- и трехмерных приложений соответственно. [ 4 ]
Основным недостатком метода FTCS является то, что для задач с большим коэффициентом диффузии , удовлетворительные размеры шага могут быть слишком малы, чтобы быть практичными.
Для гиперболических уравнений в частных производных линейной тестовой задачей является постоянный коэффициент уравнение адвекции , в отличие от уравнения теплопроводности (или уравнения диффузии ), которое является правильным выбором для параболического дифференциального уравнения . Хорошо известно что для этих гиперболических задач , любой выбор приводит к нестабильной схеме. [ 5 ]
См. также
[ редактировать ]- Уравнения в частных производных
- Метод Кранка – Николсона
- Метод конечных разностей во временной области
Ссылки
[ редактировать ]- ^ Джон К. Таннехилл; Дейл А. Андерсон ; Ричард Х. Плетчер (1997). Вычислительная механика жидкости и теплопередача (2-е изд.). Тейлор и Фрэнсис . ISBN 1-56032-046-Х .
- ^ Патрик Дж. Роуч (1972). Вычислительная гидродинамика (1-е изд.). Эрмоса . ISBN 0-913478-05-9 .
- ^ Патрик Дж. Роуч (1998). Вычислительная гидродинамика (2-е изд.). Эрмоса . ISBN 0-913478-09-1 .
- ^ Мойн, Парвиз (2010). Основы инженерного численного анализа (2-е изд.). Нью-Йорк: Издательство Кембриджского университета. ISBN 978-0-511-93263-2 . OCLC 692196974 .
- ^ Левек, Рэндалл (2002). Методы конечных объемов для решения гиперболических задач . Издательство Кембриджского университета. ISBN 0-521-00924-3 .