|
|
[求助]
udf加体反应—— DEFINE_VR_RATE宏hook到fluent中怎么知道是哪个物质的反应速率
还是说DEFINE_VR_RATE宏返回的速率要加到源项上面去???????
用udf写了宏之后,还用开fluent Volumetric(体反应)吗????
如果需要开体反应,那层流有限速率模型需要开吗???
开了后里面的Arrhenius动力学参数还需要设置吗???
这个问题是做fluent算例来的,例子中用udf写了反应速率,hook到fluent中,然而怎么知道这个反应速率是谁的反应速率????例子如下:
对于简单的双气相组分系统,组分a生成组分b的反应速率公式为:
![udf加体反应—— DEFINE_VR_RATE宏hook到fluent中怎么知道是哪个物质的反应速率]()
其中,Xa为组分a的质量分数,K1和K2为常数。
/* rate.c */
/* UDF for specifying a reaction rate in a porous media */
/**************************************************************/
#include "udf.h"
#define K1 2.0e-2
#define K2 5.
DEFINE_VR_RATE(user_rate, cell, thread, r, mole_weight, species_mf, rate, rr_t)
{
real s1 = species_mf[0];
real mw1 = mole_weight[0];
if (FLUID_THREAD_P(thread) && THREAD_VAR(thread).fluid.porous)
*rate = K1*s1/pow((1.+K2*s1),2.0)/mw1;
else
*rate = 0.;
} |
|