5西南交通大学计算流体力学

计算流体力学电子教案西南交通大学应用力学与工程系结构分析教研室喻勇2010-4-24（第一版,2004级专用）目录•第一章绪论•第二章扩散问题的有限体积法•第三章对流扩散问题的有限体积法•第四章差分格式问题•第五章压力--速度耦合问题的有限体积法•第六章有限体积法离散方程的解法•第七章非稳态流动问题的有限体积法•第八章边界条件处理第五章压力--速度耦合问题的有限体积法本章内容1压力--速度耦合问题的计算难点2交错网格技术3SIMPLE算法4SIMPLE算法的改进5-1压力--速度耦合问题的计算难点对流扩散问题的微分方程中没有直接考虑压力梯度项（可以认为压力项归入了源项中），而压力梯度是引起流体流动的最直接的动力，压力场在流场分析中是需要求解的，而它与速度场是相互耦合、相互影响的一、压力--速度耦合问题方程vvvvgradf)(12tp不可压缩流体N-S方程zvvyvvxvvtvvxpfxzxyxxxxx21MxxxxSxpvdivvdivtv)()()(gradv2维问题X方向：X方向通用形式：--见P1(1-2)式()()()()()xyxxxxxMxvvvvvvvpStxyxxyyx展开得：()()()()xyxxxxMxvvvvvvpSxyxxyyx对于稳态问题，忽略时间项：同理得y方向N-S方程：()()()()xyyyyyMyvvvvvvpSxyxxyyy()()0yxvvxy2维稳态流动的连续性方程：以上三式相互耦合，相互影响：速度场要满足连续性方程和N-S方程，压力场要满足N-S方程，同时压力场又会影响速度分布。二、计算难点难点1：同位网格有可能不能识别棋盘形压力场对于一维问题，在均匀网格中用中心差分格式求压强梯度项：222PWEPewEWpppppppppxxxx上式不能识别图5-1的锯齿关压强分布，因为由上式结果得压强梯度处处为0.同位网格：控制体中的u、v、p均存放于同一套网格的节点上。对于二维问题，有：,22EWNSppppppxxyy上式不能识别图5-2的压强分布，因为由上式结果得压强梯度处处为0.难点2：压力场难以求解在二维问题中，对于N-S方程采用分离式求解法：在u、v、p中，求其中一个变量时认为另两个为已知，经过一次计算不能得到正确结果，需要迭代求解。此时会遇到的一个问题是：压力本身没有控制方程，它是以源项形式出现在动量方程中。压力与速度的关系隐含在连续性方程中，如果压力场是正确的，则据此压力场而解得的速度场必然满足连续性方程。因此，这就要求分离式求解时能从速度场的计算结果得到改进压力场的计算式，但这一改进式无法直接得到。()()()()xyxxxxMxvvvvvvpSxyxxyyx()()()()xyyyyyMyvvvvvvpSxyxxyyy()()0yxvvxy三、解决办法难点2：压力场难以求解--SIMPLE算法难点1：同位网格有可能不能识别棋盘形压力场----采用交错网格SIMPLE算法—Semi-ImplicitMethodforPressure-LinkedEquation压力耦合方程的半隐计算格式5-2交错网格技术一、交错网格布置方式交错网格：以二维问题为例，将速度u、v及压力p（包括其它所有标量场及物性参数）分别存储于三套不同网格上的网格系统。以压力控制体为主控制体（称P控制体），速度u存于P控制体的东、西界上速度v存放于P控制体的南、北界面上。u、v各自的控制体是以速度所在位置为中心的，与P控制体在x或y方向相差半个网格步，分向前错位和向后错位。P控制体ue控制体vn控制体采用交错网格后，关于u、v、p的离散方程是通过对u、v、p控制体积分得到。例如，x方向N-S方程在uw控制体积分。()()()()xyxxxxMxvvvvvvpSxyxxyyx()()()()MxuuuvuupSxyxxyyx()()[][()()]nPnPMxsWsWuuuvuupdxdySdxdyxyxxyyxx方向N-S方程即：在uw控制体上积分得（见uw控制体）：()()()nPnPPWWsWspdxdypdyppyx上式中：将此式移抄至下页()()()ePePPSSwSwpdydxpdxppxy()()()nPnPPWWsWspdxdypdyppyx同理得：以上两式可以处理图5-1和5-2所示的压力分布。为什么？（5-6）（5-7）交错网格系统的3套网格及编号系统，使得编号及编程都比较复杂标量存于大写节点，矢量存放于小写节点。i系列编号表示u控制体的节点，j系列编号表示v控制体节点，大写的I或J表示主控制体在x、y方向的节点。u控制体i,Jv控制体I,jp控制体I,J注意编号中大小写的组合二、方程的离散包括x方向、y方向动量方程的离散以及连续性方程的离散。动量方程的离散与2~4章相同，只是多了一项压力梯度项，压力梯度项的离散见（5-6）式和（5-7）式()()()nPnPPWWsWspdxdypdyppyx()()()ePePPSSwSwpdydxpdxppxy（5-6）（5-7）1.x方向动量方程的离散unbnbPPSaa,1,,,IJIJiJiJnbnbuuuppauauVSVx,,,1,,,()iJiJnbnbIJIJiJiJauauppAb前2章通用离散方程：本章以速度代替通用变量，并考虑压强：或：,,,1,,,()iJiJnbnbIJIJiJiJauauppAb该式中的系数ai,j及anb由任意一种差分格式得到，如上风差分、混合差分、乘方格式或QUICK格式等。而系数ai,j及anb是F与D的组合，现讨论F、D的表达式。u控制体i,J(注意该图与图5-3水平错动了半格)对于均匀网格u控制体i,J的对流量F（沿流动方向线性插值）：,1,1,,,1,,11,1()2()2()2()2iJiJwwiJiJeeIjIjssIjIjnnFFFuFFFuFFFvFFFv,1,,,,,,,1,,,,()2()2IJIJiJiJiJiJIJIJIjIjIjIjFuuFvv以上四式中：其余类似。代入四式得：,1,,1,1,2,,1,1,,1,,,1,1,,,1,,,11,1,1,1,,11,()244()244()244()iJiJIJIJIJIJwwiJiJiJiJIJIJIJIJeeiJiJIjIjIJIJIJIJssIjIjIjIjnnFFFuuuFFFuuuFFFvvvFFFv1,1,1,11,,11,1244IJIJIJIJIjIjvv注意：u、v的下标区别;式中的u、v为上一迭代步结果，视为已知。1,,11,IJIJweiiiiDDxxxx对于均匀网格u控制体i,J的扩散量D：1,,1,1,111,1,11,,14()4()IJIJIJIJsJJIJIJIJIJnJJDyyDyyy方向动量方程的离散同此，即：,,,1,,,()IjIjnbnbIJIJIjIjavavppAb上式系数中包含F和D，各节点的F、D表达式为：2.y方向动量方程的离散,,1,1,1,1,1,,11,1,11,,,11,11,1,1,1,,1,2,,1,1,()244()244()244()iJiJIJIJIJIJwwiJiJiJiJIJIJIJIJeeiJiJIjIjIJIJIJIJssIjIjInnFFFuuuFFFuuuFFFvvvFFv,,1,,1,1,,,1244jIjIJIJIJIJIjIjFvv1,,1,1,11,11,1,1,14()4()IJIJIJIJwIIIJIJIJIJnIIDxxDxx,1,11,IJIJsnjjjjDDyyyy3.连续性方程的离散在主控制体中离散，方法同前。()()0yxvvxy2维稳态流动的连续性方程：()()0uvxy1,,,1,[()()][()()]0(512)iJiJIjIjuAuAvAuA即：上式在主控制体上积分并利用高斯公式得离散方程：课本中误为“-”号,,,1,,,()(510)IjIjnbnbIJIJIjIjavavppAb,,,1,,,()(58)iJiJnbnbIJIJiJiJauauppAb总结交错网格条件下，二维压力速度耦合问题的有限体积法的离散方程组为以上三式分别对应三套网格中的控制体如果压力场能满足（5-8）和（5-10）式，且由这两式解出的速度场满足（5-12）式，且该压力场和速度场为正确解。1,,,1,[()()][()()]0(512)iJiJIjIjuAuAvAuA难点2：压力场难以求解--SIMPLE算法难点1：同位网格有可能不能识别棋盘形压力场----采用交错网格备忘：压力速度耦合问题的难点5-3SIMPLE算法一、基本SIMPLE算法该算法基本思想：是一种压力预测--修正方法，通过不断修正计算结果，反复迭代，最后求出p、u、v的收敛解。,,,1,,,**(**)(514)iJiJnbnbIJIJiJiJauauppAb,,,1,,,**(**)(515)IjIjnbnbIJIJIjIjavavppAb基本思路（以二维问题为例说明）：1.假设一个压力分布p*，用它求解动量方程式，得到初始速度分布u*、v*，即：一般，由于压力场不是真实的，由此得到的速度场也不能满足连续性方程，因此需对速度场和压力场进行修正。分别设它们相应的修正量是u’、v’、p’，修正后的速度场和压力场为：u=u*+u’v=v*+v’p=p*+p’为了求得修正量，设正确的压力场为p，由(5-8)和（5-10）式（x、y向离散的动量方程，可得到正确的速度场u、v。,,,1,,,**(**)(514)iJiJnbnbIJIJiJiJauauppAb,,,1,,,()(58)iJiJnbnbIJIJiJiJauauppAb（5-8）式减（5-14）式得,,,,1,,1,,(*)(*)[()(**)](519)iJiJiJnbnbnbIJIJIJIJiJauuauuppppA即：1,,,,,,'('')'(521)nbnbIJIJiJiJiJiJauppAuaa同理得：,1,,,,,('')''(522)IJIJIjnbnbiJIjIjppAavvaa由上可见，速度修正u’、v’由两项构成。从物理意义上讲，第一项为周围节点速度引起的修正量，第二项是同一方向相邻节点压力差引起的修正量。SIMPLE算法认为p’是引起速度改变的主要因素，因此作简化处理，只保留第二