研究问题：三维非稳态、恒物性、不可压缩流体 外掠翅片管
模型如图所示：
e0e09295-5ea2-41b2-8475-af5879cebe2b-image.png
边界条件：
上下左右均为周期性边界条件；
固体壁面：无滑移耦合，恒壁温；
进出口：周期性边界条件，给定质量流量为：0.00186 kg/s
ab888039-2684-487e-a691-aae76a9f8350-image.png

采用SIMPLE算法，时间步长0.001s，初始值如下：
c2129438-a7fd-4754-8f2b-700b7f82f886-image.png
计算结果发现：进出口质量流量减小了，变成了0.00103 kg/s （做时间平均取了时均值）
请问，为啥质量流量会减少呢？少了的一部分又去了哪里？

@w251783 好的 我去看看 谢谢

在做大涡模拟（LES），但是使用的Fluent meshing划分网格，可是加密网格壁面，还是说根本实现不了。

@liuxiaodong 如果不用周期性边界就可以用Fluent自带的那些入口了，为什么一定要用周期性边界呢？

什么曲面。壁面曲面的话，法向速度都是0吧。
多相流交界面的曲面法向速度的话，就是相的体积分数梯度与速度向量的乘积。

气体多孔介质的扩散系数是一个与孔隙率和曲曲率有关的量，注意关注这两个量的设置。

不知道你这个模型的几何尺度。
但这种情况可以考虑是否是由于接触角迟滞引起的。即当液滴脱离一个表面的时候，接触角不是一个固定值，有可能小于其在静态状态时候的接触角。

需要写一个file文件，可以自定义颗粒物的状态。关于file文件怎么写，可以参照fluent帮助文件

实际上对于残差水平只是收敛的必要非充分条件，当松弛因子非常小时，残差水平很显然比较小；对于收敛判定应该主要依据：1、监测点的量基本随迭代不再变化，2、质量与热量要平衡

wall film 涉及到rebound、splash、adhesion等物理过程，必然是只对液体起作用。

@李东岳 对哦，我审题不仔细，周期边界的话，哪来的入流。

@装深沉同学 请问，最后您解决了吗？

6ce0cdba-2578-4c0f-ab51-4ceef0b045cb-image.png
配置了一台双CPU的服务器，且每个CPU是52核，Windows系统运行Flunet发现一个任务只能调用一个CPU且不能超过32核，如何解决使其能够调用两个CPU并超过32核的限制？

@tinghaik 个人觉得Gridstar是十分值得学习的，虽有不少缺点也已十分优秀，而且免费不用考虑盗版问题，使用起来与pointwise比较接近。而flowstar主要面对高速气动方向，这个就要看个人的专业了复合性了。我认为这些国产软件是国家队的主流产品，已经在逐步推广，不会轻易虎头蛇尾，有基础的学起来难度也不大，如果专业对的上，何乐而不为。

Fulent里做一个旋转域，类似于一个大圆柱套一个小圆柱，旋转域和背景网格两个计算域都涉及到流固耦合传热，且由于计算域具有周期性，只取30度的扇形区域。请问这样可以实现吗？

@Ding 你需要一句一句的看。
具体方法是：从上至下把无关的先注释掉，一句一句往里加，每一次都用fluent进行编译、初始化、模拟，看是加入那一句之后才出现的问题，这样才能看出是哪里出现了问题。
你这样单独看是很难发现问题的

请问有用Fluent做锅炉模拟的同行嘛？可以一起交流一下~

可以通过TUI实现（在控制台窗口中输入文本命令）

/mesh/modify-zones/copy-move-cell-zone

缺点是先copy再move，可以再用

/mesh/modify-zones/delete-cell-zone

删去移动前的网格域

我一直做致动线的很多年，你有问题可以交流，我Email：jrw1992@163.com。

@Marcllory 还有其他错误要更正的吗？
请提供一下速度分布，压力分布之类的信息
抱歉，根据您现在提供的东西，我并不能看出太多的东西。

@bestucan 大佬，我按照我的udf进行了模拟，结果是符合常识的，应该这种判别方法是可以的，下一步我再用论文验证一下:xiexie:

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

@guanghuilineu 那看起来发散是个很正常的结果，如果相变区速度不为0时，也变速也变密度。

如果只变速，和周围空间交换能量，即动能转换为其他形式的能量，没问题。
如果速度和密度一起变，动能是依赖于质量的。这时候，质量也变了，不守恒了，就崩了。

消失的那部分质量，带有的动能如何传递出去。

压力就是传递能量的桥梁。在传递的时候一头突变了，就断了。可不可以先使相变区域停止流动再变密度。如果这样，不会出现动能突变的状况。或者在变密度操作后立刻更新一下速度；变前计算动能；变后用动能，根据变后的质量，计算新的速度。

或者最粗鲁的方法，用极小的时间步，这样突变也可以稳住。

@知远 感谢告知

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

