24小时热门版块排行榜

返回列表

xueyahai

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 114
帖子: 12
在线: 6.2小时
虫号: 2751832
注册: 2013-10-25
性别: GG
专业: 传热传质学

[求助] DEFINE_GEOM这个宏怎么用？怎么在fluent中实现曲线拟合，求个例子

DEFINE_GEOM(nihe,domain,dt,position)这个宏怎么用？怎么在fluent中实现曲线拟合，求个例子
我写了段程序fluent总是出现负体积，程序如下：
DEFINE_GEOM(nihe,domain,dt,position)
{
real x1=0;
real x2=0;
real x3=0;
real x4=0;
real y1=0;
real xy=0;
real x2y=0;
real a=0;
real b=0;
real c=0;
real d=0;
real m=80;

real a0;
real a1;
real a2;
real x01,x02,x03,x11,x12,x13,x21,x22,x23;

x1+=position[0];
x2+=position[0]*position[0];
x3+=position[0]*position[0]*position[0];
x4+=position[0]*position[0]*position[0]*position[0];
y1+=position[1];
xy+=position[0]*position[1];
x2y+=position[0]*position[0]*position[1];
a=-1/m*x1*x1+x2;
b=-1/m*x1*x2+x3;
c=-b*b/a-1/m*x2*x2+x4;
d=-b/a/m*x1+1/m*x2;

x01=-d/c/a/m*(b*x1-a*x2)+x1*x1/a/m/m+1/m;
x02=b*d/a/c-x1/a/m;
x03=-d/c;

x11=-b/a/a/c/m*(b*x1-a*x2)-x1/a/m;
x12=b*b/a/a/c+1/a;
x13=-b/a/c;

x21=(b*x1/a/m-x2/m)/c;
x22=-b/a/c;
x23=1/c;

a0= x01*y1+x02*xy+x03*x2y;
a1= x11*y1+x12*xy+x13*x2y;
a2= x21*y1+x22*xy+x23*x2y;
Message(" a0=\n", a0 ) ;
Message(" a1=\n", a1 ) ;
Message(" a2=\n", a2) ;

position[1] = a2*position[0]*position[0]+a1*position[0]+a0;
}

回复此楼