twoPhaseEulerFoam中的IATE模型系数

[李东岳]

@upc_ngh 在 twoPhaseEulerFoam中的IATE模型系数 中说：

由于给定算例的连续相速度始终为0

这个不应该啊？为什么是0

[upc_ngh]

@李东岳 这个地方我表述错了，但是在算例中即使把Crc、Cti调整到100000000这个量级对计算结果还是毫无影响，即使存在临界数的判断，Crc这个依靠相分数判断的碰撞聚集项肯定会发生变化。还有就是对R赋值，返回却是tR，这个地方对吗，我现在是这个地方不理解，谢谢东岳老师的解答。

[李东岳]

把这一行

volScalarField R = tR();

改成

volScalarField& R = tR.ref();

再试试

[李东岳]

怎么样 有进展么

[upc_ngh]

东岳老师不好意思，我现在试试，刚看到您的回复贴

[upc_ngh]

@李东岳 东岳老师我刚刚试了一下，还是不行，会报同样的错误

[李东岳]

我测试了一下，可以编译

[upc_ngh]

@李东岳 对的，可以编译，但是运行算例还是会报错，和我以前一样，编译通过，运行出错

[李东岳]

4楼那个错误么，我测试了下运行也没有错误，你是不是动别的东西了

[upc_ngh]

那老师方便把您改动好的代码发我一下吗，我对着检查再看看，谢谢老师

[李东岳]

把这一行

volScalarField R = tR();

改成

volScalarField& R = tR.ref();

就这个 没别的了

[upc_ngh]

那最后一行还是return tR吗？

[李东岳]

[upc_ngh]

这种情况，我一会再测试一下Crc和Cti是否参与了计算，谢谢老师百忙之中抽空解答我的疑惑，谢谢

[李东岳]

有进展么

[upc_ngh]

@李东岳 目前是这样的，按照东岳老师您的意见更改完以后，是可行的，Crc和Cti均可以正常参与计算了。但是计算结果和在windows下仅仅返回R的计算结果是不一样的，不知道其中原因。目前正在使用老师您的结果在跑，观察计算结果，十分感谢东岳老师的帮助
第二种是仅仅把

return tR;
改成
return R;

[李东岳]

结果怎么样

[upc_ngh]

@李东岳 已经完全正常了，多谢东岳老师一路指点，谢谢。首先能获得结果，然后再获得精准的结果

[李东岳]

好。

IATE这个目前研究的少，你可以试试。也可以用其他模型对比下。很多模型可以选，比如矩方法、section方法、n传输方法、OPOSPM等

[upc_ngh]

好的东岳老师，我再琢磨看看，感谢感谢