求解Fluent中换热管平移类周期性边界条件和ICEM中的网格质量的困惑

浏览:165578 回答:15

参数:普通圆管,直径20mm,长度12mm,壁厚1.5mm,管内1073K的热空气流过,流速10m/s,恒壁温边界973K,管材质steel,按照流体定性温度1073K,修改热空气密度0.329kg/m3,定压比热容1.114kJ/kg K,导热系数λ=0.07176W/m•K,粘度μ=4.43*10^(-5)Pa•s,Pr=0.713.

我的问题:

使用fluent求出一根换热管壁面上Nu数,可是使用coutours of > wall fluxes >surface Nusselt Number看壁面wallNu数的分布,发现fluent算出的Nu数太大,和实际不符。如图1所示,通过Gnielinski公式计算的普通圆管的Nu数约38,而fluent计算出的Nu数小的有17.6,大的足有700+。图中Nu为负,不知道是不是表示流体被冷却?但是Nu数与计算值相差也太大了,如有大侠知道还请赐教,谢谢!

图片1.jpg

看了很多帖子,一种说法是模型管太短,流体的换热尚处于入口段,远远没到充分发展段,故严重偏离 Gnielinski 公式算出的结果。

于是,我考虑使用fluent中的周期性边界条件。先于ICEM中编辑Inlet进口和outlet出口之间的周期性(periodic)关联。生成网格文件后导入fluent分析。此时Type已经是periodic周期性边界,但是原本inletoutlet可能由于周期性设置已经合并为一个inlet,不再显示outlet了。在进出口边界条件中的periodic conditions进行设置,不知为何它只给我质量进口和压力进口2种边界条件选择,我选了质量进口。估算了下10m/s20mm管径的空气流过质量,给了0.001kg/s,改了flow direction方向和Number of Iterations的次数10次,这样是不是可以计算到原本10倍的距离上呢?

我将设置好周期性边界条件的MSH文件导入fluent后,一开始Check后没问题,没看到有负体积的网格,但在迭代第一步就出现报错,如下图。哪位大侠知道是怎么回事?

图片2.jpg

我后来怀疑是ICEM中网格质量画的太差所致,但苦于不知如何在保证周期性边界的情况下改善网格质量。下图是我的网格质量,经过global smooth之后还是提不上去,导入Fluent后迭代第一步就出错。苦于不知如何改善ICEM中的网格质量?基于周期性边界的ICEM画完网格再导入Fluent分析我还是不会,不知大家有没有这方面的例子?


图片4.jpg


图片3.jpg
邀请回答 我来回答

全部回答

(15)
默认 最新
胡坤
看你图中的网格,根本不能称之为边界层。传热问题对于边界层比较敏感。
另外你是否考虑了固体导热?如果不涉及导热的话,固体的厚度是不需要计算的。如果涉及导热,固体部分你是如何考虑的?耦合壁面?有厚度的壁面?还是其它?
2013年3月7日
评论 点赞 1
北京蓝威技术
楼主求解努塞尔数应参阅此篇博文:http://liusuanyatong.blog.163.com/blog/static/215595157201602002452913/
2016年8月4日
评论 点赞
张燕天津大学
楼主好,我也是将进出口设为周期性边界,设置好了periodic conditions,但是模拟提示没有进出来,无法计算,请问楼主怎么设置的,来壁面没有进出口的问题?
2016年1月13日
评论 点赞
xiemeiying1116
斑竹你好,我想问下你的周期性边界条件怎么设置的,我想模拟一个无限长的平板。。非常感谢你
2014年6月26日
评论 点赞
Dlasus

周期性边界在不少论文和书籍上说的很模糊呢,但是感觉周期性设置一个很重要的条件是,管子进出口的网格节点分布必须是完全相同的,如此才能一段接着一段计算迭代下去,不知道我理解的是不是对?此外,HyperMesh在这方面表现也很不错呢。
2014年5月19日
评论 点赞
Halo-Lin
所以说楼主最后没有用周期性边界条件了?我用周期性边界条件遇到的问题和你一样,弄了半天还是不知道怎么用。我查论文看到别人有用周期性模拟的,不过都写的很粗略,没什么参考价值。
2014年5月19日
评论 点赞
uuusheng

你好。我现在也是跟你算差不多一样的例子,只是我的wall边界不是恒温,二是给定heat flux ,但是计算出来的流场是充分发展的,但是温度场跟理论差距太大了。。根本不对~~ 不知道如何处理~
2013年7月5日
评论 点赞
Krs
非常感谢smart222的细心讨论,指出了我所说的普通圆管Nu=38.27的错误。因为我用错了公式,如下图。
未命名.jpg

这个Gnielinski公式据说仅适用于流体为水,而我用的是700-800摄氏度的热空气,所以计算结果错误。有关热空气在800摄氏度下通过圆管的Nu如何计算,我再去找找看。
2013年3月7日
评论 点赞
Krs

非常感谢你的提示,我自己又画了一根1.5M长的普通圆管(看了一些文献,据说入口段效应在l/d=60左右是个风水岭),管径20mm,厚度1.5mm,ICEM中已保证网格质量大于0.3,现在正在fluent里迭代,但发现每一步速度较慢,而且连续性那一项收敛得很慢,其他项收敛较快。假如计算结果符合我算出的圆管Nu理论值38.27左右,那应该就说明是我一开始没考虑到出口端效应吧?
2013年3月7日
评论 点赞
Krs
哦,对,这里输错了,是1.5mm厚度
2013年3月7日
评论 点赞

没解决?试试专家一对一服务

换一批
    App下载
    技术邻APP
    工程师必备
    • 项目客服
    • 培训客服
    • 平台客服

    TOP