# Calculating the centre of a surface

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

 October 17, 2012, 12:11 Calculating the centre of a surface #1 Member   Join Date: Sep 2012 Location: FL Posts: 76 Rep Power: 4 I need to find centroid of multiple boundaries to implement in equation in UDF. I want Fluent to calculate the coordinates (x,z) of the centroid. I have written a code for area weighted centroid. I am getting an error. Can someone tell me what am I doing wrong? ( I think its with the A(). I need to calculate area in xz planes. Is it right?) /* Calculating the centroid of the Boundary*/ begin_f_loop(f, thread) { F_CENTROID(p,f,thread); F_AREA(A,f,thread); a = (p[0]*A[1])+a; b= (p[2]*A[1])+b; } end_f_loop(f, thread) begin_f_loop(f, thread) { F_AREA(A,f,thread); d = A[1] + d; } end_f_loop(f, thread) k=a/d;l=b/d; Last edited by victoryv; October 17, 2012 at 16:59.

October 17, 2012, 19:27
#2
Senior Member

Syavash
Join Date: Apr 2010
Location: Tehran, Karaj
Posts: 216
Rep Power: 8
Quote:
 Originally Posted by victoryv I need to find centroid of multiple boundaries to implement in equation in UDF. I want Fluent to calculate the coordinates (x,z) of the centroid. I have written a code for area weighted centroid. I am getting an error. Can someone tell me what am I doing wrong? ( I think its with the A(). I need to calculate area in xz planes. Is it right?) /* Calculating the centroid of the Boundary*/ begin_f_loop(f, thread) { F_CENTROID(p,f,thread); F_AREA(A,f,thread); a = (p[0]*A[1])+a; b= (p[2]*A[1])+b; } end_f_loop(f, thread) begin_f_loop(f, thread) { F_AREA(A,f,thread); d = A[1] + d; } end_f_loop(f, thread) k=a/d;l=b/d;
The problem is right here, if you want to access Area magnitude of a boundary face, you have to implement this macro:
area=NV_MAG(A);
Goodluck

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post nishant_hull Open Source Meshers: Gmsh, Netgen, CGNS, ... 18 April 22, 2015 08:43 123catty456 Main CFD Forum 1 October 1, 2012 22:27 parkh32 OpenFOAM Pre-Processing 2 April 8, 2012 17:12 Sunil Gupta FLUENT 0 April 22, 2008 09:29 xinzhangabc CFX 0 February 24, 2005 07:29

All times are GMT -4. The time now is 17:04.