论文部分内容阅读
大坝溃决造成的洪水、以破坏方式开启非常溢洪道造成的洪水以及为保证大坝安全加大水库放流形成的洪水,传播速度快,破坏性巨大。水坝的安全问题构成了对人类生命和财产的潜在威胁。据不完全统计,在二十世纪,坝高在15m以上的溃坝事件就有200多起,造成了巨大的经济和生命损失。溃坝问题因此备受关注。溃坝引起的水流运动情况十分复杂,受到许多因素的影响,早期的研究都是通过实验来进行。随着计算技术和数值方法的迅速进步,越来越多的研究者对溃坝问题用数值计算来模拟。由于数学上浅水长波方程与可压缩气体的无粘性流方程类似,可以把浅水长波方程比拟成可压缩气体无粘性流方程,将空气动力学中所得到的许多数学概念和方法(包括数值方法)直接应用于浅水长波的运动。对于双曲型的拟线性偏微分方程组,可以在特征线上,构造相应的一组相容方程来求解,这就是特征线方法。特征线方法的物理图像清晰,数值计算精度高,在某些特殊情况下甚至可以得到封闭的分析解。本文介绍了特征线法的基本原理,推导了用以代替原浅水方程的特征线方程和相容性方程,结合溃坝水流下边界条件难以确定的特点,构造了用于溃坝计算的显格式。本文模拟了一维矩形断面低坝和梯形断面巨型坝瞬间全溃的洪水波演进过程,给出了与分析解的比较,并结合算子分裂的特征线法对二维瞬间部分溃坝时洪水波的演进过程进行了数值模拟。本文得到了激波捕捉能力强的数值模拟解,展示了溃坝波的运动特性,为溃坝引起的灾害分析及防灾减灾等工程问题提供科学依据。结果表明这种方法精度高,有较强的溃坝涌波捕捉能力,是一种有效方法。