曳力模型植入发散

[WeakForm]

谢谢李老师回复，我在调整之后能够跑算例了，但是跑一段时间后还是发散，在求解域上存在一些动能比较大的点，这种情况一般怎么去解决比较好呢，是直接限制Cd的大小还是有更好的办法呢

[李东岳]

一般曳力模型不应该植入不正确。你看那个Cd-Re那个图，天下曳力模型基本都没有太大差异。在这种情况下，可能你写了bug进去。

[WeakForm]

李老师我又整理了一下我的问题，首先现在openfoam目前在内的曳力模型都是针对牛顿流体的，然后一般在雷诺数较大的情况下，曳力系数都是随着雷诺数增加，比如说常用的：

等等。我目前计算这个是非牛顿的体系，然后非牛顿体系内的气泡曳力系数在一定范围内是随着雷诺数降低的，现存的关联式也基本上都是符合这个规律：
比如：
等等。然后这些实验一般针对的雷诺数范围都在几百以下，这样的话使用这些公式一方面在雷诺数较小的时候cd很大，我直接把最小雷诺数设置的比较大，一般牛顿的设置为1e-3,我设置成个位数，在雷诺数较大的入口等地方就会出现cd小的可怜的情况，就会导致计算发散。然后我看很多文献计算非牛顿系统并没有使用牛顿的曳力系数模型，或者使用一个下界来限制高雷诺数下的cd，dewsbury这个文章也提到了他认为雷诺数比较大的时候cd趋近0.95,
尝试类似设置效果并不好还是会发散。目前想问一下李老师，一是这种限制上下界是否合理，二是如果不存在高雷诺数下靠谱的公式，高雷诺数下再使用牛顿的公式进行融合是否合理呢。

[李东岳]

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

[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]

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