FLUENT不收敛案例+边界层分离时湍流模型的选用

今天分享一个FLUENT计算不收敛的案例及其解决办法。计算的对象是一个文丘里管,如图1所示,入口直径为28mm,出口直径为50mm。计算所用的网格如图2所示,采用四面体非结构网格,但壁面附近采用棱柱网格以更好地分辨边界层的流动注1

640.webp.jpg

                                                                                    (a)三维图

640.webp (1).jpg

                                                                                      (b)截面图

                                                                                图1  计算的文丘里管。

640.webp (2).jpg

                                                                                       图2  网格

流动的工质为空气,由于流动马赫数很低,采用不可压缩流体模型,密度=1.225kg/m3,粘性系数=1.7894×10-5Pa·s。入口边界条件为:总压=1000Pa,湍流强度=3%,水力直径=28mm。出口边界条件为:静压=0Pa。注意上面所说的压力都是表压力而不是绝对压力(绝对压力=表压力+参考压力。在FLUENT里面参考压力通过设置“operating pressure”来改变。)湍流模型采用k-ω SST。

用基于压力的求解器,压力与速度的耦合方法为“Coupled”。动量方程的空间离散采用二阶迎风格式,湍流方程采用一阶迎风格式(图3)。

 

640.webp (3).jpg

                                                                            图3  计算方法及离散格式。

采用定常(Steady)算法计算发现,随着迭代次数的增加,文丘里管入口、出口的流量总是在波动,不能收敛到确定的数值(图4)。残差曲线如图5所示,可以看出绝大部分方程的残差都远没有达到默认的收敛标准(1e-3)。修改求解过程控制参数,例如松弛因子(Relaxation Factors)、Courant Number等等,都没有明显效果。

640.webp (4).jpg

                                                                                        图4  流量随着迭代次数增加的波动。

640.webp (5).jpg

                                                                                                 图5  残差曲线。

 是什么原因导致不收敛呢?通过画出流线图(图6),可以发现文丘里管的扩张段发生了边界层分离,即流动不再依附于壁面。根据以往的经验(读者可以阅读公众号以前的文章:为何我这个流动总是算不收敛?我要砸电脑!),流场中有大规模边界层分离的时候,定常算法往往都不能收敛,只有改成非定常算法才行。所以,我们尝试改用非定常算法。从理论上来说[2],有边界层分离的时候,多数情况下分离再附区以及钝体尾迹的涡脱落区属于非平衡湍流区,即湍动能瞬时生成与湍动能瞬时耗散相差较大,因此应该使用非定常算法来计算。(相反,不分离的顺压梯度边界层属于接近平衡湍流区,即湍动能生成近似等于湍动能耗散,用定常算法就足够了)

640.webp (6).jpg

                                                                                    图6  流线图。

改成非定常算法之后,发现只要设定合适的时间步长,每一个时间步内都可以迭代收敛了。图7显示的是时间步长设为0.0001的情形,每个时间步只需迭代3~4次就能让残差下降到默认的收敛标准(1e-3)以下。问题得以解决。

640.webp (7).jpg

                                                                                     图7  用非定常算法计算时命令窗口的输出。

虽然不收敛的问题得以解决,但是对于这个题还值得多说两句。对于这种型面是光滑曲线的壁面上的逆压梯度导致的分离流动,在CFD中要准确地模拟是不容易的,计算结果对湍流模型的依赖程度很大。例如,如果我们将湍流模型换成k-ε,就会发现算出的结果中并没有明显的边界层分离。这可以从图8(a)的速度x分量的云图看出来。这里我们只显示速度的x分量小于零的区域;因为边界层分离总是伴随着回流区的,而x分量小于零意味着回流。其它几种湍流模型的计算结果也在图8中进行了对比。可以看出湍流模型对这个问题的计算结果影响很大。

640.webp (8).jpg

                                                                                          a)标准k-ε

 

640.webp (9).jpg

                                                                                              b)RNG k-ε

640.webp (10).jpg

                                                                                            c)Realizable k-ε

640.webp (11).jpg

                                                                                               d)k-ω SST

640.webp (12).jpg

                                                                                              e)v2f

640.webp (13).jpg

                                                                                               f)SA

                                             图8  几种不同的湍流模型算出的结果的比较。只显示逆向流动的区域。

