CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT > Fluent UDF and Scheme Programming

define boundary condition in just one cell face

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 16, 2010, 20:11
Default define boundary condition in just one cell face
  #1
New Member
 
omid
Join Date: Mar 2010
Posts: 17
Rep Power: 16
om1234 is on a distinguished road
Hi all
i'm going to analysis 2d centrifugal pump.in outlet line,i try to define static pressure in a single face of grid. i use this udf. but it doesn't work correctlly.
i have some questions and i wish any one could answer me.
is the define adjust updated in every iteration and is it usefull in this udf?
any one know any thing about the working of fluent memory?
is the memory values updated in each iteration,i mean is the previous values deleted and replaced by the new pressure value or not?


#include"udf.h"

DEFINE_ADJUST(pressure,domain)
{
real p;
face_t f;
int ID=7;/*this is the outlet ID*/
int k=0;

Thread *t=Lookup_Thread(domain,ID); /*this macro is used to define specific thread*/

begin_f_loop(f,t)
{
p=F_P(f,t); /*F_P read the pressure of each face and assign it to p*/
F_UDMI(f,t,k)= p; /*this macro is used to store/access the values of the user-defined memory on the face.the outlet boundary is devided to 12 subline.so i set the number of memory 12.define>user deine>memory*/
k=k+1;
}
end_f_loop(f,t)
}


DEFINE_PROFILE(pressure_outlet, t, i) /* this macro is used to define values in BCs*/
{
real pr;
real x[ND_ND];
real y;
face_t f;
int j=0;
if(!Data_Valid_P()) /*check for the availability of data*/
return;
begin_f_loop(f, t)
{
F_CENTROID(x,f,t)
y=x[0]; /* the outlet line is horizontal and i need the x coordinate*/
if (y>-0.005 && y<0.001) /* there is a face centroid in this interior*/
{
F_PROFILE(f, t, i) = 100000; /* the desired value of the defined face cell*/
j=j+1;
}
else
pr=F_UDMI(f,t,j); /* access to stored values for each face*/
F_PROFILE(f, t, i) = pr;
j=j+1;
}
end_f_loop(f, t)
}
best regard

Last edited by om1234; June 21, 2010 at 10:08.
om1234 is offline   Reply With Quote

Reply


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
[Gmsh] Import problem ARC OpenFOAM Meshing & Mesh Conversion 0 February 27, 2010 11:56
[GAMBIT] different boundary condition on one volume face slawek.mikula ANSYS Meshing & Geometry 8 May 2, 2009 12:23
Missing math.h header Travis FLUENT 4 January 15, 2009 12:48
boundary condition of face for liquid and gas Luke FLUENT 6 May 21, 2006 14:45
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Meshing & Mesh Conversion 2 July 15, 2005 05:15


All times are GMT -4. The time now is 11:49.