24小时热门版块排行榜    

查看: 375  |  回复: 1

buwangchuxin

新虫 (小有名气)

[求助] symmetry与axis情况下的udf编写是否相同

之前我使用symmetry做过二维对称条件下的udf质量源项的添加,现在改为2维的axis情况下来添加,貌似有些不同。不知道这两种情况下是一样的吗????现在我把udf附加上:#include"udf.h"
#include"mem.h"
#define D 1.82e-5
#define v1 360
#define r 1.0
DEFINE_ADJUST(adsorption_adjust, domain)
{  
        cell_t c0;
        Thread* t0;
        Thread* wall_thread;
        face_t wall =6;
    real A[ND_ND] ;
    real x[ND_ND] ;
        real At;
        real J,C;
        wall_thread= Lookup_Thread(domain,wall);
        begin_f_loop (wall, wall_thread)
        {         
    c0  = F_C0(wall,wall_thread);
        t0 = F_C0_THREAD(wall,wall_thread);
    C_CENTROID(x,c0,t0);
        C=C_YI(c0,t0,0)*1000;
    J=(-r*v1/4)/(1+r*v1/4*x[1]/D)*C;
    F_AREA(A,wall, wall_thread);
    At = NV_MAG(A);
        C_UDMI(c0, t0, 0) = J*At/C_VOLUME(c0,t0);
    }
        end_f_loop (wall,wall_thread)
}
DEFINE_SOURCE(mass_source,c,t,dS,eqn)
{   
        real source ;
        source= C_UDMI(c,t,0) ;
        dS[eqn]=0 ;
        return source;
}
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yunshu123456

新虫 (初入文坛)

这段程序貌似不能用啊
2楼2017-05-27 16:26:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 buwangchuxin 的主题更新
信息提示
请填处理意见