DPMFoam和twoPhaseEulerFoam颗粒流化

[李东岳]

DPMFoam计算二维有很大争议，你可以在这个方向上加一层网格，然后处理成symmetr边界。

可能是因为你的计算域是二维的，二维情况下颗粒是圆柱，圆柱大小和你展向的长度（一层网格）有关系，可以减少展向的长度进行计算

@zhouxu 嗯，我也一直在考虑DPMFoam计算二维的这个问题。这方面你们实测过这个影响因素？能做到跟展向无关的解么

[zhouxu]

@李东岳 没做过这方面的研究

[书生kao]

@李东岳 @星星星星晴 @zhouxu
各位大佬，问题基本解决了，需要关闭kineticCloudProperties中的cellValueSourceCorrection选项。
我在之前气固反应模型添加过程中借鉴了reactingParcelFoam和coalChemistryFoam的内容，计算中开启了cellValueSourceCorrection选项，其计算过程为

 template<class ParcelType>
 template<class TrackCloudType>
 void Foam::KinematicParcel<ParcelType>::cellValueSourceCorrection
 (
     TrackCloudType& cloud,
     trackingData& td,
     const scalar dt
 )
 {
     td.Uc() += cloud.UTrans()[this->cell()]/massCell(td);
 }

我在网上找了一个解答链接文本https://bugs.openfoam.org/view.php?id=2199
2. Normally during the particle iteration the continuous phase properties are frozen and the momentum sources from particles are applied to the next fluid solution. If cellValueSourceCorrection is true, the continuous phase properties are modified already after each particle time-step to roughly take into account the coupling from the particles to the fluid. This can be sometimes beneficial, but it can also cause numerical problems as there is no under-relaxation applied.
供各位参考。

[星星星星晴]

@书生kao 从来没有想到过这样的事情，因为我这边修改了dt的计算方式，按照理论上说我们的dt比较小，所以会把sourceterm最后计算

[vbcwl]

@星星星星晴 请问一下前辈，结果文件没有你说的这两个文件是怎么回事？

[星星星星晴]

@vbcwl 额 你说的什么文件？我不用DPMFoam的，当年这也就是提供个思路。
你要查什么field？ 有些field是不会写出来的，你得具体找你要查的什么，给他命令让他写出来，再编译，你就能查看到了。

[vbcwl]

@星星星星晴 找到原因了，是因为在移动硬盘里跑的，kingmaticCloud:UCoeff这个文件名带冒号非法

[星星星星晴]

@vbcwl 你竟然在移动硬盘里跑。。本身usb的速度限制是一方面，另外你就不怕硬盘坏了么。。牛

[vbcwl]

@星星星星晴 电脑之前师兄用的，只剩10G了，专门买来一个跑程序的硬盘，虽然读写的时候慢一点但是跑程序影响不大。

[星星星星晴]

@vbcwl 6