udf无法使边界值随着迭代增加

ya634984094

#include "udf.h"
#define Vg 1.00000e-05
#define H 0.0200000
#define Rg 0.01

DEFINE_PROFILE(density,thread,position)
{
real k,t,a;
cell_t c0;
Thread *t0;
face_t f;
begin_f_loop(f,thread)
{
c0=F_C0(f,thread);
t0 =F_C0_THREAD(f,thread);
t=RP_Get_Real("physical-time-step");
k=Vg*t/(H-Vg*t);
a=C_UDSI_M1(c0,t0,1)+C_UDSI_M1(c0,t0,1)*(1-Rg)*k;
F_PROFILE(f,thread,position)=a;

}
end_f_loop(f,thread)
}

边界条件如下图

[ 本帖最后由 ya634984094 于 2014-1-23 19:08 编辑 ]

ya634984094

能够编译，但是所获得的的边界值一直是初始值

ya634984094

已解决