求解N-S方程的压力修正.pptVIP

*t§12求解N—S方程的压力修正方法压力修正的意义：对动量方程的离散形式进行求解时，我们迭代的压力场P*计算出来的速度场，未必可以满足连续方程。所以我们要寻找一个改进压力估计值P*的方法，以使所得的带*的速度场逐渐的满足连续方程。01压力修正方法的基本思想：在对N—S方程的离散形式进行迭代求解的任一层次上，可以给定一个压力场，它可以是给定的或是上一层次计算出来的。但由它所计算出的速度场未必满足连续性方程，所以要对压力场进行修正。02修正原则：与改进后的压力场相对应的速度场能满足这一迭代层次上的连续性方程，据此来导出压力与速度修正值，并以修正后的压力与速度开始下一层次的迭代计算。03压力修正法步骤：假设一个压力场p*。利用p*，求解动量离散方程，得出相应的u*,v*(不完善的速度场)。改进压力场，要求与改进后的压力场相对应的速度场能满足连续性方程。即用p’,u’,v’分别表示压力与速度的修正量，称之为压力修正值（pressurecorrection）与速度修正量(velocitycorrection)。以正确的压力（p*+p’），及相应的（u*+u’），（v*+v’）作为本层次的解，并据此开始下一层次的迭代计算。1速度修正值的计算公式改进后的压力场和速度场满足这一迭代层次上的动量离散方程：引出两个关键问题：如何获得压力修正值p’以使与p相应的u、v满足连续性方程?获得压力修正值p’后，如何确定u*,v*?因为u*e是据p*从这一动量离散方程中导出的，因此也满足动量离散方程，虽可能不满足连续性方程。两式相减得：说明：任意点速度的改进值u’e由两部分组成：1)该速度在同一方向上的相邻两节点间压力修正值之差——产生速度修正值的直接动力。2)由邻点速度修正值所引起的，可视为是四周压力的修正值对所讨论位置上的速度改进的间接影响。为了简化，忽略间接影响，相当于假设∑anbu’e中anb=0的。得到速度修正方程：01或类似可得于是改进后的速度为：020304052求解压力修正值的代数方程压力修正值p’应满足的条件：由p’而改进的速度场u、v能满足连续性方程。将式1代入连续性方程的离散格式即可获得p’的代数方程。01在时间间隔Δt内对主控制体做积分，得：02将（1）代入并整理可得：03其中：-------压力修正值的代数方程。b项是按带星号的速度取值的离散化连续性方程左侧的负值。如b=0，说明该速度场已满足连续性方程，不必对压力进一步进行修正，p’方程迭代已收敛，即由该压力场求解的速度场定满足二守恒方程；可用各控制容积b值（不满足连续性方程的剩余质量的大小）作为判断速度场收敛的判据。收敛速度与p’的特定方程有关，如忽略项过多，会影响收敛速度。如b≠0，重新校正压力场p’后，再求b值，直至b≈0时，p’压力场可使速度分布满足连续性方程；由p’计算而得的u’，v’可使u、v能满足连续性方程，以u，v作为本层次上的速度场的解，并用它去改进动量离散化方程系数，开始进行下一层次的迭代计算；010302说明：因求解的问题是非线性的，没有必要在每一层次迭代上去获得p’真解，采用ADI方式对p’场作3-5次扫描后即可停止迭代（初始令各点p’=0）。3压力修正值方程的边界条件工程流场计算中边界条件一般有两种。边界上的压力已知：则p’=0边界上的法向速度已知，设ue已知，则u’e=0，则不必引入p’E。说明：不论是边界压力已知，还是法向速度已知，都没有必要引入关于边界上压力修正值的信息。在计算中，可令与边界相邻的主控制容积的p’方程相应的影响系数为零。014SIMPLE算法（Semi-ImplicitMethodforPressuerLinked）（解压力耦合方程的半隐方法）02SIMPLE法就是一种计算压力修正的一种方法。03半隐含义：所谓半隐式是在公式中略去了某一项的处理方法。略去间接影响保留直接影响的称为半隐式,全保留是全隐式041后部分是压力修正对u’e的直接影响；前一部分是压力修正对u’e的间接影响，忽略该项称为半隐；2计算中若保留∑anbu’nb项，则u’e方程就是全隐的代数方程。势必存在u’e的相邻点的压力修正和速度修正，相邻点又引入新的相邻点，这样使计算变复杂，无法进行逐点计算。但此项的忽略并不影响最终收敛解。3simple算法的计算步骤4假定一个速度分布，记为u0，v0，以此计算动量方程中的系数及常数项。（源项中压力梯度无法求）据p’改进速度值；（以使新速度u满足连续性方程但可能不满足动量守恒）假定一个压力场p*；求解压力