边界面循环宏下F_T(f, t)表示温度
_t f,Thread *t)2、物理量获取语句(1)面物理量获取语句F_U(f,t) x方向速度 (face_t f,Thread *t)F_V(f,t) y方向速度 (face_t f,Thread *t)F_W(f,t) z方向速度 (face_t f,Thread *t)F_T(f,t) 获取面上的温度 (face_t f,Thread *t) F_H(f,t) 获取面上的焓值 (face_t f...
face_t f;real flux1;Thread *wind;Thread *t;wind=Lookup_Thread (domain, 10);begin_f_loop(f...
注:PRINCIPAL_FACE_P宏只能用于编译型UDF中。 begin_f_loop(f,t) ifPRINCIPAL_FACE_P(f,t) { F_AREA(area,f,t); total_area +=NV_MAG(area); } end_f_loop(f,t) 4 外部Thread数据存储 每个Thread将与其网格单元或网格面相关联的数据存储在一组数组中。例如,压力存储在一个数组中,单元格c的压力是...
begin_f_loop(f,t) { real time=RP_Get_Real("flow-time"); F_CENTROID(xc,f,t); x=xc[0]; y=xc[1]; if(0<=time<=1) F_PROFILE(f,t,i)=2*a*P/(pi*r*r)*exp(-2*(pow(x-0.004-v*time,2)+pow(y-0.004,2))/(r*r)); ...
p_sum = F_P(f, t); /* 将单元中的压力值赋予 p_sum中*/ tem= F_T(f, t);r=F_P...
if (NNULLP(THREAD_STORAGE(t,SV_DENSITY))) dens = F_R(f,t); /* Set dens to face value if available */ else dens = C_R(c0,t0); /* else, set dens to cell value */ NV_DS(psi_vec, =, F_U(f,t), F_V(f,t), F_W(f,t), *, dens); ...
F_CENTROID(x,f,t); /*一个函数,它的意思是读取每个网格面质心的二维坐标,并赋值给 x。x 为名称,接收三维坐标值。f为网格面(因为这里可是取的面的二维坐标,所以 为f,若是是网格单元的话, 这里就为c)。t为指向结构体thread(这里的thread表示界线上所有的网格面的会集)的指针*/ y = x[1]; /*把质心...
F_NNODES(f,t) 面上节点的数目(int) 三、cell宏函数 1、 C_CENTROID(x,c,t) real x[ND_ND] 单元格的重心坐标x[] 注意:这句话直接使用,结果会输出到求解器中,需要使用的时候,调用x[];其中,x[0]代表x方向坐标,x[1]代表y方向坐标,x[2]代表z方向坐标。 2、 C_VOLUME(c,t) 单元格的体积(real...
end_f_loop(f, t) thread_loop_c(t, d)///遍历混合域d上的线程 { //代码// } 值得注意的是,DEFINE宏里面会包含结构宏,这样说来,DEFINE宏可以类比为C语言中的自定义函数,而结构宏则可以类比成循环语句for或while,只不过这里的循环的是网格或线程。 要...