CFD中文网

@李东岳 十分感谢！

@李东岳 谢谢东岳兄！

@李东岳 是的，东岳兄，还是挺划算的。还有招聘链接：https://mp.weixin.qq.com/s/_fZpcKRnCPfULBGYgPgQhQ

同时我们也招收博士生：

澳大利亚新南威尔士大学（University of New South Wales）化学工程系沈岩松教授课题组，现招收2022年入学的全额奖学金博士生数名。研究方向基于多相反应流体系（气-固、液-固、气-液-固），内容包括流动、传递、反应及其耦合过程的多尺度建模、离散-连续耦合模拟、高性能计算、机器学习及其工业应用。具有CFD，DEM，CFD-DEM，LB等一种或多种多相流数值建模和应用经验者优先。详情请见官方网页https://www.promo.unsw.edu.au/upcoming-events 。沈岩松教授现任UNSW反应流动实验室（ProMO Lab）主任，致力于研究多相流模拟和过程设计及其在复杂反应器中的应用，主要研究领域涉及资源和能源工业。目前主持十余项澳大利亚基金委的国家级科研项目，自有科研经费约500万澳元。2019年获澳大利亚科研基金委国家未来学者杰出青年人才计划（ARC Future Fellowship）。领衔的ProMO课题组已在Chem. Eng. J., Chem. Eng. Sci., Applied Energy, Powder Technology, Fuel, MMTB 等国际权威期刊上发表学术论文150余篇。有意者欢迎联系ys.shen@unsw.edu.au。欢迎各位老师专家推荐、转发或交流合作！另：澳大利亚现在无疫情，入境可以申请豁免。

澳大利亚新南威尔士大学（University of New South Wales）化学工程系 沈岩松教授课题组，现招收2021年入职的博士后两名，年薪A$97,368~104,143 (+17% superannuation)。招聘方向分别为CFD-DEM modelling of fluid-particle reacting flow（基于离散-连续耦合的反应流建模）和Data-driven CFD modelling of reacting flow（基于机器学习的反应流建模），详情请见官方网页https://www.promo.unsw.edu.au/upcoming-events 。沈岩松博士现任UNSW反应流动实验室（ProMO Lab)）主任。致力于研究多相流模拟和过程设计及其在复杂反应器中的应用，主要研究领域涉及资源和能源工业。2019年获得澳大利亚科研基金委未来学者杰出青年人才计划（ARC Future Fellowship）。有意者欢迎联系ys.shen@unsw.edu.au。欢迎各位老师专家推荐、转发或交流合作！另：澳大利亚现在无疫情，入境可以申请豁免。

@Ant 你好，我最近比较忙，所以可能没时间给你帮助。我建议你可以先采用CFPD Barracuda这个软件尝试一下你要计算的装置。因为这个软件上手容易，操作简单，用的方法也是MPPIC，你看看用这个软件能否得到你想要的结果。如果还是速度大，说明你的装置结构有问题，比如riser出口到旋风分离器入口的horizontal channel太细了，或者是料腿太细了，或者旋风分离器的vortex finder太细了等，原因很多，因为没有看到你的装置和模拟，所以没法给你针对性的建议，抱歉。多调试看看，MPPIC是个蛮简单的方法。

@nanxuan 你好，可以加入qq群：109526446，群里有比较详细的安装教程。

@shockm 你好，关于CFDEM的安装教程可以加一下这个qq群：109526446，群里有教程。
此外，原则上来说，Liggghts可以和OpenFOAM上很多求解器来耦合，比如，你要是做液固流动，那么可以采用VOF的求解器和DEM耦合。但是比如你要做流化床内气固燃烧，那么涉及到了颗粒的形变反应，liggghts没法完成的，所以虽然有气相的燃烧求解器，还是没法完成整个工作。所以你可以详细说一下你要做什么，有什么需求，我可以给你一些建议。

@yuan 那个CFD-Online的帖子是我当时问和写的，当时是为了做CFD-DEM的燃烧，在2.3.0/2.3.1上面做的，后来也大约调通了，但是当时有一个particletype的问题没弄太明白，后来MFIX的DEM燃烧可以用了，我就主要用MFIX了。OF的就没接着开发了。你的那个错误我倒是不太知道如何处理，@马乔 现在比较精通OF。

ZHOU, Z. Y, KUANG, et al. Discrete particle simulation of particle–fluid flow: model formulations and their applicability[J]. Journal of Fluid Mechanics, 2010, 661(661):482-510.参考这篇文章，就是Model A以及Model B的描述。

