24小时热门版块排行榜    

查看: 417  |  回复: 0

黄涛神

捐助贵宾 (小有名气)

[求助] 能量源项加载不上,求助!!!!!!!

这是一个体热源模型,求大家看看代码有没有问题。
#include"udf.h"
#define pi 3.1415926
DEFINE_SOURCE(my_source,c,t,dS,eqn)
{
        real x[ND_ND];
        real a,b,d,p,m,source,flow_time;
        C_CENTROID(x,c,t);
        flow_time=CURRENT_TIME;
        a=3.0e-3;
        b=3.0e-3;
        d=5.0e-5;
        p=15000.0;
        m=3.0*pow(x[0]-(5.0e-3)-(1.5e-2)*flow_time,2)/a/a+3.0*pow(x[1],2)/b/b+3.0*pow(x[2],2)/d/d;
        if (x[2]<=0&&m<=3)
                source=3.0*p*exp((-1.0)*m)/(pi*a*b*d);
        else if (m>3)
                source=0.0;
    dS[eqn]=0;
        return source;
}
回复此楼

» 猜你喜欢

已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 黄涛神 的主题更新
信息提示
请填处理意见