GGDH模型计算浮升力生成项

#include "udf.h"
int i,j;
#define ct 0.33333 /*C_theta*/
#define Kronecher(i,j) (i==j?1.:0.) /*define Kronecher constant*/
/*define delta u,velocity flux*/
#define delta(c,t) (C_DUDX(c,t)+C_DVDY(c,t)+\
(RP_3D?C_DWDZ(c,t):rp_axi?C_V(c,t)/C_AVE_Y(c,t):0.))
/*定义S_ij=1/2*(du_i/dx_j+du_j/dx_i),二维平面、无旋流轴对称和三维定义一样*/
#define S(i,j) ((i=0&&j=0)?C_DUDX(c,t):(i=1&&j=1)?C_DVDY(c,t):(i=2&&j=2)?C_DWDZ(c,t):\
(i=0&&j=1)||(i=1&&j=0)?(C_DUDY(c,t)+C_DVDX(c,t)):(i=0&&j=2)||(i=2&&j=0)?(C_DUDZ(c,t)+C_DWDX(c,t)):\
(i=1&&j=2)||(i=2&&j=1)?(C_DVDZ(c,t)+C_DWDY(c,t)):0.)
/*定义Taot_ij*/
#define Taot(i,j) (2*C_MU_T(c,t)/C_R(c,t)*S(i,j)-2./3.*(C_K(c,t)+C_MU_T(c,t)/C_R(c,t)*delta(c,t))*Kronecher(i,j))
/*define u^'_iT^'*/
#define ut(i) ct*C_K(c,t)/C_D(c,t)*(Taot(i,0)*C_T_G(c,t)[0]+Taot(i,1)*C_T_G(c,t)[1]