CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    icoFoam中压力泊松方程

    Algorithm
    2
    5
    1591
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • Z
      Zhoujq 最后由 李东岳 编辑

      在icoFoam中的压力泊松方程中
      f297e7ca-ac47-4411-b8ea-5df6e176af11-image.png

      可知(26)是通过(27)积分离散过来的

      下面进行了推导:
      3dc700d5-7329-4bbf-9514-2b1050a359a0-image.png

      我没弄明白的是
      051222da-7efa-479d-bb16-8c16f42dd79d-image.png

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

        \begin{equation}
        \frac{1}{\Delta V}\sum p_f\bfS_f=\frac{1}{\Delta V}\int p \rd\bfS=\frac{1}{\Delta V}\int\nabla p\rd V
        \end{equation}

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

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

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

          另外我发现方程12-26都需要除掉网格体积$\frac{1}{\Delta V}$,https://www.cfd-china.com/topic/2639/请教-icofoam解析-中的问题 当时做这个改动的时候弄混乱了

          原来的方程13这样

          \begin{equation}\notag
          \left( {\frac{{{V_\mathrm{P}}}}{{\Delta t}} + \sum { {\frac{{F_f^t}}{2}} + \sum { {\nu \frac{{\left| \bfS_f \right|}}{{\left| \mathbf{d} \right|}}} } } } \right)\mathbf{U}_\mathrm{P}^* = - \sum { {\left( {\frac{{F_f^t}}{2} - \nu \frac{{\left| \bfS_f \right|}}{{\left| \mathbf{d} \right|}}} \right)\mathbf{U}_\mathrm{N}^*} } + \frac{{{V_\rP}}}{{\Delta t}}\mathbf{U}_\mathrm{P}^t-\sum \frac{{p_\mathrm{P}^t + p_\mathrm{N}^t}}{2}\bfS_f.
          \end{equation}

          要改为

          \begin{equation}\notag
          \left( {\frac{{{1}}}{{\Delta t}} + \frac{1}{V_\mathrm{P}}\sum { {\frac{{F_f^t}}{2}} + \frac{1}{V_\mathrm{P}}\sum { {\nu \frac{{\left| \bfS_f \right|}}{{\left| \mathbf{d} \right|}}} } } } \right)\mathbf{U}_\mathrm{P}^* = - \sum \frac{1}{V_\mathrm{P}}{ {\left( {\frac{{F_f^t}}{2} - \nu \frac{{\left| \bfS_f \right|}}{{\left| \mathbf{d} \right|}}} \right)\mathbf{U}_\mathrm{N}^*} } + \frac{{{1}}}{{\Delta t}}\mathbf{U}_\mathrm{P}^t-\frac{1}{V_\mathrm{P}}\sum \frac{{p_\mathrm{P}^t + p_\mathrm{N}^t}}{2}\bfS_f.
          \end{equation}

          导致后面的都要做一些改动。应该是我写的一个bug

          感谢提醒!:xiexie: :xiexie: :xiexie: 已更新

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

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

          Z 1 条回复 最后回复 回复 引用
          • Z
            Zhoujq @李东岳 最后由 编辑

            @东岳 谢谢老师,还是我学艺不精~:xinxin:

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

              同理,下述网页中都需要处理$\sum_f p_f^{t}\bfS_f$为$\frac{1}{V_\rP}\sum_f p_f^{t}\bfS_f$,均已更正,谢谢

              http://dyfluid.com/rhoPimpleFoam.html
              http://dyfluid.com/rhoSimpleFoam.html
              http://dyfluid.com/simplefoam.html

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

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

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