水流作用下地形不稳定性计算模型 Numerical model of morphology instability under flow.pdfVIP

第32卷第6期 水道 港 口 V01．32No．6 2011年12月 Journalof andHarbor Dee．2011 Waterway 水流作用下地形不稳定性计算模型 王明明，邹志利 (大连理工大学海岸和近海工程国家重点实验室，大连116024) 摘要：河流中水流和水底泥沙之间的相互作用会产生各种复杂的地形形态，交替型沙坝是其中的一种 地貌特征，研究表明这一地貌的形成机理是水底地形的不稳定性。给出了这一不稳定性所导致的交替型 沙坝的数值模拟方法和相应结果。水流方程采用了交替方向隐式算法求鳃，水底变形方程采用了比传统 方法更加稳定的欧拉加权基本无振荡迎风格式求解，5l而建立了水流与水底地形变化耦合计算模型。运用 该模型研究了初始地形扰动产生的交替沙坝在长直河道中的演变。 关键词：地形演变；非线性；WENO；不稳定性 242．1 中图分类号：P731．23；0 文献标识码：A 文章编号：1005—8443(2011)06—0389—10 近岸区域和河流中地形演变依赖于水流、波浪、泥沙之间复杂的相互作用，水流作用下地形的不稳定性 是这一复杂相互作用的突出表现形式。这一不稳定性导致水底地形即使是在初始时刻是均匀平直的水流作 用下也会产生沿空间和时间变化的地貌特征，交替型沙坝就是这种复杂地貌特征的常见典型例子。在求解这 种地形和水动力相互作用所引起的水底不稳定性问题时，如何求解地形演变方程是地形变化数值模拟中的 一个难题，因为该方程是非线性的，会导致间断解的出现。以往多采用Lax．Wendroff格式…或者其改进后的 格式，这些格式在地形峰值后面会产生数值振荡，从而不能够模拟地形的长时问演变过程。Hundson[2]对比 二维模型。加权基本无振荡(WeightedEssentially 4]格式具有在间断区分辨率高等优点。为了进一 (EssentiallyNon．Oscillatory)构造思想的基础上提出的。ENOI 步提高该格式的计算精度，WENO格式引入了变化的加权因子，使计算结果更加稳定[5]。文章研究了地形和 水动力相互作用所引起的水底地形不稳定性中WENO格式的应用问题。对水流方程仍然采用了常用的ADI 方法…6，但对地形演变方程采用了Euler．WENO方法，在此基础上建立了水流与水底地形变化耦合计算模 y X 图1坐标系统 Coordinate Fig．1 system 收稿日期：2011-05—09；修回日期：201l一06—08 基金项目：国家自然科学基金项目；辽宁省教育厅高等学校自然科学基金 作者简介：王明明(1986一)，男，河南省禹州人，硕士研究生，主事海岸及近海工程的研究。 student． Biography：WANGMing-ming(1986-)，male，master 万方数据 390 水道港 口 第32卷第6期 型，并与传统的FFCS离散格式建立的模型的计算结果进行对比。应用该模型对交替沙坝地形的二维长直河 道的形成进行了数值模拟，给出了这种地形形成的不稳定性机理的数值模拟结果。 1水流和水底地形的耦合数学模型 考虑一有限宽度的长直河道，取菇轴和Y轴分别沿河道的长度方向和宽度方向。水底沿x方向有一个微 小的倾角i。，初始水流沿整个河道是均匀的。水流和地形变化满足以下基本方程。