udf燃面加质问题，等效源项法，请教

wind_cool2009

有没有谁涉及过燃面加质问题的，或者说等效源项法，即mass_flux*Area=mass_source*volume，
故质量源项mass_source=mass_flux*Area/volume
也就是说入口质量流量入口转化成靠近入口一层单元内加一个等效源项，添加质量源项和动量源项，但是计算的结果和质量流入口差异很大，里面肯定有问题，，，有做过的朋友，，解答一下
下面是质量流入口udf
DEFINE_PROFILE(mass_ap,t,i)
{
        real t_ap;
        face_t f_ap_massflux;
        real r_ap_massflux;
        begin_f_loop(f_ap_massflux,t)
        {
                t_ap=F_T(f_ap_massflux,t);
                r_ap_massflux=1.45e5*exp(-11000/t_ap);//unit cm/s
                F_PROFILE(f_ap_massflux,t,i)=1.95*r_ap_massflux*10;//unit kg/m2/s
        }
        end_f_loop(f_ap_massflux,t)

}
DEFINE_PROFILE(mass_b,t,i)
{
        real t_b;
        face_t f_b_massflux;
        real r_b_massflux;
        begin_f_loop(f_b_massflux,t)
        {
                t_b=F_T(f_b_massflux,t);
                r_b_massflux=1.036e3*exp(-7500/t_b);//unit cm/s
                F_PROFILE(f_b_massflux,t,i)=0.92*r_b_massflux*10;//unit kg/m2/s
        }
        end_f_loop(f_b_massflux,t)

}

但是如果把他改写成等效源项计算，，结果就有问题！
网格情况是均匀网格，，正四边形，，变长1e-6m，添加靠近界面一层单元的等效源项，其udf如下
DEFINE_SOURCE(mass_source, c, t, dS, i)
{
        real source_mass;
                real x[ND_ND];
        C_CENTROID(x,c,t);
        if((x[0]<=20e-6)&&(x[1]<1e-6))
        {
                 source_mass=0.92*1.036e3*exp(-7500/C_T(c,t))*1e7;//unit kg/m3/s
        }
        if((x[0]>20e-6)&&(x[1]<1e-6))
        {
                         source_mass=1.95*1.45e5*exp(-11000/C_T(c,t))*1e7;
        }
        else
        {
                        source_mass=0.0;
        }
        dS=0.0;
        return source_mass;
}

DEFINE_SOURCE(momentum_y,c,t,dS,i)
{
        real momen_y;
        real xx[ND_ND];
        real vg,con;
        C_CENTROID(xx,c,t);
        if((xx[0]<=20e-6)&&(xx[1]<1e-6))
        {
               con=0.92*1.036e3*exp(-7500/C_T(c,t))*10;//unit kg/m2/s
               vg=con/C_R(c,t);//unit m/s
               momen_y=con*1e-6/C_VOLUME(c,t)*vg;//unit kg/m3/s*m/s
        }
        if((xx[0]>20e-6)&&(xx[1]<1e-6))
        {
                   con=1.95*1.45e5*exp(-11000/C_T(c,t))*10;//unit kg/m2/s
                   vg=vg=con/C_R(c,t);//unit m/s
                   momen_y=con*1e-6/C_VOLUME(c,t)*vg;
        }
        else
                  momen_y=0.0;
        dS=0.0;
        return momen_y;
}

[ 本帖最后由 wind_cool2009 于 2012-4-26 09:28 编辑 ]

wind_cool2009

求指教

[ 本帖最后由 wind_cool2009 于 2012-4-26 09:28 编辑 ]

xiaoguo1216

您好，我也遇到你的问题了，能交流一下吗