Сохранение и несохранение форм уравнений сохранения

Я математически понимаю, как можно получить уравнения сохранения как в консервативном

р т + ( р ты ) знак равно 0

р ты т + ( р ты ) ты + п знак равно 0

Е т + ( ты ( Е + п ) ) знак равно 0

и неконсервативные формы. Однако я все еще в замешательстве, почему мы называем их консервативными и неконсервативными формами? кто-нибудь может объяснить с физической и математической точки зрения?

Этот вопрос рассматривается во многих сторонних темах ( здесь и здесь ), но ни одна из них не дает мне достаточно хорошего ответа!

Если кто-то может подсказать, буду очень признателен.

Ответы (4)

Что это значит?

Причина, по которой они консервативны или неконсервативны, связана с расщеплением производных. Рассмотрим консервативную производную:

р ты Икс

Когда мы дискретизируем это, используя простую числовую производную только для того, чтобы подчеркнуть суть, мы получаем:

р ты Икс ( р ты ) я ( р ты ) я 1 Δ Икс

Теперь, в неконсервативной форме, производная делится на:

р ты Икс + ты р Икс

Используя ту же численную аппроксимацию, получаем:

р ты Икс + ты р Икс знак равно р я ты я ты я 1 Δ Икс + ты я р я р я 1 Δ Икс

Итак, теперь вы видите (надеюсь!) некоторые проблемы. Хотя исходная производная математически одинакова, дискретная форма отличается. Особую трудность представляет выбор членов, умножающих производную. Вот я взял его в точку я , но это я 1 лучше? Может быть, в я 1 / 2 ? Но тогда как мы получим его в я 1 / 2 ? Простое среднее? Реконструкции более высокого порядка?

Эти аргументы просто показывают, что неконсервативная форма отличается и в некотором смысле сложнее, но почему она называется неконсервативной? Чтобы производная была консервативной, она должна образовывать телескопический ряд . Другими словами, когда вы суммируете члены по сетке, должны остаться только граничные члены, а искусственные внутренние точки должны компенсироваться.

Итак, давайте посмотрим на обе формы, чтобы увидеть, как они работают. Предположим, что сетка состоит из 4 точек, начиная от я знак равно 0 к я знак равно 3 . Консервативная форма расширяется как:

( р ты ) 1 ( р ты ) 0 Δ Икс + ( р ты ) 2 ( р ты ) 1 Δ Икс + ( р ты ) 3 ( р ты ) 2 Δ Икс

Вы можете видеть, что когда вы складываете все это, вы получаете только граничные условия ( я знак равно 0 а также я знак равно 3 ). Внутренние точки, я знак равно 1 а также я знак равно 2 отменили.

Теперь давайте посмотрим на неконсервативную форму:

р 1 ты 1 ты 0 Δ Икс + ты 1 р 1 р 0 Δ Икс + р 2 ты 2 ты 1 Δ Икс + ты 2 р 2 р 1 Δ Икс + р 3 ты 3 ты 2 Δ Икс + ты 3 р 3 р 2 Δ Икс

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

Вы можете повторить анализ, играя с изменением координаты этих терминов вне производной, например, попробовав я 1 / 2 где это просто среднее значение в я а также я 1 .

Как выбрать, что использовать?

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

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

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

В то время, когда он был представлен, он получил невероятно точные результаты. Результаты, которые были сопоставимы с полными результатами Навье-Стокса, несмотря на использование уравнений потенциала, которые не содержат вязкости. Они обнаружили свою ошибку и опубликовали новую статью, но результаты оказались намного «хуже» по сравнению с исходной схемой. Оказывается, что несохранение ввело искусственную вязкость, заставив уравнения вести себя больше как уравнения Навье-Стокса за крошечную часть стоимости.

Излишне говорить, что инженерам это понравилось. «Лучшие» результаты при значительно меньших затратах!

