@pajack 首先你定义了一个Thread *t，此时它只是一个空指针，然后你就用t去找mixture_thread，它都不知道t在哪，怎么找？thread_loop_c(t,d)只会遍历d下面的t，而相体积率和温度分属不同的domain，建议你用mp_thread_loop_c

想问下楼主最后成功实现仿真了，我目前已经完成了UDF的编写，并且在fluent中也能正常计算，但是液膜厚度一直保持在第一层网格内，不知道是什么原因。希望可以交流一下。

@bestucan 好的，多谢

有没有做欧拉颗粒流液固搅拌方向的学者，可否推荐一下这方面的经典文献

@李东岳 谢谢李老师回复，前几天关了一下，没起作用。我用其他的服务器跑了一下Fluent。

干起来不太好弄得感觉

@akatsuki 是的，但是前面已经稳定了(后面其实也稳定了)，我是读取中间保存的cas和dat文件计算的

我也会这样子，我现在也不知道是什么原因

@徐小双 应该是没有删除干净的原因。

@石油cfder 理论上越小越好，但也不一定，可以参考文献的设置或者文献中的结果来选择

好的，非常感谢

现在高版本只有Schnerr-Sauer模型和ZGB模型

最近老师让做一个垂直轴风力机的模拟，但是没做过类似的，想请问一下各位前辈，怎么在运动中设置旋转之后，增加一个旋转时的阻尼用来模拟有负载下的风机转动,麻烦前辈们了

@bestucan 多谢

@bestucan 实在是太感谢老师了，刚刚去试了下，把单核变多核就能增加新相了。应该是因为服务器是用的多核计算，所以改参数也需要换成多核。

@百善孝为先0626 您好这是用的我们课题组的超算账号，天河那边的

@李东岳 谢谢东岳老师。是的，水气系统，但含有一些其它成分，所以液相粘度会增大一些。我能理解气含率增加，小幅增加比如0.25符合预期，但0.4/0.5有点过高

@doubletrump 您好，我想向您请教，投影面积指的是？如果是一个圆柱，风垂直吹过来，指的是直径乘以高还是圆的面积？

@bestucan 是的 老师 就是直接采用C_T进行读取多孔介质年内的气体温度，然后通过努塞尔数获得对流换热系数，然后可以求解出多孔介质骨架的温度

最近遇到了类似的问题。在做高炉风口回旋区模拟，就是热风和焦炭两相流。热风从炉侧鼓入，炉顶排出；焦炭从炉顶加入，下部排出。所以，炉顶，也就是计算域的上边界，对于热风来说是出口边界，对于焦炭相来说是入口边界。在fluent中设置边界条件的时候遇到了问题。看楼上李老师说，速度可以设置，意思是不是在上部边界选择速度入口，然后热风相给负速度，焦炭相给正速度？

@maud 我上个月也遇见了和你相同的问题，确实提取不到epsilon，但是可以提起到omega，转换一下就好了

@firejet 补一下图9a7cf44e-8063-4afa-ac14-ffff490811e2-Screenshot from 2022-04-26 23-16-44.png be16cb5a-0889-4fc1-b536-d22654dae7d1-Screenshot from 2022-04-26 23-17-57.png 2e8ac1c2-54b9-4b90-84f2-b1ae0e2c515d-Screenshot from 2022-04-26 23-18-35.png 8c2ea078-a7aa-43b1-880e-149f6a54fe15-Screenshot from 2022-04-26 23-33-06.png

@qingcijiao 大佬你好，请问我也对你提出的这两个问题感兴趣，不知道你弄明白了吗？求大佬指点下

@bestucan 在 关于fluentmeshing画网格 中说：

@hitsc30 我没搞过边界层 y+ 什么的。
但是，对于库郎数大于1的状况，fluent 确实比 OpenFOAM 能算，但是能算是怎么算的就不知道了。

如果是高速流动，边界处很难出现垂直界面的速度，大纵横比就没问题，建筑上吹的风，有个回流什么的就……

当然要直接解决就是减小时间步，特别小的那种，保证在垂直界面方向一个时间步也跨不过一个网格。但是估计会很慢。

snappyhexmesh确实玩的可行性比较高，我画了一个1300万的网格，运行了4h，质量还是不太行。目前对网格有点麻木了，谢谢指导！

@cfd小白-0 好的，谢谢你

@李东岳 好的老师:xiezuoye:

@cfd小白-0 好的，非常感谢你的解答，我正在进行中

fluent中的界面只能设置全局Courant数，如果像openFOAM那样设置最大Courant是需要使用udf吗？udf手册里也没找到相关的udf

连续相：水；分散相：硅油（αd=0.01v%）。油在水中分散。速度进口，压力出口。湍流模型选SST K-w。用PBM模拟液滴直径分布。
Fluent的Mesh Quality:（网格质量）

Minimum Orthogonal Quality = 4.58757e-01 cell 63512 on zone 14 (ID: 119996 on partition: 0) at location (-1.46837e-03 -3.11654e-04 2.62275e-01)
(To improve Orthogonal quality , use "Inverse Orthogonal Quality" in Fluent Meshing,
where Inverse Orthogonal Quality = 1 - Orthogonal Quality)

