欧拉拉格朗日方法的问题

[gxf]

标准的EL方法一般指的是PCM，相当于源项只作用于颗粒质心处的网格，如果颗粒接近网格尺寸或比网格大就会造成一定误差或数值发散。
我看到很多固体颗粒燃烧的场景，会对PCM方法进行改进，比如DBM（利用扩散方程使颗粒的源项影响颗粒周围的网格）
请问各位有没有在OF上改进这种方法的经验呀？十分感谢！！！

[李东岳]

有没有sci发来瞅瞅 老铁

[gxf]

https://www.sciencedirect.com/science/article/pii/S0016236124007488
https://www.sciencedirect.com/science/article/pii/S1385894719333790
李老师，请您看下~

[李东岳]

我粗略看了一下这个，https://www.sciencedirect.com/science/article/pii/S1385894719333790
这个文章也是参考Capecelatro et al.之前发在JCP的方法。这个方法就是把源项的影响抹平到其他网格了：

1. 把源项创建volScalarFIield
2. 求解源项这个扩散方程

文章里面求解的是方程14。我在想的是另外一个思路，可以通过滤波的方法来进行。 在LES方法中的滤波函数可以看出把某个场在局部进行光顺。你也可以用滤波的方法来试一下。请参考无痛苦ns方程笔记中的LES 中的simple滤波函数以及LES 中的laplace滤波函数，这种方法不需要求解PDE。

我这个思路没见过人用过，你也可以两个方法都验证一下，发个文章是没问题的

[gxf]

好的，谢谢李老师！我去试一下

[gxf]

@李东岳 李老师，打扰一下，我计算了一个孤立液滴蒸发的case，发现网格尺寸对于蒸发速率以及液滴表面温度影响不是很大，和文献中的不太一样（蒸发速率和液滴表面温度对于网格很敏感），不知道这是因为什么，请您指导一下，可以吗？谢谢
结果如图所示：

其中x/d0代表了网格尺寸与液滴直径的比。
参考文献为Abramzon的那篇文章：https://www.sciencedirect.com/science/article/pii/0017931089900434 （Figure 2 (a) (b)）

[李东岳]

结果算出来，趋势没啥问题，但是对不准这个事很难搞。我一般也没啥头绪只能一点点debug哪里有了问题。

[gxf]

是的，趋势没啥问题，但是看起来这个case对于网格的依赖性没那么强，网格和液滴直径的比值小于1的情况也和另外的基本吻合

[Mll]

@gxf 您好，请问您计算孤立液滴蒸发的case使用的是哪个求解器呢？是否进行了重新编译？我之前试算过，结果误差很大。