请教：关于openfoam中添加源项的问题

[yingqing]

volVectorField Source_U("Source_U", -U*Source_Y);
volVectorField Source_U_linear_explicit("Source_U_linear_explicit", SourceMomentumFactor*rhoU);
solve
(
	fvm::ddt(rhoU) 
	+ fvc::div(phiUp) 
	== 
	Source_U 
	- Source_U_linear_explicit 
	+ fvm::SuSp(SourceMomentumFactor, rhoU)
);

在看别人程序的时候发现在源项的添加上多出来两项Source_U_linear_explicit 和fvm::SuSp(SourceMomentumFactor, rhoU)，这两感觉是一样的，不是很理解这样写的原因。哪位老师，能指导一下

[yingqing]

在fluent中也存在这种源项线性化，目的是增强求解稳定性。
源项通常表示为S_phi=A+B*phi的形式，其中phi表示因变量，A为源项的显式部分，B*phi为隐式部分

A=S_phi-d(S)/d(phi)*phi
B=d(S)/d(phi)

参考：《ansys fluent 二次开发指南》
rhoU做为因变量，Source_U_linear_explicit 和fvm::SuSp(SourceMomentumFactor, rhoU)大概也是为了形成A和B。

[李东岳]

2楼回帖没看懂。不过你1楼的回帖，是把方程处理成$...=(-\phi)+(\phi)$的形式，$(-\phi)$隐性离散增加对角占优