Жесткие ОДУ и устойчивость

Жесткие дифференциальные уравнения содержат процессы, развивающиеся в широко разделенных временных масштабах, поэтому явные методы вынуждены делать непрактично малые шаги для обеспечения устойчивости; их эффективное решение требует неявных методов с сильными свойствами устойчивости.

Definition

Дифференциальное уравнение называется жестким, когда оно допускает компоненты решения, затухающие в очень разных временных масштабах, так что численная устойчивость, а не точность, диктует размер шага; теория устойчивости анализирует, какие методы могут делать большие шаги без роста ошибки.

Scope

Эта тема охватывает явление и неформальное определение жесткости, линейное тестовое уравнение и область абсолютной устойчивости, понятия A-устойчивости, A(альфа)-устойчивости и L-устойчивости, причины неэффективности явных методов при решении жестких задач, а также неявные методы — неявные Рунге-Кутты и формулы обратного дифференцирования — которые их решают.

Core questions

• Что делает задачу жесткой и почему это препятствует явным методам?
• Как определяется область абсолютной устойчивости с помощью линейного тестового уравнения?
• Что требуют A-устойчивость и L-устойчивость, и почему они важны для жестких задач?
• Какие методы обеспечивают устойчивость, необходимую для жестких и дифференциально-алгебраических систем?

Key theories

Абсолютная устойчивость и тестовое уравнение

Применение метода к скалярному линейному тестовому уравнению дает коэффициент усиления; множество произведений размера шага на собственное значение, для которых этот коэффициент имеет величину не более единицы, является областью абсолютной устойчивости метода, которая должна содержать жесткие собственные значения задачи, чтобы допускать большие шаги.

A-устойчивость и L-устойчивость

Метод является A-устойчивым, если его область устойчивости содержит всю левую полуплоскость, так что он стабилен для всех затухающих мод независимо от размера шага, и L-устойчивым, если он дополнительно полностью демпфирует очень жесткие моды; эти свойства выделяют неявные методы, подходящие для жестких задач.

Mechanisms

В жесткой задаче наиболее быстро затухающая мода имеет большое отрицательное собственное значение; ограниченная область устойчивости явного метода вынуждает размер шага разрешать эту моду даже долго после того, как она физически затухла, что делает вычисления безнадежно медленными. Неявные методы, такие как метод обратного Эйлера, неявные схемы Рунге-Кутты и формулы обратного дифференцирования, имеют области устойчивости, покрывающие левую полуплоскость (или большую ее часть), поэтому они остаются стабильными при больших шагах и позволяют выбирать размер шага исключительно по точности. Каждый шаг затем требует решения (обычно нелинейной) алгебраической системы, как правило, с помощью итерации Ньютона с использованием якобиана.

Clinical relevance

Жесткость широко распространена в сетях химических реакций, горении, электрических цепях, системах управления и дискретизациях параболических дифференциальных уравнений в частных производных методом линий; распознавание жесткости и выбор соответствующего стабильного неявного решателя необходимы для получения результатов в приемлемое время, и большинство промышленных программных обеспечений для ОДУ включают автоматическое обнаружение жесткости и переключение методов.

History

Понятие жесткости было определено Кертиссом и Хиршфельдером в 1952 году, а поддерживающая теория устойчивости — A-устойчивость и барьеры порядка — была разработана Дальквистом; коды формул обратного дифференцирования Гира и более поздние неявные методы Рунге-Кутты высокого порядка сформировали практический инструментарий для жестких и дифференциально-алгебраических задач.

Key figures

• Germund Dahlquist
• C. William Gear
• Ernst Hairer
• Gerhard Wanner

Seminal works

• hairer1996
• iserles2008

Frequently asked questions

Что именно делает ОДУ жестким?

Жесткость возникает, когда система имеет компоненты, которые затухают гораздо быстрее, чем развивается интересующее решение. Единого четкого определения нет, но практический признак заключается в том, что явные методы вынуждены использовать очень малые шаги для обеспечения устойчивости, даже когда точность позволяет использовать большие шаги.

Почему жесткие задачи требуют неявных методов?

Неявные методы могут иметь области устойчивости, покрывающие всю левую полуплоскость (A-устойчивость), поэтому они остаются стабильными при больших размерах шага для быстро затухающих мод. Явные методы имеют ограниченные области устойчивости, что вынуждает использовать крошечные шаги и делает их непрактичными для жестких задач.