Вопрос по граничным условиям на шахматной сетке.

[Shadow-2005]

При проведении расчета течения в плоской каверне с движущейся верхней крышкой (скорость стенки задана) используется метод конечных объемов на шахматной сетке (давление вычисляется в центрах контрольных объемов, а компоненты скорости на гранях). Трудности возникают при аппроксимации граничного условия (движущейся верхней стенки) для касательной составляющей скорости. Ближайшая к стенке расчетная точка отстоит от стенки (движущейся с заданной скоростью) на половину шага и для расчета в ней скорости необходима некая интерполяция. Как ее провести?

[TOTAL]

Жидкость несжимаемая?
Уравнение неразрывности как аппроксимируете? Уравнения движения как?

[Shadow-2005]

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

[Shadow-2005]

up!

[TOTAL]

Shadow-2005 в сообщении #269092 писал(а):

up!

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

[Shadow-2005]

Решается система из уравнений Навье-Стокса в проекции на две оси и уравнения неразрывности. Для аппроксимации диффузионных членов используется центрально-разностная схема, а для конвективных - схема с разностями против тока с законом пятой степени Патанкара. Применение шахматной сетки необходимо для устранения нефизичных полей давления. Общий алгортим следующий
1. Задание поля давления.
2. Решение уравнений движения для нахождения компонент скорости (по осям).
3. Решение уравнения для поправки давления, получаемого из уравнения неразрывности.
4. Расчет нового поля давления путем добавления к принятому полю рассчитанной в п. 3 поправки.
5. Повтор алгоритма до достижения критерия сходимости.
Это краткое описание алгоритма. Если есть еще вопросы с удовольствием отвечу. Для меня этот вопрос крайне важен.

[TOTAL]

Shadow-2005 в сообщении #269098 писал(а):

2. Решение уравнений движения для нахождения компонент скорости (по осям).
3. Решение уравнения для поправки давления, получаемого из уравнения неразрывности.

На каком этапе требуется скорость на движущейся крышке?
Зачем граничное условие, про которое спрашиваете?
Продольная скорость не хранится в точках крышки, так зачем эту скорость находить?
Откуда в уравнении неразрывности давление?

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

[Shadow-2005]

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

[Shadow-2005]

up

[TOTAL]

Shadow-2005 в сообщении #269113 писал(а):

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

В эту точку не нужна интерполяция, в ней решается уравнение.

[Shadow-2005]

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

[TOTAL]

Всё, я прекращаю эту тему, т.к. Вы не сказали точно, в чем проблема, разговор получается непонятно о чем.