UDF求助求助急！

wangkeke

我模拟的是多个化学反应的燃烧模型，怎么用UDF定义其中一个反应的反应速率呢？我好想看到的只是统一的反应速率，没有分别单个定义的地方，还是是我没找到？请知道的告知我一下，非常感谢！

gearboy78

原帖由 wangkeke 于 2010-5-17 02:43 发表

                               
登录/注册后可看大图

我模拟的是多个化学反应的燃烧模型，怎么用UDF定义其中一个反应的反应速率呢？我好想看到的只是统一的反应速率，没有分别单个定义的地方，还是是我没找到？请知道的告知我一下，非常感谢！

DEFINE_VR_RATE的例子fluent udf帮助里面有的。

wangkeke

谢谢你回复我，我想请问可以分别定义每一个反应的速率吗？比如我有20个反应，我想对第5个反应的反应速率用UDF,在哪里可以设置呢？

zyc_8501

我看过您的贴子，用到了UDM，想向您请教一下~我想用UDM来显示uds的梯度，但是最后后处理的时候，显示的值都是零，望您指点一下是哪出问题了：
# include "udf.h"
DEFINE_ON_DEMAND(store_gradient)
{
Domain *d;
cell_t c;
Thread *t;
real gar1;
real gar2;
real gar3;
d = Get_Domain(1);
thread_loop_c(t,d)
{
begin_c_loop (c,t)
{
gar1=C_UDSI_G(c,t,0)[0];
gar2=C_UDSI_G(c,t,0)[1];
gar3=C_UDSI_G(c,t,0)[2];
C_UDMI(c,t,0) =gar1;
C_UDMI(c,t,1) =gar2;
C_UDMI(c,t,2) =gar3;
}
end_c_loop (c,t)
}
}
万分感谢~

wangkeke

你用什么后处理的？

wangkeke

谁能帮我解释一下这个例子啊，看不懂，呜呜
DEFINE_VR_RATE(vol_reac_rate, c, t, r, wk, yk, rate, rr_t)
{
real ci, prod;
int i;
/* Calculate Arrhenius reaction rate */
prod = 1.;
for(i = 0; i < r->n_reactants; i++)
{
ci = C_R(c,t) * yk[r->reactant] / wk[r->reactant];
prod *= pow(ci, r->exp_reactant);
}
*rate = r->A * exp( - r->E / (UNIVERSAL_GAS_CONSTANT * C_T(c,t))) *
pow(C_T(c,t), r->b) * prod;
*rr_rate = *rate;
/* No "return..;" value. */
}

zyc_8501

用得是fluent里面的display contour，udm显示的结果都是零啊~

wangkeke

那说明你的值没传进去啊，是不是梯度返回有问题