OF有欧拉拉格朗日求解器吗
-
铁铁们，看见东岳老师在白话多相流里面说欧拉拉格朗日法已经被广泛植入OF中，现在OF有这个求解器吗？
-
好像搜索按钮不见了。。。
-
才发现 我看看
-
欧拉拉格朗日就是DPMFoam啊
论坛最近迁移服务器，才发现搜索功能丢了，我去找找，铁铁
-
@东岳 好像还是用不了，我搜欧拉，拉格朗日等词都都不出来。
-
最近有这么一个问题，就是给定了一个稳态的流场了，然后呢我想根据这个稳态的流场驱动颗粒的运动。这个不知道DPM能做到吗
另外我想考虑的是这个流场中的颗粒的粒径是不同的。 -
icoUncoupleKinematicParcelFoam
-
-
@星星星星晴 非常感谢，太有用了，还有编译过程。。这个是哪个大学的教材吗。有没有官网呀，看看还有没有其他的教程
-
@Samuel-Tu google 大学
哈哈哈 搜的 -
@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. -
@qingdong_wang 太感谢了，麻烦贴出来学习下。我确实是想做变直径的。我是做稳态风场中雪颗粒的运动。。你是按这个教程里面开发的吗，不知道OF7能编译通过吗。。
-
@东岳 这个我不是很清楚啊， 我用laminar。关于后处理想请教您一个问题， 如何用preview将粒子显示在alpha.water 的上面。 我现在只会用extract block来显示粒子以及粒子的glyph， 但是当我把alpha.water调成visible之后，他就覆盖粒子的分布和粒子的速度矢量。 我搜了搜网上没找到怎么弄
-
@Samuel-Tu 可以通的过，我本想把我照着他做的不可压缩单相流贴在这边， 但是文件大小有限制。 那个python我也顺便贴一下。这个是配合manuelInjection model一起使用。你到时候使用的时候改一下用random函数改一下直径就好了。如果有问题， 也烦请你再联系我一下， 也让我改一下自己的错误。manualInjection.zip 在这个压缩包里，因为Python文件不能直接上传
-
@qingdong_wang 感谢了，有任何进展我都会及时和你联系的
-
稳态风场中雪颗粒的运动
@Samuel-Tu 为什么不直接用icoUncoupleKinematicParcelFoam，你不需要做任何编译
@qingdong_wang 你可以打开两个case，一个显示粒子，一个现实alpha
-
@东岳 嗯嗯，我正在看icoUncoupleKinematicParcelFoam的源代码。另外我看了下射入颗粒的直径分布，有那种随机直径分布函数，normalDistribution。但是我想做的可能是PantchInjection，然后呢，入射的颗粒是跟patch的z坐标有关的。就是说这个patch，低处颗粒直径大，高处颗粒直径小。这个好像没有这个功能。
另外一个功能是想，颗粒和某一个壁面碰撞后，反射角度和入射角度是不一样的，并且根据和壁面碰撞的能量，这个壁面还会激发颗粒从这个壁面入射。不知道OF有没有这个功能。。
感觉可能要研究下Inject的代码了，可能要自己修改 -
@qingdong_wang 非常感谢。写的代码非常清晰。我看了下tutorials的mannualInjection，但是发现PositonFile里面没有记录直径的，只记录了位置信息，不知道mannualInjection是否支持变直径这种功能。
感觉可能要深入这个源代码，看他如何初始化颗粒信息的 -
@东岳 另外，我这个是湍流算出来的流场，我看ico这个solver初始化的时候，是用的laminar的mu,不知道是否能用。。
-
@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上方。 -
@qingdong_wang 谢谢，我也发现了这个设置。说明在初始化粒子的时候，是可以对每个粒子分别设置的，我需要去研究一下粒子的构造函数是如何初始化和调用的。
-
@qingdong_wang 你加一些transparent就可以了 你试试
-
@东岳 我试过， 但是我觉得比较丑，哈哈。 我看网上有人用python的paraview module进行图层的叠加。 我再看看那个。
-
这是我们放在文章里面的图，我们觉得还行
我们把粒子显示在了U上面 -
研究了下manualInjection的构造函数，直径和速度分别有构造函数，构造了List，说明可以分别设置，有了更系统的发现再跟大家汇报
