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

第七章不可压缩流动的数值方法 §7.1基本方程 §7.2涡量-流函数方法 §73 SIMPLE方法  §7.1基本方程 单一介质不可压缩流动的控制方程 0 Ov +()=f-VP+1v v-运动黏性系数,p代表p/p 无量纲化形式 ar+v ()=f-Vp+ 压力场可以相差任一常数而对速度场无影响,为了确定全场压力值,应指定 流场中某一点的压力 控制方程的特点:不包含压力的时间导数项,表现出椭圆抛物组合型的特点。 不存在状态方程,压力不具有热力学意义;压力对速度场加以限制,使 连续性方程得到满足。压力场的波动具有无穷大的传播速度,瞬间传遍 全场,以使不可压缩条件在任何时间、任何地点满足。  目前,不可压缩流动的数值方法不如可压缩流动的数值方法成熟 求解不可压缩 Navier-Stokes方程的困难在于处理“压力-速度”耦合问题 不可压缩 Navier-Stokes方程的一种可能的求解方案 动量方程和连续方程完全耦合求解。 +D(+cin7=2Re2+“2+ 其中,D(),=D【·()+DLj()G(),=D()+D,(*)j D,(*) D,() △x (*)+() (*).,+() i+12, +1/2 △ 如果对流项采用显式处理,如:用时间方向具有二阶精度的 Adams- Bashforth格式离散 D()12=1.D(m-0.5D(mn  得到一个非常庞大、形状不规则、刚性很强的稀疏线性方程组。计算量非常大 不易收敛 如果对流项采用隐式处理,得到一个非线性方程组,其求解更为困难 因此,完全耦合求解方法在实用中很少使用。 实用中,为了处理压力速度耦合问题,需要把连续方程和动量方程在一定程度 上进行解耦。 §72涡量一流函数方法 计算二维不可压缩 Navier -Stokes方程的有效算法 基本思想:对 Navier- Stokes方程进行变换,写成涡量一流函数的形式。 7.2.1基本方程 (7.2.1) au at Ov..O_Op1(0v.av at ax Rel ox?  引入涡量方程,消去方程中的压力项 涡量m=v1-l2 (7.22)对y求偏导,(7.2.3)对x求偏导,推导出涡量满足的方程为 (7.2.4) 引入流函数ψ,使连续方程自然满足 V (72.5) 满足 (7.2.6 求解过程:通过(7.24)计算涡量 通过(7.26)计算流函数 通过(7.25)计算速度分量。 如果需要计算压力,把x方向的动量方程对x求导数,把y方向的动量方程对y求导数, 二者求和后,利用连续性方程,可以得到压力的 Poisson方程, 求解压力的 Poisson方程,可以得到压力的分布  7.2.2差分格式 数值求解涡量满足的方程(对流-扩散方程) @+uo+ vQ 1.如果采用FTCS格式 2.,+m.,-1 R 2.如果采用迎风格式 其中,() (·) 2 当∝=0时,一阶迎风格式 当=l时,二阶迎风格式  3.如果采用 Crank-Nicolson格式 R 4△x (6。+6a#) 4Ar,j+8.o+ 1(an,6,a,6:a,+8:a 2Re(△x)2(△y)2(△x)2(△y) 数值求解流函数满足的 Poisson方程 二阶精度的差分格式为  7.23边界条件 图71:二维槽道流动的计算域 (1)固体壁面边界:j=0,j=M 边界条件为 (2)进口边界:i=0 边界条件为 (y) 分量形式(假定进口流动方向是水平的)  (3)出口边界:i=M 当槽道的长度远大于其高度时,边界条件可以认为是充分发展条件, 即物理量(除压力)沿x方向的导数等于零 在涡量-流函数方法中,边界条件的处理: 1)固壁边界。固体壁面是流线, 在j=0的固壁上,可以指定v0=0 在=M的固壁上,可以指定.=a24=coms 在=固壁上,叫=-n--=y 离散形式a0=-%=2y,。+ 其中,虚拟网格点(,-1)处的y_是未知的 利用 2y 0  (2)进口边界 流函数:可以用v。4(y),通过适当的数值积分方法计算 涡量:a wl-2yo +y1. yoj +-2yo. +yo. 利用 V1.-V 41-,)v+-2y1+ (3)出口边界。假定流动是充分发展的,有 当采用一阶近似时 为了求解压力 Poisson方程,还应补充压力的边界条件。在固壁处,Vp=△a 利用法向分量作为边界条件,即vpn=Re△