浅水方程推导.doc_第1页
浅水方程推导.doc_第2页
浅水方程推导.doc_第3页
浅水方程推导.doc_第4页
浅水方程推导.doc_第5页
已阅读5页，还剩1页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1. 浅水方程推导将三维的基本方程沿水深积分平均，即可得到沿水深平均的平面二维流动基本方程。定义水深为， 、为基准面下液面水位和河床高程定义沿水深平均流速为：引用莱布尼兹公式自由表面及底部运动学条件以x方向为例三维流动的运动方程沿水深平均为非恒定项积分对流项积分首先将时均流速分解为，式中为垂线平均流速，为时均流速与垂线平均流速的差值。式中，是由于流速沿垂线分布不均匀而引入的修正系数，类似于水力学中的动量修正系数，其数值一般在1.021.05，可以近似取1.0，因此类似，可以得到上几式相加，并利用底部及自由表面运动学条件可得压力项积分（莱布尼茨公式）将代入上式后化简得：扩散项积分上式右边后两项分别为由底部创面阻力和表面风阻力引起的阻力项。式中，为无因次风应力系数；为空气密度；为风速；为风向与x方向的夹角。最后运动方程写成张量形式为2. 差分格式的稳定性与收敛性两种误差舍入误差：=有限精度的计算机上的解差分方程的精确解（稳定性）离散误差=偏微分方程的精确解差分方程的精确解（收敛性）（1） 稳定性分析误差也满足差分方程。求解稳定，则要求 Von Neumann（冯.诺依曼）稳定性分析将某一时刻分布在网格点上的误差按Fourier级数展开，然后考察下一时刻各网格点上误差的Fourier分量是衰减还是增长，以判断差分方程是否稳定。用一维热传导方程作为模型方程差分显式格式：假设误差随时间按指数函数的方式增长或衰减，即随时间按指数函数变化。将其代入差分方程得： 称为放大因子解不等式，即，得： （2） 收敛性分析收敛性是指当网格点空间趋于零时，差分方程的解无限接近于偏微分方程的解。可以证明，当网格变细，并且时，有限差分方程的解收敛于给定的扩散方程的精确解。由于收敛性问题的讨论在理论上证明较为困难（微分方程较复杂时），关于收敛性问题比稳定性问题复杂得多，所以在此给出更具实用意义的定理，Lax等价定理，即，对于一个与线性偏微分方程相容的适定的初值问题的差分格式，稳定性是差分方程解收敛于微分方程的充分必要条件。也就是说对一适定的线性初值问题，相容性加稳定性等价于收敛性3. 多维问题的常用差分格式以二维扩散方程为例（1）交替方向隐式格式(ADI)基本思想是将差分计算分成两步：第一步在一个方向是隐式的，而在另一个方向上是显式的；第二步则是两个方向交换一下，即在第一个方向上为显式，而在第二个方向上为隐式。由于只在一个方向上隐式，求解时形成的方程组是三对角方程组，所以求解大为简化。（2）时间分裂格式基本
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论