|
|
DEFINE_PROFILE(viscous_resistance, t, i) /*粘滞阻力*/
{
real x[ND_ND];
real y,dp;
cell_t c;
Thread *t;
C_CENTROID(x,c,t);
y=x[1];
if(y>0.06)
dp=0.06;
else
dp=0.03;
begin_c_loop(c, t)
{
F_PROFILE(c, t,i) =(150.*pow((1-epsi),2))/pow(epsi,3)*pow(dp,2);
}
end_c_loop(c, t)
}
这是我写的代码，没有错误，但是初始化会出现问题，我猜测是cell_t c调用了两次，求大神帮忙看看 |
|
回复此楼
