|
[Sponsors] |
March 4, 2010, 16:41 |
URGENT!!! F_UDSI at symmetry boundary
|
#1 | |
Member
Join Date: Mar 2009
Location: Istanbul, Turkiye
Posts: 47
Rep Power: 17 |
Hi,
In FLUENT udf manual following statement is given, Quote:
Then, I tried to test this in a simple 3D uds diffusion problem where the domain has symmetry boundary condition. After convergence, I applied following function to printout the face values of uds at all boundaries: Code:
face_t f; Thread *tf; thread_loop_f(tf, domain)/* loops over all face threads in a domain*/ { begin_f_loop(f, tf) /* loops over faces in a face thread */ { if (BOUNDARY_FACE_THREAD_P(tf)) { Message("face id = %d, UDS0 value = %12.5e",f, F_UDSI(f, tf, 0)); } } end_f_loop(f, f_thread) } thanks [/COLOR] Last edited by gemini; March 5, 2010 at 03:34. |
||
March 8, 2010, 05:51 |
|
#2 | |
Senior Member
Join Date: Feb 2010
Posts: 164
Rep Power: 17 |
Quote:
#include "udf.h" DEFINE_ON_DEMAND(test) { face_t f; cell_t c0; Thread *tf,*t0; Domain*domain=Get_Domain(1); real NV_VEC(dist_vector); real f_centoid[ND_ND],c_centroid[ND_ND]; real face_value,cell_value; thread_loop_f(tf, domain)/* loops over all face threads in a domain*/ { if (THREAD_F_AXIS==THREAD_TYPE(tf)) /* if the thread is an axis thread*/ { begin_f_loop(f, tf) /* loops over faces in a face thread */ { F_CENTROID(f_centoid,f,tf); //centroid of a face c0=F_C0(f,tf); t0=THREAD_T0(tf); C_CENTROID(c_centroid,c0,t0); //centroid of the neighbour cell NV_DD(dist_vector,=,f_centoid[0],f_centoid[1],0,-,c_centroid[0],c_centroid[1],0); face_value=C_UDSI(c0,t0,0)+NV_DOT(C_UDSI_G(c0,t0,0 ),dist_vector); cell_value=C_UDSI(c0,t0,0); Message0("x coord.=%g,face value=%g,cell value=%g\n",f_centoid[0],face_value,cell_value); } end_f_loop(f, f_thread) } } } |
||
March 10, 2010, 04:20 |
|
#3 | |
Member
Join Date: Mar 2009
Location: Istanbul, Turkiye
Posts: 47
Rep Power: 17 |
Quote:
Can your script be used for calculation of face value at the interior faces by just removing the following line? if (THREAD_F_AXIS==THREAD_TYPE(tf)) /* if the thread is an axis thread*/ Also, THREAD_F_AXIS didn't worked for SYMMETRY bc's. Instead "THREAD_F_SYMMETRIC" works well... Thanks. Last edited by gemini; March 10, 2010 at 05:00. |
||
March 10, 2010, 22:50 |
|
#4 | |
Senior Member
Join Date: Feb 2010
Posts: 164
Rep Power: 17 |
Quote:
|
||
March 11, 2010, 02:23 |
|
#5 | |
Member
Join Date: Mar 2009
Location: Istanbul, Turkiye
Posts: 47
Rep Power: 17 |
Quote:
Your script was very useful, and as you said, it can be used for calculating face fluxes of UDS and then they can be used anywhere in your model. thanks again. |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Implementation of boundary conditions for FVM | Tom | Main CFD Forum | 7 | August 26, 2014 06:58 |
Symmetry and boundary conditions | eric | FLUENT | 3 | August 30, 2012 10:09 |
inlet velocity boundary condition | murali | CFX | 5 | August 3, 2012 09:56 |
Normal velocity at the symmetry boundary | Seventy | FLUENT | 1 | October 26, 2009 07:52 |
Boundary conditions? | Tom | Main CFD Forum | 0 | November 5, 2002 02:54 |