2D欧拉拉格朗日计算
-
个人经验来看,我对2D欧拉拉格朗日模拟理解存在一些疑惑,考虑相间动量耦合,在欧拉框架下,可以写为:
方程中不存在网格体积,因此个人认为2D欧拉欧拉方法是适用的。但是对于2D欧拉拉格朗日,曳力导致的动量交换计算大体如:
其中 为每个粒子的曳力, 为网格体积。对于2D网格,如果还是按照通常的计算标准,网格厚度应该不对结果有影响,但是很明显,网格越厚,力越小,导致2D欧拉拉格朗日的模拟结果依赖于网格厚度,看起来不是很合理。在Kadyrov et al.
Impacts of solid stress model on MP-PIC simulation of a CFB riser with EMMS drag
的这篇文章中,用的就是2D欧拉拉格朗日模拟,同时,文中还专门讨论了网格厚度对结果的影响,结论是没什么影响。Summing up we can say that, change of domain thickness value does not significantly impact the accuracy of 2D simulations
但是我个人很疑惑,不知道哪位老师可以讨论下
-
本来第一个想猜的是程序里的“体积”“面积”都是广义的,就是“网格尺寸”和“网格边界”,
当三维时“网格尺寸”和“网格边界”分别是体积和面积,
当二维时“网格尺寸”和“网格边界”分别是面积和长度。但是太没有依据了
对着 empty 类型边界条件翻了翻找依据,翻到这翻不下去了 153行读取边界条件类型 ,以后 C++ 知识补上来了再翻。如果全是薄3D假装2D,那 empty 这个类型有点多余了。估摸抠出来 empty 会影响什么就知道它2D是怎么被假装出来的了。
算法上不太懂,代码上挖大概可以这样挖。
公众号里说 axisymmetric 会丢失信息,我记得 fluent 里也有,同时也有 sxisymmetric swirl。后者好像能捕获前者丢失的信息。OpenFOAM 里也有,不知道一样不。
-
@李东岳 在 2D欧拉拉格朗日计算 中说:
几天要推送一个2D/3D模拟的文章,可以关注下
@李东岳 李老师,公众号找了一圈,好像没找到,这个公众号文章链接可以发一下么?最近同学投了一篇文章,是基于欧拉-拉格朗日方法的,做的沙尘障碍物绕流2D研究。一个审稿人意见是他不接受2D Model,并且明确表明如果下一稿还是2D Model的话不用给他审直接拒绝。我就想到之前公众号好像有一篇推文。但是搜不到了。
5/10