一个简单而完整的UDS例子.doc

--------------------------校验:_____________-----------------------日期:_____________一个简单而完整的UDS例子一个简单而完整的UDS例子gearboy标准方程形式:(默认)Fluent中各项的定义如下:名称表达式各项在UDS中需要定义的量在UDS中对应的宏非稳态项离散形式的非稳态项(详细见UDF帮助)DEFINE_UDS_UNSTEADY对流项通量DEFINE_UDS_FLUX扩散项扩散率DEFINE_DIFFUSIVITY假定需要求解方程为当,即如果计算达到稳定时,,此方程解析解为,其中c为常数(假定其为0),可由边界条件确定。因此,当解稳定时,,等值线为一系列45度的斜直线。将此方程与标准方程对比发现:名称表达式各项在UDS中需要定义的量需要在UDS中定义的宏非稳态项离散形式的非稳态项DEFINE_UDS_UNSTEADY对流项通量(单位向量)DEFINE_UDS_FLUX扩散项扩散率DEFINE_DIFFUSIVITY源项22DEFINE_SOURCE边界条件(b代表边界处)DEFINE_PROFILE*注:源项和扩散率为常数,可以不用定义宏,但为了说明问题,仍采用宏来定义。假定如图所示的计算域,其中正方形边长为10,原点在正方形中心,四边均为wall。DEFINE_UDS_UNSTEADY(MyUnsteady,c,t,i,apu,su){realphysical_dt,vol,phi_old;physical_dt=RP_Get_Real("physical-time-step");vol=C_VOLUME(c,t);*apu=-vol/physical_dt;/*implicitpart*/phi_old=C_STORAGE_R(c,t,SV_UDSI_M1(i));*su=vol*phi_old/physical_dt;/*explicitpart*/}DEFINE_UDS_FLUX(MyFlux,f,t,i){realNV_VEC(unit_vec),NV_VEC(A);//声明矢量变量F_AREA(A,f,t);NV_DS(unit_vec,=,1,1,1,*,1);//单位矢量赋值returnNV_DOT(unit_vec,A);//矢量点积}DEFINE_DIFFUSIVITY(MyDiff,c,t,i){}DEFINE_SOURCE(MySource,c,t,dS,eqn){dS[eqn]=0;}DEFINE_PROFILE(MyProfile,thread,index){realx[ND_ND