04-第四章-速度场计算-SIMPLE算法

第四章速度场计算－SIMPLE算法•第一节速度场求解的困难•第二节交错网格及动量方程的离散•第三节压力修正算法•第四节SIMPLE算法及其发展与改进•第五节开口系统的流场计算•第六节封闭系统的流场计算•第七节设定计算区域值的方法第一节速度场求解的困难•一、方程的通用型式•二、速度场求解的困难一、方程的通用型式•连续性方程：•动量方程（x方向）：•动量方程（y方向）：•能量方程：xpyuxuyuvxuutu1)(2222ypyvxvyvvxvutv1)(22220yvxutppcQyTxTckyTvxTutT)(2222一、方程的通用型式•写成统一型式Syyxxyvxut)()()()()(方程S连续性方程100动量方程（x方向）uxp/动量方程（y方向）vyp/能量方程Tpck/pcQ/二、速度场求解的困难二、速度场求解的困难•1.一阶导数的存在•考虑一维常密度问题，有：•可离散为xp22dxuddxdpdxduu2111111)(222xuuuxppxuuuiiiiiiiii二、速度场求解的困难•考虑连续性方程：•由于•故•所以•亦有•可得到不合理的锯齿形压力场。constuconst11iiiuuu11iipp22iiippp二、速度场求解的困难二、速度场求解的困难•2.求解方法的困难•压力本身并没有控制方程，它是以源项的型式出项在动量方程中。压力与速度的关系隐含在连续性方程中，如果压力场是正确的，则根据此压力场而求得的速度场必满足连续性方程。如何构造求解压力场的方程，或者说在假定初始压力分布后如何构造计算压力改进的方法，是一个关键问题，即所谓的速度和压力的耦合问题。二、速度场求解的困难•3.解决方案的选择：•（1）涡量－流函数法；•（2）交错网格。√第二节交错网格及动量方程的离散•一、交错网格•二、动量方程的离散特点•三、建立离散化方程编制程序一、交错网格•所谓交错网格就是把速度u、v及压力p（包括其它所有标量及物性参数）分别存储在三套不同网格上的网格系统。u控制容积和主控制容积（即压力控制容积）之间在x方向有半个网格步长的错位，而v控制容积与主控制容积之间则在y方向上有半个步长的错位。一、交错网格一、交错网格一、交错网格二、动量方程的离散特点•1.积分用的控制容积不是主控制容积而是u、v各自的控制容积；•2.压力梯度项从源项中分离出来。对ue的控制容积，该项积分为nsEPEPnsEPyppdypdxdyxp)()()(二、动量方程的离散特点•关于ue的离散方程有以下形式：•其中•类似可得到关于vn的离散方程：•3.各控制容积界面上的流量、物性参数需要用插值的方法进行处理。eEPnbnbeeAppbuaua)(00eeCuavSbyAenNPnbnbnnAppbvaua)(三、建立离散化方程编制程序•1.三类变量的节点编号方法•对主节点（i，j），其控制容积西界面上的流速为ui,j，南界面上的流速为vi,j。•对主节点Φi，i＝1－L1，j＝1－M1。对ui,j，i＝2－L1，j＝1－M1；对vi,j，i＝1－L1，j＝2－M1三、建立离散化方程编制程序三、建立离散化方程编制程序三、建立离散化方程编制程序•2.与边界相邻接的速度控制容积与内部速度控制容积的不同。三、建立离散化方程编制程序•3.边界压力的递推计算21,3,2,3,1)(llppppjjjj第三节压力修正算法•一、速度修正方程•二、压力修正方程一、速度修正方程•先考虑如何修正速度方程。设原来的压力为p*，与此对应的速度分别为u*，v*。压力修正值为p′，相应的速度修正值为u′，v′。则改进后的速度与压力分别为u=u*+u′,v=v*+v′,p=p*+p′,代入动量方程有：eEEPPnbnbnbeeeAppppbuuauua)]()()()(''''一、速度修正方程•由于u*，v*是根据p*求解出来的，故：•两式相减有：•可以认为任一点上的速度改进值由两部分组成：一部分是与该速度在同一方向的上的相邻两节点之间的压力修正值之差，这是产生速度修正值的直接动力；另一部分是由邻近速度的修正值所引起的。eEPnbnbeeAppbuaua)(eEPnbnbeeAppuaua)(''''一、速度修正方程•为计算简便考虑，后一项可忽略不计，即anb=0，于是可得到速度修正方程：•或•类似可得到•其中•改进后的速度：eEPeeAppua)(''')())(('''''EPeEPeeeppdppaAu)('''NPnnppdvnnnaAd)(''*EPeeeppduu)(''*NPnnnppdvv二、压力修正方程•对连续性方程在时间间隔内对主控制容积进行积分，采用全隐格式，可得到•代入速度修正方程并整理成p′的代数方程：0])()[(])()[(0xvvyuuyxtsnwePPbpapapapapaSSNNWWEEPP'''''二、压力修正方程•其中：•b的数值代表了一个控制容积不满足连续性的剩余质量的大小，可以用各控制容积的剩余质量的绝对值的大小作为速度场迭代是否收敛的一个判据或指标。ydaeeEydaxdannNxdassESNWEPaaaaaxvvyuuyxtbnsewPP])()[(])()[(0二、压力修正方程•压力修正方程的边界条件：•（1）若边界压力已知则p′=0，故aB=0。•（2）若法向速度已知，则u′=0，p′=0，故aB=0第四节SIMPLE算法及其发展与改进•一、SIMPLE算法的计算步骤•二、SIMPLE算法的讨论•三、SIMPLER算法•四、SIMPLEST算法•五、SIMPLEC算法一、SIMPLE算法的计算步骤•SIMPLE（Semi-Implicitmethodforpressure-LinkedEquations）算法，即求解压力耦合方程的半隐方法，是Patankar和Spalding在1972年提出来的。所谓“半隐”是指在计算速度的修正值时，忽略了邻近速度修正值的影响，否则即称为“全隐”。一、SIMPLE算法的计算步骤•1.假定一个速度分布，记为u0，v0，以此计算动量离散方程中的系数及常数项；•2.假定一个压力场p*；•3.依次求解两个动量方程，得到u*，v*;•4.求解压力修正方程，得到p′；一、SIMPLE算法的计算步骤•5.据p′改进速度场；•6.利用改进后的速度场求解与速度场耦合的变量；•7.利用改进后的速度场重新计算动量离散方程的系数，用改进后的压力场作为下一层次迭代计算的初值，重复计算。直到获得收敛的解。二、SIMPLE算法的讨论•1.速度修正中忽略了邻近速度修正值的影响，不影响最后收敛的解，但加重了修正值p′的负担，使得整个速度场的迭代速度减慢，故对p′作亚松弛，•其中一般可取0.8左右。•类似速度修正值也可考虑亚松弛，一般松弛因子可取0.5左右。'pppPP二、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算法（SIMPLERevised）。三、SIMPLER算法•2.压力方程的推导•动量离散方程可写成：•其中前一项称为假拟速度，速度可记为：)(EPeuenbnbeppdabuaue)(EPeeeppduu)(NPnnnppdvv三、SIMPLER算法•将上两式代入连续性方程的离散形式，可得到与压力修正方程形式相同的压力方程：bpapapapapaSSNNWWEEPPxvvyuuyxtbnsewPP])()[(])()[(0三、SIMPLER算法•3.SIMPLER算法的计算步骤•（1）假定一个速度场u0，v0，计算动量方程的系数。•（2）据已知的速度计算假拟速度，。•（3）求解压力方程。•（4）把求解的压力作为p0，求解动量方程，得u*，v*。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算法中，为求解方便，略去了速度修正值中的和，从而犯了速度和压力不协调一致的错误。为此在速度修正方程两端同时减去，可得到：•可略去前一项，于是有：'nbnbua'nbnbva'enbuaeEPenbnbenbeAppuuauaa)()()(''''')('''EPnbeeeppaaAu)('''NPnbnnnppaaAv五、SIMPLEC算法•2.SIMPLEC（协调一致的SIMPLE算法）的特点•（1）以代替。•（2）在SIMPLEC算法中，p′不再亚松弛，即取。)/(nbeeaaAeeaA/1P第五节开口系统的流场计算•一、开口系统流场计算的关键•二、出口边界条件的处理方法一——充分发展•三、出口边界条件的处理方法二——取均匀的流场•四、出口边界条件的处理方法三——从内点的速度大分布来获得出口截面上的速度分布。第五节开口系统的流场计算•一、开口系统流场计算的关键•开口系统流场计算的关键在于出口截面位置的确定和出口截面上法向流速的确定。第五节开口系统的流场计算•二、出口边界条件的处理方法一——充分发展•一般只有当出口区域有一平直