fluent总能比of得到更满意的结果？

[悬浮的猪]

同样的网格，也近乎相同的边界条件，最后都达到了相应的参差需求，然而对比实验数据，of的结果惨不忍睹，连大致的趋势都没有正确，相反fluent却得到了更令人满意的结果。至于湍流模型是我自己添加的，尽量仿照文献公式编译，与fluent中自带的湍流模型保持一致，具体来说式低雷诺数k-e系列的，AB模型，研究的现象是冲击射流的流向努塞尔数变化，话不多说下面式对比图

[悬浮的猪]

[dzw05]

建议用排除法比较。（1）首先计算层流，对比fluent和OpenFOAM结果，这一步基本不会有问题，否则OpenFOAM不用混了。（2）再计算湍流，进行对比。个人觉得可能你的湍流模型可能有问题。

[cccrrryyy]

如果FLUENT和OPENFOAM都用同样的湍流模型呢？我比较过k-epsilon，区别不大

[悬浮的猪]

@cccrrryyy 你好，请问您验证的是什么样的算例呢？我现在反思了一下或许我这边的问题还要出在边界条件上，亦或是温度方程上，湍流模型基本上与fluent保持一致。

[悬浮的猪]

@dzw05 谢谢您的建议，我的湍流模型是拿来验证的，基本按照其出处文献的说明来写，我会再看看比较一下的，谢谢您的建议。

[cccrrryyy]

@悬浮的猪 我的是高温气体喷入室温的环境中，用理想气体状态方程，在RANS下做了验证。当时只是大致算了一下，我关注的是速度和温度，所以没比较努塞尔数。速度和温度分布，包括k和epsilon，FLUENT和FOAM差别是有的，但并不大，基本上可以忽略。光从结果图上看，FOAM算出来的努塞尔数快到100了，明显太大。你的工况很复杂吗，为什么需要自己弄湍流模型呢？RANS的东西现在每个软件应该都是非常成熟了，足以处理绝大部分问题。

[李东岳]

如果楼主可以分享一下openfoam算例，有空的话，并且网格数量比较小的话.. 我可以尝试调节一下，

[悬浮的猪]

@cccrrryyy 老师您好！我准备在原有的rans低雷诺数模型上写一个自己稍加改动的模型，加上一些自认为算是有所改善的尺度或者说是一些阻尼函数，那么我的工作意义也就是在于对这些新兴的尺度或者说阻尼函数在特定的模型上的适用性，那么第一步就是现在of上还原这个模型，并且能够与一些权威的软件上的湍流模型相验证，然后再在其上做出调整，做一些数值试验。我现在的这个算例算出来的k与epsilon的分布差异还是蛮大的，所以现在我认为这与我自己的某些细节设置存在关系，最后，谢谢老师抽出时间回答。

[悬浮的猪]

@东岳 好的谢谢李老师，我尽快将我的算例以及模型贴上来，网格的数量在十几万左右，比较小是一个二维的网格

[李东岳]

图片中对比的哪一条线的什么物理量？
有实验结果的文章么发来看看

[李东岳]

我发现你用3D的网格，缺把前后面设置为了empty边界条件进行2D计算？

[cccrrryyy]

@悬浮的猪 恩，我的意思是你可以先用两个软件都有的现有模型，确定你的设置中某些方面是正确的，这样可以缩小你检查问题的范围。

[悬浮的猪]

@东岳 李老师这样做是不正确的吗？of进行2d计算我的师兄师姐们是这样教我的，将不用的那个面empty。小白问题还请李老师指点！见笑了老师

[悬浮的猪]

@cccrrryyy 谢谢老师的建议，我准备用都有的Laundershmark做一下试验

[李东岳]

@悬浮的猪 做2D计算就把你的Z方向处理成一层网格并给empty边界，做3D计算需要把你Z方向厚度同时加密并测试LaunderSharmaKE，同时更改壁面边界条件

[悬浮的猪]

@东岳 李老师必须是严格的一层网格嘛？那我之前是有问题的，我做的是2d数值试验

[李东岳]

@悬浮的猪 那是严重错误的计算

[悬浮的猪]

李老师给出的建议一针见血，刚刚看到李老师的建议我就去进行修改了，我把改正后的对比图放上来，细节决定成败，李老师今天的给的建议让我感觉很多细节上我都没有到位还需要和老师们多多学习

[悬浮的猪]

李老师这是我修改后的结果

[悬浮的猪]

这也说明我往里面写的模型大体正确

[李东岳]

你这个挺好啊新结果

[李东岳]

目前有什么进展么

[linhan.ge]

这个标题容易吓到要入of的新人。

[悬浮的猪]

@东岳 抱歉东岳老师最近才看到回复您，我整理了一下数据内容，疫情期间写了一篇关于冲击射流传热的文章，主要是探讨不同湍流热扩散系数求解方案对预测结果的影响，最近正在修改文章

[悬浮的猪]

@东岳 文章正在给导师以及相关的老师看，目前我在研究相关磁流体领域类似的传热现象，期望能再研究生生涯的末期再写一篇文章，恐怕到时候又是一连串问题要请教东岳老师

[悬浮的猪]

