[Gmsh] listing elements according to which physical volume they belong to

[Abiti Adili]

Thank you very much. Now I have the lists of the elements with the physical volume tags they belong to.

I have attached a png file of the 3D view of the mesh generated, geomtry.geo and omnibus.geo.

It seems to me that
 when I clicked on 3D in Mesh option  it seems like it started generating more densely around those two balls.

when clicked on "refine by splitting", it still gets dense around those two balls.

You can see it from the attached file of the 3D view of the mesh.

Is there any reason for this mesh to be denser around the ball and not so dense in  the complement of the balls? Or I am just judging this by my eyes while it is not actually like that? Is there any way I can confirm that it is indeed denser around the balls? and If that is true, am I doing something wrong in the geo files?

Thank you very much.

Adil


________________________________________
From: G. D. McBain <gdmcbain at protonmail.com>
Sent: Sunday, November 25, 2018 11:43:56 PM
To: Abiti Adili
Cc: gmsh at onelab.info
Subject: Re: [Gmsh] listing elements according to which physical volume they belong to

©\©\©\©\©\©\©\ Original Message ©\©\©\©\©\©\©\
Le lundi, novembre 26, 2018 12:43 PM, Abiti Adili <aadili1 at lsu.edu> a ¨¦crit :

> Hello,
>
> I have a geometry.geo file that creates a box with two small balls in it. I used omnibus.geo file that utilize the geometry.geo. Using this file, I created a 3D mesh over the box(as well as the two balls in it). I need to be able to have a mesh( created over the box) so that the gmsh file lists the tetrahedral elements according which physical volume( ball1, ball2, or the complement of them in the box) they belong to.
>
> The files I currently have are not quite generating what I wanted, it is only listing the elements belong to the complement. I would greatly appreciate if someone can take a look at them and point out the mistakes I made.

The last lines of geometry.geo have:

    Volume(1) = {2};
    Volume(2) = {3};
    Volume(3) = {1, 2};

The last should read

    Volume(3) = {1, 2, 3};

if the second sphere is also to internally bound it.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: geometry.geo
Type: application/octet-stream
Size: 3006 bytes
Desc: geometry.geo
URL: <http://onelab.info/pipermail/gmsh/attachments/20181126/117ade82/attachment-0002.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: omnibus.geo
Type: application/octet-stream
Size: 124 bytes
Desc: omnibus.geo
URL: <http://onelab.info/pipermail/gmsh/attachments/20181126/117ade82/attachment-0003.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: The_box_and_2_spheres.pdf
Type: application/pdf
Size: 1208173 bytes
Desc: The_box_and_2_spheres.pdf
URL: <http://onelab.info/pipermail/gmsh/attachments/20181126/117ade82/attachment-0001.pdf>