COMSOL如何模拟这样的粉体燃烧过程

送江

尝试设计一个金属粉末混合物燃烧过程，被压成圆柱体的金属粉末材料，放置于一高压反应器内，内冲入高压气体（此气体可以和金属混合物材料发生反应）。当圆柱体的一端由点火器点燃，圆柱体开始从一端自发燃烧至另一端（利用金属燃烧时的放热维持这样的燃烧过程）。

这个模拟已经做了有快半年了，但还是没有做出来，设ETA 为转化率，模拟的主要过程是，高压气体扩散进入金属粉末圆柱体的过程，金属粉末间的反应，粉末圆柱体的传热过程，用的是PDE，热传插件，做过很多关于边界条件和反应区域的划分，但还是很不理想，计算结果不耦合，算到一半就断掉了。

希望各位朋友帮忙看看，帮着参谋参谋，非常非常感谢大家

送江

PDE中的source term 这项是用来表示什么的？

kinggun

似乎可以考虑：多孔介质流体流动+多孔介质传热，耦合燃烧反应，其中孔隙率是反应的函数。

送江

您说的这个以孔隙率为反应的函数确实没想到过，我会好好考虑。另外，现在在PDE中设置eta(反应速率)为变量，是否可以再设置一个PDE用来控制孔隙率呢？反应物的量(浓度)是否可以也用PDE来控制，即 dC_metal / dt, 由于高压气体是反应物的一种，压力和反应速率如何控制？同样适用PDE吗？非常感谢您

kinggun

呵呵，这些都可以用PDE，反应的话，就是用反应动力学方程来描述，即动力学常数和浓度的函数。

kinggun

顺便提一句，如果想研究燃烧过程中的结构力学，不妨参考地球科学中的Biot固结的例子，即孔隙弹性。

送江

请问，有没有一些实例是关于这个模拟过程的，非常感谢

送江

另外，圆柱体的金属粉末材料的边界条件如何定义呢？

ansys/ls-dyna

也很关注这个圆柱体的燃烧模拟过程，LS是怎么做的？有没有相关的粒子和文献？