CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

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

    OpenFOAM
    2
    3
    959
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • D
      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,但是我这里有点迷糊呢:mihu:

      我希望我得到的少一点,少一点,再少一点......

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

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

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

        线上CFD课程开始报名:http://www.dyfluid.com/class.html

        CFD高性能服务器 http://dyfluid.com/servers.html

        1 条回复 最后回复 回复 引用
        • 李东岳
          李东岳 管理员 最后由 编辑

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

          线上CFD课程开始报名:http://www.dyfluid.com/class.html

          CFD高性能服务器 http://dyfluid.com/servers.html

          1 条回复 最后回复 回复 引用
          • First post
            Last post