Просто и по делу @ tpg2114!... Большое спасибо, что нашли время написать такое объяснение. Вы точно дали мне то, что я искал. Ваше здоровье
@ user2536125 Рад, что смог помочь. Это был вопрос на моих квалификационных экзаменах на докторскую степень несколько лет назад :) Не забудьте принять этот ответ, если вы удовлетворены тем, что он отвечает на ваш вопрос. Вы просто нажимаете галочку слева под стрелками голосования!
Согласен, хороший и простой ответ на то, что быстро становится очень сложным.

Вы показываете уравнения Эйлера, сокращенные формы уравнений Навье-Стокса, которые являются законами сохранения массы, импульса и энергии наряду с гипотезой Стокса. Итак, консервативный/неконсервативный не имеет ничего общего с законами сохранения.

В консервативной форме вы можете сразу интегрировать производные один раз на контрольном объеме и сохранить величины потока через контрольные поверхности (т. е. метод конечных объемов). Это консервативная и неконсервативная форма (поток).

В алгоритмах CFD форма имеет важное значение для размещения ударных волн и скоростей распространения. Для сжимаемых потоков вы хотите использовать консервативные формы (где у вас есть толчки).

Спасибо за ответ. Было бы очень полезно, если бы вы могли быть более щедрыми в своем последнем комментарии об алгоритмах CFD. Не могли бы вы указать, почему можно использовать консервативные формы для сжимаемых течений?
Да, но мне нужно будет сделать это, когда у меня будет больше времени, чтобы сделать это правильно.
@user2536125 user2536125 Переменные примитивного потока прерывисто изменяются в ударной волне, поэтому оператор градиента, используемый в форме без сохранения, теряет всякое физическое (и вычислительное) значение. Проще говоря, мы не можем решить матричное уравнение, в котором есть бесконечные значения. Форма сохранения прекрасно решает эту проблему, потому что переменные сохранения непрерывны на всем протяжении ударной волны, даже если давление, температура и плотность непрерывны.

Одна важная причина (возможно, не единственная) для обозначения «консервативной» и «неконсервативной» форм уравнений сохранения связана с их численным решением. При записи в консервативной форме после дискретизации дифференциального уравнения численным методом, таким как метод конечных объемов, результирующее алгебраическое уравнение по-прежнему придерживается принципа сохранения, независимо от размера конечного объема. Другими словами, в алгебраическом уравнении все еще можно увидеть баланс сохраняющейся переменной (массы, импульса, концентрации частиц). С другой стороны, когда уравнение в неконсервативной форме дискретизируется, скажем, методом конечных разностей или методом конечных элементов, сохранение гарантируется только при уточнении числовой сетки. В этом случае баланс нельзя увидеть в алгебраическом, дискретное уравнение. Следовательно, эта метка имеет смысл только при работе с уравнениями сохранения. Бессмысленно пытаться написать дифференциальное уравнение, не имеющее отношения к принципу сохранения, в консервативной или неконсервативной форме.

Разве это не то, что tpq2114 написал в своем ответе?

Все ответы хороши, но те не отвечают на вопрос. Вопрос очень простой и гласит:

Вопрос: когда математическая модель называется консервативной, а когда неконсервативной?

О: математическая модель называется консервативной, когда она непосредственно вытекает из сохранения математических величин в пределах контрольного объема, т.е. из закона сохранения.

Более подробно, если помимо простого сохранения (т. е. сохранения массы или импульса) мы также применяем цепное правило , то мы получаем неконсервативную математическую модель, названную так потому, что она больше не описывает сохранение величины в контрольном объеме, поскольку она не вытекает непосредственно из закона сохранения. Но разве эти две формулировки не эквивалентны? НЕТ!!! Они эквивалентны, только если решение непрерывно. Мы не можем применить цепное правило, если решение не является непрерывным! Следовательно, неконсервативная форма не может правильно описать распространение разрывов (т.е. скачков). См . книгу Торо для более подробной информации.

Резюмируя: речь идет о математической модели, а не о численной дискретизации. Очевидно, что существует сильная связь между числовой и математической моделью, и все остальные ответы великолепны. Но они не отвечают на вопрос.