【水力机械水动力学】5-弱可压缩流体理论及应用.docxVIP

旗开得胜 PAGE 1 PAGE 1 读万卷书 行万里路 5 弱可压缩流体理论及应用 5.1控制方程 Navier-Stokes 方程是流体力学的基本方程。 5.1. 连续方程： （5-1） 运动方程： （5-2） 上述二式与坐标无关，为矢量，未知数为 。运动方程中有 ，而连续方程中仅有 没有，无法耦合求解。若先假设压力场，据此可解出速度场，但解出的速度场一般都不满足连续方程，这是因为压力场假设不正确。实际上运动方程中的是压力梯度，或者说是压力差，这就使得假定的压力场很难准确。 5.1. 连续方程： （5-3） 运动方程： （5-4） 上述方程中比不可压缩流体方程多了一个未知量，此时方程组不封闭，求解时需要加上状态方程或，即密度与压力或体积的关系。 对水体可写出如下状态方程： （5-5） 式中：为水的密度，K为水的弹性模量。对上式稍做变化可得到：，式中：C为波速。将其代入可压缩流体连续方程（1-3）有： （5-6） 若流体不可压缩，则，即，方程变成不可压缩流体连续方程。 5.1 引入人为压缩性和人为状态方程对不可压缩流体的连续方程施加干扰，将连续方程写为： （5-7） 式中：为极小的参数，为人工密度。人为状态方程可写为： （5-8） 式中：为人为压缩系数。经过如此处理，原则上可由方程（1-7）解出各时刻的密度，再由式（1-8）解出各时刻的压力。由于要求的值很小才能在每个时刻都得到比较精确的结果，因而计算过程中时间步长必须很小，这就限制了这种方法的广泛应用。 5.1 因为在水动力学中，主要关心的是压力而不是密度，所以需要将以粘性为独立变量的原始连续方程转换为以压力为独立变量的连续方程。对于正压流体，状态方程可写为： （5-9） 式中a 是声速。 应用方程（5-9）消去原始连续方程中的密度，可得到下列包含压力项的连续方程： （5-10） 上述方程是压力输运方程，式中最后一项表示由于体积变化而导致的压力改变。 当马赫数很小时密度变化也很小，状态方程可由下列直线方程代替： （5-11） 上面的方程中，下标表示参考点的参数。考虑线性状态方程，连续方程变为： （5-12） 通过无量纲分析可知，对流项的阶数为马赫数的平方，因此可以忽略不计。然而可压缩水动力学的关键是保留了第一项，当水流迅速变化时该项的数值会很大。对于流动突然开始的问题，不管压缩性多小，在接近t = 0 时该项的数值非常大。因此对于可压缩水动力学，最终采用的连续方程为： （5-13） 式中 是体积弹性模量。 如果流体是不可压缩的，K为无穷大，可得到熟知的不可压缩流体的连续方程： （5-14） 由于缺少时间导数项，方程（5-14）只能用于外部区域，即时间t 很大的情况亦是如此。但方程（5-13）既可用于外部也可用于内部区域，这是因为当时间t增大时，方程（5-13）中的第一项趋于零并自动消失，这时方程（5-13）连续地趋向于方程（5-14）。宋教授称满足方程（5-13）的流动为弱可压缩流动。但很明显，称其为可压缩水动力学连续方程或可压缩边界层方程更合适。 不可压缩流体的Navier-Stokes 运动方程或相应的Euler 方程在接近零马赫数时没有奇点问题。在内部以及外部区域，可压缩水动力学Navier-Stokes 运动方程的表达式为： （5-15） 方程（5-13）和方程（5-15）组成包含四个未知数p，u，v，w的四个方程组。这些方程组可写成下列守恒形式： （5-16） 式中：