@noodles 0.准备几何模型和网格模型。 1.考虑核心物理过程 确定是否考虑冷凝。 不考虑冷凝的话，按照传质过程模拟，可以定义UDM和扩散系数。 2.考虑边界条件和模型细节。 3.计算和结果分析。

@bart 这部分可以从湍流模型的假设和封闭性开始入手分析，我也不是特别懂，而且内容很多。0_1544074389927_第三章，湍流模型.pdf 附件是中科大fluent 讲稿，和fluent 帮助文件类似，可以参考。其他部分可以参考陶文铨数值传热学第九章湍流模型部分。

这个问题提得非常好！图2 中的速度量级上已经到了300m/s 显然是不合理的。 我遇到过类似的情况。当时我模拟2D轴对称情况下的颗粒运动也出现了局部速度过大。 我将对称轴修改为壁面之后发现错误消失。
我推测和维度有关，你可以测试一下3D的情况。 fluent 在多相的情况中还有很多不稳定的情况，需要再探索一下。

@搬运工不好当 这种问题不是特别想回答了，监控你关注的量即可。 设置动画和自动保存。

@书生kao 好哒，谢谢大神，我去试试看。

@东岳 十分感谢，看了一下我的设置还是以物理核心为准吧，而且测试我感觉，多几个核心的速度几乎没有什么差距

@murphy 我主要考虑的是过程时间非常短，传热可能并没有那么迅速，所以不施加两相传热模型反而更加符合实际情况

@anubis 另一种并行方式也用过，无效

不求解动量方程的话，将相变区间比热设置为温度随时间变化即可。开融化凝固模型的话给恒定比热以及相变潜热就行，不用给比热随温度变化的函数

进口的速度不是superficial velocity，进口的速度应该是physical velocity。这俩个之间要换算。换算公式是：

#Superfacial velocity Ug = 0.1 #Volume flux Vflux = Ug*Aout #Inlet disperse phase fraction alpha = 0.01 #Inlet velocity Uin = Vflux/Ain/alpha

这是我用python下的计算器，应该会明白了。你那面的physical velocity就是代码中的Uin

你两相间的传热项选的什么？两相的温度是怎样的？

欧拉多相流的话不是应该是网格尺寸大于颗粒尺寸吗，和气泡也会有具体的关系吗？

老师们，这种关于一个涡的分析图怎么画，涡核能识别出来，但这种有坐标的局部图不知道咋弄

好的好的，支持支持

fluent dpm 中有parcel 这个概念，比如一个parcel 中有10^21个粒子，这样的话，单位体积中有100 个parcel 这样的话，就可以用了。
另外，如果droplet 的属性一样的话，不考虑droplet之间相互作用的话，就可以直接用典型粒子代替。
如果考虑相间作用的话，可能需要其他方案。

在SCI里面用云图和矢量图解释不能做正统的解释，并且云图很容易造成视觉假象。一般是用数据图。

0_1540861155976_捕获.JPG

比如这个图，就是分析不同的阻力模型系数预测的Cd系数，可以看出他们有一些不同。接下来在文章中就可以说

由于Ma and Ahmadi模型没有考虑液体粘度的作用，其预测的阻力系数相对于其他模型预测的要低

单纯给一个云图，比如下面这个图，可以分析下涡的存在、涡的特征、相分数分布的均一程度等。

重要的是要能从理论上解释CFD预测的现象

如果不能从CFD理论解释plot以及各种云图的形状以及差异。那应该是CFD模型理解的不够深。比如有些传输模型，湍流耗散率比较大的时候，会导致产生一个较大的源项。这样在看云图的时候，就可以说

从图中可以看出，A区域湍流耗散率比较大，进而导致比较大的源项，从而变量A在此区域浓度较大。和xxx模型是一致的

比如一些模型会增加一个径向的力，那么在考虑这些模型的时候，你的浓度分布在颈向上应该分布不均：

从图中可以看出，由于考虑了径向的xxx力，其将变量A从壁面向管壁中心push，导致一个径向的浓度分布，和实验值更加温和。相比较没有考虑xxx力的情况，浓度分布差异较大。因此建议在通过CFD分析的时候开启xxx力模型获得更好的结果

不知道你用什么网格生成软件。例如ICEM，interFace界面里外一般是单独做2套网格

@astremon 诶，还没有朝这个方向想过，以前一直udf里面只有只有函数和常量......

嗯嗯，是我的公式有问题，已经解决了，谢谢

你说的对。
QMOM这面重构不出来NDF，对于矩方法如果要重构NDF，需要用EQMOM。

顶一下

@东岳 谢谢李老师的解答，我好好研究下你发给我的资料。

@两月三年 定义的flux可以是与速度无关的。参考这个例子https://wenku.baidu.com/view/af5875254b35eefdc8d3338f.html，这里例子里的flux就是单位向量乘以面积的向量。单位向量是无散的，肯定没有问题。我想知道的是如果不是无散的，为什么帮助里说会导致overshoot或undershoot。

