坡面降雨入渗产流规律的数值模拟研究2行政.docVIP

坡面降雨入渗产流规律的数值模拟研究2行政论文范文大全 坡面降雨入渗产流规律的数值模拟研究 3 模型求解 　　模型的求解也分为两部分。运动波方程是一个非线性的对流型方程，求解采用了简单的一阶显式迎风格式。经过比较，一阶显式迎风格式对此问题的计算结果在幅值、相位、守恒性几个方面的综合效果较其他一些格式为优。我们曾使用过preissmann四点偏心格式，它在计算此类坡面薄层水流运动波方程时耗散过大，表现为达到平衡产流的时间偏快。我们也曾试过中心差分格式和高阶迎风格式，结果发现标准的中心差分格式仍然和在普通的线性对流型方程中的表现一样是绝对不稳定的。一些变形的中心差分格式能得到较好的结果，但在守恒性和稳定性上不如一阶迎风格式。高阶迎风格式的结果也与此类似，其中二阶迎风格式实际计算时稳定的currant数比一阶迎风小一半。除此之外，一阶迎风格式的计算是最为简单的，边界条件也易于处理。 　　将式(1)第一式写为 (8) 其中f代表通量，f=f(h)=q=1/nh5/3s01/2,sr为源项，sr=pcosθ-i。 　　计算格式形式按有限体积法写为 hin+1-hni/δt+fi+1/2-fi-1/2/δx=srni (9) 　　一阶迎风格式的通量写为 fi+1/2(h)=1/hhi5/3s01/2 (10) 　　[