Fluent的计算结果与理论公式值不匹配

[李东岳]

网格实际上是因为想获得完全发展湍流，把通道弄很长所以才会有这么多。一个光滑通道的网格也不怕泄密啥的。

另外，你这个就是矩形通道么？这样的话RANS不会有涡出来。

[CarbonO2]

@李东岳 project1.blk project1.tin

emmmmmmm我不太了解没有涡指的是什么，这一步也只是为了下一步大涡模拟获得完全发展的速度入口，如果没有影响的话倒也无所谓。

另外，刚刚发现我模拟出的结果好像是层流（也可能是顶端速度变化很小的湍流？），Re=20000的流动应该是湍流，这也许就是问题所在？

[李东岳]

如果就是矩形通道，如果就是RANS，那一定是层流。都没必要用这么长的网格，直接画个一维网格高度上计算就行了。我觉得现在你对这个文章里面的一些概念，或者一些CFD的概念可能有点歧义。我建议你深层次理解一下，然后我再给你对比。目前你给我的这套网格，是完全没有必要的。

[CarbonO2]

@李东岳 研一，这几个月才开始接触CFD，目前属于赶鸭子上架，有些概念确实来不及深入了解 可能之前我话没说清楚有点歧义，我知道光滑通道的RANS模拟不会有涡，但是我在帖子里也没提到涡所以我不太清楚大佬你提涡是为了什么。RANS模拟不会有湍流的那种状态我也清楚，我主要想表达的是速度分布的状态是官方教程里左边图的样子，但是从雷诺数来看应该是右边的样子所以会不会是这里的问题，关于这方面的知识我确实需要深入了解一下。

一维网格的概念没接触过，最开始倒是有考虑过二维网格，不过因为涉及到雷诺数参考长度的选取以及profile文件导入三维网格可能出现的问题就转做三维了。模型搞这么长主要还是为了保险，毕竟官方教程里的湍流公式按水力直径算得4.6m才能完全发展。目前对我来说网格多倒不是问题，反正还得上课，课上完了差不多也算完了。

现在我最想深入了解的还是条件一样为啥算出来的东西不一样而且结果还是错的的问题

[李东岳]

条件一样为啥算出来的东西不一样而且结果还是错的的问题

如果条件一样，算出来的肯定是一样的。肯定是某些细节方面没弄好。fluent这面我不太熟悉所以帮不上你扣细节。

[CarbonO2]

@李东岳 唉，还能有哪些方面的细节呢，我又重新看了遍文献，湍流模型，边界条件都一样，改物性参数和求解方法也都试过了，就剩一个改模型参数了，早上试了一下也没啥效果。

把湍流模型换成了sstkw，这次倒是能有右边图的速度曲线了，但是压降还是只有理论值的5/6,范宁摩擦系数还是和理论值接近并且大于总摩擦系数。总不能△P不是拿静压算的吧 大佬你要是有时间的话就帮我算算这个模型用RANS模拟到完全发展之后到底能不能用压降算出和布拉休斯公式一样的数，没有时间的话就算了。总之感谢大佬这么多天以来和我的交流。

[李东岳]

用RANS模拟到完全发展之后

这种几何，可以用1维模型来模拟，这就是充分发展的。所以我就没想清楚为什么要用几百万网格算呢。你可以看一下这个算例，就是用RANS模拟的充分发展管道 http://dyfluid.com/openfoam.html#y

你也可以联系下文章的通讯作者

[CarbonO2]

@李东岳 额，做这个光滑通道也不是为了验证什么就是为了拿它出口参数当LES的入口参数，三维的直接导出出口profile文件再导入LES的case里就行了，不用三维的可能要搞UDF但是组里也没人会所以多一事不如少一事直接算三维的得了 。算例里好像没给压力的数据。

[李东岳]

光滑通道这个基本就类似个抛物线的形状。这个相对很光滑。你用来做LES，如果壁面没有障碍物的话，出不来湍流。另外LES那面进口也是个深坑，你看看这个贴 https://www.cfd-china.com/topic/4317/les定义入口速度的问题-dsrfg方法

[CarbonO2]

@李东岳 这点我清楚，实际上这么多参数里也就总摩擦系数对不上，其他的基本上全对上了。LES的话得在文献里是说在入口用vortex method施加扰动，而且模拟的是有肋的，所以基本上不担心涡的问题。流向周期性的模型也可以在RANS模拟之后用TUI命令施加扰动，之前试过了涡能长期存在，不过那个case没算完结果对不对就不清楚了。