请教个X方向源项的编程问题

yt19870303

X方向的源项表达式为 Sx=-U * Vx / B 其中U为流体的动力粘度，Vx为X方向的速度，B为材料的渗透率。B=1e-12。请教一下，这个程序怎么编，主要是流体的动力粘度，不知该怎么办。

注：
    模型为层流，气体为理想气体

以下为本人自己编的，但是感觉动力粘度的处理上不怎么正确，请高手赐教！
#include "udf.h"

#define B  1e-12

DEFINE_SOURCE(cell_x_source, cell, thread, dS, eqn)
{
  real source;

  real u;

  u=C_MU_EFF(cell,thread);

      /* source term */
      source = -u*C_U(cell,thread)/B;

      /* derivative of source term w.r.t. x-velocity. */
     ds=-u/B

  return source;
}

[ 本帖最后由 yt19870303 于 2009-4-27 14:13 编辑 ]

peaker007

我不是高手， 。
你的程序中未包含源项的偏导数项的语句。

yt19870303

忘写了，源项的偏导数项已经添加上去了