[Gmsh] Not getting a volume mesh

[Ruth V. Sabariego]

Hi Wayne,
For that you have to reclasify the elements of your geometry and create the volume with the new surfaces.

Try this:
**********************
Mesh.RemeshAlgorithm=1; //(0) nosplit (1) automatic (2) split only with metis

Merge "try1-geom.gmsh";
RefineMesh;
CreateTopology;

ll[] = Line "*";
For j In {0 : #ll[]-1}
  Compound Line(newl) = ll[j];
EndFor

ss[] = Surface "*";
For j In {0 : #ss[]-1}
  ssn[] += news ;
  Compound Surface(news) = ss[j];
EndFor


Surface Loop(10) = { ssn[] };
Volume(11) = {10};

Physical Volume("Cylinder.2") = {11};
**********************


Regards,
Ruth

--
Dr. Ir. Ruth V. Sabariego
University of Liege, Electrical Engineering & Computer Science, 
Applied & Computational Electromagnetics (ACE),
phone: +32-4-3663737 - fax: +32-4-3662910 - http://ace.montefiore.ulg.ac.be/



On 24 Mar 2012, at 23:41, Wayne Christopher wrote:

> Hi Ruth,
> 
> Thanks - that fixes my volume problem.
> 
> My other problem remains, though.  I want to define geometry using a msh file, and then create a mesh with new surface elements, not the ones I am using for geometry surfaces.    Is there a way to do this?
> 
>     Wayne
> 
> On 3/23/12 2:09 PM, Ruth V. Sabariego wrote:
>> 
>> Hi Wayne, 
>> 
>> You need to define the geometrical Volume:
>> Volume(11) = {10};
>> the Physical Volume is used just for tagging the elements and their further use in a e.g. FE solver.
>> 
>> Regards,
>> Ruth
>> 
>> --
>> Dr. Ir. Ruth V. Sabariego
>> University of Liege, Electrical Engineering & Computer Science, 
>> Applied & Computational Electromagnetics (ACE),
>> phone: +32-4-3663737 - fax: +32-4-3662910 - http://ace.montefiore.ulg.ac.be/
>> 
>> 
>> 
>> On 23 Mar 2012, at 18:41, Wayne Christopher wrote:
>> 
>>> I'm trying to mesh a volume defined by surfaces that are defined by tessellations.  I don't want the mesh to be based on the surface triangles - those are purely geometry and in most cases finer than the volume mesh that I want.
>>> 
>>> When I run gmsh like this:
>>> 
>>> gmsh try1.geo -3 -format msh -saveall -smooth 3 -o try1-out.gmsh
>>> 
>>> (files attached) I do not get any tetras in the output, and it appears that the output mesh is a refinement of the input mesh.
>>> 
>>> What am I doing wrong?  Thanks,
>>> 
>>>    Wayne
>>> 
>>> <try1.geo><try1-geom.gmsh><try1-out.gmsh>_______________________________________________
>>> gmsh mailing list
>>> gmsh at geuz.org
>>> http://www.geuz.org/mailman/listinfo/gmsh
>> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20120325/384fca97/attachment.html>