应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » fluent中添加氧气源项，编译能通过，一计算就出错，附程序，求大牛解答
71/1返回列表
查看: 984  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

longma135

金虫 (小有名气)

[求助] fluent中添加氧气源项，编译能通过，一计算就出错，附程序，求大牛解答

UDM已经定义，因为程序较多，只贴了氧气的源项，因为其他的源项加入的时候没有错误。
DEFINE_SOURCE(source_o2_g,c,t,dS,eqn)                                   /*氧气的质量源项，定义在气相中*/
{
        Thread *mixture_t=THREAD_SUPER_THREAD(t);                           /*混合相的指针*/
        Thread *ts=THREAD_SUB_THREAD(mixture_t,1);                                                                                          
        real Ro13,Ro14,Ro15,Ro16,Ro17,Ro18,Ro19;
        int n;
        face_t f;
        Thread *tf;
        Thread *t1;
        cell_t c1;
        Ro17=0.0;
        if(C_UDSI(c,ts,4)>1170.0)
                Ro13=1170.0;
        else if(C_UDSI(c,ts,4)<730.0)
                Ro13=730.0;
        else
                Ro13=C_UDSI(c,ts,4);
        Ro14=2500.0*exp(-6420.0/Ro13);                                                                                                     
        Ro15=(Ro14+2.0)/(2.0+2.0*Ro14);                                            /*氧气的化学反应当量系数*/
        Ro16=-C_UDMI(c,mixture_t,8)*16.0-C_UDMI(c,mixture_t,9)*16.0-C_UDMI(c,mixture_t,10)*16.0-C_UDMI(c,mixture_t,7)/12.0*Ro15*32.0;                                          
        c_face_loop(c,t,n)                                                                                                           
        {
                f=C_FACE(c,t,n);
                tf=C_FACE_THREAD(c,t,n);
                if(F_FLUX(f,tf)>0.0)
                {
                        Ro17+=F_FLUX(f,tf)*C_YI(c,t,5);
                }
                else
                {
                        t1=F_C1_THREAD(f,tf);
                        c1=F_C1(f,tf);
                        Ro17+=-F_FLUX(f,tf)*C_YI(c,t,5);
                }

        }
        Ro18=Ro16+Ro17/C_VOLUME(c,t);                                                                                                
        return Ro18;
        dS[eqn]=0;
}
添加源项计算后就出现这种问题
fluent中添加氧气源项，编译能通过，一计算就出错，附程序，求大牛解答

[ Last edited by longma135 on 2014-3-13 at 16:12 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-03-13 16:10:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longma135

金虫 (小有名气)
我删除了c_face_loop循环，就可以运行了，我想问下是什么问题，f_flux不能再两相流中使用么
赞一下 回复此楼
2楼2014-03-15 18:14:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lianghtbuct

木虫 (小有名气)
楼主你好，能请教你一个问题吗？我看到你程序中有C_YI(c,t,5)，这个5应该是指第5个组分吧？能问一下5是在哪里定义的吗？怎么定义？
      我刚学UDF，很多都不懂，谢谢楼主了！
赞一下 回复此楼
努力！加油！
3楼2014-11-20 21:39:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longma135

金虫 (小有名气)
引用回帖:
3楼: Originally posted by lianghtbuct at 2014-11-20 21:39:34
楼主你好，能请教你一个问题吗？我看到你程序中有C_YI(c,t,5)，这个5应该是指第5个组分吧？能问一下5是在哪里定义的吗？怎么定义？
      我刚学UDF，很多都不懂，谢谢楼主了！

5是第六个组分，从0开始的，定义组分需要打开组分输运方程，多看下fluent的英文帮助。
赞一下 回复此楼
4楼2014-11-21 10:04:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lianghtbuct

木虫 (小有名气)
引用回帖:
4楼: Originally posted by longma135 at 2014-11-21 10:04:04
5是第六个组分，从0开始的，定义组分需要打开组分输运方程，多看下fluent的英文帮助。...

就是在组分列表里按顺序排的吗？
赞一下 回复此楼
努力！加油！
5楼2014-11-21 10:18:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longma135

金虫 (小有名气)
引用回帖:
5楼: Originally posted by lianghtbuct at 2014-11-21 10:18:29
就是在组分列表里按顺序排的吗？...

恩啊，含量最大的放在最下面，组分标号最大，其他的随便排，还是建议看看帮助，更详细一点
赞一下 回复此楼
6楼2014-11-21 10:56:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lianghtbuct

木虫 (小有名气)
引用回帖:
6楼: Originally posted by longma135 at 2014-11-21 10:56:27
恩啊，含量最大的放在最下面，组分标号最大，其他的随便排，还是建议看看帮助，更详细一点...

好的，非常感谢！
回复此楼
努力！加油！
7楼2014-11-21 11:06:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 longma135 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