两相流动量方程求解的一点疑惑

[ZY-GONG]

请教各位大佬一个问题，两相流，在动量方程这块想处理相间作用力，OF我目前了解到的求解器是采用显式处理，

，但是我想做成隐式处理，弄成两个类似于a1Ul+b1Ug=c1,a2Ul+b2Ug=c2的矩阵,消元后再解线性方程组，这种有实现的可能吗，还是说这种应该做成迭代求解的模式呢？非常感谢!

[ZY-GONG]

这个是我想求解的方程形式，这段时间看了点相关东西，请问各位前辈，foam-extend里的blockMatrix能实现我想要的求解方式吗，然后就是曳力项应该OF里应该怎么写呢？是fvm : : Sp ( K, Ug ) - fvm : : Sp ( K, Ul )吗？希望各位大佬不吝赐教啊