Численное моделирование гидродинамики в каналах. Программирование

По коду code_1.zip нужно дополнить следующее:

1) Нужно проверить метод блокировки контрольных объемов. Например, добавить в расчетную область какое-либо препятствие и посчитать течение. Поcтроить рисунки, например. в ParaView. Нужно убедиться, что в блокированных контрольных объемах скорость держится равной 0.

2) Далее нужно реализовать граничные условия на давление с постоянной мощностью потока. Здесь ориентирутесь на статью Зубкова, а также детально изучите расчет поля давления по Панаткару. На правой границе давление равно нулю, а на левой раcсчитивается из условия постоянной мощности.

3) Реализовать граничные условия для скорости на левой и правой границах. Скорость на левой и правой границах считается неизвестной, поэтому ее определять будем, исходя из уравнения неразрывности. Этот алгоритм можете взять в книге Versteeg. An Introduction to Computational Fluid Dynamics. Раздел 9.5 The constant pressure boundary condition.