rhoPimpleFoam中pEqn.H中有关rho的疑问
-
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
,但是我这里有点迷糊呢 -
在这里讨论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
-