MacCormack法应用于输沙计算研究.docVIP

MacCormack法应用于输沙计算研究 　　摘 要：本研究探讨将MacCormack方法应用于河流的推移质计算，为流域输沙计算提供科学的数值方法。研究基于MacCormack法的差分方案，将水流计算分为预测计算和修正阶段，而后将推移质输沙计算公式与MacCormack计算方案耦合，并应用于设定条件的事例计算。计算结果形象地再现了河流在改变流量条件下输沙量随时间的变化，该方法适用于计算常流和变流量条件下的河道输沙量、也可用于河床演变，河道疏浚的计算。 　　关键词：MacCormack法，有限差分，推移质，输沙 　　中图分类号：TV148文献标识码： A 文章编号： 　　1前言 　　 流域的水沙过程包括降雨产流、汇流的全过程, 以及泥沙侵蚀、输移、沉积的不同运动形态. 对这一物理过程的全面描述有其必要性和现实意义[1]。输水的同时输沙, 对河流起着泄洪排沙、维持河道正常演变的作用[2]。泥沙运动既应满足一定的力学规律, 同时也是一种随机现象[3]。输沙计算，一般先求解水流条件, 得到水深、流速、流量等水力因子, 将水力因子代入输沙方程中, 得到含沙量或输沙率过程[4]。 　　 在现代河流泥沙物理模型及数学模型中,经常遇见输沙计算问题, 但由于河流动力学所面临的自然现象, 是边界复杂多变的天然河流及水流泥沙运动的多相体, 许多问题的运动规律尚无定论, 特别是水流挟沙力的计算, 人们还不能全面地把握它的力学机制, 因此, 多少年来, 泥沙界的学者们对水力输沙特性在不懈地探索[5]。崔侠等（1987）从冲淤过程的机理出发, 推导出了二度恒定非均匀流中泥沙扩敬方程, 求出了含沙量沿程变化规律, 着重讨论了水流挟沙力与含沙量的关系[6]。张洪武等（1992）从水流能量消耗和泥沙悬浮力之间的关系出发, 考虑了泥沙存在对卡门常数和泥沙沉速等的影响, 给出了半经验半理论的水流挟沙力公式并在黄河的输沙计算中应用[7]。胡海明等（1996）较为系统地研究了非均匀沙运动机理及输沙率的计算方法, 在建立泥沙运动交换模式的基础上, 指出现有非均匀泥沙起动流速公式存在的不足并对此作了修正[8]。黄永健等（1997）结合灌渠的结构特点及水流泥沙特性, 对一维非耦合、非均匀、不平衡输沙基本方程组进行了合理的简化, 根据实体模型试验及实测资料, 确定了自流与提灌时干、支渠分流的不同分沙分水比并应用于黄河下游的输沙计算[9]。迄今为止的输沙计算方法普遍存在两个不足，一是经验型公式（模型）较多，其应用范围常被限制在研究开展的流域或河段；二是计算精度普遍较低，对输沙过程的模拟精度不高。本研究探讨将MacCormack方法用于推移质的输沙计算，为河流输沙，河床演变以及河道疏浚等提供科学的计算方法。 　　2 MacCormack方法 　　 MacCormack方法由Robert W. MacCormack在1969年提出，是计算流体力学中用来求解双曲形偏微分方程数值解的一个普遍离散方法，为二阶有限差分方法。该计算方法非常简洁，易于理解和实现编成计算，应用MacCormack方法计算时包括两个阶段：预测阶段和修正阶段，该方法非常适合线性方程的求解。由于采用二阶差分对预测阶段和修正阶段进行计算，数值计算结果的精度较高。对于非线性方程，应用MacCormack方法也可得到较好的结果。 　　3 算法建立 　　3.1 基础方程式 　　 用于径流输沙过程计算的基础方程式包括地表径流连续方程式，运动方程式以及输沙量计算公式[10]。 　　·地表径流连续方程式（一维稳定流） 　　 （1） 　　·运动方程式 　　（2） 　　·输沙量计算公式（Ashida / Michiwue 经验公式） 　　（3） 　　式中，Δt为计算的时间步长，1s；Δx为河流纵向计算的单位步长，50 m；b为河道宽度，m；Q为流量；g为重力加速度，m/s2；h为水深，m；ib为河道平均坡度；ie为水力坡度；qbi为单位时间内单宽河道输沙量，m3/s；s为沙在水中的比重，kg/m3；d为流沙的平均粒径，0.5 mm，为推移质；τ*c为界限推移力，kg；τ*为推移力，kg；u*c为界限摩擦速度，m/s；u*为摩擦速度，m/s。 　　3.2 有限差分 　　 计算时，需要依据地表径流连续方程式及运动方程式耦合输沙计算在时间上离散化，即有限差分。因为计算依赖于一定的初始条件和边界条件，所以预测阶段采用后退差分法，修正阶段采用中间差分法[11,12]。 　　·预测阶段差分式 　　 （4） 　　·修正阶段差分式 　　 （5） 　　式中，i为计算的栅格编号，n为计算的时间编号，A为过流面积，m2; kv为人工粘性系数，0.00001；其它因子同上。 　　3.3 计算程序编译