|
|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
/* UDF for computing the ic force on a charged particle */
#include "udf.h"
#define Q 1.0 /* particle electric charge */
#define BZ 3.0 /* z component of ic field */
#define TSTART 18.0 /* field applied at t = tstart */
/* Calculate ic force on charged particle. ic */
/* force is particle charge times cross product of particle */
/* velocity with ic field: Fx= q*bz*Vy, Fy= -q*bz*Vx */
DEFINE_DPM_BODY_FORCE(particle_body_force,p,i)
{
real bforce;
if(P_TIME(p)>=TSTART)
{
if(i==0) bforce=Q*BZ*P_VEL(p)[1];
else if(i==1) bforce=-Q*BZ*P_VEL(p)[0];
}
else
bforce=0.0;
/* an acceleration should be returned */
return (bforce/P_MASS(p));
}
麻烦各位大侠看看。这是从fluent中的help文件中的, if(i==0) 这个命令是什么意思?i指的是什么? |
|