13 SIMPLE 算法的发展.ppt

§13 SIMPLE 算法的发展 针对SIMPLE 的算法的三大假设问题： 1、假定一个初始的速度场及与其相互独立的压力场，两者无任何关系。 2、在导出速度修正值计算式时没有计及邻点速度修正值的影响。 3、在一个层次的计算中，动量离散方程中的各个系数及源项均为定值。 发展的算法：其中针对1、2二种假设的算法是SIMPLER，对第2种假设的算法是SIMPLEC和SIMPLEX算法，针对后两种假设的算法是SIMPLE的date修正。 13.1 SIMPLER算法（1980） 1、SIMPLE算法中的问题： 1）u0，v0与p*相互独立不协调，从而影响迭代速度； 2）虽然p’采用了亚松弛处理,但改进的速度场和压力场仍不能很好的同步，从而影响了整个流场的迭代速度。 2、SIMPLER算法采取的措施: 1）在假设速度场后，计算出与其相适应的压力场，而不是另假设压力场； 2）p’只用于修正速度,而另寻压力场的改进方法。 5) 根据u*、v*，求解压力修正值p’； 6) 利用p’修正速度,但不修正压力； 7) 利用改进后的速度u，v,计算动量方程系数。重复(2)-(7),直到收敛为止。 13.2 SIMPLEC算法 1、引入原因：忽略周边速度对速度修正值的影响，犯了一个不协调的错误。 2、 SIMPLEC算法采取的措施：为了能略去周边速度的影响，同时又能使方程基本协调，试在u’e方程的等号两端同时减去∑ anbu’e，设u’nb与u’e具有相同的数量级,故忽略此项造成的影响要比忽略∑anbunb’项造成的影响小得多。 13.3 同位网格上的SIMPLE算法 1、 交叉网格特点： 1）优点：成功地解决了速度与压力的藕和关系问题； 2） 缺点：随着数值计算的问题由二维发展到三维，由规则发展到不规则区域，由单重发展到多重网格，交叉网格的程序编制的复杂与不便日益突出。 2、同位网格特点：把各个变量均置于同一套网格上而且能保证压力与速度不失藕的方法。 思考题 1、不同算法与SIMPLE比较有什么进展？ 2、SIMPLER算法中第三步与第五步有什么关系？第六步如何进一步进行？ 3、每一层次中系数及源项不变如何解释？ 4、SIMPLE算法假设是否影响收敛解？ 5、同位网格与交叉网格有何不同？如何解决1-压差？ * 3、初始压力场的改进: (如何由已知速度场计算压力场) 动量离散方程写成： 为假拟速度（在已知速度分布的时候可以算出上式中的第一项，则该项 有速度的量纲，称为假拟速度） 则： 同理： 代入连续方程的离散格式得： (1) 注意：与p’公式中b不同，其余相同。 其中： 4、具体步骤: 1) 估计一个速度场u0，v0,计算动量方程系数； 2) 根据u0，v0计算假拟速度； 3) 求解压力方程(1),得到压力场p*； 4) 根据解出的压力场p* , （与 u0，v0是协调的）求解动量方程，得u*、v*； 5、讨论： 1）在SIMPLER算法中，压力修正值p’的代数方程与压力改进 方程及边界条件与SIMPLE算法的类似性，均给程序编制带来方便； 2）初始压力场与速度场相协调，压力场不必亚松弛； 3）迭代层次数可以减少，但每一层次计算中SIMPLER所花的时间要多些，但总的计算时间要少。 上式可写成： （2） 同理： 3、讨论： 1）这是协调一致的SIMPLE算法，简称SIMPLEC算法。 2）SIMPLEC与SIMPLE的计算步骤基本相同，只有两点微小的 区别： ①以Ae/( ae -∑ anb)代替SIMPLE中的Ae/ ae； ②在SIMPLEC算法中，p’不再亚松弛，取ap=1。 *