SST K-w湍流模型入口边界条件的UDF编写问题
本帖最后由 zf392222982 于 2013-6-28 11:25 编辑以下是我用于RNG K-E模型的速度,湍动能,耗散率的UDF
#include "udf.h"
DEFINE_PROFILE(velocity_inlet, thread, position)
{
real x;
real z;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z=x/10;
F_PROFILE(f, thread, position) =26.83*pow(z,0.15);
end_f_loop(f, thread)
}
}
DEFINE_PROFILE(k_inlet, thread, position)
{
real x;
real z;
real i;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z=x;
if(z<5)
{i=0.23;}
if(z>=5,z<=350)
{i=0.1*pow(z/350,-0.21);}
if(z>350)
{i=0.1;}
F_PROFILE(f, thread, position) =1.5*(26.83*pow(z/10,0.16)*i)*(26.83*pow(z/10,0.16)*i);
end_f_loop(f, thread)
}
}
DEFINE_PROFILE(e_inlet, thread, position)
{
real x;
real z;
real i;
real l;
real k;
face_t f;
begin_f_loop(f, thread)
{
F_CENTROID(x,f,thread);
z=x;
if(z<5)
{i=0.23;}
if(z>=5,z<=350)
{i=0.1*pow(z/350,-0.21);}
if(z>350)
{i=0.1;}
l=100*pow(z/30,0.5);
k=1.5*(26.83*pow(z/10,0.16)*i)*(26.83*pow(z/10,0.16)*i);
F_PROFILE(f, thread, position) =pow(0.09,0.75)*pow(k,1.5)/l;
end_f_loop(f, thread)
}
}
我现在要选用SST K-w湍流模型,那么我对入口边界条件的UDF编写有几个问题
1、上面的UDF是不是只有速度入口velocity_inlet是可以直接使用的,k和e在SST里是无法使用的?
2、SST k-w模型如果要编写适合我过新荷载规范的边界条件还需要其他什么参数?
3、SST模型里的湍动能k和比耗散率w是不是得重新编写?如果要重新编写,谁能帮发对上面的UDF进行修改后发我QQ392222982@qq.com一份?
页:
[1]