269有进出口的容器内流体流动数值分析.docVIP

有进出口的容器内流体流动数值分析 济南市同圆建筑设计院有限公司 潘雷 摘要：本文对带有进出口的容器内的流体流动问题，建立二维稳态流场模型，并运用SIMPLE算法程序进行数值求解，对出口边界的处理采用总体质量守恒的方法，从而得出带有进出口的容器内的二维流速矢量图，并对不同流速、不同流体的流场进行了分析比较。 关键词 ：流体 SIMPLE算法 稳态流场 矢量 0．引言： 计算带有进出口的容器内的流场，是流体力学和传热学中常见的问题。工程实际中，像此类问题比较常见，例如房间内的送排风，带有下进口上出口的储水箱等。可将流场简化为一个形状规则的空间内的流体流动问题，如下图所示，流体以一定的速度流进一空间，尔后流出。空间长度为2，高度为1。求解此空间的流速分布。 1．数学模型 假定为常物性的不可压缩流体流过的空间，流动的守恒型控制方程及边界条件。 控制方程： 质量守恒方程 动量守恒方程 边界条件： 进口截面 固体壁面 出口边界 从数学的角度应给出的分布，实际上常难以实现。本文采用使出口法向流速分布满足总体质量守恒的方法。 2. SIMPLE数值解法： 采用SIMPLE数值解法求解本问题。在离散过程中采用内节点法确定所求解的物理量的几何位置，节点位于子区域的中心，这时子区域就是控制容积，划分子区域的曲线簇就是控制体的界面线，边界节点的控制容积为零。采用均匀网格及错列网格。采用控制容积法有限差分公式求解问题时如果将速度分量的位置与控制容积的节点位置定义在一处，则会出现不真实的速度场与压力场，为了解决这个问题，采用所有的速度分量都不再定义在结点处，而是定义在控制容积的交界面处，也就是错开了半个控制容积。其他的变量（包括压力变量）仍定义在控制容积的结点位置。因此，速度分量的位置和控制容积的结点错开了。按这个原则划分的网格称为错列网格。其好处在于可以直接应用速度分量计算控制容积交界面上的质量流量，得到合理的速度场与压力场。 最后利用SIMPLE算法程序求解本问题。SIMPLE算法原理：一般情况下，压力场往往未知，而假设的压力场又不可能很精确，因此按对流扩散方程求出的速度场必定不能精确满足连续性方程。这时，如果用连续方程来校正压力场，就可用校正后的压力场作为改进值，根据动量方程的离散化形式重复计算速度场，直到求得的解（u,　v,　w）同时满足动量方程和连续方程为止。其采用乘方格式与错列网格。 SIMPLE算法计算步骤： (1) 假定一个速度分布，记为u0,v0,以次计算动量离散方程中的系数及常数项； (2) 假定一个压力场； (3) 依次求解两个动量方程，得； (4) 求解压力修正值方程，得； (5) 根据改进速度值； (6) 利用改进后的速度场求解那些通过源项物性等与速度场耦合的变量，如果并不影响流场，则应在速度场收敛后再求解； (7) 利用改进后的速度场重新计算动量离散方程的系数，并用改进后的压力场作为下一层次迭代计算的初值，重复上述步骤，直到获得收敛的解。 物理参数的确定。本文中的流体为动力粘度系数的流体，密度为Simple主程序中的默认值。 松弛因子的选取。在实施SIMPLE算法的过程中，速度与压力的修正值都应作亚松弛处理。对速度，为限制相邻两层次之间的变化，以利于非线性问题迭代收敛，要求亚松弛，本文取速度亚松弛因子值为0.3，速度亚松弛因子值为0.7。对压力，由于在速度修正值公式中略去了邻点的影响，所解得的修正速度是合适的，但对压力修正值本身，则是被夸大了，因而要亚松弛，本文取压力亚松弛因子值为0.3。 3. 边界条件处理 入口水平速度，其中为给定值。固体壁面上水平速度、垂直速度都为零。因为给定的问题中只有速度条件，没有温度条件，所以在分析时不考虑温度的影响。处理出口边界是本文中最关键的，目前广泛采用的一种处理方法是假定出口截面上的节点对第一个内节点已无影响，因而可以令边界节点的影响系数为零。这样出口截面上的信息对内部节点的计算就不起作用，也就无需知道出口边界上的值了。这种处理的物理实质相当于假定出口截面上流动方向的坐标是局部单向的。为了在数值计算中应用这一简化处理方法而又不致引起过大误差，应做到：（1）在出口截面上无回流；（2）出口截面应离开感兴趣的计算区域比较远。 确定出口法向速度的方法有两种，按出口截面前一站的流速来确定出口截面的流速并要求其满足总体质量守恒的方法，这适用于出口截面采用局部单向化假设的情形。 第一种方法：假定出口截面上各点的发向速度的相对变化率为一常数，对图1所示的流动，设出口截面上不同的处都有以下关系成立： 由此得： 确定系数的条件示满足总体质量守恒： 其中为入口质量流量，为流动截面面积。于是有： 第二