不可压缩流动数值方法初步.pptVIP

第七章不可压缩流动的数值方法 §7.1基本方程 §7.2涡量-流函数方法 §7.3 SIMPLE方法  7.1基本方程 单一介质不可压缩流动的控制方程 v-运动黏性系数,p代表p/p 无量纲化形式 Ou.v(uu)=f-vp* Re 压力场可以相差任一常数而对速度场无影响,为了确定全场压力值,应指定 流场中某一点的压力 控制方程的特点:不包含压力的时问导数项,表现出椭圆-抛物组合型的特点 不存在状态方程,压力不具有热力学意义;压力对速度场加以限制,使 连续性方程得到满足。压力场的波动具有无穷大的传播速度,瞬间传遍 全场,以使不可压缩条件在任何时间、任何地点满足  日前,不可压缩流动的数值方法不如可压缩流动的数值方法成熟。 求解不可压缩 Navier-Stokes方程的困难在于处理“压力-速度”耦合问题 不可压缩 Navier-Stokes方程的一种可能的求解方案 动量方程和连续方程完全耦合求解 其中,D(*)=D[(+D,[j(*);(*)=D(*)+D() D(). -12 2(+,=(*)r2 (*) (*),+() 1/2, +12 (),-2(*),+();() 如果对流项采用显式处理,如:用时间方向具有二阶精度的 Adams- Bashforth格式离散 D()212=1.5D(  得到一个非常庞大、形状不规则、刚性很强的稀疏线性方程组。计算量非常大, 不易收敛 如果对流项采用隐式处理,得到一个非线性方程组,其求解更为困难 因此,完全耦合求解方法在实用中很少使用。 实用屮,为了处理压力速度耦合问题,需要把连续方程和动量方程在一定程度 上进行解耦。 §72涡量一流函数方法 计算二维不可压缩 Navier-Stokes方程的有效算法 基本思想:对 Navier-Stokes方程进行变换,写成涡量一流函数的形式 7.2.1基本方程 (72.1) a+t2o+,a=-p2+1 ay ax rel oxto (7.22 ap 1(av. a2v l=+1  引入涡量方程,消去方程中的压力项 涡量m 72.2)对y求偏导,(72.3)对x求偏导,推导出涡量满足的方程为 ①.+0+1O (724) 引入流函数ψ,使连续方程自然满足 aψ= (7.25 y满足 求解过程:通过(7.24)计算涡量 通过(726计算流函 通过(7.2.5)计算速度分量 如果需要计算压力,把x方向的动量方程对x求导数,把ν方向的动量方程对y求导数 者求和后,利用连续性方程,可以得到压力的 Poisson方程, 求解压力的 Poisson方程,可以得到压力的分布