【求助】大家帮忙修改一下udf程序

gzxcqu

fluent模拟一个矩形池内的二元溶液流动传热过程（二维），矩形池两侧有恒定的温度和浓度，且左侧大于右侧，矩形池上界面为与空气接触的自由表面，底部绝热，各固壁满足无滑移边界条件，udf使用来处理自由表面的上的力平衡的即 μ*du/dy=-Rt*dT/dx-Rs*dS/dx,T,S分别为温度和浓度UDF通过编译，也可运行，计算结果不正确，设置监测变量后发现，我设定的边界条件没有运行，麻烦大家帮忙改一下，非常谢谢
#include "udf.h"
#include "sg.h "

#define           Rt              -8.5e-3
#define           Rs              1.7e-4  
#define   density          995.077
DEFINE_PROFILE(p1_bc, thread, position)

{
  face_t f;
  real A[ND_ND];
  real dG[ND_ND],dr0[ND_ND],es[ND_ND],ds,A_by_es;
  real Gsource;
  Thread *t0=thread->t0;

  if (!Data_Valid_P()||!FLUID_THREAD_P(t0))
          return;

  begin_f_loop (f,thread)
    {
      cell_t c0 = F_C0(f,thread);
      BOUNDARY_FACE_GEOMETRY(f,thread,A,ds,es,A_by_es,dr0);
       if (NULLP(T_STORAGE_R_NV(t0,SV_DENSITY_G)))
        Gsource = 0.;        
      else     
        Gsource =-(Rt*C_T_G(c0,t0)[0]+Rs*density*C_YI_G(c0,t0,0)[0]);
     F_PROFILE(f,thread,position) =Gsource;
    }
  end_f_loop (f,thread)
}

[ 本帖最后由 gzxcqu 于 2011-5-23 10:20 编辑 ]

gzxcqu

自己顶！！！

gzxcqu

gzxcqu

置顶！！！！

gzxcqu

怎么没人发表看法呢？？？？？

gbd20022002

先不讲你的ＵＤＦ是否正确，你的所谓力平衡，是以什么物理量设定到边界上的？是在边界上设定的吗？
你所谓有自由表面，让它以流场边界形式出现？有没有考虑采用多相流模型？