楼: Originally posted by 47425804 at 2012-03-10 22:18:47:
除了楼上两位说的问题外，我觉得还有点问题，供参考。
face变量同时用在了出口和入口，如果实际上出口和入口的网格情况完全一致，那还是可以的。
但是，如果出口网格的编号方向和入口网格的编号方向刚好是相反的 ...

楼: Originally posted by kylafree at 2012-03-10 14:08:02:
你上面那个写法是错误的，连编译都通不过，你没试着编译一下吗？

#include "udf.h"
#define INT_YI 0.03
#define UDF_FILENAME "udf_loc_velo"
real AVG_YI;


DEFINE_ADJUST(outlet_avg_yi, d)
{
FILE *fp;
real sum_yi_A=0.0;
real sum_yi=0.0;
int j=0;
face_t f;
Thread *thread_out=Lookup_Thread(d,4);
  begin_f_loop(f,thread_out)
{
  sum_yi_A+=F_YI(f,thread_out,0);
  j++;      
}
end_f_loop(f,thread_out)

  AVG_YI=sum_yi_A/j;
fp=fopen(UDF_FILENAME,"w");

fprintf(fp, "%f",sum_yi_A); /* write out kount to data file */
fprintf(fp, "%f",AVG_YI);
fprintf(fp, "%d",j);
fclose(fp);

printf("%f\n",sum_yi_A);/* write out kount to data file */
printf("%f\n",AVG_YI);
printf("%d\n",j);
}


DEFINE_PROFILE(inlet,thread_in,index)
{
    real inlet_yi=0.0;
    face_t inlet_face;
  
    real time=RP_Get_Real("flow-time");   
    begin_f_loop(inlet_face,thread_in)
    {
          if(time==0)
        {
         inlet_yi=INT_YI;
          F_PROFILE(inlet_face,thread_in,index)=inlet_yi;
        }
        
          else
        {
          inlet_yi=AVG_YI+0.2*0.5;
          F_PROFILE(inlet_face,thread_in,index)=inlet_yi;
        }
    }
        end_f_loop(inlet_face,thread_in)   

printf("%f\n",inlet_yi);                     
}

楼: Originally posted by jesseliuxl at 2012-03-13 20:06:49:
以前做的udf，出口值调节后赋给入口，可供参考下

#include "udf.h"
#define INT_YI 0.03
#define UDF_FILENAME "udf_loc_velo"
real AVG_YI;


DEFINE_ADJUST(outlet_avg_yi, d)
{ ...