CFD中文网

想开个贴讨论下多相流界面追踪类方法目前的前沿水平，例如能算到的最高参数，各类方法都有什么限制等等。目前多相流算法里，计算界面追踪的多相动力学，比如droplet collision, droplet splashing, droplet impact on solid wall，这几个典型案例，目前有用VOF，level-set，phase field，LBM, SPH, 欧拉拉格朗日耦合其它模型的。这几个案例在不同Re和We下会产生许多种演化形态，实验都做出来了，但感觉计算还没跟上。

调研过不下百篇的文献，基本上关键词相关的都翻了一遍，在文献中看到droplet collision大概最好的也就是Re=6500,We=1500，密度比DR=1000左右的，大致是这个级别，像droplet impact有篇08年的用TransAT算到了RE=11000，而最近16年左右的用levelset这些也算到Re=8000，We=1500左右。大多数都是采用了一些改进的数值技巧来稳定模拟，而且网格也都比较密，计算资源消耗也都比较大，估计都得上超算，当然和实验也都对的上，基本都是到了prompt splattering这样的形态。

想讨论下大家都做过或者见过的这类模拟最高参数是多少? 或者往上算各类方法的限制在哪？ 因为多相流在算到高RE和We的时候会数值发散，稳定性应该是主要限制，还有例如VOF和level set会在薄液层处由于数值格式产生人工破碎的非物理情况，需要单独使用人工识别薄液层再稳定的方式，对不规则碰撞并不实用。

wedge几何只能用wedge，wedge不能用symmetry。wedge与symmetry的流动特征不一样

deb046f8-a6ad-47bb-8824-432519675318-9110c9a264e0aa2a0bccdf0b156cb74.jpg

请问还可以用什么命令定义颗粒的直径和密度呢？
create_atoms和set吗？
是的话怎么用呢，求大佬举例，不胜感激

@李东岳 好嘞，感谢信任，我是很喜欢咱们的笔记本，希望做得更好:146:

@东岳 原来是连续方程的引用啊，谢谢李老师啦

请教一个比较初级的问题，查了帮助文档没有找到想得到的答案。
ICEM采用Octree方法划分非结构网格的时候，帮助文档和网上的一些教程中都写到：“建议全局网格参数设置中的Max element值要是2的幂”
这里2的幂是指：2的正整数幂？2的正数幂？2的实数幂？
实际测试过感觉比如填0.05啥的也可以划分，网格也确实要比2的密很多。所以想问下各位高手这里这个数值填写应该有什么经验或者规律呢？尤其是模型尺寸很小的话，有时候确实需要很小的尺寸啊~

用于求解不可压的pseudo-compressible method, 被认为有更好的收敛特性，见A. J. Chorin, JCP, 1967, pp. 12-26

http://www.cfd-china.com/topic/2744

我重新推出来一个带偏导数的压力方程..

ef44221d-062f-4315-b9c1-50e96688d88e-image.png

一个公式真是折腾死了。

我看了好几个人讲关于von Neumann稳定性分析。

12步入门Python这个女老师讲的。 MIT qiqi wang讲的PDE课程

还是有些地方很模糊。

@sxz0823 你好，请问有关concentration和parcelconcentration您是怎么计算的，求指导。

如果能在CFD界重新发一下，那真是再好不过了。万分感谢李老师。

确实和概率论有关，谢谢，记录一下

\begin{equation}
\left( \begin{matrix}
k\\
x
\end{matrix}\right )=\frac{k!}{(k-x)!x!}
\end{equation}

@hhj0219 和y+定义一样。其实就是：流向网格间距除以壁面第一层距离，然后再乘以y+。

是用八叉树算法吗？设置part上的最大尺寸也需要是2的幂吗？对这方面有了解吗？我也有这方面的一个疑问呢，之前在论坛发帖子还没人回。帮助文档里面也没写清楚，网上查了一些也没查到

@李东岳 在 无线鼠标掉帧 中说：

最近10年一直用mx anywhere系列，用了3个，都有无线鼠标掉帧问题（我怀疑是跟无线网卡一起用有信号干扰，然后每次插硬盘必掉帧成狗，拔了就好，插上就掉帧）。也懒得处理。有时候卡顿到不行。不知道你们用其他的无线鼠标有没有这个问题。anywhere主要是那个滚轮比较好，pdf可以一刷刷一大片。

东岳老师遇到这个问题极大可能是主机不同的信号互相干扰了，买一个usb延长线，把2.4G接收器插到延长到桌面的usb口就能解决这个问题

@l-j刘侃 目前这个现象只是刚刚发现，暂未做深入研究。李老师对流动本身就是三维的论述点醒了我，流动本身就可以分解为平均速度和脉动速度，所以，二维发文章还是慎重点好。

感谢分享！

@suzuhaz 大师求助帮忙看一下简单的一个模型块的划分0_1532943202163_05.zip

@东岳 我最近用reactingMultiphaseEulerFoam做算例，也是设置的时候某一相为0，湍流模型一直出现divide by zero，找了半天没发现错误，更改之后搞定了，谢谢东岳老师