udf--有爱帮助资料.doc

解释型速度进口 /*********************************************************************** myexample.c UDF for specifying steady-state velocity profile boundary condition ************************************************************************ #include udf.h DEFINE_PROFILE(inlet_x_velocity, thread, position) { real x[ND_ND]; /* this will hold the position vector */ real y, h; face_t f; h = 0.016; /* inlet height in m */ begin_f_loop(f,thread) { F_CENTROID(x, f, thread); y = 2.*(x[1]-0.5*h)/h; /* non-dimensional y coordinate */ F_PROFILE(f, thread, position) = 0.1*(1.0-y*y); } end_f_loop(f, thread) } 解释型周期压力出口 /********************************************************************** unsteady.c UDF for specifying a transient pressure profile boundary condition ***********************************************************************/ #include udf.h DEFINE_PROFILE(unsteady_pressure, thread, position) { face_t f; real t = CURRENT_TIME; begin_f_loop(f, thread) { F_PROFILE(f, thread, position) = 101325.0 + 5.0*sin(10.*t); } end_f_loop(f, thread) } 动量方程添加源相 /****************************************************************** UDF that adds momentum source term and derivative to duct flow *******************************************************************/ #include udf.h #define CON 20.0 DEFINE_SOURCE(cell_x_source, cell, thread, dS, eqn) { real source; if (C_T(cell,thread) = 288.) { source = -CON*C_U(cell,thread); dS[eqn] = -CON; } else { source = dS[eqn] = 0.; } return source; } 定义材料属性 /********************************************************************* UDF for specifying a temperature-dependent viscosity property **********************************************************************/ #include udf.h DEFINE_PROPERTY(cell_viscosity, cell, thread) { real mu_lam; real temp = C_T(cell, thread); if (temp 288.) mu_lam = 5.5e-3; else if (temp 286.) mu_lam = 143.2135 - 0.49725 * temp; else mu_lam = 1.;