这是我的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));
}

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

@bestucan 多谢多谢，非常感谢！！！！ 受益良多！

@bestucan 大佬你好啊，我好像发现我的反应机理是需要使用缩核模型来描述的，那这种是不是需要编写一个UDF来限制他的反应，也想请问一下去哪里学习化学反应的UDF呀，感觉找不太到教程？多谢啦

@bestucan 多谢大佬这么详细的回答啊，非常非常感谢！我在用多相流模型的时候温度没有上升也没有产物，但是我在使用组分运输模型的有些情况温度有上升，我也不敢确定反应是否发生。不知道是不是这种比较新的化学反应不能使用fluent自带的化学机理呀，是不是要自己导入他的化学机理呢？

@bestucan 多谢大佬提醒，我的计算没有问题，只是查看结果的时候发现多相流模型中反应选项卡里设置的反应并没有发生，经过我这一段时间的学习后我现在选择用组分运输模型里的颗粒表面化学反应，设置了一个混合物，但是我现在不太明白这个混合物我是需要将他当成一种材料，还是说单纯地下方通入氢气，进入催化剂床层就能自动发生反应？然后现在计算很难收敛，我准备将网格进行加密，不知道这样正不正确？

@bestucan 不好意思呀大佬，我想问的就是fluent里面的颗粒表面化学反应模型能不能实现在流化床层的同时，气体和床层颗粒进行化学反应，不好意思不好意思，我这样算了好几次都不收敛

请问一下各位大佬，我的模型是一个带有固体床层的反应器，当下方来气时，流化床层而且气体和床层颗粒发生反应，但是好像2021版本的多相流模型里面带的反应选项不能起作用了，而用组分运输混合物又好像无法实现我所预想的过程，想请教一下各位大佬有什么解决方法吗？不好意思麻烦啦，多谢各位

我想在fluent三维计算时设置两不同化学反应分布在两不同区域发生，就像这个图中里面圆柱区域和外面圆环柱区域，

不知道如何设置，请各位大佬指教，不好意思多多麻烦

@hustzzw
就像这个图中里面圆柱区域和外面圆环柱区域