@linhan-ge 哈哈，老师你好，这个只不过是我个人的一点误解，当不得真的

[李东岳]

过去很久了，这个是用什么求解器算的我有点忘了？可以提供下算例么？我放在dyfluid.com里的[验证算例]一栏，参考 http://dyfluid.com/interfoam.html#testcase 会标注你的贡献

不方便就当我没说哈 都可以理解

[悬浮的猪]

@东岳 李老师，怎么样提供算例呢，我非常乐意提供的，只是前段时间有事没有来得及上网，对此很抱歉

[李东岳]

@悬浮的猪 非常感谢，可以直接把算例发我邮箱 li.dy@dyfluid.com，然后我会把相关东西放在dyfluid，然后在文中标注你的贡献

[悬浮的猪]

@东岳 在 fluent总能比of得到更满意的结果？ 中说：

li.dy@dyfluid.com

李老师我已经将相关的算例，求解器，湍流模型都发过去了，您查收一下

[李东岳]

@悬浮的猪 谢谢，收到。问一下，这面必须使用你自己的求解器才能获得比较好的结果么？另外，下图数据的位置信息是否有？我想做个类似的图放在网站。因此需要自己后处理一下。

[悬浮的猪]

@东岳 李老师我这个求解器没什么特别，只是在simpleFoam中加了一个温度场，相关的文献我给老师发过来

[李东岳]

@悬浮的猪 感谢，我先看看，然后在本帖更新

[李东岳]

大佬，这面后处理努塞尔数特征长度取得是0.037么

[悬浮的猪]

@东岳 李老师我取的是雷诺数10400下的射流工况，w取0.01423为特征长度

[李东岳]

http://dyfluid.com/buoyantSimpleFoam.html#testcase

目前用buoyantSimpleFoam + LaunderSharmaKE算出来的结果就是这样了，中间那块跟实验值有差异

[悬浮的猪]

@东岳 李老师我后来也用LS模型，计算了一下，对于湍动能的设置可能李老师还需要改进一下，我计算的结果要比这个稍好。这个结果我用原来设置的参数计算得到过，但我觉得不太完美，就修改了湍动能和湍流耗散的设置得到了更高的结果

[李东岳]

@悬浮的猪 谢谢，我试试 感谢分享！

[jinlinna]

@东岳 李老师，你好，我看了您的算例，想请教您湍动能为0.54的初始值是如何设置的呢？我用湍动能经验公式算出来不是这个值呢？谢谢！

[jinlinna]

@悬浮的猪 如果可以的话，我可以看下您的求解器吗？在往simpleFoam里面添加温度输运方程时我遇到了问题，感激不尽！

[李东岳]

@jinlinna

老师我后来也用LS模型，计算了一下，对于湍动能的设置可能李老师还需要改进一下，我计算的结果要比这个稍好。这个结果我用原来设置的参数计算得到过，但我觉得不太完美，就修改了湍动能和湍流耗散的设置得到了更高的结果

你看他说湍流场对结果有影响，你试试调整一下边界条件

[jinlinna]

@东岳 谢谢李老师的回复！我确实修改了湍动能和湍流耗散率的初始值，得到了不同的结果，但我疑惑的是湍动能和湍流耗散率的初始值设置都是调试出来的吗？
目前我找到的大家都说用湍动能经验公式估算这个值，但用这个公式计算出来的往往与文献中的结果不符。那如果是自己做研究，而不是验证文献，又该如何给定这个值呢？

[悬浮的猪]

@李东岳 东岳老师，关于这个主题帖子，4月份会有一篇相关用OpenFoam完成计算的论文刊出，有兴趣的小伙伴可以参考着看一下，里面加入了一个直接求解热通量的双方程标量方程，也许传热相关的小伙伴会感兴趣，如果有小伙伴需要，我可以选择性提供模型代码，以及调试好的计算算例。谢谢东岳老师，因为马上毕业，不想自己弄出来的一点东西就没有了，还希望广大网友完善，谢谢！

[李东岳]

@悬浮的猪 恭喜恭喜 静候大作！

[悬浮的猪]

关于本专题讨论已发表了一篇文章，在知网可以查询到，小论文是基于低雷诺数κ-ε模型冲击射流传热的数值模拟，大论文的第五章有更详细的讨论，大论文题目是三种典型贴壁射流的数值模拟，过段时间应该也可以检索到。有错误之处还请大家指出，欢迎讨论谢谢大家，希望能帮助有需要的伙伴。

[悬浮的猪]

有需要相关代码的朋友可以私聊我，我会酌情提供

[李东岳]

@悬浮的猪 有pdf不，不在学校下载不下来

[悬浮的猪]

@李东岳 抱歉啊李老师，失踪了一段时间，主要是前段时间在谋未来，刚刚有了还算满意的着落，我现在可以把硕士写的三篇小论文都发上来，希望还来得及

1. 基于低雷诺数κ-ε模型冲击射流传热的数值模拟_朱泽辉.pdf

2. A Comparative Study of the Buoyancy-Opposed Wall Jet.pdf

3. A comparative study of heat transfer_68D53.pdf

[悬浮的猪]

哈哈哈哈好了，总算善始善终