CFD中文网

谢谢大佬，还有一个问题，就是下图右半部分该怎么关联，想要得出图2的效果，但是目前的关联方式关联出来全是负体积

各位大佬们，连续的圆弧对其中一个圆弧关联之后旁边的圆弧就没有边能关联了，该怎么解决

各位大佬好，模型是一个0.125x0.5x0.6m的光滑通道，上下壁面均为无滑移的300K恒温壁面，进口温度320K，进出口及左右壁面均采用周期性边界条件，模拟步骤是先用RANS模拟，之后使用TUI指令添加扰动然后切换到LES模拟。按理说流体温度应该小于等于320K而且随时间逐渐减小，但是云图画出来最高温度有323K，又模拟了0.5s后最高温度没什么变化其他地方的温度似乎也没什么变化。周期性设置如图（压力梯度最开始为0）
我想知道是为啥

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

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

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

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

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

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

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

@李东岳 project1.blk project1.tin

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

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

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

网格链接（两个网址是同一个文件）：
https://pan.baidu.com/s/1qTxoeUEgPGzIB2kSDf0YVA?pwd=1ruw
https://www.aliyundrive.com/s/XHWpd4c8FfU 提取码: 2n2l

网格是用icem画的，导出的solver里没找到openfoam所以用的是fluent

@李东岳 RANS

又发现了一个新的问题，表面摩擦系数算出来比总摩擦系数还大

@李东岳 补充下具体细节：
网格：

X方向流向，总长5m，网格长约6mm，均匀分布
Y方向展向，总长0.5m，网格长约3mm，均匀分布
Z方向高度方向，总长0.125m，底面和顶面第一层网格高度0.3mm，增长率1.1，网格长度最大值3mm，根据算出的结果完全发展后的y+值约为0.8
方法：

SIMPLEC，全二阶迎风格式

物性参数：

320K空气的参数，文献里没给物性常数，但是可以根据h+的公式算出的√ρ/μ的值，进而推断出文献里貌似用的是默认参数，改成默认参数之后压降还是不对而且应力更小了
收敛判定：

残差能量1e-9其余1e-5

@李东岳 感谢大佬的回复，能对上，热导率也差不多能对上，就是摩擦系数小了很多。另外根据文献图16的结果他的切应力是0.016pa左右，我计算出的结果只有0.012pa。

求助各位大佬，我正在模拟光滑通道的流动以复现一篇文献，但是计算结果比理论公式值低了约15%
文献中使用相同设置能得出正确结果，网格无关性已验证，增加到2000W网格也无法得出正确结果。我想知道问题出在哪里？以下是我的模拟细节：

那篇文献的网址：https://doi.org/10.1016/j.ijheatmasstransfer.2021.121573