Методы конечных разностей

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

Definition

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

Scope

Эта тема охватывает построение разностных аппроксимаций из разложений Тейлора, дискретизацию эллиптических, параболических и гиперболических уравнений в частных производных (УЧП), явные и неявные схемы интегрирования по времени (такие как прямой метод Эйлера, обратный метод Эйлера и метод Кранка-Николсона), анализ устойчивости по фон Нейману и концепцию согласованности-устойчивости-сходимости, специализированную для разностных схем.

Core questions

• Как выводятся точные разностные аппроксимации производных и как количественно оценивается их погрешность аппроксимации?
• В чем разница между явными и неявными схемами интегрирования по времени с точки зрения устойчивости и стоимости?
• Как анализ фон Неймана определяет устойчивость разностной схемы?
• Как тип уравнения определяет подходящую схему и любые ограничения на размер шага?

Key theories

Согласованность, устойчивость и сходимость

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

Анализ устойчивости по фон Нейману

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

Mechanisms

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

Clinical relevance

Методы конечных разностей широко используются для задач на регулярных областях и структурированных сетках: теплопроводность и диффузия, распространение волн и сейсмическое моделирование, вычислительная электродинамика (метод конечных разностей во временной области) и оценка опционов с помощью уравнения Блэка-Шоулза; их простота и легкость расширения до высоких порядков делают их предпочтительным выбором, когда геометрия проста.

History

Математические основы были заложены в работе Куранта-Фридрихса-Леви 1928 года по разностным уравнениям для УЧП; анализ устойчивости фон Неймана во время войны и теорема эквивалентности Лакса 1950-х годов заложили современную теорию, и разностные методы остаются основным инструментом вычислительной физики и инженерии.

Key figures

• Richard Courant
• Kurt Friedrichs
• Hans Lewy
• John von Neumann
• Randall J. LeVeque

Seminal works

• leveque2007
• morton2005

Frequently asked questions

В чем разница между явной и неявной схемой?

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

Почему методам конечных разностей может быть отдано предпочтение перед методами конечных элементов?

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