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怎么来的肯定要涉及动能啊。

好的，谢谢

@zhangshengchang 谢谢！:146:

@aaaachenjian 纳米尺度是不能用Fluent来模拟的。首先不一定满足连续性假设，其次要考虑分子间作用力等情况了，可以考虑用MD 分子动力学模拟，一个工具就是上面@zhangshengchang 提到的Lammps

@zhangshengchang 在 请问在fluent2020R2中怎么设置v2f湍流模型，我在viscous中没有找到？ 中说：

allow-v2f-model

捎带着多说两句，老版本中fluent 命令行实际上是Lisp 语言的一种方言 Scheme， 其一个特点是（干 事） （）里面就是命令。
激活v2f命令用的是 （allow-v2f-model）

@zhangshengchang 这个真不知道 哈哈哈 谢谢前辈 我去试一下

这已经是比较成熟的方法了 应该aixprocess在15年前就已经成型了 做木材洗涤工艺的除尘优化 如果说无脑的办法 还是直接写碰撞列表做加减法 也可以直接设置某一方质量半径都为0 然后另外一方做简单的质量累加 但是实际上真实的情况比这复杂得多 液固的碰撞的结局以及最后的补集情况是多样性的 这有很多相关文献 但是对于UDF编写来说 因为是不同injection 还是有一定难度的

射流 假设你是超音速流体 那么这是一个非常复杂的问题 这个问题是直接和初始条件相关的 对于瞬态 就是每一步的上部计算结尾的状态 因为它是对于初始条件非常敏感的 所以一个微小的积累可以导致很大不同的后续结果 实际上就是一句话 超音速这块实际上还是很难控制的 最好还是辅以实验 一般优化效果需要达到15%以上 才能说是设计真的有优势 如果要做严肃精确研究 fluent目前对于这种还不是太适合 但是fluent的优势是可以综合多场一起做出一个工程上可以大致量化的结论

vof模型只有一个混合相的温度 所以你这个方式其实难以实现 现在新出来了一个欧拉-欧拉-vof 可以尝试一下 实际上更好的办法是直接采用欧拉-欧拉同样算一遍 这样你可以方便的获取出口温度 实际上这样的方法也不会有太大的误差

除开物理上的不守恒性的考虑 还有计算时的不收敛考虑 你可以通过fixed values的方式来操作 这种方式对于某些燃烧反应器还是可以用用的47207c3e-7358-418a-99bc-b76a77800a0a-image.png

连续方程 如果你是在模拟一个工程实际设备 这种不能达到-4的情况是比较常见 通常在-2左右 这代表在你的设备设计中间 真正的运行的时候 设备并不能保证一个进气和出气均衡的稳态（比如内部有一部分涡旋致使一部分进气永远在这里堆积） 而其他方面的原因上述几位已经讲解过了 对于工程模拟来说 连续性不达标的问题是比较常见的 你可以看看你的设备内部是不是存在涡旋流 或者是因为本身你的边界条件的设定是与实际不符合的 比如你在一个明显不均匀的出口面取了一个平均压力出口值 这样有的时候你需要扩充出口 直到真正符合工况的外界条件

TP_N 输出的是在一个颗粒簇中间（一般dpm都是用质量流量来算的）的真实颗粒数量 但是这些东西在17.0版本之前好用一点

@caizhao https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1002/2016JD025316
可以参考一下这篇文章
上面链接打不开的话，可以用下面这个
https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2016JD025316

有偿求助！看了些论文，只看懂在多孔介质区域设置为非平衡项。，具体的参数又是如何推导出来呢？

@发芽的土豆 :
感谢赐教，多谢了

@guanghuilineu 请问最后解决这个问题吗？你这个液滴是三维的吧？我最近也在测试VOF计算表面张力的效果，误差也不小，但没你这么大

在方柱计算中，对升力系数进行FFT分析，取样数据点在10000左右，得到如下结果，主频与实验较为吻合，在高频处出现谐振，是由于方柱前缘分离剪切层中K-H不稳定性导致的升力系数高频振荡吗？谐振又是怎么产生的呢？
诚心向各位前辈请教，求指导！
FFT Power.jpg

大家好，我正在用fluent做一个表面反应的模拟，因为还在测试阶段，就使用了一个二维的平板来做，选用的体系为乙酸甲酯水解反应，选用了完全多组分扩散，因为液体的扩散系数一般在1E-9的量级，所以就设置了1E-9的扩散系数，发现输出的反应速率和我手算的相差很大，当我将扩散系数设小比如1E-5时，发现输出的反应速率和手算的数据就吻合的比较好，请问大家知道我这种情况是为什么吗？
很简单的二维平板。1.png
这是我的反应设置。2.png
这是加入了正反应和逆反应的情况，误差在30%double.png
这是只有一个正反应的情况，结果直接相差两个数量级。 dan.png

非常感谢您的回复与帮助，我尝试一下。

@楼道男孩 SA用来射流中不多，一般用 k-w之类的 或者v2f

空气龄这个概念在通风领域应用比较广泛，我对此不大了解，但觉得这个物理概念上和停留时间相似。
因此扩散系数本质上于空气或者示踪气体扩散系数相似。
相关理论可以参考一下 https://age-of-air.forumotion.com/

openfoam 里面的空气龄有个类似的solver， 可以参考一下 （19楼）
https://www.cfd-online.com/Forums/openfoam-post-processing/68300-mean-age-air.html

@vllst 你好请问你的问题解决了吗？

@zhijie 感谢您提供的相关宏和帮助，我再继续寻找一下相关资料

尝试一下SPH或者MPS方法？

