Численное решение дифференциальных уравнений в частных производных
Эта область разрабатывает методы, которые дискретизируют дифференциальные уравнения в частных производных в пространстве и времени, заменяя непрерывные операторы алгебраическими системами, решения которых аппроксимируют поведение полей, управляемых физическими законами.
Definition
Численное решение дифференциальных уравнений в частных производных — это построение и анализ методов, которые аппроксимируют решения ДУЧП путем дискретизации пространственной области (и времени), что приводит к конечным системам алгебраических уравнений.
Scope
Она охватывает три основных подхода к дискретизации — методы конечных разностей, конечных элементов и конечных объемов — применяемые к эллиптическим, параболическим и гиперболическим уравнениям; анализ согласованности, устойчивости и сходимости (включая теорему эквивалентности Лакса и условие Куранта-Фридрихса-Леви); а также большие разреженные линейные и нелинейные системы, которые возникают в результате дискретизации.
Sub-topics
Core questions
- Как дифференциальные операторы в пространстве и времени дискретизируются в устойчивые, сходящиеся алгебраические системы?
- Как согласованность и устойчивость объединяются, чтобы гарантировать сходимость, как в теореме эквивалентности Лакса?
- Как тип ДУЧП — эллиптический, параболический или гиперболический — определяет подходящий метод и ограничения устойчивости?
- Как эффективно решаются возникающие большие разреженные системы?
Key theories
- Теорема эквивалентности Лакса
- Для согласованной конечно-разностной аппроксимации корректно поставленной линейной задачи с начальными условиями устойчивость является необходимым и достаточным условием сходимости; эта теорема является краеугольным камнем, который сводит доказательство сходимости к проверке согласованности и устойчивости.
- Условия устойчивости и число Куранта
- Явные схемы для нестационарных ДУЧП устойчивы только при определенных ограничениях на размеры шагов; для гиперболических задач условие Куранта-Фридрихса-Леви требует, чтобы численная область зависимости содержала физическую, ограничивая временной шаг относительно пространственной сетки.
- Вариационные принципы и принципы сохранения
- Методы конечных элементов основаны на слабых (вариационных) формулировках и проекции Галеркина, в то время как методы конечных объемов обеспечивают дискретные законы сохранения; каждый подход предоставляет путь к согласованным дискретизациям с доказуемыми аппроксимационными свойствами.
Clinical relevance
Численные методы решения ДУЧП являются вычислительной основой моделирования в инженерии и физических науках — структурной и твердотельной механике, гидро- и аэродинамике, теплопередаче, электромагнетизме, геофизике, моделировании погоды и климата, а также реконструкции медицинских изображений — везде, где непрерывные полевые уравнения должны быть решены на сложных геометриях, исключающих аналитические решения.
History
Конечно-разностный анализ ДУЧП начался с работы Куранта-Фридрихса-Леви 1928 года; метод конечных элементов возник из строительной механики и вариационной математики в 1940-60-х годах, а методы конечных объемов развивались параллельно с вычислительной гидродинамикой, при этом теорема эквивалентности Лакса обеспечила объединяющую основу сходимости в 1950-х годах.
Key figures
- Richard Courant
- Peter Lax
- Olga Ladyzhenskaya
- Randall J. LeVeque
Related topics
Seminal works
- morton2005
- leveque2007
Frequently asked questions
- Почему существует три различных подхода к дискретизации?
- Конечные разности наиболее просты на регулярных сетках, конечные элементы естественным образом обрабатывают сложные геометрии и вариационные задачи, а конечные объемы обеспечивают локальное сохранение, что делает их идеальными для потоков жидкости. Выбор зависит от геометрии, типа уравнения и того, какие свойства должны быть сохранены.
- Что означает условие Куранта-Фридрихса-Леви?
- Для явных схем в нестационарных гиперболических задачах условие Куранта-Фридрихса-Леви ограничивает максимально возможный временной шаг относительно пространственного шага сетки, гарантируя, что информация не распространяется более чем на одну ячейку сетки за один шаг. Нарушение этого условия приводит к неустойчивости.