Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

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

    OpenFOAM
    4
    12
    1504
    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.
    • F
      fubianhanshu last edited by

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

      F 1 Reply Last reply Reply Quote
      • F
        fubianhanshu @fubianhanshu last edited by

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

        F 1 Reply Last reply Reply Quote
        • F
          fubianhanshu @fubianhanshu last edited by

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

          U 1 Reply Last reply Reply Quote
          • U
            upc_ngh last edited by

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

            F 1 Reply Last reply Reply Quote
            • U
              upc_ngh @fubianhanshu last edited by

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

              1 Reply Last reply Reply Quote
              • 星
                星星星星晴 讲师 last edited by

                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;*/
                

                m.sui20@foxmail.com

                F 1 Reply Last reply Reply Quote
                • U
                  upc_ngh last edited by

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

                  1 Reply Last reply Reply Quote
                  • bestucan
                    bestucan 版主 副教授 last edited by

                    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 Reply Last reply Reply Quote
                    • F
                      fubianhanshu @upc_ngh last edited by

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

                      U 1 Reply Last reply Reply Quote
                      • F
                        fubianhanshu @星星星星晴 last edited by

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

                        星 1 Reply Last reply Reply Quote
                        • 星
                          星星星星晴 讲师 @fubianhanshu last edited by

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

                          m.sui20@foxmail.com

                          1 Reply Last reply Reply Quote
                          • U
                            upc_ngh @fubianhanshu last edited by

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

                            1 Reply Last reply Reply Quote
                            • First post
                              Last post

                            CFD中文网 | 东岳流体 | 京ICP备15017992号-2