CFD中文网

0_1516841984913_捕获.PNG
动图1：河床沙石冲刷运动

河流的流动会改变河床的形状。底部的沙子以及大石头通常会分为两个区：

顶层的bed-load zone（有研究河海的么，这个怎么翻译？），其中的颗粒会随之移动； 底层：沙子以及石头会被困住，但是会缓慢的爬行；

在碎石较多的河床中，较大的沙石倾向于堆积在顶层，这种现象被称之为armoring。个人感觉之所以称之为armoring是因为较大的沙石趋向于在表面，类似一种武装。

实验表明，在这个区域，大沙石有一个净得向上的运动速度，小颗粒具有一个小的向下的移动速度。 目前对于这种现象的原因尚不清楚。不过现存了一些潜在的理论。例如支持程度比较高的巴西果效应（Brazil Nut Effect）。

0_1516841984913_捕获.PNG
动图2：有人晃荡大米里面的密度差不多的小石头

替代文字
动图3：进行的计算模拟

巴西果效应是指如果把两种颗粒的混合物置于容器中，然后施加外加的振荡，体积比较大的颗粒会上升到表层，而较小的颗粒会沉降到底部。关于这种古老的效应的动力学机制至今仍众说纷纭。

DEM模拟：文中作者采用LIGGGHTS对这个现象进行了模拟，下图可以看出，随着时间的推移，大颗粒漂浮出来。

0_1516841984913_捕获.PNG
图4：Nature Communication中的LIGGGHTS模拟

看起来很有意思 :cheeky:

本文来自Nature Communication

可以加载多个

雷诺数计算和其定义有关。比如槽道流动中特征长度采用的是槽径，平板流动就是板长，而颗粒雷诺数就是采用颗粒粒径。
你搜过其他论文里面的选取原则吗？

前期用sprayFoam计算了喷雾燃烧，审稿人提问计算过程中是否限制了parcel的最小数目和最小粒径，请问大家知道OpenFOAM中怎么实现的么？
当时并没有考虑到这方面的因素，只是根据文献（Srinivasan et al., Quality and reliability of LES, Springer, 2011）计算了下每个网格对应的parcel数目小于8，便开启流场计算了。目前的喷嘴设置模型如下图，破碎模型选用默认的ReitzDiwakar Model
720bc43b-0cf4-4576-abd2-28f1f77be5b9-图片.png

@cloud1990 在 VOF求解自由水面问题，求解空气区域吗？ 中说：

Image

Left: Correct jet shape predicted by TruVOF technique used in FLOW-3D. Right: Incorrect jet shape predicted by pseudo-VOF technique used by other CFD codes.

Image
Left: FLOW-3D‘s TruVOF technique predicts jet impingement on wall and some outflow. Right: Pseudo-VOF methods don’t predict realistic jetting of fluid on side walls.

从结果上看，flow3D 所谓的TruVOF是有问题的，air velocity 和　air resistance = 0?　这个可能比较适合粗网格计算，作为一个tricky的存在．但理论上，这个TruVOF是无法解决wind wave的问题．

请问大家有没有遇到过多项流计算案例中内部场的速度量被zero覆盖的？谢谢

已经收到了

https://www.cfd-china.com/topic/5198

参考这个帖子，有人可以做个STAR-CCM、COMSOL、CFX的算例么？

@bestucan 好的 谢谢提醒 我试试

@李东岳 谢谢李老师

看起来确实没辙了

@luomuliunian 大佬这个专业了

近壁面网格长宽比本来就很大，如果用RANS计算，有的模型对近壁面y+要求很小，比如SST k-omega要求y+≈1，但是对x+或者z+基本没有明确的限制，为了减少网格数量，一般流向和展向的网格无量纲参数可以在百千量级，这样的话近壁面网格长宽比就会很大，也是可以正常计算的。

@李东岳 李老师，有关boxTurb中需要初始给定的Ea和k0这个值就是取默认值10和5吗？还是说有什么计算方法呢？

@李东岳 @lsprxd @赵一铭 对于标准壁面函数我是这样理解的：在求解流场的时候还是用原来的k方程和epsilon方程来求解，只不过是将在壁面第一层网格（有可能是网格面上的值，有可能是网格中心的值）上的信息用壁面函数求解的结果来修正，而且壁面函数其实只求解壁面第一层网格的信息，第二层网格之后就是输运方程求解的结果了。不知我这样理解是否正确？

请自己看手册。。。
祝早日脱离苦海

@李东岳 李老师，我使用的是指数率的风速，参考风速11m/s，高度0.4m。之前用LES算过yPlus是符合要求的，第一次用RANS算。

@李东岳 http://dx.doi.org/10.1016/j.compchemeng.2015.09.008
这篇文章对三种重构方法做了对比。