请教高手如何编写炸药自热反应源项的UDF，急急急！

leory

如题，请教高手，如何进行炸药自热反应源项的UDF的编写，即Arrhenius反应律
S=pQZexp（-E/RT）
p为炸药密度，Q为反应热，Z为指前因子，E为活化能，R为普适气体常数，T为温度
多谢高手相助啊，辛苦了

gearboy78

原帖由 leory 于 2010-2-4 07:08 发表

                               
登录/注册后可看大图

如题，请教高手，如何进行炸药自热反应源项的UDF的编写，即Arrhenius反应律
S=pQZexp（-E/RT）
p为炸药密度，Q为反应热，Z为指前因子，E为活化能，R为普适气体常数，T为温度
多谢高手相助啊，辛苦了

fluent里面的层流反应模型，反应速率A*pow(T,b)*exp(-E/RT)，b为温度指数，A为指前因子，E为活化能，R为通用气体常数。对比一下可知，只要采用层流反应模型，然后温度指数取零，指前因子取PQZ就可以了。

030301029

原帖由 gearboy78 于 2010-2-4 07:56 发表

                               
登录/注册后可看大图

fluent里面的层流反应模型，反应速率A*pow(T,b)*exp(-E/RT)，b为温度指数，A为指前因子，E为活化能，R为通用气体常数。对比一下可知，只要采用层流反应模型，然后温度指数取零，指前因子取PQZ就可以了。

道理是这样，但是具体操作却不行。
温度是变量，与计算的所在区域，反应程度等有关。
因此，一定需要采用udf来实现。
具体可以参见udf的帮助文档例题，稍微做一些修改。