浅水方程推导word版本.docVIP

浅水方程推导 将三维的基本方程沿水深积分平均，即可得到沿水深平均的平面二维流动基本方程。 定义水深为， 、为基准面下液面水位和河床高程 定义沿水深平均流速为： 引用莱布尼兹公式 自由表面及底部运动学条件 以x方向为例三维流动的运动方程沿水深平均为 非恒定项积分 对流项积分 首先将时均流速分解为，式中为垂线平均流速，为时均流速与垂线平均流速的差值。 式中，，是由于流速沿垂线分布不均匀而引入的修正系数，类似于水力学中的动量修正系数，其数值一般在1.02—1.05，可以近似取1.0，因此 类似，可以得到 上几式相加，并利用底部及自由表面运动学条件可得 压力项积分 （莱布尼茨公式） 将代入上式后化简得： 扩散项积分 上式右边后两项分别为由底部创面阻力和表面风阻力引起的阻力项。式中，为无因次风应力系数；为空气密度；为风速；为风向与x方向的夹角。 最后运动方程写成张量形式为 2. 差分格式的稳定性与收敛性 两种误差 舍入误差：=有限精度的计算机上的解差分方程的精确解（稳定性） 离散误差=偏微分方程的精确解差分方程的精确解（收敛性） 稳定性分析 误差也满足差分方程。 求解稳定，则要求 Von Neumann（冯.诺依曼）稳定性分析 将某一时刻分布在网格点上的误差按Fourier级数展开，然后考察下一时刻各网格点上误差的Fourier分量是衰减还是增长，以判断差分方程是否稳定。 用一维热传导方程