[Gmsh] outward normal of surface triangles
Negi, Ashish
ashish.negi at honeywell.com
Thu May 3 10:56:03 CEST 2012
Yes, this will solve the problem provided I have tetrahedral elements available for given Triangle. I wanted to know if GMSH provides some way to get outward normal without me finding its corresponding volume element. I haven't found anything like that in their code till now. I should better use method given by you.
Thanks a bunch,
Ashish
From: Ashish Garg [mailto:ashish.garg.iisc at gmail.com]
Sent: Thursday, May 03, 2012 1:59 PM
To: Mark Starnes
Cc: Negi, Ashish; gmsh at geuz.org
Subject: Re: [Gmsh] outward normal of surface triangles
Hello Ashish,
Ya got your problem.
See you know all tetrahedra at the boundary surfaces. It means u know all the vertex of the tetrahedra ( out of four vertex three r makin boundary faces right ) .
try this
v0 = face[i].vertex[0];
v1 = face[i].vertex[1];
v2 = face[i].vertex[2];
r01 = vertex[v1] - vertex[v0];
r12 = vertex[v2] - vertex[v1];
face[i].normal = (r01 X r12);
// Check orientation of normal
fourth_vertex = face[i].vertex[3]
check_normal_orientation = face[i].normal * ( fourth_vertex- vertex[v0]);
if ( check_normal_orientation > 0.0 )
face[i].normal *= -1.0 ;
this will resolve your problem for sure.
regards
Ashish Garg
TIFR-CAM, Bangalore
On Thu, May 3, 2012 at 1:22 PM, Mark Starnes <mark.starnes at fttech.co.uk<mailto:mark.starnes at fttech.co.uk>> wrote:
Hi Ashish,
Does it help if you change the sign of every line entry number
in the Plane Surface definition that is a problem?
Best regards,
Mark.
On Thu, May 03, 2012 at 06:23:34AM +0000, Negi, Ashish wrote:
> Hi Ashish,
>
>
>
> Thanks for replying.
>
>
>
> Let me explain my problem again. I have attached 3 pictures below. First, shows
> a Tetrahedral mesh of square duct. So, all surface elements are Triangle.
> Second, shows normal vectors computed using cross product if Triangle vertices
> are saved in clockwise sense (in the picture). Third, shows normal vector
> computed if Triangle vertices are saved in anti-clockwise sense. So, second
> picture gives me outward normal which I need.
>
>
>
> I found that all Triangles do not give me outward normals if I assume clockwise
> sense (in following picture). So, I wanted to know if GMSH has some other way
> to provide me outward normal for each surface Triangle.
>
> [cid] [cid]
>
>
>
> Thanks,
>
> Ashish
>
>
>
> From: Ashish Garg [mailto:ashish.garg.iisc at gmail.com<mailto:ashish.garg.iisc at gmail.com>]
> Sent: Wednesday, May 02, 2012 1:09 PM
> To: Negi, Ashish
> Cc: gmsh at geuz.org<mailto:gmsh at geuz.org>
> Subject: Re: [Gmsh] outward normal of surface triangles
>
>
>
> Hi Ashish
>
> give the condition that the cross product should be greater than zero. I guess
> it will resolve the problem.
>
> Ashish Garg
>
> On Wed, May 2, 2012 at 12:38 PM, Negi, Ashish <ashish.negi at honeywell.com<mailto:ashish.negi at honeywell.com>>
> wrote:
>
> Hello Everyone,
>
>
>
> I?ve been using GMSH library to open Nastran BDF file and extract mesh data for
> my software. I generally load a mesh consisting of only Tetrahedra. Hence,
> surface of model have only Triangle elements.
>
>
>
> I want to find outward normal at each Triangle. I get coordinates of each
> Triangle and then take cross product to get normal to the element. But, I?m not
> sure if it will always be outward. Can someone suggest a solution to this
> problem?
>
>
>
> Thanks,
>
> Ashish
>
>
>
>
>
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org<mailto:gmsh at geuz.org>
> http://www.geuz.org/mailman/listinfo/gmsh
>
>
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org<mailto:gmsh at geuz.org>
> http://www.geuz.org/mailman/listinfo/gmsh
--
Dr. Mark Starnes
Principal Physicist
FT Technologies
Church Lane
Teddington
TW11 8PA
web: www.fttech.co.uk<http://www.fttech.co.uk>
DDI: 0208 614 2730
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20120503/908d5eb4/attachment.html>