|July 27, 2009, 15:24||
UDF for inhomogenous Neumann boundary condition
Akhilesh Kumar Sahu
Join Date: Jul 2009
Posts: 1Rep Power: 0
I am new to FLUENT
I have to set du/dy=G, G=0.01,0.02. (d is partial derivative) @ some boundary.
discretized form of the above eq would be
here u_f is u velocity at face center and u_c is same at cell center of corresponding cell and del_y is distance between face center and cell center.
I have written the following udf to set this..
#define dy 0.08754886724
#define G 0.01
#define uavg 0.1
#define H 10.0
c0 = F_C0(f,thread);
t0 = THREAD_T0(thread);
no problem in interpretation..
but results are not coming fine.
I have checked the face center and cell center value of u at that boundary, they are coming same, like for symmetry boundary condition.
Please help me in finding the bug in the udf.
|June 3, 2012, 09:01||
Join Date: Mar 2009
Blog Entries: 14Rep Power: 18
DEFINE_PROFILE(wall_x_velocity, t, i)
real s, ds, A_by_es;
s=2.0; /*Specify your own value*/
F_PROFILE(f, t, i) = C_U(F_C0(f,t),THREAD_T0(t))+s*ds;
You may also want to check BOUNDARY_FACE_GEOMETRY in the manual to improve the approximations in the UDF above; indeed it is assumed that the wall parallel velocity is along the x direction. A more correct implementation should first obtain the wall parallel velocity from the full 3D one, use this one to define the velocity at the wall from the specified wall normal gradient (also accounting for grid non-orthogonalities) and then project this velocity back to the three coordinates. This UDF provides you all the ingredients to accomplish the task, you just need to put in some vector algebra.
However, in the wall panel you also have the option to directly specify the shear at the wall, making the above UDF actually not required.
|Thread||Thread Starter||Forum||Replies||Last Post|
|CFX Solver : Sudden crash||Herv√©||CFX||2||June 16, 2008 06:40|
|UDF for pressure oulet boundary condition||Jubs||FLUENT||0||February 13, 2007 01:26|
|NEUMANN BOUNDARY CONDITION IN MESHLESS METHODS||K√©melli||Main CFD Forum||2||July 15, 2006 01:41|
|udf boundary condition||stefan||FLUENT||0||December 4, 2002 10:10|
|UDF for boundary condition||Shankar||FLUENT||4||February 26, 2002 21:00|