7.3_Navier-Stokes方程的解.ppt

高等流体力学 7 Navier-Stokes 方程的解 (第3部分) 7.6 缓慢流动 只有在一些特殊的流动条件下，才能在定解条件下求得N-S方程的精确解。除了N-S方程的精确解外，在流体力学的发展过程中，人们曾不断地寻求N-S方程的近似解。 这些近似求解有： 7.6 缓慢流动 ①小Re数条件下的近似解： 全部略去惯性项，全部保留粘性项； 全部略去惯性项，部分略去粘性项； 部分略去惯性项，全部保留粘性项。 ②大Re数条件下的近似解，如边界层的近似计算，将在第6章讲述。 7.6.1 Stokes流动 这是在N-S方程中全部略去惯性项(即非线性项)的流动问题。 当流体质点流动的流线(在定常流动时，与迹线重合)形状是直线(或同心圆)时，N-S方程中的惯性项(指位变惯性项)全部消失。如果不可压缩流体流动的流线形状与直线相差很小(或与同心圆周相差很小)时，则位变惯性项也很小，在近似处理时可以把它全部略去。 7.6.1 Stokes流动 可以认为，由于惯性力与速度平方成比例，而粘性力与速度的一次方成比例，当流速很小时，粘性力在合外力中占主导地位，而且粘性力比惯性力大得多(即Re1)，近似处理时可将惯性力全部略去。这种Re1的流动在流体力学中通常称为缓慢流动或蠕动。 7.6.1 Stokes流动 略去惯性项的N-S方程通常称为Stokes近似方程。在直角坐标系中，有 7.6.1 Stokes流动 对N-S方程式分别求x、y、z的偏导数，得 三式相加，得 7.6.1 Stokes流动 由于 所以 当忽略质量力的作用或者质量力可视作常数(如重力场中，fx＝fy＝0，fz＝-g)时，有 因此，缓慢流动中的压强p(x,y,z)是调和函数。 在给定边界条件时，由上述方程，可以求得流动的压强场，进而由Stokes近似方程求出速度场。 7.6.1 Stokes流动 现按照Stokes近似方程探讨球的缓慢运动问题 物体在流体中作等速直线运动引起的流体运动，如果相对于固定在物体上的坐标系来说，则流体运动是定常的，在求解这类流动问题时，往往采用固定在物体上的坐标系来进行研究。这样就把物体在静止流体中作等速直线运动的问题化为无界流体对静止物体的定常绕流问题，而这无界流体的自由来流速度与物体的运动速度大小相同，方向相反。以下即用物体的绕流问题来研究球的缓慢运动。 7.6.1 Stokes流动 设有半径为r0的物体，在充满不可压缩粘性流体的无界空间中以缓慢的速度U0作等速运动。由于球体半径r0甚小，运动速度也甚小，而运动粘性系数ν较大，因此Red0=U0d0/ν很小。Stokes认为，在小Re数的情况下，惯性力与粘性力相比有可能完全略去；此外，他假定绕流中不发生分离现象，且忽略质量力，又考虑到运动的轴对称性，即流动与φ无关，如下图所示。 7.6.1 Stokes流动 7.6.1 Stokes流动 这样，球坐标下的连续方程及N-S方程(Stokes近似方程)可简化成 (1) (2) (3) 7.6.1 Stokes流动 方程组的边界条件为 球面上(即r＝r0)：ur＝0，uθ＝0 无穷远处(即r→∞)：即ur＝U0cosθ， u?＝-U0sinθ，p＝p0 采用分离变量法求解上述方程组。 7.6.1 Stokes流动 由于流动的轴对称性，ur(r,θ)，u?(r,θ)，所以假设方程组的解具有如下形式，以便于分离变量。 (4) 将此形式的解代入方程组，先求出f1(r)、f2(r)及f3(r)，再获得ur(r,θ)、u?(r,θ)及p(r,θ)的具体表达式。 7.6.1 Stokes流动 7.6.1 Stokes流动 代入连续方程及N-S方程，得 7.6.1 Stokes流动 即 (5) (6) (7) 根据待求方程组的边界条件，有 (8) 7.6.1 Stokes流动 由连续方程(5)，得 (9) (10) (11) 将式(9)、(10)及(11)代入式(7)，得 (12) (13) 7.6.1 Stokes流动 将式(9)及(13)代入式(6)，得