CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    MPPICFoam中如何输出某个颗粒在各个时间步的位置信息?

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

      欧拉拉格朗日方法的MPPICFoam求解器计算的颗粒输运,该怎么输出或处理粒子的位置信息,比如:给每个粒子编号,追踪某个固定编号(id=100)的粒子在每个时间步的空间位置、速度?纠结纠结······QQ截图20200926134725.png 跪求懂得前辈指点

      F 1 条回复 最后回复 回复 引用
      • F
        fubianhanshu @fubianhanshu 最后由 编辑

        @fubianhanshu 自己先回复下吧,在paraview中可以处理,我先贴出来一个图,某个固定ID粒子的Z坐标随时间的变化曲线图,具体的还在摸索,不是很清楚,希望得到大家更好的建议1.png

        F 1 条回复 最后回复 回复 引用
        • F
          fubianhanshu @fubianhanshu 最后由 编辑

          @fubianhanshu 上面的图没输出好,贴个清晰点的1.png

          U 1 条回复 最后回复 回复 引用
          • U
            upc_ngh 最后由 编辑

            您好,请问上面那个固定ID粒子的纵坐标图怎么做出来的?可以说的在详细一点吗?

            F 1 条回复 最后回复 回复 引用
            • U
              upc_ngh @fubianhanshu 最后由 编辑

              @fubianhanshu 前辈您好请问这个粒子高度变化图可以给指点一下吗?跪求:chouchou: QQ:1514969704希望得到前辈的指点

              1 条回复 最后回复 回复 引用
              • 星
                星星星星晴 最后由 编辑

                hi 关于粒子编号什么的 在of中会有origid 和origProcID两个东西来显示该粒子的编号,但是问题是,如果用并行计算,这个东西是怎么走的,目前我是不太了解。 如果你有特定的parcel需要追踪,在kinematicparcel 的move中 单独写一个输出,选择好特定的parcel,然后输出就好了,但是并行这个真的很头大。。希望给你一点思路

                 /*   std::ofstream outfile;
                    outfile.setf (ios_base::fixed, ios_base::floatfield);
                    outfile.precision (7);
                    outfile.open ("Injectior_Magnitude_kinematic_parcel", ios_base::app);
                            << p.position().component (0)           << tab
                            << p.position().component (1)           << tab
                            << p.position().component (2)           << tab
                            << mag(p.U())                  << tab
                            << p.U().component (0)                          << tab
                            << p.U().component (1)                               << tab
                            << p.U().component (2)                              << tab
                            << p.d()*1e6                   << tab            
                
                            << nl;*/
                

                在写博士毕业论文,寻找postdoc或者职位中~
                m.sui20@foxmail.com

                F 1 条回复 最后回复 回复 引用
                • U
                  upc_ngh 最后由 编辑

                  好的得到数据之后再自己作图就好了,差不多明白了,谢谢前辈!

                  1 条回复 最后回复 回复 引用
                  • bestucan
                    bestucan 版主 副教授 最后由 编辑

                    paraview 工具栏 edit 里有个 find data,可以用”与或“逻辑添加条件过滤目标信息。

                    滚来滚去……~(~o ̄▽ ̄)~o 滚来滚去都不能让大家看出来我不是老师么 O_o

                    异步沟通方式(《posting style》from wiki)(下载后打开):
                    https://www.jianguoyun.com/p/Dc52X2sQsLv2BRiqnKYD
                    提问的智慧(github在gitee的镜像):
                    https://gitee.com/bestucan/How-To-Ask-Questions-The-Smart-Way/blob/master/README-zh_CN.md

                    1 条回复 最后回复 回复 引用
                    • F
                      fubianhanshu @upc_ngh 最后由 编辑

                      @upc_ngh 6bbe091e-840b-4d00-9664-21490f4862aa-image.png 选完后显示全域,装上Python后,用find找也很方便,就像楼下道友说的

                      U 1 条回复 最后回复 回复 引用
                      • F
                        fubianhanshu @星星星星晴 最后由 编辑

                        @星星星星晴 我都是合成后处理的,选定某个粒子,整个时间序列ID号应该是不变的,多核分割后怎么处理界面的也不太清楚,我注意到分割后并行计算时,每个块内的粒子是动态变化的

                        星 1 条回复 最后回复 回复 引用
                        • 星
                          星星星星晴 @fubianhanshu 最后由 编辑

                          @fubianhanshu so, 这就是个问题,因为不太懂具体是怎么追踪的,我的想法是在injection model中,直接给parcel 添加一个maker, 然后在paraview中 单独输出这一个粒子就能得到你想要的信息了,就是费劲点,或者就是在move中添加一个输出这样做。

                          在写博士毕业论文,寻找postdoc或者职位中~
                          m.sui20@foxmail.com

                          1 条回复 最后回复 回复 引用
                          • U
                            upc_ngh @fubianhanshu 最后由 编辑

                            @fubianhanshu 好的我自己试一下,谢谢各位前辈 @星星星星晴 @fubianhanshu

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