Rheointerfoam求解不稳定

[WeakForm]

目前在用open foam9版本的rheotool求解粘弹性单气泡上升问题，
使用PPTLog模型或者FENE-CRLog模型稳定性都很差，Co数缩小到0.05依然发散，GiesekusLog模型能够收敛但是气液界面看起来很怪，不是光滑的水滴型。

仿真主要设置参考的官方case，即theta的散度离散采用高阶形式，
PPT模型和fene模型主要发散都发生在theta计算特征值或者计算应力时候发散，使用clipping办法限制相关量的范围也没有取得很好的效果
各位前辈有没有好的建议