CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

How do I compute the normal vector to a surface with FEM

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 30, 2009, 09:00
Default How do I compute the normal vector to a surface with FEM
  #1
New Member
 
Join Date: Apr 2009
Posts: 17
Rep Power: 17
Vasilis is on a distinguished road
I am trying to compute the normal vector to an inclined plate with the FEM.
Based on theory, the normal vector is given by

n=(-dz/dx,0,1)

This is the normal to the top surface. There is no variation in the y-direction. Consider it as a 2D geometry. The plate has length L, and height H, and it is inclined with a slope equal to PHI, so that it's bottom right corner is located at (a,0) and it's top left corner is located at (b,1).

If I use the FEM to compute the dz/dx derivative, I will find that it is equal to zero. Does this make sense?
The only way to compute a non-zero value for the dz/dx derivative is to solve the problem where the domain was originally a rectangular which deformed to the one I have now. Is this the only way to compute the normal vector with FEM?
Vasilis is offline   Reply With Quote

Old   July 31, 2009, 06:10
Default
  #2
New Member
 
Join Date: Apr 2009
Posts: 17
Rep Power: 17
Vasilis is on a distinguished road
I am trying to calculate the normal on the outer surface of a circle. The circle lies in the y-z plane.

I compute the value of the derivative dz/dy (which is equal to (y-y0)/(z-z0), if the equation for the circle is (z-z0)^2+(y-y0)^2=R^2 ), and there is no good agreement between FEM and theory.
Any ideas why this is happening?
Vasilis is offline   Reply With Quote

Old   July 31, 2009, 07:37
Default
  #3
New Member
 
Join Date: Apr 2009
Posts: 17
Rep Power: 17
Vasilis is on a distinguished road
My mistake, I can not compute the normal vectors because both x and z are independent variables.
Vasilis is offline   Reply With Quote

Old   July 31, 2009, 08:59
Default
  #4
New Member
 
Join Date: Apr 2009
Posts: 17
Rep Power: 17
Vasilis is on a distinguished road
So, how do I calculate the normal vectors in a moving boundary domain.
If we assume that x=f(X,Z) and z=h(X,Z) taking the derivative of dx/dZ will give me the rate of deformation. But, we need to compute quantities such as dx/dz.

Does anybody has any idea about how to compute the normal vectors with FEM?
Vasilis is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to get the surface normal vector ohrmond CFX 3 February 1, 2018 06:45
[Gmsh] Error : Self intersecting surface mesh, computing intersections & Error : Impossible velan OpenFOAM Meshing & Mesh Conversion 3 October 22, 2015 12:05
[Gmsh] Problem with Gmsh nishant_hull OpenFOAM Meshing & Mesh Conversion 23 August 5, 2015 03:09
NACA0012 geometry/design software needed Franny Main CFD Forum 13 July 7, 2007 16:57
CFX4.3 -build analysis form Chie Min CFX 5 July 13, 2001 00:19


All times are GMT -4. The time now is 12:02.