我想用fluent自带的模型(未加udf)做经典案例Sn-Pb金属凝固宏观偏析模拟,做出的结果图虽然跟文献相似,但偏析超标。我一度想是不是动量方程源项boussinesq假设的参考值(温度,浓度)没起作用,于是自己又扁了一个boussinesq源项UDF,S=−ρg[1+β_T(T−T_ref)+β_s(C_l−C_ref)]。但运行出来的结果,别说对照文献了,跟本来的图都差远了。下面是它的编程;
#include"udf.h"
#define CON 0.0005886
DEFINE_SOURCE(cell_y_source,cell,thread,dS,eqn)
{
real source;
source=9.81*C_R(cell,thread)+C_R(cell,thread)*CON*(C_R(cell,thread)-499.);
dS[eqn]=CON*C_R(cell,thread);
return source;
}
#include"udf.h"
#define Cr -0.051993
DEFINE_SOURCE(y_source,cell,thread,dS,eqn)
{
real source;
real arc;
arc=1./(1.0656-0.0656*C_LIQF(cell,thread));
source=Cr*C_R(cell,thread)*(arc*C_YI(cell,thread,0)-5.);
dS[eqn]=Cr*C_R(cell,thread)*arc;
return source;
}
请各位大佬给予帮助!!!
![fluen凝固熔化及宏观偏析-经典案例Sn-Pb合金模拟求助]()
w.png
![fluen凝固熔化及宏观偏析-经典案例Sn-Pb合金模拟求助-1]()
文献.png |