10流场的计算.pptVIP

流场的计算 第十讲 流场的计算 流场的计算 流场的计算 流场的计算 流场的计算 流场的计算 流场的计算 流场的计算 流场的计算 SIMPLE算法 SIMPLE算法 SIMPLE算法 SIMPLE算法 SIMPLE算法 SIMPLE算法 SIMPLE算法 SIMPLE算法 SIMPLE算法 SIMPLE算法 SIMPLER算法 SIMPLER算法 SIMPLER算法 SIMPLER算法 SIMPLER算法 SIMPLEC算法 SIMPLEC算法 SIMPLEC算法 PISO算法 PISO算法 1986年Issa提出了PISO（Pressure Implicit with Splitting of Operators）算法，即压力的隐式算子分割算法，它源于非稳态可压缩流体的无迭代计算所建立的一种压力速度计算程序，后来在稳态流动中也较广采用。 SIMPLE，SIMPLER，SIMPLEC算法为一步预测，一步修正；PISO算法则是一步预测，二步修正。PISO算法的预测步与SIMPLE算法相同；第一步修正也与SIMPLE法相同，采用压力修正方程，在完成第一步修正后，再寻求第二步的修正，以便更好地同时满足运动方程和连续性方程，并加快每单个迭代步的收敛速度。 由于PISO算法经过2次的压力的修正，需要单独对二次压力修正方程的源项设立储存空间；同时，在每一次迭代中，PISO算法涉及较多的计算，相对复杂。尽管如此，也正是通过了2次压力修正，迭代过程更易收敛，计算速度更快。特别对于非稳态问题，PISO算法有明显的优势。相对地，在稳态问题中，SIMPLER与SIMPLEC则更合适。 * * * * 流体仿真与应用 ◆分析前面反映流场运动规律的控制方程，将会发现如下问题 ①运动方程中的对流项包含非线性量 对于第一个问题，解决的办法是迭代法。迭代法是处理非线性问题经常采用的方法，它是从一个估计的速度场开始，通过迭代逐步逼近速度的收敛值。对于第二个问题，如果压力已知，求解速度不会特别困难，只需用第5章介绍的方法，导出运动方程所对应的速度分量的离散方程，求解速度。 ②每个速度分量既出现在运动方程中，又出现在连续方程中，方程错综复杂地耦合在一起。更为复杂的是压力项的处理，它现在运动方程中，但却没有可用以直接求解压力的方程。 压力也是待求的未知量，在求解速度场之前，压力场是未知的 ◆解决因压力所带来的流场求解难题的方法 ▼非原始变量法和原始变量法。 ●非原始变量法是从控制方程中消去压力的方法。 非原始变量法存在明显的问题，如有些壁面上的边界条件很难给定，计算量及存储空间很大，因而，其应用不普遍。 ●原始变量法是直接以原始变量u，v，w，p作为因变量进行流场求解，该类方法也称基本变量法。 目前，广泛使用的是这类方法中的SIMPLE算法，以及在SIMPLE算法基础上改进的SIMPLER算法、SIMPLEC算法和PSIO算法等。 ◆交错网格 ▼运动方程中压力梯度离散所遇到的困难 运动方程对这样一个波形压力场的“感受”竟然与均匀的压力场的“感受”一样，因为相间压力值处处相等，显然，这是不能接受的结果。 ◆交错网格 ▼连续方程离散所遇到的困难 这就意味着离散化的连续方程将包含2个相间节点的速度差，而不是相邻节点的速度差；同样，锯齿波形的速度场完全不合乎实际的速度场，却满足离散化的连续性方程，对于二维和三维问题的数值计算，即使满足连续方程，也同样可能存在不合理的解。 压力和速度出现的问题主要来源于压力或速度的一阶导数项；相反，二阶导数则一般不出现此问题。解决这一离散困难的方法是采用交错网格。交错网格（Staggered Grid）又称为移动网格（Displaced Grid），是F.H.Harlow等人在提出著名的MAC法时首先使用的。 ◆交错网格 交错网格是将标量型变量（如压强、温度、浓度）的网格与矢量型变量—速度的网格系统错开。 ◆交错网格 标量型变量的控制容积称为主控制容积，相应的网格节点称为主节点。 ◆交错网格 由于交错网格，也避免了连续性方程所遇到的困难。另外，在主控制容积中，速度节点的位置正好是在标量输运计算时所需要的位置。因此，不需要任何插值就可得到主控制容积界面上的速度。 在计算过程中，所有存储于主节点的物性值在求解u，v，w方程时，必须通过插值才能得到所需位置上的值 。 其次，由于u，v，w，p及其它变量的网格系统不同，在求解离散方程时，往往需要一些相应的插值。 由于三套网格系统，节点编号必须仔细处理方可协调一致。 ◆运动方程的离散 ●对u，v，w方向上的运动方程积分所用的控制容积不是主控制容积，而是各自的控制容积 。 ◆运动方程的离散 ●运动方程中的压力梯度项从源项中分离出来。 ◆压力与速度的修正 对于离散的运