在这里讨论buoyantPimpleFoam,rhoPimpleFoam和buoyantPimpleFoam大同小异
首先,buoyantPimpleFoam解析中的方程(3)中的$A_{\mathrm{P}}$包含密度$\rho$;
1/Ap也是除掉密度$\rho$的
不可压缩算法和可压缩算法的压力方程不一样,可压缩算法的压力方程需要乘以$\rho$,不可压缩不需要,buoyantPimpleFoam解析中的方程(15)都需要乘以$\rho$,你看那个$\frac{\rho^{n}}{A}$,所以,
你会发现surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU))中出现了rho*rAU