Maximum Aspect Ratio = 3.17826e+01 cell 63834 on zone 14 (ID: 881495 on partition: 12) at location ( 3.37841e-05 -1.15102e-02 2.24530e-01)

PBM设置如下：离散法求解，7个bins，比例指数为2.4，给定的最小分散相粒径为0.4mm,破碎和聚并和均为luo模型。
12.png
Bin Sizes for oil (7) :-
Bin Number Bin Size (mm)

Bin-0 11.143047
Bin-1 6.4
Bin-2 3.6758347
Bin-3 2.1112127
Bin-4 1.2125733
Bin-5 0.69644045
Bin-6 0.4

计算过程中发散，提示如下：
Error at Node 14:

Divergence in Luo Breakup kernel.
Please stop iterating and check solution - parent dia = 1.114305e-02, daughter dia = 3.258247e-03!!

有路过的大佬帮忙看看不？卡了好久，不知怎么解决？提前谢谢了~

往一个体积为30平方米的水池里喷射带有泥沙的水流，开始阶段水池内没有水。计算域中包含
gas（空气）
liquid （水）
solid （泥沙颗粒）
目的：观察泥沙颗粒对水池底部的冲蚀程度。

应该采用 Mixture+DPM模型
还是采用 Euler+DPM模型 ？

各位老师好，我做的模拟是叶轮在水中旋转，用的是滑移网格法，我的问题是我在cfdpost后处理中发现我的流线无法通过我的旋转区域的交界面，我在交界面附近截几个面看各个面的流速矢量图，发现越靠近交界面，流速矢量就越小!以下是我上面所说情况的配图，希望老师能帮我看看我是哪里出现问题的，怎么改正
*

列表-1cdb9ca2093b85da.png -3b3d5b7f16bb0108.png 156efe9894e3d845.png 3ac00c8398c2c821.png -5acc652ec5c2f874.png 602c7722dca5eefa.png -6b94f2b2e4a3d40f.png -3653bc60cd994076.png -732c3ebbc1cd125.png

大家知道如何查找或修改软件中焦-汤系数吗？找了好久也没找到，不知道是不是背后程序里？

@ruibin_tongji 看来这就是算法的区别导致结果的差异，我觉得应该根据和现实情况相符的现象去选择模型和算法。

@jpzhao 有monitor和report的功能，具体看看相应版本的手册哈
捕获.JPG

好歹给个工况、示意图之类的瞅瞅:143:

77811b41-69c8-4b20-8956-a2cbf7487022-image.png
确实，多相流模型和EDC模型是不兼容的，我也遇到了这个问题，持续关注下

@ustbwenwu 并没有

@jamescheng 我研究的问题，颗粒通过圆柱物体与管道间隙的时候，应该认为是稠密颗粒流。我看到过fluent 动网格 加 欧拉 的文章，只有一篇。加dpm的有挺多。目前基本确认加edem不可行。目前我还没开始做，在做别的，后续肯定要做了。谢谢您的指导。

@bestucan 好的，我再去看一看帮助文件，谢谢老师啦:146:

对，那你按道理最多用24个核心

@李东岳 回复晚了抱歉，谢谢大佬提携指导

@李东岳 求老师解答

@yfleung 设置x,y的起止范围就行

@bestucan 您好，谢谢您的回复，我也从CFX的theory guide里面找到了这句话，我原以为这个“available”是指CFX将DDES等封装在软件内部，而不是像Fluent那样可以通过GUI页面来选择而已，另外我从已有的文献中看到有研究利用CFX进行DES仿真也成为DDES。还有690465b4-739e-4bfc-b548-e0a08fd1a3e2-1645959803(1).png 这张图片也是从CFX-solver theory guide里面截取的，我感觉里面的fp就是实现DDES的，只是我们不能进行选择操作而已。
但是，我仍然感觉您说得是对的，可既然DDES要在控制grid induced separation中起到更好的效果，为什么不耦合呢？

因为……要确定加速度吧，比如典型的物理题，得知道物体当前的速度和受力，才能分析接下里的运动。以超声速进去之后得减速，减速的加速度得用静压和流域压的差算吧。

出口，都出去了，也不用管他出去之后静压和外部压差什么的，反正都流域外了。难道出口静压大，所以压出去更快，然后加快出口速度？

我猜的:mianmo:

@cccrrryyy 是的，谢谢。

@bestucan 感谢，我先试试。

fluent 的可以看《ANSYS Fluent Advanced Add-On Modules》，Part V 讲的就是 MHD，其中 chapter 3: Implementation 中开头提到，MHD 计算依托于 UDS 和 UDM。而 UDS 和 UDM 的计算流程《ANSYS Fluent UDF Manual》中有，在其中搜 “solution procedure for” 就能搜到不同求解器的求解流程图。UDS 和 UDM 都在其中。

OpenFOAM 的可以看这个的 3.3.6 节。也可以搜类似的文章。

可以试试matlab后处理:chouchou:

@bestucan 感谢您的回复，我也一直尝试修改模型（尝试了不同网格类型，计算步长、加注速度、入口温度、模型尺寸等）。目前发现影响最大的还是模型尺寸，但是温度先升后降的结果还没有解决