CFD编程例子.doc

一、速度入口#include udf.hDEFINE_PROFILE(inlet_velocity, /* function name */ th , /* thread */ nv) /* variable number */face_t f;real xND_ND;real f_time = RP_Get_Real(flow-time;begin_f_loop (f,th) F_CENTROID(x,f,th); if (f_time=T_SAT) m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;/*如果液相单元的温度高于蒸发温度，液相向气相的质量转移*/ dSeqn = -0.1*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*定义源项对质量转移偏导*/ else m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th)/T_SAT;三、有关流沙的udfDEFINE_PROFILE(inlet_velocity,mix_thread,position) real xND_ND; /* this will hold the position vector */ real y,m; cell_t c; /water (primary 0),sand (secondary phases 1) Thread *pt; pt0=THREAD_SUB_THREAD(mix_thread,s_col);/water pt1=THREAD_SUB_THREAD(mix_thread,f_col);/sand m=C_VOF(c,pt0)+C_VOF(c,pt1); begin_c_loop(c,thread)/* loops over all faces in the thread passed in the DEFINE macro argument*/ C_CENTROID(x,c,thread); y = x1; for(y=0.5;y0) C_PROFILE(f,thread,position) = 0.5; else C_PROFILE(f,thread,position) = 0.0; end_c_loop(c,thread)四、粒径分布include udf.h DEFINE_PROFILE(particle_distribution, thread, position) real xND_ND; real y; cell_t c; begin_c_loop(c,thread) C_CENTROID(x,c,thread); y=x1; F_PROFI