Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    使用MPPICFoam计算时,颗粒初始位置此文件如何建立?

    OpenFOAM
    3
    13
    1407
    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.
    • 李
      李梦瑶 last edited by

      各位前辈好,我想用MPPICFoam来模拟流场中颗粒的运动,需要设定一个特定颗粒浓度的初始流场,换而言之,将一定量的颗粒均匀分布在计算域内,我目前打算用mannualInjection的注入模型人工设定颗粒位置,但是不知道如何设置(已知颗粒以初始浓度为0.016均匀分布在高1.5米的立管(三维圆柱体)中;(计算颗粒总数为40万)),恳请各位老师指导一下,谢谢!

      hurricane007 zhe 2 Replies Last reply Reply Quote
      • hurricane007
        hurricane007 @李梦瑶 last edited by

        @李梦瑶 Goldschimidt 算例里面, constant/kinematicCloudPositions,里面你能看一行一行的括号,里面有三个数字,就是每个颗粒的坐标,把这堆东西都删掉,换成你自己需要的就可以了。当然,40万颗粒的话,找个别的语言, Julia/Python/Matlab等来写个小代码生成、输出一个颗粒的位置信息的文本比较好。

        李 zhe 2 Replies Last reply Reply Quote
        • 李
          李梦瑶 @hurricane007 last edited by

          @hurricane007 好的,谢谢您!

          1 Reply Last reply Reply Quote
          • zhe
            zhe @李梦瑶 last edited by

            @李梦瑶 你好,我看到关于MPPIC,所以想给你一点儿提示,希望能帮助你。MPPIC Foam也是可以有很多种粒子注射方式的。像你所说的手动注入,是最直观的一种,也可能算是最费力不讨好的一种。因为MPPIC并不是很好的追踪粒子位置,它是通过对粒子速度修正才模拟粒子的路径的。那么也就是说当你排列出40万粒子的时候,下一个时间步长,就都打乱了。它会根据修正后的效果重新分配粒子的位置。我建议你那么大数量级的粒子堆,用zoneInjection好一些。调试好你的初始浓度,自动分配粒子位置。希望能帮到你。

            李 2 Replies Last reply Reply Quote
            • 李
              李梦瑶 @zhe last edited by

              @zhe 好的,非常感谢您的建议,我去试一下!

              1 Reply Last reply Reply Quote
              • 李
                李梦瑶 @zhe last edited by

                @zhe 前辈,麻烦再问一下,初始浓度在哪设置?

                zhe 1 Reply Last reply Reply Quote
                • zhe
                  zhe @李梦瑶 last edited by

                  @李梦瑶 这个不是直接输入的量。只是有一个最大不能超过的值可以规定,至于你所想要的特定的浓度值,我的经验是你需要知道你每个网格在划分时的大小。然后通过在注入粒子时粒子size的确定,大致就可以知道每个网格的浓度。因为它是系统平均分配的。你可以尝试几个步长的case试试,看看注入网格后粒子浓度的值。如果有更好的方法,希望你能分享。

                  李 1 Reply Last reply Reply Quote
                  • 李
                    李梦瑶 @zhe last edited by

                    @zhe 非常感谢您!

                    hurricane007 1 Reply Last reply Reply Quote
                    • hurricane007
                      hurricane007 @李梦瑶 last edited by

                      @李梦瑶 话说你最后有没有搞定这个cellZoneInjection啊,能不能给个模板啊

                      李 1 Reply Last reply Reply Quote
                      • 李
                        李梦瑶 @hurricane007 last edited by

                        @hurricane007 很遗憾没有 还是用的手动入射模型

                        hurricane007 1 Reply Last reply Reply Quote
                        • hurricane007
                          hurricane007 @李梦瑶 last edited by

                          @李梦瑶 看来只能at大神了。 @zhe 老师能不能给个zoneInjection那部分的模板啊。。。

                          1 Reply Last reply Reply Quote
                          • zhe
                            zhe @hurricane007 last edited by 李东岳

                            @hurricane007

                                    model1              
                                    {
                                        type            cellZoneInjection; //for injection at the patch
                                        SOI             14.15;   // Start time of injection
                                        parcelBasisType fixed;// parcel basis enumeration
                                        //patch       top;     
                                        cellZone       insandZone;
                                        //duration        100;           
                                        nParticle       1;   // number of particle per parcel
                                        numberDensity    5e7;//numberDensity multiplied by cellZoneVolume gives the number of particles to be injected
                                        U0              (0 0 0);
                                        sizeDistribution
                                        {
                                            type        normal;
                                            normalDistribution
                                            {
                                                expectation 5e-4;
                                                variance 1.25e-4;
                                                minValue 2e-4;
                                                maxValue 8e-4;
                                            }
                                        }
                                    }
                            
                            

                            你可以用这个,是可以的。

                            hurricane007 1 Reply Last reply Reply Quote
                            • hurricane007
                              hurricane007 @zhe last edited by

                              @zhe 好久没来了才看到回复,感谢!!

                              1 Reply Last reply Reply Quote
                              • First post
                                Last post

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