高手给看看这个udf

dawey

对计算域中某一点（x,y）给出导热系数的计算公式
conduct_r(x,y)=4*sigma*[(T1**2+T2**2)*(T1+T2)]/(3*beta*L)
T1=T(x-L,y)
T2=T(x+L,y)
L=0.002
beta=100.0
sigma=5.67*10**8
请问应该如何写成property的udf？
高手能不能给指点一下
关键是从温度场中取出T（x-L,y）和T（x+L,y）

qixin2000

可能要用uds才能解决

jieremy

C_T()
define_command....

沙海飞

  real a[ND_ND];
  real x;
  real y;
  face_t f;
  begin_f_loop(f, thread)
    {
      F_CENTROID(x,f,thread);
      x = a[0];
      y = a[1];
      ………………
      x和y都有了，下面就好写了
    }
  end_f_loop(f, thread)
}