第二章流体流动的数学模式2012.pptVIP

定义沿水深平均流速和时均流速的关系 引用莱布尼公式 定义水深 自由表面及底部运动学条件 略去流速等水力参数沿垂直方向的变化 沿水深方向的动水强度分布符合静水强度分布 沿水深平均的连续方程 最终形式 沿水深平均的动量方程 表面风阻力 表面风阻力 注意以下几方面问题 方程在推导过程中引用了牛顿流体所满足的本构关系，对类似高含沙水流的非牛顿流体不适用。 方程在推导过程中对流体作了均质不可压的假设，因此只能在含沙量较小的情况下近似使用，含沙量较大时，应考虑密度的变化。 在垂向积分过程中，略去流速等水力参数沿垂直方向的变化，并假定沿水深方向的动水强度分布符合静水强度分布，因此所研究问题的水平尺度应远大于垂向尺度，流速等水力参数沿垂直方向的变化较之沿水平方向的变化要小得多 流动控制方程的分类 二维流动，控制方程在数学上的一般形式为 对于求解区域中的任一点 双曲型方程 ： 过该点有两条实的特征线 抛物型方程： 过该点有一条实的特征线 椭圆型方程： 过该点无实的特征线 不同类型微分方程在特性上的主要区别是它们的依赖区（domain of dependence）和影响区(domain of influence)不同 椭圆型方程 椭圆型方程描述的是稳态问题，这种问题的变量与时间无关，而需要在空间的一个闭区域上求解 任一点P的依赖区域是包围该点的求解区域边界的封闭曲线，P点的影响区域则是整个求解区域R 又称为边值问题（boundary value problem）。如：稳态传热、有回流的流动、对流换热 抛物型方程 抛物型方程描述的是步进问题，因变量与时间或类似于时间的变量有关 求解区域是一个开区间，计算时从已知的初值出发，称为步进法（marching method） 分布与瞬时以前的情况和边界条件相关，称时间步进，如1D非稳态导热；下游的分布仅与上游的变化相关，称主流步进，如2D稳态边界层流动 双曲型方程 P点的依赖区是上游位于特征线之间的区域，影响区是下游特征线之间的区域 双曲型方程数值求解也是一个步进过程。如：无粘性流体的稳态超音速问题，无粘性流体的非稳态问题 守恒与非守恒形式的控制方程 守恒与非守恒控制方程的形式 对流项写作散度形式 物理量都在微分符号内 守恒型方程（conservative form） 非守恒形式(non-conservative form) 将瞬态项和对流项中的物理量从微分符号中移出 最终形式 课后练习：推导非守恒型控制方程 守恒与非守恒型控制方程的区别 以能量方程为例 对空间任意，有限大小的容积V做积分 利用Gauss定律 含义：该体积内单位时间内能量的增加=由于流动进入该容积的能量+由于热传导进入该容积的能量+内热源的生成热，有限大小容积的能量守恒表达式 从微元体的角度看，控制方程的守恒型与非守恒型是等价的，都是物理守恒定律的数学表示。但对有限大小的计算体积，两个形式是有区别的。 几点说明： 不论节点布置的疏密程度如何，根据守恒型控制方程导出的离散方程也具有对任意大小容积守恒的特性。 非守恒型方程便于对离散方程进行理论分析，而守恒型控制方程能保持物理量守恒的性质，便于克服对流项非线性引起的问题。 第二章 流体流动的数学模式 流体流动的基本控制方程 三维直角坐标系及微元体 质量守恒方程（mass conservation equation） [单位时间内微元体中流体质量的增加]=[同一时间间隔内流入该微元体的净质量] 引入矢量符号 写成 瞬态三维可压流体的质量守恒方程，又称连续方程 若流动处于稳态，则密度 不随时间变化 若流体不可压，则密度 为常数 动量守恒方程（momentum conservation equation） [微元体中流体动量对时间的变化率]=[外界作用在该微元体上的各种力之和] 该定律实际上是牛顿第二定律 流体微元体上的压力 分子粘性作用而产生的微元体表面上的粘性应力 F 微元体上的体力 若只有重力 其它情况举例 对任何类型的流体（包括非牛顿流体）均成立 对于牛顿流体，粘性应力 与流体的变形率成比例 动力粘度（dynamic viscosity） 第二粘度（Second viscosity） 一般取 牛顿流体动量守恒方程展开形式 动量守恒方程简称动量方程（momentum equation） 也称运动方程，还称为Navier-Stokes方程 动量守恒方程的广义源项 表达式如下 一般情况下，