【求助】大家帮忙修改一下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

紧急求之，

xrs333

不用发明轮子了。用VOF模型就解决了。

gzxcqu

愿闻其详，麻烦版主再说明白点。

gzxcqu

这个不涉及多项流动问题啊？？？？

飞翔鸟

对啊，这么简单个模型直接VOF就行了，还要什么UDF啊？就算你不想用VOF模型，假设自由液面是个平面直接给滑移墙面就好了。

gzxcqu

使用vof怎样设定边界上的力呢？边界上是个温度梯度和浓度梯度引起的毛细力，这个怎样加进去呢？？？

飞翔鸟

哪里的边界？自由液面？壁面的接触角？表面张力和接触角啊。你可以自定义表面张力和接触角。

gzxcqu

我那个力的边界条件是在自由液面上的，自由液面作为边界，呵呵，并且表面张力不是定值，而是随温度和浓度梯度变化的。

gzxcqu

您方便加我qq吗 1049283760，咱们qq上再讨论。

gzxcqu

置顶

gzxcqu

置顶

飞翔鸟

你那个力不是外力，是内力，不需要你单独加上去。严格的讲，自由液面不是你的边界。正因为你的张力不是定值，跟温度和浓度有关系，所以用UDF定义张力即可。
我没有QQ啊，很老土。

gzxcqu

老兄，我这个还没弄出来，麻烦您能不能帮个忙，写个udf程序，小弟不胜感激，谢谢。