平面二维数学模型在变率研究中的应用.pdfVIP

十、河流模拟技术 。1481。 平面二维数学模型在变率研究中的应用 陈凤玉1昊祥华2 f1．长江科学院 武汉430010； 2．长江水利委员会设计院 武汉430010) 摘要：本文在有限元算法基础上，建立了平面二维有限元水沙数学模型。该模型对控制方程采用显式 有限单元法进行离散，运用二阶预测一校正多步法、龙格一库塔四阶单步法对有限元方程进行求解。推 导了由于实体模型变态引起平面=维水沙数学方程的变化，并在一突然放宽矩形断面河道中，利用本文 所建立的数学模型研究了不同变率对回流区水流结构与泥沙淤积的影响。 关键词：平面二维 水流泥沙 有限元法 变率研究 1前言 目前在泥沙工程领域研究具体工程问题时，特别在预测将来可能出现的情况时，主 要运用实体模型和数学模型两种研究手段。由于水流泥沙运动的复杂性，在建立数学模 型时必须要做若干假设和简化，使得计算结果与实际情况不可能完全一致。实体模型的 优势是可以显示出三维条件下水流和泥沙的运动特性。但实体模型还存在许多问题尚待 解决，如变态模型的适用范围和变率限制，成为长期探讨的重要问题之一，以往多用概 化实体模型进行研究，取得不少研究成果，得到定性上的认识。但针对具体的工程问题， 其河道地形和边界条件都极为复杂，许多客观因素都不允许去研究多种变率下的模型试 验成果。基于此点，本文试图采用平面二维水沙数学模型对某些变率情况作初步探讨。 回流和回流区的泥沙淤积问题是河流工程和水利工程规划、设计中常遇到的疑难问题之 一。如三峡水库变动回水区河段河道宽窄相间、浅滩交错，存在许多港池回流区；盲肠 河段口门回流区的泥沙淤积；两江汇合和分汊河道汇流区也存在大小及强度不等的回流 区。工程上通常采用变态物理模型试验进行研究。本文选用单边突然展宽河道，进行变 态数学模型的模拟计算。 2平面二维有限元水沙数学模型 2．1 基本方程 平面二维水沙数学模型是以垂线平均的水流和泥沙因素为研究对象，给出它们在平 面上的分布及变化情况。基于控制体概念并遵循质量守恒定律和牛顿第二定律，平面二 作者简介：陈凤玉，女，1968年出生，高级工程师。 ．1482． 第六届全国泥沙基本理论研究学术讨论会论文集 维水流泥沙数学模型由以下方程控制： 水流连续方程 一0Zs+型塑+—O(H—V)：0 (1) at ax 却 水流运动方程 c2， 西 叙 办 。 p 一舐2 咖∥ 詈+u豢+y券一g竺ax一号¨c害+等， iOV+【，娑+y娑：一g堡一2M(害+等) O)。x3( 西 缸 加 P 一叙2 却∥ (rx,Vy)=等拓2+矿2(u，y) (4) 式中：C为谢才系数。 悬移质泥沙扩散方程(泥沙连续方程) 瓦OS+嘻+矿等=Ss(窘+害)一詈cs一剐 c5， 河床变形方程 ∥竿=aco(S—S。) (6) 以上各方程式中，X、Y为笛卡尔纵、横坐标分量；U、V分别为垂线平均流速在x、Y方 向上的分量；磊、Z6分别为水位和河底高程；f为时间变量；日为相对某一基准面的垂线 水深，H=zs—z6；p为水密度；M、岛分别为紊流黏滞性系数和泥沙紊动扩散系数；炙、 元分别为底部切应力在x、Y方向上的分量；S、＆分别为垂线平均含沙量和挟沙力；P’、 口分别为泥沙千密度和恢复饱和系数；国为泥沙的静水沉速。 2．2有限元离散模式的建立 对平面二维水流泥沙控制方程采用Galerkin加权余量法逼