谢谢~应该就是版本bug了
感觉fluent实在玄妙 用了18.0 就算起来了 sigh

@lee459317530 加油了

我用17.2版本的话，图形显示区 title 显示的字体特别小，还挤在一起，模拟过程中显示多个windows窗口的话还不好找，都给挤没见了，用户体验越来越不好了

应该不会吧，我之前用的dell工作站，CPU是e5-2697v3，双路，windows系统，开了超线程，两颗CPU都可以调用，

可能是显示的问题，保存后，关闭，重新打开，查看一下是否正确。

Patch-filed function

来看看

目前就是这个疑问，2D以及2.5D（轴对称网格）的拉格朗日算法，怎么处理体积分数的计算问题？

用Vof 是算的气液两相吗？

@渤文 您好，问一下如果做流固耦合结构会发生很大的柔性变形可以使用overset计算吗？我找到的overset的算例基本上结构都是刚体运动。

@hegui 用瞬态吧，就是计算量太大了

@东岳 按照教程安装ansys 17.2 后，fluent可以正常打开，但是单独打开OCEM时，报如下错误：

ICEM_ACN is "/opt/ansys_inc/v172/icemcfd/linux64_amd". Hallo root. ANSYSLIC_DIR is set to /opt/ansys_inc/shared_files/licensing LD_LIBRARY_PATH is "/opt/ansys_inc/v172/icemcfd/linux64_amd/lib:/opt/ansys_inc/v172/icemcfd/linux64_amd/bin:/opt/ansys_inc/v172/icemcfd/linux64_amd/dif/iges:/opt/ansys_inc/v172/icemcfd/linux64_amd/../../Framework/bin/Linux64". args = libGL error: unable to load driver: swrast_dri.so libGL error: failed to load driver: swrast Window information: depth = 24 class = 4 visual depth = 24 class = 4 X Error of failed request: BadValue (integer parameter out of range for operation) Major opcode of failed request: 154 (GLX) Minor opcode of failed request: 3 (X_GLXCreateContext) Value in failed request: 0x0 Serial number of failed request: 4785 Current serial number in output stream: 4789 Exit from ICEM CFD

我的机子没有独显，用的是集显，会是这个原因吗，还请各位大神帮忙解答一下。

你好，请问你在做EULERIAN和VOF的时候是耦合计算的还是先后计算的呢？在多相流的模型选择上应该是只能选一种吧，耦合的话是需要自编程序来说实现么。我也是在做反应多相流，希望可以互相交流，谢谢！

@东岳 :huahua: got it

有一个东西 我知道 就叫做真空吸尘器
实在看不出真空的论点在哪里

节点如果是一个监测坐标 那么他一般会固有在一个网格内 那么你只要编写一个算法 根据这个坐标来获取这个cell的id就可以了
但是在并行下 这个cell_id是要包含在你的thread_id 还有你的cpu_id之下的 相信你对这块是比较了解的
最后 我们所说的压力 其实都是通过插值法来计算的 你如果想要得出比较精准的压力值 需要采用结构网格 并且在压力离散上 采用QUICK格式
确保上述之后 最后只要把你的所落cell的压力值 也就是C_P（c,t）输出就行 这些方法就很多了 可以用UDMI作为一个后处理输出 也可以导入成为一个自定义的标量
最一般的方法 可以通过一个并行的节点传递给node0 然后node0传递到host 由host打印输出

可以尝试利用ON_DEMAND宏 直接编辑一个dpm的初始化.inj 然后利用dpm的入射方法下拉菜单中的file导入实现一个比较好的初始化

在完成你的求解前 务必详细阅读手册 有任何不明的地方 可以联系我的QQ 18743641
我相信这不是一个一帆风顺的工作 比如ntdll.dll的报错

请问楼主解决问题了吗？
我是占用停在70%后上不去了，不管开多少核都一样

@硫酸亚铜 请问16.0以后的版本是不是没有线程限制了？

@东岳 这一年纯粹UDF排错了 唉

fluent的内部是用C语言写的 到了16.0之后也开始兼容C++
控制台直接兼容shell windows下可以直接用控制台写CMD命令
文本命令采用scheme 早期GUI采用scheme 后期采用C++

这个问题比较复杂，有兴趣的话，可以加QQ18743641聊一下。
轴向，径向，切向还是按照解析几何的方法换算一下就行，只是坐标系的切换而已。这个宏在每个颗粒的轨迹可能与壁面碰撞的时候都会运算一次，所以很容易输出所有的碰撞相关的物理量，唯一的缺陷是这个碰撞是不计时间，只考虑碰撞那一刻的。由于fluent内部trap算法保密，所以很多问题也是很难处理，你不得不自己编个碰撞算法。

@Mellonie
Fluent那面必须指定x轴为轴对称？

@分子在运动 705718754@qq.com,谢大神。 我也在做UDF， 但是老是报错。。。很头疼

是的，可以。

pressure-outlet 水中还会进入空气