在fluent中使用密度基隐式求解器进行孔板的DES模拟，出口后的马赫数最高可达1.5。由于流速很快，如果要保证CFL小于1，时间步长将会特别小，现有的计算条件无法满足。我困惑的点在于以下几点
（1）对于隐式求解，在控制面版里可以设置CFL，增大CFL在定常计算中可以加快收敛，但在非定常中起着什么作用呢，而且似乎无法对CFL数进行监测。我测试了一下在保证时间步长、时间步不变的情况下，改变CFL似乎并没有对总的模拟时间有影响。
（2）对于隐式求解，时间步长的选择一定需要CFL小于1吗？理论上不是说隐式是无条件稳定的吗可不可以取的大一点？最大迭代数/时间步怎样选择合适呢？
5595a209-0180-4ae3-a17d-3d1c76c838f8-image.png

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

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

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

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

@李东岳 直接除掉普朗特是什么意思？是直接令VT=D吗？那我的VT是怎么得到呢，是通过数值模拟数值模拟得到的吗？如果是的话，扩散系数还是一个变化的值呀。

simplefoam如果单一流体，不存在重力与浮力，这部分重力与压力互相制衡了。因此不会出现你想要的效果。

在3D正交笛卡尔内场网格中一个cell分别有几个face neighbours，edge neighbours，point neighbours呢？只共享一个vertex的是point neighbours，还是说只要共享顶点就是point neighbours呢？如图的说法有点违背了我以前的认知，希望大家给个正确的说法，谢谢:xinxin: :xinxin:
0678cdbe-c1b6-45b6-8464-6b7078b89e6e-image.png

自问自答一下
C-AMG和SA-AMG的区别在于粗化方案不同以及插值矩阵的构造是不一样的，C-AMG基于细网格生成粗网格，可以理解为更精细，因此计算收敛性会更好；SA-AMG根据聚类的思想来的，所以在setup阶段会大大缩减时间，但是收敛性不如C-AMG；而两者都比较鲁棒性，因此并不是只针对某种方法有效。
而后一个问题目前调研到的是，AMG在商软中均可以独自使用或者作为预条件使用；而Fluent、COMSOL、STAR-CCM+这三个软件的官网都有对应的介绍，基本上作为预条件都是与CG/BiCGSTAB联用，其他的还有GMRES。
找了好久也算是把资料找到了

可以这么理解

@vonboett

if I may implement two kinematic clouds in DPMFoam, one that accounts for particle collisions (larger stones but smaller than gridsize) and one that only accounts particle-wall collisions (many gravel particles). This way, I could model erosion..."

Dear Albrecht,

Thanks for your comments. The E-L algorithmn is expensive since it needs to track the lagrangian particles. So it depends on your particles' number. From what you've written, it looks you are on your right way doing this.

D

@匡之舒云 我在做射流燃烧，可以加个联系方式聊一下:xiezuoye:

感谢分享感谢分享

@lyc 我也是菜鸡，这玩意我师兄去年也想过，一直没成功，今年我做课题时候要添加自己颗粒属性方程可能接触一点。但是如何实现我也在研究中，你可以看看我之前和其他大佬的一些疑问，他们指出过在哪些文件里面进行属性添加修改，不过都不全，导致我目前现在想要修改内容也还差一点

如果我10月份有时间的话，或许会写个解析，更新到 http://dyfluid.com/cfd.html
就是我这个时间比较紧

@李东岳 在 关于LES模拟中采用大时间步长的问题 中说：

你可以放大库朗数，或者你可以用URANS更快。这些都没问题。真正分析的时候库朗数放小就行。并且库朗数一般是局部库朗数的最大值。如果1000个网格里面有1个网格特别小。会出现999个网格的库朗数都在1以下，但这1个网格在10左右。都不是问题。

好的，谢谢李老师！我明白了

@李东岳 感谢回复。我是想先写个简单的特定算例学习一下算法，格式边界条件啥的都直接写死，然后再慢慢做的通用一点，能搞出来肯定开源。

我今天看一篇论文，一个简单的方形腔，通过添加金属颗粒来改变腔内相变材料的热性质。文中一个研究点是改变金属颗粒的体积分数，按理说相同的纳米颗粒体积分数下混合流体的热性质都是一样的，根据Gr数计算公式，Gr数不应该变化呀，但是相同的颗粒体积分数下文中又设置了10000与100000两个量级的Gr数，不知道怎么弄的，求大佬解惑。
对了，我找遍那篇文章都没发现方形腔的物理尺寸，不知道这个与Gr的控制有没有关系

居然还发公众号了

三维的类似情况，我最后一个打错字了。

@李东岳 明白了，谢谢李老师。排版排了好久，清楚一点方便讨论:papa: