|
|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
我用define-on-demand宏,将计算区域所有的单元和面上的负压全部置零,为什么边界面上的压力没有改变?
程序是这样的:
#include "udf.h"
DEFINE_ON_DEMAND(Static_Pressure)
{
Thread *t;
cell_t c;
face_t f;
Domain *d;
d=Get_Domain(1);
thread_loop_c(t,d)
{
begin_c_loop(c,t)
{
if (C_P(c,t)<103000)
C_P(c,t)=103000;
else
C_P(c,t)=C_P(c,t);
}
end_c_loop(c,t)
begin_f_loop(f,t)
{
if (F_P(f,t)<103000)
F_P(f,t)=103000;
else
F_P(f,t)=F_P(f,t);
}
end_f_loop(f,t)
}
}
[ 本帖最后由 l55366 于 2011-11-12 14:37 编辑 ] |
|