2D треугольная сетка в Gmsh с использованием двоичного изображения в качестве входных данных - PullRequest
0 голосов
/ 11 января 2019

Мне нужно создать двумерный файл треугольной сетки из пористой среды для моделирования CFD. У меня есть двоичное PNG-изображение распределения зерен, и я хочу получить сетку для сети пор и не иметь элемента сетки в зернах (ссылка на image.png ) Это означает разбиение только подключенного домена пикселей «1» и наличие пустот вместо доменов пикселей «0».

Я успешно загрузил изображение в Gmsh, используя следующую геометрию image_mesh.geo (вдохновлено этой страницей http://www.montefiore.ulg.ac.be/~geuzaine/)

// Merge the image (this will create a new post-processing view,
// View[0])
Merge "image.png";

// Apply the view as the current background mesh
Background Mesh View[0];

// Build a simple geometry on top of the background mesh
w = View[0].MaxX;
h = View[0].MaxY;
Point(1)={0,0,0,w};
Point(2)={w,0,0,w};
Point(3)={w,h,0,w};
Point(4)={0,h,0,w};
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,1};
Line Loop(5) = {3,4,1,2};
Plane Surface(6) = {5};

Оттуда я не могу найти какой-либо вариант, чтобы превратить изображение в геометрию (двухмерную связанную плоскую область вокруг круглых пустот) и соответствующую двумерную треугольную сетку. Есть ли у вас какие-либо предложения о том, как это сделать с Gmsh? Большое спасибо !!

...