有一个密度方程(类似标量输运方程)
其中表示附加速度,表示密度通量.
那么
所以得到OF代码如下
fvScalarMatrix myRhoEqn
(
fvm::ddt(rho)
+ fvm::div(ustphi,rho)
+ fvm::div(phi,rho)
+ fvc::div(rhoPhi)
==
fvOptions(rho)
);
fvOptions.constrain(myRhoEqn);
myRhoEqn.solve();
fvOptions.correct(rho);
这个推导过程正确吗?