CFD中文网

Liu et al., 2016 "A Comprehensive Benchmark Suite for Simulation of Particle Laden Flows Using the Discrete Element Method with Performance Profiles from the Multiphase Flow with Interface eXchanges (MFiX) Code" 一文中对MFIX-DEM的并行能力进行了验证，采用了Weak Scaling和Strong Scaling两种评估方法，结果证明在核数达到1000左右，MFIX-DEM code的并行能力还很不错。但是注意到，在做自由堆积等的时候，文章回避了负载均衡的问题，也就是如果流化床下半部分非常密集，而上半部分非常稀疏，这是对就需要一种新的分区方式，使得每个核都能计算差不多的颗粒数。现在的MFIX好像还达不到这个功能，但我知道有人开发出来了。文章还认为，当颗粒数达到10(8)，即使用1000核，每个核也有10(5)个颗粒，还是非常有挑战性的。气固两相流的计算时间主要花费在了DEM颗粒彭碰撞计算，颗粒和流体之间的数据交换，曳力的求解。
看到文章的一点心得，和大家分享一下。

同问,我也正在做类似的东西,但是算到最后连续性方程都是收敛的,但湍流方程总不收敛

同问，目前在用PBM处理气液两相流，用discrete方法求解PBM，也想试试qmom方法，上面几个矩代表什么意义？要模拟气液流中的气泡，该怎样设置粒径，数量还是。。。？

并没有用过pyfoam

@winsway_zero 前辈，请问gmesh怎么划分block

明白了，感谢各位的回复

这几天备课OKS，今年大幅度增加湍流的内容。前几天也在朋友圈发了个图。简单说一下。LES处理壁面几种方法：

LES-RANS混合方程：用RANS求解壁面区域，LES求解主流区。也就是你需要调用DES、DDES等之类

Wall-modelled LES：对壁面区域进行模化，主流区采用LES。这个进一步还区分不同的模型，比如代数类、PDE类、ODE类等。

http://www.dyfluid.cn/theory.pdf 混合方程我在这里更新了一点内容。

另外没理解你说的壁面雷诺数什么意思。

捕获.JPG

Ref: Journal of Computational Physics 305 (2016) 589–603

fvm::xxx(u)，其中的u是需要求出来的，fvc::xxx(u)，其中的u为当前时间步的值，其返回一个场。所以，需要你要求某个场u，就用fvm。

@ibelief 谢谢，我再看一看。

@搬运工不好当 说：

OpenFOAM2.3.0中有buoyantPressure这个边界条件吗?没有的话要用哪个边界条件来代替 buoyantPressure呢？

OpenFOAM 2.3.0 中没有buoyantPressure，替代的边界条件为fixedFluxPressure，详细可参考cfd-online的这个帖子。

把define_source里你想用的量给赋给uds，然后在define_adjust里调用，这样的量还能在后处理里显示。或者粗鲁点，直接定义个全局变量，这里赋值那里读。

你说的pri_th这是应该是这些源的名称，第一相第二相混合相。名称怎么变都行。他们起作用的时候是在计算的时候把他们选成谁的源。

谢大婶答。
一开始看到底下三个按钮直接崩溃了（原谅我两天两宿基本没睡觉，实在冷静不下来）。。后来发现了是给定的波浪参数有问题。。。
PS：原来都说C语言亦或是Python啥的从入门到放弃，自从我看到了Fortran里面的各种变量定义之后还没有注释，要求硬生生理解看懂程序，于是。。。
嗯，Fortran这种东西，哇，真是美妙。
想退退不了呢，(^o^)/~

雷诺数从物理意义上说是惯性力比上粘性力，想问为什么雷诺数的定义会涉及到惯性力这么一个虚拟力，是怎么考虑的，对于管道流动来说，如果是匀速运动，则流体的压力和粘性摩擦力平衡，对于所研究的流体单元来说是不存在惯性力的，如果管道流动存在加速度，则流体的压力和粘性力的合力导致流体的加速，可以采用一个虚拟的惯性力来等效流体的加速，则可以认为流体在压力、粘性力和惯性力三者作用下平衡。 所以，雷诺数为什么没有压力项，而用一个虚拟的惯性力来定义？

@linhan-ge 在 多相流VOF 中说：

VOF-DEM Study of Solid Distribution Characteristics in Slurry Taylor Flow-Based Multiphase Microreactors
强啊前辈，我这进展还缓慢着呢
刚看了您“VOF-DEM Study of Solid Distribution Characteristics in Slurry Taylor Flow-Based Multiphase Microreactors”的模型部分，请问您是在interFoam上加了Laplacian filter就可以抑制假速度、算出液膜了吗？因为我现在就做气液流动，速度也跟您差不多（1cm/s~10cm/s），所以想请教一下这个方案可行性。我直接用interFoam算的液膜厚度非常厚，不知道问题在哪里。
另外我看到您是设置接触角90°，那最开始设置气相区域的时候需要留出一个液膜厚度吗？还是直接把整个管径宽度都设置成气体都可以？我比较担心如果把整段管都设置成气体会让气体“粘”在壁面上。
谢谢！

又研究了一下明白了：如果手动添加各项异性的全隐式源项，
需要拷贝一份A()并且构造成为张量，然后添加隐式源项，之后再构造HbyA，这样就能将各向异性的源线用全隐的方式添加到面通量上面了。

@evensun 能在高校待时间长的，必然都有外快。各有各的路子，有些大家能看见，有些大家看不见。否则待个5 6年必然就走了。

@wwzhao 嗯嗯，明白了，非常感谢您的解答！:xiexie: :xiexie:

什么是粘性边界条件，粘度的边界条件？粘度不都一样嘛 :quwan:

@李东岳 同一个计算域的网格，有单元的节点信息，单元相连的单元信息，面信息。但是两个不同计算域网格之间没有这种拓扑信息。这种情况下如何并行剖分网格同时保证负载平衡，减少通信量

@veen 从源码编译of是个挺痛苦的事情，建议不把of安装在root，可以考虑试一试easybuild，它是基于module的环境管理，应该比较成熟了。我从0开始，最开始在PC上面编of-3.0，到后来在实验室的机群上home编译of-4.0, foam-extend等等，多版本的管理用module还是方便些。印象中，gcc（如果你用gcc的话）和mpi的版本很重要。