或者说是，在使用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）。

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

在你的pipeline brower里面打开两个.foam算例，针对每个生成plot over line滤镜，然后同时显示两个plot over line。

@xpqiu 谢谢，就是时间步问题

@xpqiu ，谢谢解答

simpleForm与icoFoam这些求解器在运行时如何通过MPE监测API的调用频次与时间呢？

相似问题：http://www.cfd-china.com/topic/734/采用什么工具可以评估mpi调用过程中每个api调用的频次还有api执行所花费的时间呢