Метод конечных объемов

[larts]

Решаю двумерное уравнение Навье-Стокса, на равномерной структурированной сетке, обе проекции скорости и давление беру в одной точке, то есть не использую отдельную сетку для давления. Все вроде работает, но результаты выдает с небольшими колебаниями, вроде, 50-52-51-55-53-58 т.д. В шахмотном порядке. Кто разбирается, подскажите, пожалуйста, это особенность такой разностной схемы или я где-то накосячил?

[Pphantom]

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

[larts]

А, так нет, колебания не по времени, а по пространству..

[Pphantom]

Тогда хуже. Потоки при этом оказываются нулевыми?

[larts]

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

[Pphantom]

МКО основан на идеологии вычисления изменения расчетных величин в объемах (ячейках) за счет источников/стоков и потоков из соседних ячеек. Соответственно, если источниковых членов нет (что для большинства ячеек обычно верно), то постоянная разница чего-нибудь в соседних ячейках должна приводить к появлению потоков этого "чего-нибудь" через границы ячеек. То, что разница не пропадает, означает, что либо Вы неверно считаете потоки, либо проблемы появляются где-то на стадии вывода. Отсюда и вопрос - чему равны значения потоковых членов?

[larts]

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

[Pphantom]

Но такое "шахматное" распределение точно не меняется со временем? Потому что если максимумы сдвигаются туда-сюда, то это нормально, а вот если получается статическая картина...

[larts]

Нет, не меняется. Если получается статическая картина, то что?
Колебания возникают только когда появляются стенки. Если просто взять область, поставить гран. условие "скорость справа 10 и слева 10", решение найдется абсолютно правильно.
Я подозреваю, у Вас нет экстрасенсорных способностей, но, возможно, Вы уже сталкивались с такой ситуацией и знаете, в чем может быть проблема?

[Pphantom]

larts в сообщении #851035 писал(а):

Нет, не меняется. Если получается статическая картина, то что?

То это в общем случае означает какой-то ляп в программировании. Ненулевые потоки означают ненулевые производные расчетных величин по времени.

larts в сообщении #851035 писал(а):

Колебания возникают только когда появляются стенки.

А, это уже интереснее. Как Вы их ставите?

larts в сообщении #851035 писал(а):

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

Один вариант я уже предложил, другие не попадались.

[TOTAL]

Pphantom в сообщении #851228 писал(а):

larts в сообщении #851035 писал(а):

Нет, не меняется. Если получается статическая картина, то что?

То это в общем случае означает какой-то ляп в программировании. Ненулевые потоки означают ненулевые производные расчетных величин по времени.

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

[larts]

Pphantom в сообщении #851228 писал(а):

А, это уже интереснее. Как Вы их ставите?

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

[Pphantom]

TOTAL в сообщении #851234 писал(а):

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

В стационарном в общем случае - да, но тут, судя по описанию, все проще.

larts в сообщении #851337 писал(а):

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

Да, тогда такое я когда-то видел, правда, не у себя. Автор, кажется, объяснял это получающимися артефактными "стоячими волнами", поскольку при этом стенки фактически работают как отражатели. Вязкости у Вас ведь нет?

[larts]

Вязкость есть. Сейчас попробую переделать гран. условия, посмотрю, что получится.

[larts]

Pphantom в сообщении #850963 писал(а):

Тогда хуже. Потоки при этом оказываются нулевыми?

Всем привет, еще раз. Столкнулся с еще одной проблемой, ситуация похожая. В случае, где скорость должна быть одинаковой во всех ячейках, 1.0, она распределяется полосами, 2.0-0.0-2.0 и т.д. Потоки нулевые. Подскажите, пожалуйста, как это можно исправить.