|
|
发表于 2011-9-2 16:06:22
|
显示全部楼层
dS[eqn就是求导,例题如下:
#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]; con=50×ρ×y
source = -con*fabs(C_U(c, t))*C_U(c, t); source=-con×|Vx|×Vx
dS[eqn] = -2.*con*fabs(C_U(c,t)); dS/dVx=-2×con×|Vx|
return source; 返回source的值
}
可见,
dS[eqn] 表示:dS/dVx=-2×con×|Vx|,还有不明白的吗? |
|