欧拉2拉格朗日方法在三维水流模型中的时间阻力ρ0∫H-水科学进展.PDF

第 20 卷 第 6 期 水 科 学 进 展 Vol20 ,No6 2009 年 11 月 　　 ADVANCES IN WATER SCIENCE 　 Nov. , 2009 　 欧拉拉格朗日方法在三维水流模型中的时间阻力 1 ,2 3 1 1 胡德超 , 赵维阳 , 钟德钰 , 张红武 ( 1 清华大学水沙科学与水利水电工程国家重点实验室 , 北京 　100084 ; 2 长江科学院河流研究所 , 湖北 武汉 　430010 ; 3 长江航道规划设计研究院 , 湖北 武汉 　430011) ( ) 摘要 : 在三维数模的实践中 , 发现使用欧拉拉格朗日方法 ELM 离散对流项会引发计算水位随时间步长减小而增 大的异常现象 , 本文结合数值实验和数学分析对其进行了阐明和研究 。数学推导分析表明 , 这些异常是由于 ELM 中的线性插值与水平流速的非线性垂线分布之间的不适应所致 , 插值误差随着插值点离单元侧面距离的减小而增 大 。在此基础上 , 对异常现象的本质进行了解释 。此外 , 通过使用垂向 σ坐标网格上的三维水动力学模型进行计 算 , 得出 ELM 的时间阻力现象具有普遍性 , 但对边界适应能力强的网格系统可以减轻其影响 。 关 　键 　词 : 水力学 ; 数值方法 ; 三维 ; 欧拉拉格朗日; 时间阻力 ( ) 中图分类号 : TV1332 　　　文献标识码 : A 　　　文章编号 : 2009 ( ) 欧拉拉格朗日方法 ELM 兼有Lagrangian 方法的稳定性和 Eulerian 方法的固定网格 , 其 自然迎风性可消除 欧拉方法固有的 Courant 数稳定限制 , 是一种无条件稳定的显式对流项离散格式[ 12 ] , 在 目前水动力学模型领域 获得了越来越多的应用[37 ] 。另外 , ELM 具有非守恒的本质 , 前人对此也提出一些改进的方法[8 ] , 但这些工作 大多局限于一维 、二维模型 。在三维模型领域 , Dimou [3 ] 曾在数值实验中观察到 ELM 的计算误差随着时间步长 减小而增大的现象 , 验证 了 Baptista[ 1] 在一维模型中用傅立叶分析得 出的结论 , 但没有给 出分析和解释 ; J ankowski [4 ] 曾提到 ELM 具有较强的数值粘性 , 但也未给出具体说明; 文献[ 57 ] 则仅在模型中应用了 ELM 。由 此可见 , 在三维模型领域 , 由于问题的复杂性 , ELM 特性的相关研究还