(经典好帖!!!!!!!!)udf编译的经验总结

hakuzlj

#include "udf.h"
#include "math.h"

real const Cu=0.09;   
real const Sct=0.7;
real const Cm=0.5;
real const A0=353.71;  //组分A初始浓度
real const B0=2599.7685;  //组分B初始浓度
real const D0=2475.97;  //组分D初始浓度

real const k2=0.67313;  //反应2的反应常数
real Gama; //扩散系数 Γt
real r;  //微混合系数r
real Co;//在完全湍流区，Co=2,但Fox给出了起回归公式。
int   flag1;//确定是否计算标量1和标量2的梯度平方项（修正项），0:不计算；1:计算。
int   flag2;//确定标量1和标量2Cφ的值，1:利用Ying Liu和Fox回归的公式计算；0: 取2.0(适于

完全湍流，Rel>100)；
int   flag3;//确定是否计算标量3和标量4的梯度平方项（修正项），0:不计算；其它:计算。
int   flag4;//确定标量3和标量4Cφ的值，1:利用Ying Liu和Fox回归的公式计算；其它: 取2.0(适

于完全湍流，Rel>100)
real F_Co(real Re)
{
        real C0=0.0;
        if(Re>=100||Re<=0.2) C0=2.0;
        else
      {C0=C0+0.4093+0.6015*log10(Re)+0.5851*log10(Re)*log10(Re)+0.09472*log10(Re)*log10

(Re)*log10(Re);
      C0=C0-0.3903*log10(Re)*log10(Re)*log10(Re)*log10(Re);
      C0=C0+0.1461*log10(Re)*log10(Re)*log10(Re)*log10(Re)*log10(Re);
      C0=C0-0.01604*log10(Re)*log10(Re)*log10(Re)*log10(Re)*log10(Re)*log10(Re);
        }
        if(C0>2.0) C0=2.0;
        return C0;
}
编译以后，Fluent里出现如下错误

zhangtao2030

py19870531

嗯 正好能用上。 谢谢楼主~