| 查看: 1498 | 回复: 0 | |||
[求助]
新手自己写了个UDF，欲替代FLUENT里的S-A湍流模型，总是出现问题，求高手指点
|
|
#include "udf.h" /*Turbulence model constants*/ #define C_B1 0.1355 #define C_B2 0.622 #define SIG_V 0.66667 #define C_V1 7.1 #define C_W2 0.3 #define C_W3 2.0 #define K 0.4187 #define C_W1 3.206 /*Modeling the turbulent viscosity*/ real X(cell_t c, Thread *t) { return C_UDSI(c,t,0)/C_MU_L(c,t)*C_R(c,t);} real f_v1(cell_t c, Thread *t) { return X(c,t)*X(c,t)*X(c,t)/(X(c,t)*X(c,t)*X(c,t)+C_V1*C_V1*C_V1);} DEFINE_TURBULENT_VISCOSITY(sa_mut,c,t) { return C_R(c,t)*f_v1(c,t)*C_UDSI(c,t,0); } /*Modeling the turbulent production*/ real f_v2(cell_t c, Thread *t) { return 1.0-X(c,t)/(1.0+X(c,t)*f_v1(c,t));} real ome(cell_t c, Thread *t) { return sqrt(SQR(C_DUDY(c,t)-C_DVDX(c,t))+SQR(C_DUDZ(c,t)-C_DWDX(c,t))+SQR(C_DVDZ(c,t)-C_DWDY(c,t))); } real S(cell_t c, Thread *t) {return ome(c,t)+2.0*MIN(0.0,(C_STRAIN_RATE_MAG(c,t)-ome(c,t)));} real S_W(cell_t c, Thread *t) { return S(c,t)+C_UDSI(c,t,0)/K/K/C_WALL_DIST(c,t)/C_WALL_DIST(c,t)*f_v2(c,t);} real G_V(cell_t c,Thread *t) {return C_B1*C_R(c,t)*S_W(c,t)*C_UDSI(c,t,0);} /*Modeling the turbulent destruction*/ real r(cell_t c, Thread *t) {return C_UDSI(c,t,0)/S_W(c,t)/K(c,t)/K(c,t)/C_WALL_DIST(c,t)/C_WALL_DIST(c,t);} real g(cell_t c, Thread *t) {return r(c,t)+C_W2*(pow(r(c,t),6)-r(c,t));} real f_w(cell_t c, Thread *t) {return g(c,t)*pow((1.0+pow(C_W3,6)/(pow(r(c,t),6)+pow(C_W3,6))),0.166667);} real Y_V(cell_t c, Thread *t) {return C_W1*C_R(c,t)*f_w(c,t)*SQR(C_UDSI(c,t,0)/C_WALL_DIST(c,t));} /*Source term*/ DEFINE_SOURCE(sa_source,c,t,dS,eqn) { dS[eqn]=0; return G_V(c,t)-Y_V(c,t)+C_B2*C_R(c,t)/SIG_V*(SQR(C_UDSI_G(c,t,0)[0])+SQR(C_UDSI_G(c,t,0)[1])+SQR(C_UDSI_G(c,t,0)[2])); } /*Diffusivity*/ DEFINE_DIFFUSIVITY(sa_diffusivity,c,t,eqn) { return (C_MU_L(c,t)+C_R(c,t)*C_UDSI(c,t，0))/SIG_V; } 基本是按照斯坦福的一个ppt里写KE模型的udf改写的，参考FLUENT帮助文件里SA模型的方程。但是编译时有很多问题。另外，dS[eqn]没法确切表达，师兄说可以设成0，不知道是否合理。不胜感激。 |
回复此楼» 收录本帖的淘帖专辑推荐
 Fluent技巧 |
» 猜你喜欢
中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培
已经有4人回复
-
调剂推荐
已经有4人回复
-
总分322求生物学/生化与分子/生物信息学相关调剂
已经有5人回复
-
0703化学338求调剂！
已经有3人回复
-
284求调剂
已经有7人回复
-
271求调剂
已经有6人回复
-
312求调剂
已经有9人回复
-
352求调剂
已经有4人回复
-
085602化学工程求调剂。
已经有3人回复
-
085600材料与化工306
已经有6人回复
点击这里搜索更多相关资源
6