CFD中文网

想要计算的物理问题是跨临界流体（二氧化碳，或者是氟利昂制冷剂）在大范围压降条件下的相变过程，说人话就是比如入口是超临界的CO2，经过Laval nozzle膨胀到亚临界；
对于这种情况必须考虑真实流体热力学方程了，目前我用的实现方法是通过调用NIST REFPROP生成的real gas property lookup table的办法来计算气液相流体的热物性，这个方法在计算亚临界范围内，大压降的工况也可以计算，收敛情况相对满意，但是一旦计算跨临界状态就发散了，在外国的CFD论坛上找过，类似的答复是NIST流体数据库使用的方程在临界点附近似乎存在不收敛的问题，如图所示:

1c5761fa-95be-44a0-a527-4e2f418b4098-9bacb2a3045cb6950be3e5620394144.png
现在很苦恼，一直卡在这个问题上面，请教一下各位大神，在fluent里面计算跨临界流体热物性还有什么别的实现方法吗

你好，请问您解决这个问题了吗

@李东岳 老师，最近我也遇到了这个问题，单相流计算的时候会提醒，但是结果也可以收敛，但是计算完初场的时候，改为多相流空化计算，马上就发散了，不知道怎么可以解决这个问题，已经重新画了好几次网格了

应该采用laplace 光顺而不是弹簧光顺

在WallSpringSliderDashpot.C中，有这样几句：

scalar etaN = alpha_*sqrt(p.mass()*kN)*pow025(normalOverlapMag); vector fN_PW = rHat_PW *(kN*pow(normalOverlapMag, b_) - etaN*(U_PW & rHat_PW)); ... scalar etaT = etaN; ... (if条件：颗粒不滑移的条件下) fT_PW = - kT*tangentialOverlap_PW - etaT*USlip_PW;

这里的颗粒碰撞中的切向阻尼系数，etaN，等于法向值，然后分别应用到求碰撞力的切向和法向分量，fN, fT，中。是不是意味着，我在颗粒配置文件中通过alpha设置的回弹系数值e, 就是同时设定了切向和法向的回弹系数，eT, eN, 都等于e？

请问这个理解对不对？

好的，谢谢岳哥

@李东岳 确实

我觉得STARCCM+的网格功能挺6的，笛卡尔/多面体等，控制的参数有限，基本上零基础跟着教程做几遍就差不多了，网格生成速度也很快，比较复杂的几何还是很容易解决的；
其他的用point wise多一点，复杂的几何自己划分块有点麻烦，如果不能分层很容易混乱，简单几何还是非常便利的；ICEM复杂几何切起来还是比较舒心的；^_^

@星星星星晴 谢谢，我再看看

可以考虑ParaView反向代理，就是你的显卡性能弱的服务器连接你的带显示屏的电脑，把渲染和显示工作放在你的电脑上处理

估计网格质量不太好

可以相信这一家：https://www.phoronix.com/review/amd-ryzen9-7950x3d-linux/

简单的点线面操作撤销是没问题的，检查拓扑或者是操作包含的数量太多ICEM就没办法了，解决这个问题的最好办法就是你在做包含多步操作的问题时候先另存为一下是最保险的，问了一下学长他也是用这个方法的

IMG_20231127_193223.jpg
在计算单相流作为初场，可以收敛，但是湍流粘度比超限，继续设置为多相流空化计算的时候，开始发散，不知道有哪些解决办法，已经重新画过网格。
边界条件是速度入口，压力出口，湍流sst k-w，稳态计算

网格流动性是什么..

如果能在CFD界重新发一下，那真是再好不过了。万分感谢李老师。

您好，请问您这个问题解决了吗？能分享一下解决办法吗？

顶一下，如果有同学想组队团报，可以联系QQ:876440251
我问了客服小姐姐，说优惠可以享受hhh

文章题目为《The inviscid impingement of a jet with arbitrary velocity profile》，主要是分析了法向二维冲击射流流场的解析解，但是个人比较疑惑的是为什么控制方程只包括了压力的泊松方程（如下）而没有包括二维的欧拉方程？
d64e75d4-5884-42a4-b208-20c90325cc29-image.png

全文见附件，希望有朋友可以解答。The inviscid impingement of a jet with arbitrary velocity profile.pdf

@李东岳 在 学习各种算法的用处在哪 中说：

我那个笔记非常basic，才60多页，后续还要学习非常多非常多的东西。我最近打算增加点内容，循序渐进。

2024年，4年过去了，现在333页了。