非边界区域 外加速度的udf

zhouliang125

试图在区域1、2、3分别加入速度方程 Vx=C1*X 、  Vx=C2*X 、 Vx=C3*X  


我写的好像有问题， 算出来和不外加速度方程结果一样

zhouliang125

#include "udf.h"
#define CON100.0
DEFINE_SOURCE(cell_x_source, cell,thread,dS,eqn)
{
real x[2];
real source;

if(x[1]>=0.35&&x[1]<=0.55)
{C_U(cell,thread)=0.0926*4*65.35*x[0]; source=CON*C_U(cell,thread);dS[eqn]=CON;}
else if(x[1]>0.25&&x[1]<0.35)
{C_U(cell,thread)=0.0569*4*65.35*x[0];source=CON*C_U(cell,thread);dS[eqn]=CON;}
else if(x[1]>0.55&&x[1]<0.65)
{C_U(cell,thread)=0.0145*4*65.35*x[0];source=CON*C_U(cell,thread);dS[eqn]=CON;}
return source;
}

zhouliang125

望高手指点

invader

这udf最起码也得对空间坐标赋值才行啊，
没有赋值的向量一点意义都没有，都是随机数，
#include "udf.h"
#define CON 100.0
DEFINE_SOURCE(cell_x_source, cell,thread,dS,eqn)
{
        real x[ND_ND];  此处虽然定义了，空间坐标变量，但是后面一直未赋值，
        real source;
                      *****************         此处应用某相关宏获取当前赋值区域的坐标
        if(x[1]>=0.35&&x[1]<=0.55)
        {
                C_U(cell,thread)=0.0926*4*65.35*x[0];
                source=CON*C_U(cell,thread);
                dS[eqn]=CON;
}
        else if(x[1]>0.25&&x[1]<0.35)
        {
                C_U(cell,thread)=0.0569*4*65.35*x[0];
                source=CON*C_U(cell,thread);
                dS[eqn]=CON;
}
        else if(x[1]>0.55&&x[1]<0.65)
        {
                C_U(cell,thread)=0.0145*4*65.35*x[0];
                source=CON*C_U(cell,thread);
                dS[eqn]=CON;
}
        return source;
}