Постройте результаты моделирования конечного объема на полярной сетке (python и/или paraview)

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

Моя проблема заключается в следующем: у меня есть численный симулятор, который решает задачу потока жидкости в 2D или 3D с использованием радиальной (цилиндрической) сетки. Симулятор использует подход конечного объема, при котором одно значение каждой интересующей переменной (например, давление, концентрации) рассчитывается для каждого блока сетки.

Я хочу визуализировать профили в 2D (для случая 3D я хочу построить срезы в z-направлении) с помощью Python (matplotlib). Однако мне не нужен контурный график с гладким профилем, а скорее одно значение с цветовой кодировкой внутри каждой ячейки сетки. Я также хотел бы иметь возможность строить линии сетки.

В качестве альтернативы или в качестве дополнительного подхода, интересно, можно ли легко решить ту же проблему с помощью Paraview? То есть, как можно написать необходимые файлы .vtk, основываясь на информации об ограничивающих радиусах, азимутальных углах и координатах z ячеек моделирующей сетки?

Если кто-то может помочь мне с любым из этих вопросов, я был бы очень благодарен!


person Axiom    schedule 26.10.2016    source источник


Ответы (1)


Для ParaView самым простым решением, которое я могу придумать, является запись файла VTK. Вы можете загрузить этот файл, разрезать его в любом направлении и довольно легко просматривать данные по центру ячейки с линиями сетки и без них.

Ваш файл VTK будет содержать координаты xyz, которые вы преобразуете из полярной системы координат, в которой, как я предполагаю, они представлены. Ваш тип ячейки будет VTK_HEXAHEDRON ячеек. Похоже, вы хотите хранить свои данные как то, что VTK называет массивом данных «ячейки» (в отличие от массива данных «точек», который связан с точками сетки). Дополнительные сведения об устаревшем формате файлов VTK см. на странице http://www.vtk.org/wp-content/uploads/2015/04/file-formats.pdf.

person Cory Quammen    schedule 27.10.2016