Непонятный вывод

Question

Я был бы признателен за помощь, чтобы прояснить природу моих результатов моделирования. Это моделирование рассматривает обтекание куба внутри цилиндрического канала в 3D. Математической моделью являются уравнения Навье-Стокса и уравнение неразрывности (div U = 0). Я разработал код для двумерного численного моделирования, и результат верен, это наводит меня на мысль, что вариационная проблема, которая одинакова для двумерных и трехмерных задач (поправьте меня, если я ошибаюсь), хорошо сформулирована в FEniCS. . Когда я пытаюсь запустить моделирование трехмерной геометрии, я получаю следующее () как в терминах полей давления и скорости (u,p).

Я хочу объяснить, как я управлял геометрией в несколько шагов:

1. Создайте его в FreeCAD

2. Импортировал в ГМШ

3. Определите границы с помощью команды add-> поверхности и пометьте их цифрами от 1 до 4
   3.1) Стена, которая соответствует тегу '1', представляет поверхность, связанную с цилиндром
   3.2) Obstalce, который соответствует тегу '2', представляет внутренний куб поверхностей
   3.3) входа (выхода) , которые соответствуют тегу «3» ( «4») представляют собой впускной и выпускной поверхностях , связанных с цилиндром

4. Определите объем с помощью add-> volume и пометьте его 5

5. Преобразуйте файл .msh в файл .xml, который создает два дополнительных файла: physical_region.xml и facet_region.xml. Эти два файла должны автоматически содержать информацию о границах. Затем реализуется следующая строка для разработки геометрии в python3 перед вычислением моделирования:

   % Прочитать сетку
   mesh = Mesh ('file.xml')
   cd=MeshFunction('size_t', mesh, 'prova_physical_region.xml')
   fd=MeshFunction('size_t', mesh, 'prova_facet_region.xml')
   bcu_inflow = DirichletBC(V , Константа ((2, 2, 0)), '3')
   bcu_wall = DirichletBC(V, Constant((0, 0, 0)), '1')
   bcu_cube = DirichletBC(V, Constant((0, 0, 0)), '2')
   bcp_outflow = Дирихле BC(Q, Константа (0), '4')

Моделирование не дает никаких ошибок и вычисляет результат, показанный в (выводвыводе), который не имеет смысла в нашем физическом контексте (гидродинамика), тогда я хотел бы знать, в чем проблема. Подчеркну, что в FEniCS сетка заряжена хорошо.