曳力模型植入发散

[李东岳]

你的算例，非牛顿流体，粘度场怎么样，最大值能有多少

[WeakForm]

粘度最大280mpas，用幂律定律描述的话流变指数n0.8。从我实验测试的气泡速度来看剪切变稀的粘度最低能到160多，仿真的粘度场也差不多这个范围

[李东岳]

你用SN曳力模型，针对这个非牛顿流体，能计算正确的流场结果么

[WeakForm]

SN模型和ishiizuber模型算出来大差不差，感觉误差算是可接受

[李东岳]

看起来还是曳力模型的问题。你可以把你的模型Cd跟Re的线，以及SN的Cd跟Re的线，放到一个图里么？

[WeakForm]

，我画的ishiizuber是看HZDR他们baselinemodel那个模型，表达式如下
openfoam自带的ishiizuber趋势应该是和sn模型差不多

[李东岳]

在Re比较小的时候，Cd差的很大。在算例初始的时候，速度的滑移非常小，一般Re非常小。所以你的模型可能会有稳定性问题。我建议你把你的模型，在Re比较小的时候做个clip，而不是原生的模型，来提高数值稳定性。

[WeakForm]

好的李老师我去试试

[李东岳]

我在2018年左右的时候研究过曳力对稳定性的影响，当时做了一个数值稳定性的改进，但具体怎么搞的我已经记不清了。clip是个非常好处理的方法。但是我记得Re比较小的时候存在一个问题。当时用拉格朗日粒子验证的。

这个是需要写sci的么，还是你做着玩的

[WeakForm]

是的老师，打算评估一下各种非牛顿曳力模型作为小论文一部分内容