利用reactingFoam计算正癸烷在水平圆管中的裂解反应（18组分生成物）、热物理模型采用
0_1492507996639_upload-7ce04876-1271-4e49-b10f-e26cf46d3989
计算结果中出口速度偏小、正癸烷在圆管中的停留时间过长，导致出口处的正癸烷质量分数偏小，考虑到与其密度有关；因为正癸烷裂解反应是在超临界条件下进行的，算例中采用的热物理模型不能满足计算要求，故考虑采用rhoReactingFoam求解器，热物理模型采用
0_1492508023620_upload-b90bd2f7-b8c7-4863-8558-4089770b0ad8
各个热物性参数通过与温度之间的多项式关系进行计算。例：
0_1492507617167_upload-2786160a-cf82-4436-b6a6-f40fb8d70e2f
但计算一步就出现错误
0_1492507781160_upload-d1e18a85-91ef-448b-a634-cacd5cda8fc0
目前不知道什么原因导致的，请问哪位大神能给指点一下？

利用reactingFoam计算正癸烷在水平圆管中的裂解反应（18组分生成物）、热物理模型采用0_1492507194722_upload-7d5ad516-4737-4153-b8ec-01e682a983c3 ，计算结果中出口速度偏小、正癸烷在圆管中的停留时间过长，导致出口处的正癸烷质量分数偏小，考虑到与其密度有关；因为正癸烷裂解反应是在超临界条件下进行的，算例中采用的热物理模型不能满足计算要求，故考虑采用rhoReactingFoam求解器，热物理模型采用0_1492507500590_upload-d139e8d7-727d-4acd-8d0b-71cc7303f484
各个热物性参数通过与温度之间的多项式关系进行计算。例：
0_1492507617167_upload-2786160a-cf82-4436-b6a6-f40fb8d70e2f
但计算一步就出现错误
0_1492507781160_upload-d1e18a85-91ef-448b-a634-cacd5cda8fc0
目前不知道什么原因导致的，请问哪位大神能给指点一下？

在用tecplot360画速度矢量图时，箭头分布比较密，想把它调稀疏点，调point只能改变I方向上的点数，J,K方向的无法调，请问有遇到类似的吗？

楼主将cantera移植到最新OpenFOAM的努力成功了吗？还望分享

或者说是，在使用backwardDual physDeltaT CoDeltaT这个ddtScheme的时候，如何去确定我想要的deltaT 和求解的subCycle 的数量呢？

@wwzhao 这个是我自己写的边界，没写别的造波方式。。。。

@温宝霖 恭喜。
你模拟的是什么问题？能否发个网格截图出来看看？听起来挺有意思。

@piteqiu 使用laplacianFoam，添加#include "fvOptions.H"，方程中添加源项（参考scalarTransportFoam），在fvOptions中定义codedSource。

src/fvOptions/sources/general/codedSource/CodedSource.H

你需要依据时间确定你的区域添加源项就可以了。可能需要一些网格的东西：mesh.C()返回网格位置，mesh().time()返回当前时间。

OF有标准版，什么都不需要改，如果仔细看程序的话。
Wall shear在有fv3的情况下不是特别好而已。

在https://www.cfd-online.com/Forums/openfoam-paraview/122425-error-calculators-equation.html#post447832
中有提到用Paraview计算平均温度的方法，

0_1491533074551_upload-3aaf8874-b0d3-40b7-822c-db8ddb6b5aa6
不知道这种面积平均得到的温度是不是就是平均温度。

@xpqiu 边界条件刚上来的那个图是什么意思呢？

@xpqiu :big_mouth:

发散的问题很常见，并不一定要换求解器。跨音速用密度基确实合理。提供更多信息看看是如何发三的。

可以，你试试mapFields，具体如何使用你看OpenFOAM的用户指南里面有使用方法，icoFoam下面有个算例可以专门练习使用mapFields。你先试试，不成功再联系！

9个月过去了，仔细看了一下一楼的结果，依然没看出什么区别..

可能是温度边界条件有点问题。你的E发散了。

@程迪 如果你详细看源代码，你会发现UEqn.D()是无量纲的，而UEqn.A()是强行把量纲附上去的，我觉得OF这些地方是有偷懒或者不一致的地方的。

这句话我也有所感觉。

一开始上手确实有点蒙。非常感谢！

@cfd-china 求解的方程是（19）。

请问如何在算例中监测释热值，求大神赐教