计算流体CFD2004-5.pptVIP

第五章　不可压缩粘性流的有限差分法 回顾 2.2 Navier-Stokes 定常N-S方程组(二维标量式) ②守恒形式： 守恒变量 三维流体力学控制方程组的性质列表如下： 5.2人工压缩性方法 5.3 压力修正法 压力修正方程推导 例1：一维稳态动量方程 例2 1、SIMPLE 2、SIMPLER 3、SIMPLEC 4、SIMPLEX 5、SIMPLET 6、MAPLE算法 SIMPLE系列的比较 7、PISO算法 SIMPLE系列算法的加速迭代 5.4 涡方法 ?和?的边界条件： 流函数?，涡矢量 ? 补充压力方程 * 华中科技大学能源与动力工程学院 5.1 概况 不可压缩流 M0.3 不可压缩 N-S方程组，两类变量 原始变量（p,u,T） 守恒变量（ρ, ρu, ρE） 不可压缩流的处理，两种方法 原始变量：ACM，SIMPLE，PISO 涡方法： κ-ω 由基本方程(质量、动量、能量三大守恒) ①N-S方程组(三维，矩阵式，原始变量，非守恒) 得N-S方程组： 或(2.2.9a,b,c) (2.2.9b) (2.2.9c) (2.2.9a) （2.2.11） U U守恒流动变量，F对流通量变量，G扩散通量变量,B源项 U三维展开 椭圆一双曲 椭圆一抛物 抛物一双曲 抛物一双曲 非定常 椭圆一双曲 椭圆 亚声速： 椭圆一抛物 超声速： 双曲一抛物 椭圆一双曲 定常 不可压 Euler 不可压 N．S 可压 Euler 可压 N．S 椭圆型与双曲或抛物型混合的方程组处理起来很困难，必须在数学上进行讨论与分析。（定解条件提法，论述求解方法、分析求解时应注意的问题。） 单一类型方程组——平衡型（椭圆型）或发展型（双曲、抛物型）比较容易处理，分别为纯边值、初边值问题 在不可压缩流场的控制方程的动量方程中，压力梯度是以源项形式出现，压力没有独立的方程(在可压缩流体中压力与密度问的关系由状态方程所规定)。在以速度、压力为求解变量的原始变量法中，为了解决压力没有独立的方程的困难，先后发展了 ●人工压缩性方法 ●压力修正算法(pressure—correction method)等方法。其中以压力修正算法——如SIMPLE系列的应用最广 。 ● κ-ω方法不是原始变量法，它避开了压力项。 不可压定常N-S方程是椭圆型方程，其他不可压情况都是混合型方程。动量方程中有压力项，使压强信息瞬间传遍全域，不能按时间推进求解 其中　是人工密度，与人工压缩性参数β的关系为： ，在稳态条件下 。 是伪时间。 SIMPLE(Semi-Implicit Method for Pressure-Linked Equations)算法虽然最初是对不可压缩流体发展起来的，但近年来已成功地推广到了可压缩流体。 5.3.1压力连接方程的半隐含法（SIMPLE） 已知压力场，可由动量方程求得速度分布，而压力应该由质量守恒条件约束。但无压力信息，迭代无法进行。为此，引入压力修正方程 两方程之间压力无关联 设压力修正为p’,相应速度修正为u’,v’ 忽略有关项后 ut’+ px’/ρ=0 u’= -τpx’/ ρ 同理 vt’+ py’/ρ=0 得压力修正方程 τ(pxx’+pyy’)/ρ= ūx+ ūy 将ū, , p*与u=ū+u’,v= +v’,p=p*+p’分别代人动量方程后相减 ut’+ūux’+u’ūx+u’ux’+ uy’+v’ūy+v’uy’+px’/ρ=ν(uxx’+u yy’) 由连续方程ux+vy=ūx+ +ux’+vy’=ūx+ +[-τ(pxx’+pyy’)/ρ]=0 如果ui = ui-1 =ui+1；则Pi-1 =Pi+1; p的锯齿波不能检测出 p x 下角标nb表示四个邻点，Ae为压力作用面积，系数a e,a nb之值取决于格式，b为源项。 图5.3.1 交错网格的计算区域 即节点上p , T ； 边界上u , v ； 例如二维 在交错网格上离散: 有限容积的边界e、n e n (as，aw)可仿此写出 b为控制容积P的剩余质量流量，对稳态流动，有 压力修正方程 采用SIMPLE算法的步骤： (1)假定一个速度场，记为u(0)，v(0)，由此计算ae, an，anb，及b； (2)假定一个压力场，记为p*； (3)求解动量离散方程，得ū ， ； (4)计算压力的修正值p’：要求与(p*+p’)相对应的(ū +u’)、( +v’)能满足连续性方程。这样的压力修正方