第四章速度场计算simple算法.pptVIP

一、SIMPLE算法的计算步骤 5.据p′改进速度场； 6.利用改进后的速度场求解与速度场耦合的变量； 7.利用改进后的速度场重新计算动量离散方程的系数，用改进后 的压力场作为下一层次迭代计算的初值，重复计算。直到获得收敛的解。 二、SIMPLE算法的讨论 1.速度修正中忽略了邻近速度修正值的影响，不影响最后收敛的解，但加重了修正值p′的负担，使得整个速度场的迭代速度减慢，故对p′作亚松弛， 其中 一般可取0.8左右。 类似速度修正值也可考虑亚松弛，一般松弛因子可取0.5左右。 二、SIMPLE算法的讨论 2.压力修正方程是椭圆形方程，即压力是向各个方向传播的，只有在可压缩流体的超音速流动中，压力的传递才会有单向的特性，这时应采用可压缩流体的p′方程。 二、SIMPLE算法的讨论 3.压力的参考点选取问题。一般只考虑压力的相对值，无须指定某一压力参考点，但若指定某一特点的参考点，其迭代收敛速度会变慢。 二、SIMPLE算法的讨论 二、SIMPLE算法的讨论 4. p′方程的求解方法。可采用Patankar提出的交替方向线迭代（ADI）加块修正的方法。 二、SIMPLE算法的讨论 5. p′方程迭代收敛准则的选取。 （1）简单的规定实施交替方向线迭代与块修正运算的轮数。 （2）规定p′方程余量的范数小于某一数值。 （3）规定终止迭代时的范数与初始范数之比小于允许值。 二、SIMPLE算法的讨论 6.终止整个问题的迭代准则。 （1）各节点上前后两次解偏差的绝对值或相对偏差的绝对值小于允许值。 （2）要求在内点上连续性方程余量的代数和及节点余量的最大绝对值小于一定的数值。 （3）要求连续性方程余量的范数小于一定的数值。 （4）要求在整个区域内动量方程余量之和与入口动能的比值小于一定的数值。 三、SIMPLER算法 1. SIMPLER算法的原理：SIMPLE算法得出的p′对速度的修正是相当好的，对压力的修正则过分了。虽然对p′采用了亚松弛处理，但未必恰到好处。由此产生了下列想法： p′只用来修正速度，压力场的改进则采用更合适的方法，此即Patankar提出的SIMPLER算法（SIMPLE Revised）。 三、SIMPLER算法 2.压力方程的推导 动量离散方程可写成： 其中前一项称为假拟速度，速度可记为： 三、SIMPLER算法 将上两式代入连续性方程的离散形式，可得到与压力修正方程形式相同的压力方程： 三、SIMPLER算法 3.SIMPLER算法的计算步骤 （1）假定一个速度场u0，v0，计算动量方程的系数。 （2）据已知的速度计算假拟速度 ， 。 （3）求解压力方程。 （4）把求解的压力作为p0，求解动量方程，得u*，v*。 三、SIMPLER算法 （5）据u*，v*求解修正压力值p′。 （6）利用p′修正速度，但不修正压力。 （7）利用改进后的速度，计算动量方程的系数，重复第2步到第7步的计算，直到收敛。 三、SIMPLER算法 4.SIMPLER算法的特点 在SIMPLER算法中，初始的压力场是与速度场是协调的，不必亚松弛，使SIMPLER方法的迭代层次数可以减少；但每一层次的计算中所花费的时间比SIMPLE算法要多。总的说来，SIMPLER算法所花费的时间比SIMPLE算法少。 四、SIMPLEST算法 由Spalding 提出，在PHOENICS软件中得到应用，其特点是： （1）对流项采用迎风格式； （2）把邻点的影响系数表示成对流分量及扩散分量之和，并把对流部分全部归入源项。 五、SIMPLEC算法 １.SIMPLEＣ算法的原理： 在SIMPLE算法中，为求解方便，略去了速度修正值中的　　　和　　　　　，从而犯了速度和压力不协调一致的错误。为此在速度修正方程两端同时减去　　　，可得到： 可略去前一项，于是有： 第一节 速度场求解的困难 一、方程的通用型式 二、速度场求解的困难 一、方程的通用型式 连续性方程： 动量方程（x方向）： 动量方程（y方向）： 能量方程： 一、方程的通用型式 写成统一型式 二、速度场求解的困难 二、速度场求解的困难 1.一阶导数 的存在 考虑一维常密度问题，有： 可离散为 二、速度场求解的困难 考虑连续性方程： 由于 故 所以 亦有 可得到不合理的锯齿形压力场。 二、速度场求解的困难 二、速度场求解的困难 2.求解方法的困难 压力本身并没有控制方程，它是以源项的型式出项在动量方程中。压力与速度的关系隐含在连续性方程中，如果压力场是正确的，则根据此压力场而求得的速度场必满足连续性方程。如何构造求解压力场的方程，或者说在假定初始压力分布后如何构造计算压力改进的方法，是一个关键