FLUENT学习udf编程实例

{ /*Argument types face_t f Thread *d int nband int n_a int n_b real ray_direction real ray_direction real en int *internal_reflection real *specular_reflectivity real *specular_transmissivity */ real angle,cos_theta; real PI=3.141592; cos_theta=Nv_DOT(ray_direction,en) if(angle>45&&angle<60) { *specular_reflectivity=0.3; *specular_transmissivity=0.7; }
real time_step; real flow_time=CURRENT_TIME; if(flow_time<0.5)
time_step=0.1; else
time_step=0.2; return time_step; }
/********************************************************* 使用一个自定义标量计算空气平均龄的 UDF *******************************************************/ #include"udf.h"
{real t=RP_Get_Real("flow-time"); int n; float c; n=t/25; c=t-n*25; omega[0]=0; omega[1]=0; omega[2]=0; if (c>=0&&c<=10)
vel[0]=0; vel[1]=0; vel[2]=0.2; else if (c>10.0&&c<=15.0)
if(temp>288.) mu_lam=5.5e-3; esle if(temp>286.) mu_lam=143.2135-0.49725*temp; else mu_lam=1.; return mu_lam; }
#include"udf.h" #include"math.h" DEFINE_CG_MOTION(velocity,dt,vel,omega,time,dtime)
/*散发增加 5%*/ *emission*=1.05;
}
/*UDF，输出半透明壁面上的扩散的反射率和扩散的传播率 */
#include"udf.h"
DEFINE_DOM_DEFFUSE_REFLECTIVITY(user_dom_diff_refl,t,nband,n_a,n_b, diff_ref_a,diff_tran_a,diff_ref_b,diff_tran_b)
/*寻找关于气体的脉络（初步） */ /*寻找关于固体的脉络（第二相） */ thread_g=THREAD_SUB_THREAD(mix_thread,s_col);/*气体相*/ thread_s=THREAD_SUB_THREAD(mix_thread,f_col);/*固体相*/
DEFINE_DIFFUSIVITY(mean_age_diff,c,t,i) {
return C_R(c,t)*2.88e-05+C_MU_EFF(c,t)/0.7; }
/*UDF，改变 DO 模型中的散发项*/
#include"duf.h"
DEFINE_DOM_SOURCE(user_dom_source,c,t,ni,nb,emission,in_scattering,abs_coeff,sc at_coeff) {
}
#include"udf.h" DEFINE_EXECUTE_AT_END(execute_at_end) { Domain *d; Thread *t; /*对耗散进行积分。*/ real sum_diss=0; cell_t c; d=Get_Domain(1);/*若是多相，则是混合域*/
thread_loop_c(t,d) { if(FLUID_THREAD_P(t)) {begin_c_loop(c,t) sum_diss+=C_D(c,t)*C_VOLUME(C,T); end_c_loop(c,t) } } printf("Volume integral of turbulent dissipation:%g\n",sum_diss); fflush(stdout);
}
/************************************************************** UDF：在 FLUENT 中，用户化缺省的 Syamlal 阻力法则 **************************************************************/ #include"udf.h"
static real v_prev=0.0;
DEFINE_CG_MOTION(piston,dt,vel,omega,time,dtime) {
Thread *t; face_t=f; real NV_VEC(A); real force,dv;
/*重新设置速度*/ NV_S(vel,=,0.0); NV_s(omega,=,0.0);
vel[0]=-0.1; vel[1]=0; vel[2]=0; else vel[0]=0; vel[1]=0; vel[2]=-0.2; }
/******************************************************** 计算气浊率的 UDF *******************************************************/
#define Pi4.*atan(1.) #define diam2 3.e-4 DEFINE_EXCHANGE_PROPERTY(custom_drag,cell,mix_thread,s_col,f_col)
{ Thread *thread_g,*thread_s; real x_vel_g,x_vel_s,y_vel_g,y_vel_s,abs_v,slip_x,slip_y, rho_g,rho_s,mu_g,reyp,afac, bfac,void,vfac,fdrgs,taup,k_g_s;
/*设置速度分量*/ vel[0]=v_prev;
} /********************************************************* 改变关于一个依赖于时间的解的时间步长的 UDF *********************************************************/ #include"udf.h" DEFINE_DELTAT(mydeltat,domain) {
if(!Data_Valid_p()) return;
/*获得脉络指示器，这个运动关于这个指示器而被定义 */ t=DT_THREAD(dt);0.0; begin_f_loop(f,t) {
F_AREA(A,f,t); Force+=F_P(f,t)*NV_MAG(A); } end_f_loop(f,t) /*计算速度中的变化，如 dv=F*dt/mass,使用显得的 Euler 公式更新速度*/ dv=dtime*force/50.0; v_prev+=dv; message("time=%f,x_vel=%f,force=%f\n"time,v_prev,Force);
{ /*Argument types
Thread *t int nband int n_a int n_b real *diff_ref_a //边 a 上得扩散的反射率 real *diff_tran_a //边 a 上得扩散的传播率 real *diff_ref_b //边 b 上得扩散的反射率 real *diff_tran_b //边 b 上得扩散的传播率 */ printf("diff_ref_a=%f diff_tran_a=%f\n,"*diff_ref_a,*diff_tran_a); printf("diff_ref_b=%f diff_tran_b=%f\n,"*diff_ref_b,*diff_tran_b);
}
/*UDF，在半透明的壁面上，沿面 f 的方向改变镜子的反射率和传播率 */ #include"udf.h"
DEFINE_DOM_SPECULAR_REFLECTIVITY(user_dom_spec_refl,f,t,nband,n_a,n_b, ray_direction,en,internal_reflection, specular_reflectivity,specular_transmissivity)
DEFINE_CHEM_STEP(user_chem_step,ifail,n,dt,p,temp,yk) {
int i; real c=1./n_spe,decay=exp(-dt); for(i=0;i<spe;i++)
yk[i]=(yk[i]-c)*decay+c; }
/************************************************************************ *运动方程的一个自由度（X 方向） *汇编的 UDF ***********************************************************************/ #include"udf.h"