那究竟哪个结果最符合实际呢?这只能通过实验验证来决定。不过,有一点还是比较明确的,那就是对于这类型面是光滑曲线的壁面上的逆压梯度导致的分离流动,用k-ε模型模拟是不合适的;k-ε模型往往会显著低估边界层分离的程度。这一点在FLUENT的User’s Guide中已经明确地指出:

“The draw-back of some k-ε models is their insensitivity to adverse pressure gradients and boundary layer separation. They typically predict a delayed and reduced separation relative to observations. This can result in overly optimistic design evaluations for flows that separate from smooth surfaces (for example, aerodynamic bodies, diffusers).”

为了说明k-ε模型的这个缺点,我们再展示一个例子。这是来自于斯坦福大学Gianluca Iaccarino教授的湍流模拟讲义[1]中的例子,计算的是一个扩压器中的流动。可以看出,通过与实验数据的对比发现k-ε模型大大地低估了流动分离

640.webp (14).jpg
640.webp (15).jpg
640.webp (16).jpg

 

640.webp (17).jpg
640.webp (18).jpg
640.webp (19).jpg

如果壁面形状是不光滑的,比如有突然向外转折,那么反而简单一些,因为边界层分离通常会在转折处发生。

西华大学的学生罗亮提供了文中的算例;The University of Melbourne的研究生刘丽媛阅读了本文的初稿并提出了很好的修改建议;另外北航航空科学与工程学院的研究生靳一超和李健与作者进行了有益的讨论;在此一并表示感谢。

参考文献

[1] https://web.stanford.edu/class/me469b/handouts/turbulence.pdf

[2] 张兆顺, 崔桂香, 许春晓. 湍流大涡数值模拟的理论与应用. 北京: 清华大学出版社, 2008

流体力学及仿真

FLUENT不收敛案例+边界层分离时湍流模型的选用的评论2条

FLUENT不收敛案例+边界层分离时湍流模型的选用的相关案例教程

1.Q:Fluent Meshing划分体网格时出现重叠节点错误 问题描述:面网格调整好后,划分体网格出现这种错误,如图所示: A:尝试合并这些节点。菜单栏boundary->merge boundary或者输入命令/boundary/merge-duplicates 。 2. Q:Fluent Meshing如何提高poly网格质量? 问题描述:在画网格的过程中,遇到这样的问题,skewness
Fluent UDF 许多学习了一段时间Fluent的朋友们会开始接触到一个叫UDF的东东,这个东东的全称叫User Define Function,翻译成中文可以是“用户自定义功能”或“用户自定义函数”。 许多想入门UDF的朋友在学习UDF前面临的第一个问题应该是UDF环境变量设置的问题,这个问题可能已经成为新手入门UDF的一个门槛了。UDF环境变量设置现在网上百度也能搜到一些教程,这里我不操作
自然和工程中多数流动现象都是多相的混合流动。物理上物质的相分为气相、液相和固相。但在多相流系统中相的概念意义更广泛。在多相流中,相被定义为一种对浸没其中的流体及势场有特定的惯性响应及相互作用的可分辨的物质。例如,同一种物质的不同尺寸颗粒都可以被看作不同的相,因为相同尺寸的颗粒集合对流场具有相似的动力学响应。 多相流具有多种存在方式,以两相体系为例,可分为:气液多相体系;气固多相体系;液固多相体系;
来源:安世亚太 自从ANSYS 2020R1版本开始,Fluent软件就支持中文界面。很多人问怎么启动这个中文界面。 不需要卸载重装,因为即使在重装时候选择中文语言,也只是安装程序界面是中文,但是安装后,Fluent软件默认打开依然是英文界面。 首先,按正常操作,在电脑菜单栏中启动Fluent软件界面。此处可以是单独启动Fluent模块,或者在Workbench平台启动Fluent模块都可以,如下
本教程的目的是准确地模拟在300千瓦BERL燃烧室的燃烧过程。这类问题可以通过物质输运模型或非预混燃烧模型来模拟。在本教程中,将使用非预混燃烧模型来建立和解决天然气燃烧问题。 1 启动FLUENT并导入网格 (1)在Windows系统下执行“开始”→“所有程序”→ANSYS 19.1→Fluid Dynamics→FLUENT 19.1命令,启动FLUENT 19.1。 (2)在FLUENT La
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    2 0