编写udf遇到的问题,质量分数梯度C_YI_G(c,t,i)的应用?
浏览:1494
请教一下大家质量分数梯度的用法,C_YI_G(c,t,i)怎么使用,我编写程序的时候老是会报错 ,这里我想调用组分的质量分数梯度,用到了组分输运模型。
.\..\src\0813(1).c(18) : error C2297: “*”: 非法,右操作数包含“real [2]”类型,就是这个错误,是不是和C_YI_G(c,t,0)是矢量有关,代码:
#include "stdio.h"
#include "udf.h"
#include "math.h"
#include "mem.h"
#include "dpm.h"
#define pi 3.14
#define mu 1.8*pow(10,-5)
#define d pow(10,-6)
#define cgm -0.26
#define Dv 0.027
DEFINE_DPM_BODY_FORCE(partical_diffuse_force,p,i)
{
real beforce;
int j;
cell_t c=P_CELL(p);
Thread *t=P_CELL_THREAD(p);
beforce=-3*pi*mu*P_DIAM(p)*((18/29)/C_YI(c,t,1)+cgm)*C_YI_G(c,t,0)*Dv*29/18;
return(beforce/P_MASS(p));
}麻烦大家帮忙看看