倒阶梯流动数值模拟.docVIP

倒阶梯流动的数值模拟 一、课题介绍 问题的定义 上图所示为一个倒阶梯流道，流动的介质为常温的水，入口速度为20m/s，出口为常压，流道尺寸如图中所表示。 边界面定义 流体进口命名为inlet，出口命名为outlet，上壁面为upwall，下壁面从左到右分成三段，命名为downwall-1、doanwall-2和downwall-3。 二、倒阶梯流道网格划分 采用Gambit软件对倒阶梯流道进行网格划分。在划分过程中分别采用： 1．均匀的结构化网格； 在Gambit中划分节点时，所有线节点的Ratio选取1，得到如图1所示的均匀的结构化网格，局部放大图见图2。 图1 均匀的结构化网格 图2 均匀结构化网格的阶梯附近放大图 2．加密的结构化网格 根据要求对流体流动的边界层区域进行网格加密，再通过调整所有线节点的Ratio，对阶梯附近进行局部加密，示意图如下： 图3 加密后的结构化网格 图4 加密后结构化网格的阶梯附近放大图 划分完网格后，设上壁面upwall以及下壁面downwall-1、doanwall-2和downwall-3为壁面边界条件（wall）。入口边界Inlet为速度入口边界条件（velocity-inlet）。另外，考虑到管道较长，出口处流体已经充分发展，无回流，因此设置出口边界Outlet为质量出口边界条件（Outflow）。 三 Fluent 数值模拟 描述流动、传热问题的微分方程常常是一组非常复杂的非线性偏微分方程。除对某些简单问题外，很难获得这些偏微分方程的精确解。对于大多数实际工程流动、传热问题，必须采用试验研究或SIMPLE算法包括各种改进方案是不可压缩流体的Navier-Stokes方程数值求解中应用最为广泛的算法，并且已经成功地应用于可压缩流体的数值计算中。随着数值计算技术的发展，对SIMPLE算法的改进又产生了多种方法，如SIMPLER,SIMPLEC,SIMPLEX,SIMPLET等。发展新算法的目的主要是进一步提高收敛特性以及收敛速度。目前，SIMPLE算法最为成熟也具有较好的精度，所以本文是基于SIMPLE算法计算的。标准化k-ε模型k-ε 模型。 在直角坐标系下，的控制方程通用形式如下式： （3-1） 其中：Φ分别代表速度u ，v和w, k-湍动能，ε-湍流耗散率，SΦ-气相引起的源项，SpΦ-固体颗粒引起的源项。式中各项具体见下表-1： 数学模型控制方程 其中，=从临近单元中心到壁面的距离；=自由流速度；=流体的动力学粘性系数；=从边界层起始点开始沿壁面的距离。上面的方程基于零攻角层流流动的Blasius解。 当然，近壁面处的网格间距对于湍流也十分重要。理论上讲，模型，RSM模型以及LES模型都仅适用于湍流核心区域（一般都远离壁面）。如果近壁面的网格划分足够好，Spalart-Allmaras模型和模型也可以用来解决边界层的流动。 由于平均流和湍流的剧烈的交互作用，使得湍流流动的数值结果比层流流动的数值结果对网格的依赖更加敏感。关于这一点，从本文的算例中也可以看出。因此，成功的湍流流动计算在生成网格时必须要做一些考虑，尤其是当所要求解的区域平均流动变化很快并且存在有大比率应变的剪切层时，建议使用足够好的网格。具体而言，可以采用对壁面附近进行加密等方法。 85cm m 15cm m u＝20m/s 3cm 5cm