CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    拉格朗日粒子在平行運算下的問題

    OpenFOAM
    3
    8
    1223
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • R
      r07522121 最后由 编辑

      大家好,先簡單說一下我的狀況:
      最近在仿真超音速下sprayFoam的問題,固定Equivalence ratio情況下, 粒子在每個timeStep下的數量都非常巨量。但平行運算下即使我使用64核計算機計算情況下運算速度覺得不理想。後來發現如果parcel injection並非均勻的話(i.e. 某位置injection), 即使Mesh 被decompose成n個domain,大量什至全部的拉格朗日粒子仍然只會被一個processor所處理,Eularian-phase早就算完,處理器都在等某processorN計算粒子。

      目前想到是手動decompose 整個computational domain(例如injection的附多切好幾個subdomain),看到decomposeParDict有MANUAL的type,但需手動一個一個cell的label放在input file,而且對於cell number非常多的問題非常不實際。 想請問各位大大有解決方法或者想法嗎?

      Regards,
      Rick

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

        http://www.cfd-china.com/topic/231/openfoam并行计算如何手动划分区域?_=1585089700462

        2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        1 条回复 最后回复 回复 引用
        • 刘
          刘雄国 最后由 编辑

          可neng你的粒子源就在很小区域位置释放,所以导致进程不平均。建议在每个时间步网格重组,每个进程都有global mesh,在global mesh进行粒子输运, 就能解决这个问题啦!However, it is a little hard.

          R 1 条回复 最后回复 回复 引用
          • 刘
            刘雄国 最后由 编辑

            如果是静态网格 global mesh 可以只重组一次,就能节省通信时间。动态网格这个方法也不是特别好,需要测试一下。OpenFOAM这个粒子输运库我想改改,太慢了

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

              OpenFOAM这个粒子输运库我想改改,太慢了

              是觉得粒子跟踪算法慢还是并行通信慢?

              2023年线下CFD课 报名正式启动 http://dyfluid.com/class.html
              CFD高性能服务器 http://dyfluid.com/servers.html

              刘 1 条回复 最后回复 回复 引用
              • R
                r07522121 @刘雄国 最后由 编辑

                @刘雄国 您好,

                是的因為粒子源在某很小區域噴發出來,所以通常只會被一個或兩個processor domain所計算,拖慢整體計算效率。

                請問一下你說的是reconstructMesh嗎? 可是我沒有使用Dynamic mesh。請問該如可implenment?

                謝謝

                刘 1 条回复 最后回复 回复 引用
                • 刘
                  刘雄国 @李东岳 最后由 编辑

                  @东岳 老师不好意思,我回答不了,和具体应用有关。

                  1 条回复 最后回复 回复 引用
                  • 刘
                    刘雄国 @r07522121 最后由 编辑

                    @r07522121 google MPI load balance

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