Inked1_LI.jpg
这个是fluent自带的udf帮助文件里面的一个例子，也用到了NV_DS，可以对于比一下

我去试一哈，感谢大佬:146:

感觉DDPM比较适合你的例子，就是计算慢了些，或者用双欧拉

@李东岳 谢谢李老师。这个问题我当时解决了，忘记回论坛答复了。我的处理办法是，在并行的时候，分节点node_n把所有壁面信息如壁面摩擦速度、壁面摩擦速度传递给node0，由node0统一处理把全部壁面信息传递给所有分节点node_n，然后分节点分别计算壁距离、y+。思路是这样，写在这里供大家参考吧

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

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

@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°，那最开始设置气相区域的时候需要留出一个液膜厚度吗？还是直接把整个管径宽度都设置成气体都可以？我比较担心如果把整段管都设置成气体会让气体“粘”在壁面上。
谢谢！

话说什么叫较大体积力？考虑重力的气液两相流动算吗？反正我看到文献力一般处理这种都是PRESTO！

我被这个问题困扰了很久...我发现局部气含率特别高的时候（云图上出现红色区域），计算不守恒，液相质量减少。有遇到过类似问题的cfder吗？欢迎交流!

报的错都没贴出来，人眼debug啊:135:

从m_lg=0开始试，一项一项往里加东西，看加到哪一项出的问题

能，模型是否适用和这个没太大关系，和状态有关系。
如果流动特别慢都没有气泡，就不适用了。

@freeazma 如果我想在 fluent 里面同时查看/获取/输出固体和流体温度，能实现吗？

感谢各位老师和同行的鼎力相助，目前入口最高压力能够达到30MPa，和实验匹配效果较好。
0.2271-200steps-improve2.jpg
0.2271-200steps-improve2-速度.jpg
最后再一次感谢大家！

@cccrrryyy 有可能, 但是确定流体域中大部分为湍流，只有细缝附近不是很确定，这种在Fluent里面如果一部分为层流一部分为湍流，感觉就不怎么好处理了。

@b占余文乐 CFD-POST里面倒是能创建一个区域，但只能是球状区域，也是有限制。

@cccrrryyy 好的，谢谢:140:

@wseric 嗯嗯，fluent自带后处理是可以做的，但出的图不太好看，tecplot出图漂亮一点。目前，基本可以解决这个问题，就是速度三角形换算问题，公式推导出来编公式计算就行。

@peak_han 建议使用质量入口条件，由于入口处的压力分布不好描述，可以延长流体入口长度，保证流体充分发展。

@luoluo 感谢回复，我有一个疑问，由于脉动速度的长时间相关总是等于零，经过足够长的时间后初始场的随机状态对湍流场的发展几乎没有影响，因此对初始场添加扰动的原因还有些疑问，为什么要添加这个扰动，如果不加扰动，直接用RANS计算得到的时间结果作为初场，进出口也给周期边条，采用大涡模拟能否自发生成湍流？之前算过，好像是产生不了湍流的，只是其中的原因一直不清楚。

求解颗粒流有很多方法，成熟的方法fluent里有现成的就不用udf了。如果写udf当然看用什么方法了。

围观中

边界设置好气温不就行了。传热就三种途径，辐射、对流、传导，空气的传导热作用比较弱。鸭绒就是通过绒毛阻止空气流动形成空气墙保暖的。边界处的空气除非进入流场参与热交换，也不会有辐射（假设太阳在头顶……），传导又弱。

没发现解释就自己解释啊，解释清楚了往论文里写:chigua2:

Ansys19.1 student，网格超过51.2万时的报错信息：

Error at host: This is an educational executable, and can only be used
with cases containing less than 512000 cells.
Please exit this fluent session and start another session to continue.

Error: This is an educational executable, and can only be used \n with cases containing less than 512000 cells. \n Please exit this fluent session and start another session to continue.
Error Object: #f

e09cef93-36a6-411c-8259-623e46a90966-image.png

你这个是三维的LES单相流模拟，应该不难分析，常见的不收敛的原因有很多，上面大神们说的时间步长、网格质量、松弛因子都是很不错也很常见的调试首选项，还有就是边界类型也挺关键，再就是几何域合不合适，上下游延伸段够不够，比如存在很狭小的流动域之类的、动网格更新策略，都可以作为改进收敛性的切入点。

@bestucan 谢谢您的分享，这些我之前看过一点。但是这个的定义就还是在初场中定义的湍流积分尺度，用到的是跟湍流强度等有关的计算。但是这样的公司似乎不能用在模拟结束之后的该参数的计算，具体为啥我感觉我也说不上来:135:

https://www.dataphysics-instruments.com/us/knowledge/understanding-interfaces/surface-tension/
这里有一句

Thus, for liquids the surface tension and the surface energy are equal

至于动能，k就是湍动能啊，如果不只是要湍动能，那k怎么来的肯定要涉及动能啊。