[Gmsh] Problem creating mesh

[Cian Wilson]

Dear Penny,

gmsh2triangle is a script that comes with fluidity so I'm cc'ing the 
fludiity e-mail list on this reply.

In this case gmsh2triangle is failing because of a problem in your .geo 
file (so it's not really a problem with gmsh2triangle or fluidity).  You 
should replace:

Physical Surface(20) = {13,14,15};

with:

Physical Surface(20) = {16,17,18};

as the ids in the expression-list on the lhs should be to the elementary 
surfaces not the line loops (see 
http://www.manpagez.com/info/gmsh/gmsh-2.4.0/gmsh_34.php).

I also believe your:

Line Loop(19) = {1,2,3,-12,8,9,10,6,-11,5};

is superfluous.  Removing it and making the changes above both gmsh and 
gmsh2triangle work for me locally.

Hope this helps,
Cian

--
Cian Wilson
Postdoctoral Research Scientist
Lamont-Doherty Earth Observatory, Columbia University

On 06/25/2012 05:53 AM, Jackson Penelope A wrote:
>
> I am trying to do a simple example for Fluidity which first requires 
> me to set up a mesh in Gmsh, but it's not working and I can't see why. 
> I am using the following .geo file:
>
> Point(1) = {5.0, 0.0, 0.0, 0.1};
>
> Point(2) = {35.0, 0.0, 0.0, 0.1};
>
> Point(3) = {35.0, 2.0, 0.0, 0.1};
>
> Point(4) = {30.0, 2.0, 0.0, 0.05};
>
> Point(5) = {30.0, 4.05, 0.0, 0.1};
>
> Point(6) = {0.0, 4.05, 0.0, 0.1};
>
> Point(7) = {0.0, 2.05, 0.0, 0.1};
>
> Point(8) = {5.0, 2.0, 0.0, 0.05};
>
> Point(9) = {5.0, 2.05, 0.0, 0.05};
>
> Point(10) = {30.0, 2.05, 0.0, 0.05};
>
> Line(1) = {1,2};
>
> Line(2) = {2,3};
>
> Line(3) = {3,4};
>
> Line(4) = {4,8};
>
> Line(5) = {8,1};
>
> Line(6) = {7,9};
>
> Line(7) = {9,10};
>
> Line(8) = {10,5};
>
> Line(9) = {5,6};
>
> Line(10) = {6,7};
>
> Line(11) = {8,9};
>
> Line(12) = {10,4};
>
> Line Loop(13) = {1,2,3,4,5};
>
> Line Loop(14) = {6,7,8,9,10};
>
> Line Loop(15) = {4,11,7,12};
>
> Plane Surface(16) = {13};
>
> Plane Surface(17) = {14};
>
> Plane Surface(18) = {15};
>
> Line Loop(19) = {1,2,3,-12,8,9,10,6,-11,5};
>
> Physical Surface(20) = {13,14,15};
>
> //Inflows
>
> Physical Line(21) = {2};
>
> Physical Line(22) = {10};
>
> //Outlets
>
> Physical Line(23) = {5};
>
> Physical Line(24) = {8};
>
> //Top and bottom
>
> Physical Line(25) = {9};
>
> Physical Line(26) = {1};
>
> //Less sure about these
>
> Physical Line(27) = {6,11};
>
> Physical Line(28) = {3,12};
>
> The commands "gmsh -2 -v 0 a2wayflow.geo" and "gmsh -2 -v 0 
> a2wayflow.geo" both seem to work fine. The problem occurs when I do 
> "gmsh2triangle -2 a2wayflow.msh" for use in Fluidity, I get the error 
> "Unable to determine dimension of problem". Using the GUI, pretty much 
> the same thing happens. Mesh generates no problem, looks OK, saves 
> without complaining, then running the gmsh2triangle command gets the 
> same error message.
>
> I've thoroughly looked at it several times and a colleague who is much 
> more experienced with it than me can't spot the problem either, can 
> you tell me where I'm going wrong? Or is the problem at the Fluidity end?
>
> Penny Jackson.
>
> Penelope Jackson
>
> Scientist
>
> *[**dstl**]*Sensors & Countermeasures Department
> Rm 101, iSAT-C, Bldg 5
> Dstl Porton Down
>
> T 9 6801 8413
>
> "This e-mail and any attachment(s) is intended for the recipient 
> only.   Its unauthorised use, disclosure, storage or copying is not 
> permitted.  Communications with Dstl are monitored and/or recorded for 
> system efficiency and other lawful purposes, including business 
> intelligence, business metrics and training.  Any views or opinions 
> expressed in this e-mail do not necessarily reflect Dstl policy."
> "If you are not the intended recipient, please remove it from your 
> system and notify the author of the email and centralenq at dstl.gov.uk 
> <mailto:centralenq at dstl.gov.uk>"
>
>
> _______________________________________________
> 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/20120625/41bde189/attachment.html>