#include "udf.h"
DEFINE_UDS_UNSTEADY(MyUnsteady,c,thread,i,apu,su)
{
real physical_dt, vol, phi_old;
physical_dt = RP_Get_Real("physical-time-step");
vol = C_VOLUME(c,thread);
*apu = -vol / physical_dt;         /*implicit part*/
phi_old = C_STORAGE_R(c,thread,SV_UDSI_M1(i));
*su = vol*phi_old/physical_dt;           /*explicit part*/
}
DEFINE_UDS_FLUX(MyFlux,f,thread,i)
{
real NV_VEC(unit_vec), NV_VEC(A);  //声明矢量变量
F_AREA(A, f, thread);
NV_DS(unit_vec, =, 1, 1, 1, *, 1);      //单位矢量赋值
return NV_DOT(unit_vec, A);        //矢量点积
}
DEFINE_DIFFUSIVITY(MyDiff,c,thread,i)
{
return 1.6667e-3;
}

DEFINE_SOURCE(uds_source,c,thread,dS,eqn)
{
  
  dS[eqn]=0.0;
  return  0.0;
}
DEFINE_PROFILE(density,thread,position)
{
real x[ND_ND];
real y,t,d;
cell_t c;
begin_c_loop(c,thread)
{
C_CENTROID(x,c,thread);
y=2*x[1];
t=RP_Get_Real("physical-time-step");
  d=C_YI(c,thread,1);
C_PROFILE(c,thread,position)=(d*9.0*10E-8)/(y-10E-7*t);
}
end_c_loop(c,thread)
}
DEFINE_PROFILE(xshear_stress,thread,position)
{
face_t f;
cell_t c;
Thread *tc;
begin_f_loop(f,thread)
{
c=F_C0(f,thread);
tc=THREAD_T0(thread);
F_PROFILE(f,thread,position)=-10.0*C_UDMI(c,tc,0);
}
end_f_loop(f,thread)
}
DEFINE_PROFILE(yshear_stress,thread,position)
{
face_t f;
cell_t c;
Thread *tc;
begin_f_loop(f,thread)
{
c=F_C0(f,thread);
tc=THREAD_T0(thread);
F_PROFILE(f,thread,position)=-10.0*C_UDMI(c,tc,1);
}
end_f_loop(f,thread)
}
DEFINE_SOURCE(ymomentum_source,c,thread,dS,eqn)
{
  real source,tem;
C_UDMI(c,thread,0)=C_T_G(c,thread)[0];
C_UDMI(c,thread,1)=C_T_G(c,thread)[1];

  tem=C_T(c,thread);                                         
  source=-0.4*tem;
  dS[eqn]=0.0;
  return  source;
}

20楼: Originally posted by gaojinglong at 2012-03-30 10:47:56:
我来是想请问大牛们一个关于UDF的问题的，请问大侠们，fluent的UDF能不能导入自己的本构模型？谢谢啊

3楼: Originally posted by kylafree at 2012-03-16 16:26:36:
C_UDMI(c,thread,0)=C_T_G(c,thread);
C_UDMI(c,thread,1)=C_T_G(c,thread);
注意这两行代码，里面用到了梯度宏，请查有关梯度宏的HELP

4楼: Originally posted by ya634984094 at 2012-03-16 16:36:53:
help里写的不是很清楚，问了师姐，说可以用在定义源项里~求指教

5楼: Originally posted by kylafree at 2012-03-16 16:43:49:
不好意思，我看错了。UDM的个数开足了没？

6楼: Originally posted by kylafree at 2012-03-16 16:51:54:
你的DEFINE_PROFILE宏里面用了这个语句： F_PROFILE(f,thread,position)=-10.0*C_UDMI(c,tc,0); 而C_UDMI(c,tc,0)是在DEFINE_SOURCE里面定义的。如果DEFINE_PROFILE宏先于DEFINE_SOURCE宏被调用的话，那么C_UDMI( ...