OF有欧拉拉格朗日求解器吗

[Samuel-Tu]

铁铁们，看见东岳老师在白话多相流里面说欧拉拉格朗日法已经被广泛植入OF中，现在OF有这个求解器吗？

[Samuel-Tu]

好像搜索按钮不见了。。。

[李东岳]

才发现 我看看

[李东岳]

欧拉拉格朗日就是DPMFoam啊

论坛最近迁移服务器，才发现搜索功能丢了，我去找找，铁铁

[Samuel-Tu]

@东岳 好像还是用不了，我搜欧拉，拉格朗日等词都都不出来。

[Samuel-Tu]

最近有这么一个问题，就是给定了一个稳态的流场了，然后呢我想根据这个稳态的流场驱动颗粒的运动。这个不知道DPM能做到吗
另外我想考虑的是这个流场中的颗粒的粒径是不同的。

[李东岳]

icoUncoupleKinematicParcelFoam

[星星星星晴]

@Samuel-Tu
可以的，具体教程见下
https://www.foamacademy.com/wp-content/uploads/2018/03/particles_slides.pdf

[Samuel-Tu]

@星星星星晴 非常感谢，太有用了，还有编译过程。。这个是哪个大学的教材吗。有没有官网呀，看看还有没有其他的教程

[星星星星晴]

@Samuel-Tu google 大学
哈哈哈 搜的

[qingdong_wang]

@Samuel-Tu 这是德国罗斯托克大学的， 我也找这个做的， 不过我做的是气泡流， 加上粒子的影响。 你的颗粒是一开始静止的吗？ 我写了一个python脚本， 来注射粒子。你需要的话， 我可以贴上来，但是我写的是均匀直径的， 只要用random做一下就可以变直径了。

[李东岳]

@qingdong_wang 气泡流欧拉拉格朗日模拟？结果怎么样？我现在也在做，发现dispersion模型需要改进。

[星星星星晴]

@东岳

我们这边同事用的这个dispersion模型

J. Lipowsky, M. Sommerfeld, “LES-Simulation of the formation of particle
strands in swirling flows using an unsteady Euler-Lagrange approach,” in
6th International Conference on Multiphase Flow ICMF 2007, Leipzig, Germany,
2007.

[Samuel-Tu]

@qingdong_wang 太感谢了，麻烦贴出来学习下。我确实是想做变直径的。我是做稳态风场中雪颗粒的运动。。你是按这个教程里面开发的吗，不知道OF7能编译通过吗。。

[qingdong_wang]

@东岳 这个我不是很清楚啊， 我用laminar。关于后处理想请教您一个问题， 如何用preview将粒子显示在alpha.water 的上面。 我现在只会用extract block来显示粒子以及粒子的glyph， 但是当我把alpha.water调成visible之后，他就覆盖粒子的分布和粒子的速度矢量。 我搜了搜网上没找到怎么弄

[qingdong_wang]

@Samuel-Tu 可以通的过，我本想把我照着他做的不可压缩单相流贴在这边， 但是文件大小有限制。 那个python我也顺便贴一下。这个是配合manuelInjection model一起使用。你到时候使用的时候改一下用random函数改一下直径就好了。如果有问题， 也烦请你再联系我一下， 也让我改一下自己的错误。manualInjection.zip 在这个压缩包里，因为Python文件不能直接上传

[Samuel-Tu]

@qingdong_wang 感谢了，有任何进展我都会及时和你联系的

[李东岳]

稳态风场中雪颗粒的运动

@Samuel-Tu 为什么不直接用icoUncoupleKinematicParcelFoam，你不需要做任何编译

@qingdong_wang 你可以打开两个case，一个显示粒子，一个现实alpha

[Samuel-Tu]

@东岳 嗯嗯，我正在看icoUncoupleKinematicParcelFoam的源代码。另外我看了下射入颗粒的直径分布，有那种随机直径分布函数，normalDistribution。但是我想做的可能是PantchInjection，然后呢，入射的颗粒是跟patch的z坐标有关的。就是说这个patch，低处颗粒直径大，高处颗粒直径小。这个好像没有这个功能。
另外一个功能是想，颗粒和某一个壁面碰撞后，反射角度和入射角度是不一样的，并且根据和壁面碰撞的能量，这个壁面还会激发颗粒从这个壁面入射。不知道OF有没有这个功能。。
感觉可能要研究下Inject的代码了，可能要自己修改

[Samuel-Tu]

@qingdong_wang 非常感谢。写的代码非常清晰。我看了下tutorials的mannualInjection，但是发现PositonFile里面没有记录直径的，只记录了位置信息，不知道mannualInjection是否支持变直径这种功能。
感觉可能要深入这个源代码，看他如何初始化颗粒信息的

[Samuel-Tu]

@东岳 另外，我这个是湍流算出来的流场，我看ico这个solver初始化的时候，是用的laminar的mu,不知道是否能用。。

[qingdong_wang]

@Samuel-Tu 你可以看一下这个kinematicproperties

injectionModels
   {
       model1
       {
           type            patchInjection;
           massTotal       140e-3;
           SOI             0;
           parcelBasisType fixed;//mass;
           nParticle       1;
           patch           inlet;
           duration        1;
           parcelsPerSecond 1e5;
           U0              (0 0 0.1);
           flowRateProfile constant 1;
           sizeDistribution
           {
               type        normal;
               normalDistribution
               {
                   expectation 550e-6;
                   variance 50e-6;
                   minValue 400e-6;
                   maxValue 800e-6;
               }
           }
       }
   }

这个sizedistribution就是设置粒子直径， 这个是正态分布， 你可以看看其他设置方法， 看看有没有这种达到你想要的方法。
@东岳 这个方法我试验过， 他的图层是按照次序自然叠加的吗？ 比如说第一个case显示粒子， 第二个case显示alpha， 我就是这么做的， 但是alpha还是覆盖粒子。 是不是我应该第一个case显示alpha， 然后第二个case显示粒子， 这样就可以显示？因为我就是想要粒子显示在alpha上方。

[Samuel-Tu]

@qingdong_wang 谢谢，我也发现了这个设置。说明在初始化粒子的时候，是可以对每个粒子分别设置的，我需要去研究一下粒子的构造函数是如何初始化和调用的。

[李东岳]

@qingdong_wang 你加一些transparent就可以了 你试试

[qingdong_wang]

@东岳 我试过， 但是我觉得比较丑，哈哈。 我看网上有人用python的paraview module进行图层的叠加。 我再看看那个。

[李东岳]

这是我们放在文章里面的图，我们觉得还行 我们把粒子显示在了U上面

[Samuel-Tu]

研究了下manualInjection的构造函数，直径和速度分别有构造函数，构造了List，说明可以分别设置，有了更系统的发现再跟大家汇报