UDF 设置出口压力边界条件

ababyfriend

出口压力设为p(x,t)=z*q(x,t),其中z暂设为常值，q是体积流量
/**************************************************************************************/
/***   unsteady flow                                                                                       ***/
/***   UDF for specifying a transient pressure profile boundary condition    ***/
/**************************************************************************************/
# include "udf.h"

DEFINE_PROFILE(unsteady_pressure2,thread,index)

{
  face_t f;

  begin_f_loop(f,thread)
    {
       F_PROFILE(f,thread,index)=127323960*F_FLUX(f,thread)/1060;        /*z=127323960*/
     }
  end_f_loop(f,thread)

}
监控显示此出口压力一直是零，请教各位大神是在怎么回事？

fty0083

你加一行在F_PROFILE下面

Message("FLUX=%g\n",F_FLUX(f,thread));

然后只算一步，你看下输出的数值是不是全是0，如果全是0，那么计算没问题。如果不是0，那么你监视的压力可能不是你要的压力。而且，压力分为静压、动压、全压。你区分下。

ababyfriend

F_FLUX直接引用好像没有响应，我按照论坛里的例子重新计算了一下，现在的问题是想输出流量，不知用什么宏？
DEFINE_PROFILE(stenosis_pressure,thread,position)
{
   /*real sum_volumeflowrate=0;          /*要计算的出口面体积流量*/
   real NV_VEC(A);                     /*面积向量*/
   real NV_VEC(psi_vec);               /*速度向量*/   
   face_t f;

   begin_f_loop(f,thread)
   {

    F_AREA(A,f,thread);

    NV_D(psi_vec,=,F_U(f,thread),F_V(f,thread),F_W(f,thread));

    sum_volumeflowrate+=NV_DOT(psi_vec,A);

    F_PROFILE(f,thread,position)=800000000*sum_volumeflowrate;

    }
   end_f_loop(f,thread)


}
DEFINE_RW_FILE(writer,fp)
{
printf("Writing UDF data to data file sv.dat\n");
fprintf(fp,"%d",sum_volumeflowrate); /* write out kount to data file */
}
这样可以输出sum_volumeflowrate吗？

迷糊xiao妞儿

你好！我是fluent初学者，请问你UDF用的什么参考书？

ababyfriend

也没有什么参考书，就是网上流传的帮助文档

fty0083

不行，sum_xxx没定义过，没法输出。

fty0083

UDF Manual + C语言基本教材（个人看的谭浩强的）