CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    向各位大佬请教一下气固反应的UDF

    Fluent
    2
    3
    139
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • H
      hustzzw 最后由 编辑

      请问一下各位大佬,我像如图定义了一个多相流的化学反应,并且编写了一个UDF来定义他的反应速率,可是遇到问题,计算过程中化学反应不能发生并且总是发散,想请问一下各位是否我的设置有错误还是UDF编写有错误?感谢各位大佬的指教。56507974-7eea-413f-aef2-b5cd60fd3526-image.png

      1 条回复 最后回复 回复 引用
      • H
        hustzzw 最后由 编辑

        这是我的UDF
        /Heterogeneous net reaction/

        include "udf.h"

        //# include "mem.h"
        //# include "sg_mphase.h"

        define FLUID_ZONE_ID 3 /Zone ID of Reaction/

        define R 8.31434 /Gas Constant J/mol-K/

        define Pre 2.65 /Pre-exponential factor, 1/s-kPa/

        define E1 538.2e5 /Activation Energy, J/mol/

        define rho_ca2fe2o5 14.375 /Molar density of Ca2Fe2O5, kgmol/m3/

        DEFINE_HET_RXN_RATE(het_rxn_rate, c, t, r, mw, yi, rr, rr_t)
        {
        //int zone_id;
        Thread **pt = THREAD_SUB_THREADS(t);
        Thread *prim_t = pt[0]; /Thread for primary Phase/
        Thread *sec_t = pt[1]; /Thread for secondary Phase/
        real T_SEC = C_T(c,sec_t); /Phase secondary temperature, K/
        real P = C_P(c,prim_t)/1000; /Static pressure of the primary phase, kPa/
        real y_h2 = C_YI(c,prim_t,0); /Mass fraction of gas species in primary phase/
        real Nsum;
        y_h2 *= 1/mw[0][0];
        Nsum = y_h2 ;
        y_h2 *= 1/Nsum;

        rr = rho_ca2fe2o5Py_h2Preexp(-E1/(RT_SEC));
        }

        G 1 条回复 最后回复 回复 引用
        • G
          gemini @hustzzw 最后由 编辑

          @hustzzw 你这是*号都被格式化吃掉了吧。。看起来奇奇怪怪的

          可以加些输出语句,看看y_h2这些量算的对不对。(小心除零)

          1 条回复 最后回复 回复 引用
          • First post
            Last post