rhoPimpleFoam中pEqn.H中有关rho的疑问

Dingcy

OpenFOAM中的动量方程公式两边都同时除了密度rho，最近想实现把密度rho作为一个变量加在每一项中，虽然我研究的是不可压缩问题，但是参考了可压缩的rhoPimpleFoam求解器，在其pEqn.H文件中发现其动量预测的过程中代码是这样写的

volScalarField rAU(1.0/UEqn().A());
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU));

也就是rhorAUf=rho/Ap呗，但是在方程离散的过程中不应该是rho*Ap吗？rhorAUf不应该等于1/rho*Ap吗？
这里我看了东岳流体的解析，是先不带rho半离散化得到一个与icoFoam不可压缩的动量预测表达式相同的式子，然后又把速度表达式带入到可压缩带有rho的动量方程表达式中，得到了rho/Ap,但是我这里有点迷糊呢

东岳

你说的这个有点意思，我之前没详细写因为写起来比较麻烦，这是不可压缩向可压缩转变的一个流程。我有空更新一下。

你参考这个 http://dyfluid.com/buoyantPimpleFoam.html rhoPimpleFoam那个我需要大修一下，风格没更新

东岳

在这里讨论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