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,但是我这里有点迷糊呢:mihu:



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

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



  • 在这里讨论buoyantPimpleFoam,rhoPimpleFoam和buoyantPimpleFoam大同小异

    1. 首先,buoyantPimpleFoam解析中的方程(3)中的$A_{\mathrm{P}}$包含密度$\rho$;

    2. 1/Ap也是除掉密度$\rho$的

    3. 不可压缩算法和可压缩算法的压力方程不一样,可压缩算法的压力方程需要乘以$\rho$,不可压缩不需要,buoyantPimpleFoam解析中的方程(15)都需要乘以$\rho$,你看那个$\frac{\rho^{n}}{A}$,所以,

    你会发现surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU))中出现了rho*rAU


Log in to reply