悬沙运动的快速特征有限元模拟方法分析讲解.docVIP

悬沙运动的快速特征有限元模型 柳海涛1 王小松1 夏庆福1 中国水利水电科学研究院，北京，100038，htliou@163.com) 摘要：本文在已有的研究基础上，发展了一种快捷的水沙特征有限元计算模型，来模拟以对流作用为主的悬沙运动及河床变形。该模型采用一种预测-校正方法来对特征点进行搜索，提高了特征含沙量的计算精度，对粘土质河床的处理引入粘性土冲刷平衡流速作为判别指标。另外对长系列计算下时间步长的取值、含沙量的数值耗散、大型矩阵的存储方法以及动边界处理也给出了具体方法。通过典型算例表明，模型具有较好的适应性和可靠性。 关键词：特征有限元模型；特征点；对流为主的悬沙运动；冲淤平衡流速 1 前言 当前，二维泥沙数学模型在基本理论和计算方法上已相当成熟，国内外出现了许多模块化和商业化的计算软件，但这些软件在使用时一是没有源代码，用户无法根据实际需要进行修改；二是其采用的数值方法在一些特殊的情况下并不适用。因此，为解决实际问题，在现有的模型的基础上进行二次开发是一种必要的手段，在文献[1]中，作者就提出了一种简洁的求解悬移质泥沙运动的特征有限元方法，该方法在求解精度和速度上基本可以胜任河道泥沙计算，其理论精度较以往的集中质量法以及许多常用的隐格式都要高。然而，该方法在泥沙特征点的搜索、单元干湿判别、粘性土质的冲刷、大型稀疏矩阵的存贮和求解方法等方面的研究还不够深入。本文在已有水流模型的基础上，发展了一种基于特征有限元方法的水沙计算模型，并对含沙量特征点的搜索方法、长系列水沙计算时间步长的取值、粘性土质和局部刚性河床的处理等方面进行研究，同时对模型中含沙量的数值耗散、单元的干湿变化的处理、大型矩阵的存储技术和计算方法等问题进行改进，使其更加实用。最后通过典型算例来验证该模型的适应性和可靠性。 基本方程 本文模型的水流控制方程如下 (1) (2) (3) 其中，为流速（m/s）；和为水深和床面高程(m)；为综合紊动粘滞系数(N·s/m2)；为谢才系数，。 泥沙输移控制方程 (4) 河床变形方程 (5) 其中为泥沙扩散系数(m2/s)，为淤积物干容重(kg/m3)，为颗粒沉速(m/s)，和分别为垂线平均含沙量和水流挟沙力(kg/m3)，为恢复饱和系数，为冲刷判别系数，其作用将在后面介绍。 根据两相流理论，对于泥沙输移方程，因其忽略惯性力、升力、泥沙应力梯度等的作用，扩散仅适用于含细颗粒的低浓度挟沙水流驰豫时间水流微元自身的拉格朗日积分时间尺度颗粒能充分跟随水流紊动 限于作者的水平，上述问题还有待下一步研究，本文模型暂时采用如下假定： (1) 悬移质为较均质细沙，床沙质为粘性土或相同的均质细沙。 (2) 挟沙力采用一维水流挟沙力公式：，其中和取常数。 (3) 泥沙扩散系数取常数；恢复饱和系数在冲刷时为1.0，淤积时为0.25。 模型的数值方法 3.1 特征有限元方法 对于对流占优问题，Galerkin法等价于中心差分格式，因缺乏足够的耗散，往往导致数值振荡采用高分辨率格式对流项进行重构，即通过引入几乎相等的扩散与反扩散以保证格式的高精度，同时利用限制因子保证影响系数的非负性解的保单调性质量集中方法预报-校正-迭代 （6） 其中为时刻的特征点K处的含沙量值，K点的坐标为 (7) 这样泥沙输运方程作如下离散 (8) 将式(6)代入上式，经整理得到如下有限元离散方程 (9) 系数矩阵 (10) 右端向量 (11) 基函数 单元节点含沙量 单元节点的特征点处的含沙量 单元节点的水流挟沙力 这样，对每一个单元e求出单元系数矩阵A和右端向量B，然后将其合成总体矩阵，即可求得时刻的含沙量S。 对于河床高程的修正，仍采用式(5)中的河床变形方程，在求解Z之前，先应用如下的隐式差分格式：