笛卡尔坐标系下不可压NS方程在方形通道中的奇点及-兰州交通大学.PDF

28 1 Vol. 28 No. 1 2009 2 Journal of Lanzhou Jiaotong University Feb. 2009 : 1001- 373(2009) 01-0125-05 N-S * 林志敏, 杨丽芳, 王 涛, 王良璧 ( , 730070) :在笛卡尔坐标系下, 用来描述方形通道中不可压缩 动的N-S 方程存在奇点. 以方形通道中充分发展层 为研究对象, 采用数值方法对动量方程的守恒性和奇点的影响区域进行了分析. 数值结果表明: 假定沿主 方向的 压力梯度为常数时, 求解简化的 N-S 方程而获得的主 方向的速度分布在横截面上的 个角点的附近区域内不满 足动量方程; 而采用连续性方程构造的压力修正方程的三维数值方法可以获得主 方向的压力梯度恒为常数, 但 速度分布在横截面上的 个角点的附近区域内也不满足动量方程. 得出奇点的影响区域主要集中在角点附近区 域, 在该区域的壁面上其影响较为显著. 为了获得物理意义上的解, 在数值计算中需要对奇点作 一些特殊处理. : N-S 方程; 方形通道; 奇点; 压力梯度 : O357. 1 : A [ 1] , . N-S , . 1 , , , , ¨# v = 0 , N-S 2 ( 1) 2 2 2 2 Q( v # ¨) v = - ¨p + L¨v L( 5 u/ 5y + 5 u/ 5z ) = dp / dx ( 2) ( v # ¨) v : , u = v = w = 0. , , N-S . N-S , . , ) ) ) . , . , N-S 1/ r , r = 0, . , . 1 , . Fig. 1 Schematic of internal flow in square duct , / 0 A E BF CG D H , , / 0 . : 2 2 2 2 , N-S 5 u/ 5y + 5 u/ 5z = 0 ( 3) , , ( 2, 3) , , N-S dp / dx = 0. , A E BF CG DH , * : 2008- 10-06 : ( 1982- ) , , , . 126 28 , . .