http://www.cfd-china.com/topic/81/sprayfoam-只喷固体该如何设置
可以看一下这个帖子，就是稀相的。此外，你的时间步是不是设置的有点过大？

其实我感觉你的这个使用稀相来处理更好一点，因为扩散时候一般不考虑颗粒之间的碰撞。DPMFoam是考虑空隙率的密相求解器，你虽然关掉了曳力以及软球，硬球模型，但是空隙率这一项还是会出现在连续性和动量方程里面的。体积分数没有下限，我说的粒径50微米是考虑到这个求解器是用来求解流化床等密相问题的。如果把它改成稀相的求解器，也就是只有颗粒和流体之间有相互作用力，那么颗粒粒径再小很多也没关系的。OpenFOAM千人群里面有一个昵称是“Champion”的，qq号是：252978797，邮箱：chpjz0391@vip.qq.com的澳大利亚同学和你做的挺像的，他是做通风的，细小颗粒在大管道内通过，会在边角富集。他使用的是稀相求解器，好像根据ParcelFoam改的。你可以和他交流一下。

不能这样扩大啊，因为你的颗粒粒径扩大1000倍，但是流体粘度，颗粒物性这些没有改变，那么计算出来的颗粒最小流化风速就和原来不一样了，同样的风速也许就吹不起来扩大1000倍的颗粒了。所以，我们做scale-up，也就是放大时候，要遵循一定的规律的。比如你可以先把这个几何模型放大2倍，甚至4倍，其他保持不变；下一步几何模型不变，颗粒粒径放大2倍或者4倍，不要一下子上千倍，容易出问题。
此外，manual injection是采用初始颗粒堆积分方式，你计算的是什么背景？流化床吗？颗粒跑出边界是跑出哪个边界，顶部的边界吗？你可以上传一下图片，或者你的设置，这样大家会比较清楚。rebound边界跟Fluent等软件设置一样，亦即颗粒和壁面有粘附、反弹以及穿透三种作用。
DPMFOAM是处理密相气固流动，我认为处理颗粒粒径的下限到50um左右，因为再小的颗粒需要考虑更多的力了，比如静电力等；没有上限的，可以非常大，因为颗粒处理用软球模型，但是注意网格尺寸要随时保证大概是颗粒粒径的3-5倍。
如果对于基础知识掌握不够，那么可以看CFD中国CFDEM部分有个论文推荐，把那几篇论文好好看看，了解一下CFD-DEM的背景知识。

这个问题有好多种原因。首先，你颗粒数减小了，那么你将床体几何尺寸随之改了吗？比如原来堆积是16 mm， 你颗粒减少了，是不是就只能堆积4 mm，那么相同气速吹进去，有可能就无法形成稳定鼓泡状态了；其次，你颗粒数减少到1/10，那么nparticle你设置为100 是什么意思？nparticle的意思是一个代表颗粒内包含的颗粒数，要是传统的CFD-DEM一般把它设置为1，要是设置为100的话，那么颗粒的初始位置文件kinematicCloudPositions里面就要保证颗粒不重合。总之这个问题你要弄清楚物理问题，参考文章Numerical investigation and comparison of coarse grain CFD – DEM and TFM in the case of a 1 MWth fluidized bed carbonator simulation。我建议刚开始时候设置nparticle=1，然后颗粒数减少为原来的1/2，同时把床体厚度减少为原来的1/2。就是一次先改变一个变量，慢慢调试。

@.J. 讲解的很透彻，应该就是这样的。

@hurricane007 OpenFOAM的DPMFoam我添加过化学反应，不是很难的。你好好关注一下lagrangian/intermediate库。

@.J. 感谢你的回复，那么请问在这个里面考虑空隙率和不考虑影响大吗？

是的，可以这么理解，CFDEM的子模型更多一些，比如只针对空隙率的求解，它就有两三种方法，而DPMFoam就是最原始的中心法求空隙率。所以感觉计算特别简单的CFD-DEM耦合问题，以及要对CFD-DEM进行扩展，比如加能量方程啥的用DPMFoam比较好，因为相对简单一些。如果几何模型复杂，或者不对CFD-DEM进行修改，用CFDEM好一些，因为它程序整体复杂一些。

ICEM挺好用的。

这个几何结构我研究较多，有几种方法：1.改变几何结构，不要相切，往里面移动一些；2. 网格不要画那么密，这一块想切的地方粗一点，质量反而提升；3.使用混合网格，其余的地方结构，那个入口到相切的地方用混合。
总之，不改结构，用结构化网格绝对提高不到100%。