非稳态计算向压力方程中添加源项
-
@Tong 在 非稳态计算向压力方程中添加源项 中说:
非稳态计算过程中需要连续的向模型中注入流体,
具体是什么过程?
由于尺度的问题无法使用“显式”的入口,而是使用一个“隐式”的方法:添加一个体积源项到压力方程中实现注入;
虽然看的不是很懂但是感觉很厉害
但是,你的代码看起来并不“老练”,你压力方程的源项,出现在动量方程中,倒是可以理解。但为何不直接加在动量方程中?即使加在压力方程中,直接一减
fvc::div(phiHbyA))-(wellSource/mesh.V())
看起来很粗暴的样子,最好换成通量话说回来,你这个应该是方程植入有问题。并且,你把源项放在压力方程中得到phi之后,速度需要进行减去或者加上
不是大问题,好好检查下你的方程和代码的一致性