各位老师w71h11131532_1601213641_696.png ，我想模拟一下油滴在超临界水里下落过程中的对流和扩散问题：
初步打算用VOF模型来捕捉相界面，但有一个问题VOF并不能解决扩散的问题。请问有没有什么比较好的解决方案？
如果用UDF/UDS来解决扩散项的问题，请问能不能分享下类似的代码？

在利用per phase求解湍流方程时，可以用以下语句获取各相指针。
Domain *dp = Get_Domain(2); /primary phase domain pointer/
Domain *ds = Get_Domain(3) ;/secondary phase domain pointer/
Thread *tp = Lookup_Thread(dp, 7);
Thread *ts = Lookup_Thread(ds, 7);
2,3分别是各相domain_id，7为区域ID，可以从用户图形界面获得

@纵浪大化 好的，常常交流联系

@西湖冷月
nice兄弟，d=====(￣▽￣*)b

kOmegaSST和kEpsilon印象中是比较好收敛的两个模型，其他的用起来似乎都不太顺手，当然最重要的是结果正确。

@ibelief 在 求教fluent后处理-导入曲面（线）-输出曲面（线）参数的实现方法 中说：

或者是分割或者什么，探索一下。 祝

好的 谢谢老师

@bestucan 开的是双精度的，之前用pbm都是自己写UDF的，这次用的自带的，第一次遇到这个问题。。。回头有空搞搞，这感觉是软件自身的问题，我装个低版本的c试试，现在装的是 2012，改明儿装个2010，我还是感觉2010比较稳定。

各位老师，最近在研究双欧拉模型模拟液固两相流动，遇到以下两个疑惑，希望老师给予帮助:xiexie:
①想利用UDF对颗粒相剪切粘性系数和固相压力等参数进行修正（如26fbd023-2fb6-48f6-842a-733347fa08d1-图片.png ），公式中包含颗粒拟温度项，在编写UDF的过程中不知道如何获取该参数？？:134:
是否可以利用颗粒脉动速度C表示的颗粒拟温度7bf0c5a9-9020-484b-ba1e-2ba848de0117-图片.png 该方法中如何获取颗粒脉动速度？？？曾尝试用速度获取宏 C_U(c,t); C_V(c,t); C_W(c,t)得到颗粒速度，所得结果明显错误
②湍流模型中的湍流黏度是否要定义:141: 378f2ea3-f977-4a04-904a-05db21608467-图片.png
k-e湍流模型中湍流黏度定义为e94ab208-99c9-49db-9e65-664f5af4574b-图片.png
在输出湍流黏度时，是默认按照该公式进行计算的吗？？？只有有别于改公式时才需要写UDF定义吗:141:

@veen 谢谢，我学习一下

这种情况宁可相信Tecplot中Calculate variables和自己编程计算的结果，Fluent导出的速度场是没有问题的。

@关伯兰在笑 逐行排查吧，所有分母都添加一个很小的量。

这个提醒就是有负质量了，或许是负网格造成的，check mesh 一下

@东岳 因为时间成本问题吧，现阶段还没有这方面的打算，之后应该会学。。。其实只是想问下这个求解器，如果简单方便入手的话想要了解一下。。。

@关伯兰在笑 我也在做PBM模型的相关模拟，希望可以私下交流，QQ：1327010187

@百善孝为先0626 可以的话还是想麻烦您把我的联系方式告知他，这边看能不能合作一下

1.个人思路是在DPM中把雾滴设置成一个injection，粉尘颗粒设置成另一个injection，然后运算，想看到这个喷雾开启前后粉尘浓度的变化，但是冥冥之中感觉没有这么简单，希望各位老师能指点一下如何实现这个过程?
2.DEM可以模拟两种不同颗粒的碰撞过程吗？比如水滴和粉尘碰撞。

@东岳 非常感谢东岳老师的答复。 我想要求的是阻力，但是目前来看输出的Fx和输出的drag大小完全相等，而且瞬时值有正有负，这个平板一直往x正方向移动，按理其受到的阻力应该是一直朝x负方向的一个力。我目前的理解是所得的Fx和drag均是平板运动过程所受的总的合力，但是如何分离其受到的阻力和摆动产生的推力，分别得到这两个力是我目前想要解决的问题。

@楼道男孩 谢谢东岳大神的回复，我已经找到了问题的源头，文献中对应的是低雷诺数下ke模型，并不是标准的ke模型，是我之前知识盲区，导致一直没有找到解决的办法，现在正在尝试进行模拟验证，谢谢大家的关注！

只有在可以设置 flux 的边界条件上，才可以指定梯度边界。
比如，温度的边界条件可以指定 heat flux，这时你就可以用 DEFINE_PROFILEDEFINE_PROFILE 定义一个值，挂载在 flux 的边界条件上。
UDS 也可以指定 flux，同理可以用 UDF 指定梯度边界条件。

在帮助文档找到了模型。:xinlei:

用 DEFINE_SOURCE 宏，根据需求挂载在 k 或者 epsilon 的方程后面，试试

我想算一个高粘原油的沿程压降，参数：速度3m/s、密度940kg/m3、粘度0.2p.s、管径0.2731m、雷诺数3850.71。我不管用什么湍流模型，什么网格精度算出来的结果都和达西公式算的结果差20%多.但是同样的工况只是介质换成了水，结果就只差5%了。我去fluent的帮助文件找了关于沿程摩阻系数但是没找到，所以我想问问有没有知道fluent里面是怎么计算压降的，沿程摩阻系数是怎么取的?有没有算过类似东西的大佬，用的什么湍流模型、什么计算方法，可不可以分享一下设置细节。谢谢

直接的宏应该是没有的 你可以自己用UDF写个类似于level set的函数计算每个位置到界面的距离