求助一个阻力源项的问题，不好意思啊各位大侠，图片不知道怎么插入，有点问题。

passflame

The following UDF, named xmom_source, is used to add source terms in FLUENT. The source code can be interpreted or compiled. The function generates an  -momentum source term that varies with  position as



Suppose




where




Then




The source term returned is




and the derivative of the source term with respect to  (true for both positive and negative values of  ) is




/*******************************************************************/
/* UDF for specifying an x-momentum source term in a spatially     */
/* dependent porous media                                          */
/*******************************************************************/

#include "udf.h"

#define C2 100.0

DEFINE_SOURCE(xmom_source,c,t,dS,eqn)
{
  real x[ND_ND];
  real con, source;

  C_CENTROID(x,c,t);
  con = C2*0.5*C_R(c,t)*x[1];

  source = -con*fabs(C_U(c, t))*C_U(c,t);
  dS[eqn] = -2.*con*fabs(C_U(c,t));

  return source;
}
以上这个程序是从fluent中help文件中复制的，
这个程序定义的x方向的动量源项
更据上述程序中的公式可知道source的单位为：
C_R(c,t)是密度,单位是
x[1]是y轴坐标，单位m
fabs(C_U(c, t))是x方向上速度的绝对值，单位m/s
C_U(c,t) 是x方向上速度，单位m/s
那么x方向上的动量源项的单位是：
而fluent面板中规定x方向上的动量源项单位是n/m3，如下图：

那这两个地方岂不是矛盾了？，

[ 本帖最后由 passflame 于 2012-11-28 16:00 编辑 ]