24小时热门版块排行榜    

查看: 476  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 05:54 by qL1cU46YIdK0
[基金申请] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 05:43 by qL1cU46YIdK0
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-06 05:22 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 05:21 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 04:58 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ncTVBM6AOJHF 2026-11-04 16/800 2026-11-06 04:58 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 13/650 2026-11-06 04:49 by qL1cU46YIdK0
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 04:47 by qL1cU46YIdK0
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +7 ncTVBM6AOJHF 2026-11-04 19/950 2026-11-06 04:38 by qL1cU46YIdK0
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 4/200 2026-11-06 02:37 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 01:20 by qL1cU46YIdK0
信息提示
请填处理意见