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

[李梦瑶]

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

[hurricane007]

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

[李梦瑶]

@hurricane007 好的，谢谢您!

[zhe]

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

[李梦瑶]

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

[李梦瑶]

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

[zhe]

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

[李梦瑶]

@zhe 非常感谢您！

[hurricane007]

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

[李梦瑶]

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

[hurricane007]

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

[zhe]

@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]

@zhe 好久没来了才看到回复，感谢！！