Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    icoFoam中压力泊松方程

    Algorithm
    2
    5
    1428
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • Z
      Zhoujq last edited by 李东岳

      在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 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by 李东岳

        \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://dyfluid.com/servers.html

        1 Reply Last reply Reply Quote
        • 李东岳
          李东岳 管理员 last edited by 李东岳

          另外我发现方程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://dyfluid.com/servers.html

          Z 1 Reply Last reply Reply Quote
          • Z
            Zhoujq @李东岳 last edited by

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

            1 Reply Last reply Reply Quote
            • 李东岳
              李东岳 管理员 last edited by

              同理,下述网页中都需要处理$\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://dyfluid.com/servers.html

              1 Reply Last reply Reply Quote
              • First post
                Last post

              CFD中文网 | 东岳流体 | 京ICP备15017992号-2
              论坛登录问题反馈可联系 li.dy@dyfluid.com