求浓度边界uds出错，高手帮帮忙~

ya634984094

DEFINE_PROFILE(density,thread,position)
{
real x[ND_ND];
real y,t,d0;
cell_t c;
face_t f;
Thread *t0;
begin_f_loop(f,thread)
{
c=F_C0(f,thread);
t0 = F_C0_THREAD(f,thread);
F_CENTROID(x,f,thread);
y=2.0*x[0];
t=RP_Get_Real("physical-time-step");

  d0=C_UDSI_M1(c,t0,0);
F_PROFILE(f,thread,position)=d0+(d0*9.0e-08*t)/(y-1.0e-07*t);
}
end_f_loop(f,thread)
}

jntan

从你贴的方程来看，你的C浓度是流场时间 t 的函数，而你的UDF里对字母 t 的定义是时间步长 delta t，并非流场的绝对时间。建议你把RP_Get_Real("physical-time-step")改成RP_Get_Real("flow-time")试试。