CFD中文网

@allanZhong 另外曾在cfd-online上看到过一个帖子，里面众多大佬对youtube上xflow的视频口诛笔伐，主要质疑的点是网格不能局部细化，导致远离避免的流场也是被划分的超密，对计算资源、存储资源是一种浪费

@tens 目前来看确实如你所说没有这个功能（VOF turn to Larangian）。

@mohui 您好，请问您解决这个问题了吗？我想给同种液滴的两个液滴染色，便于观察融合过程中的各个形态变化，求指教。

而且对于稍微复杂一点的结构，边界层经常出错。反正人家重心也不在这上面

@mohui 如果使用湍流模型，那么就是模化（或部分模化）了流动的脉动。真实流动的脉动越强，湍动能也更大；但是你计算的时候是不会算出来流动的脉动的，只会得到较大的k值。
你在k方程加扰动，对结果影响应该不会很大。

of5.0里面已经有fvOption了，

tmp<fvScalarMatrix> epsEqn ( fvm::ddt(alpha, rho, epsilon_) + fvm::div(alphaRhoPhi, epsilon_) - fvm::laplacian(alpha*rho*DepsilonEff(), epsilon_) == C1_*alpha()*rho()*G*epsilon_()/k_() - fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_) - fvm::Sp(C2_*alpha()*rho()*epsilon_()/k_(), epsilon_) + epsilonSource() + fvOptions(alpha, rho, epsilon_) );

或者复制模仿realizableKE的格式，赋值于epsilonSource() 或者 kSource()．

@mohui 那我感觉应该可以用，我只是知道，没用过，期待你的结果

最近也基于pisoFoam植入了一下RKprojection foam，单个计算步不限制残差，projection method的确快很多，但是限制残差之后速度就变得特别慢，并且该算法内存占用量是pisoFoam的20倍。
不知道有没有大佬看到过把这个算法植入到多相流求解器的呢？想尝试一下把这个方法放到多相流计算里边去。

好的 谢谢 我去看一下

这个问题使用fluent算貌似同样会出现这个问题，其实这时候已经出现了部分空化，处理不当产生负压。

还有这种神操作 :duang:

可以啊，Why not？

我插一句：openfoam这个compressibleInterFoam的温度方程，应该是多相流里面最后一个尚未完善的地方了。温度并不是守恒变量。理想情况应该求解e方程。

你要清楚你想实现的语义是什么，再去找具体的语言机制。

虚函数的语义是多态，也就是你有一个基类的指针引用可以自动调用到实际子类的函数。

纯虚函数的语义是一个接口。

有纯虚函数的类是虚基类起的是定义调用接口的作用，同时虚基类不能实例化，对编程也是一种限制性保护。

p.s.简单一点儿，都整成虚的也可以。调用开销通常没那么大。

我也忘了..

@mohui 谢谢哈，

@mohui 感谢你的解释

如果单纯的mules求解每一相是无法保证alpha1+alpha2=1的

现在我就想确定，alpha2通过对应方程推导，用MULES求解不可行的吗？

MULES的出发点是保证变量的有界。所以如果你用MULES求解alpha1，那么理论上alpha1是有界的。alpha2的求解没有调用MULES，他是通过1-alpha1算出来的，如果alpha1有界，alpha2也有界。

如果你用算法同时求解alpha1和alpha2，如何处理耦合？这俩个变量是耦合在一起的，就像速度和压力。你不能单独的去分离求解，如果分离求解就就需要迭代。迭代就导致计算速度变慢。因此现存大厂据我所知都是只求解alpha1，然后alpha2=1-alpha1。

我在去你年底验证了MULES，理论上是可以保证有界，但是真实计算的时候，还是会越界。尤其是物理模型比较复杂的时候。个人觉得这方面内容搞出来，绝对是个好文。