|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
udf只能作用在第一相上,无法加到第二相上,请高手指点一下
代码如下:
#include "udf.h"
#define pi 3.1415927
DEFINE_SOURCE(laser,c,t,dS,eqn)
{
real x[ND_ND];
real source;
real eta, P, R, x0, x1;
eta=0.7;
P=30*5;
R=0.19e-3*2;
x0=0;
x1=0.002;
C_CENTROID(x,c,t);
source = eta*P*3/pi/pow(R,2)*exp(-3*(pow(x[0]-x0,2)+pow(x[1]-x1,2))/pow(R,2));
dS[eqn] = 0;
return source;
}
|
|