论文部分内容阅读
Savage-Hutter( SH)方程是描述Coulomb颗粒材料在地表面上流动的薄层近似方程，是滑坡、泥石流等地质灾害研究中的一个重要的力学模型。本文对SH方程的一维和二维形式，利用Godunov有限体积法进行积分离散，用一阶常值重构或者二阶MUSCL重构网格界面两侧的守恒变量，并采用Lax-Friedrichs数值通量逼近Riemann问题的通量，或者采用无需解Riemann问题的中心型格式(NOC)。作为初步研究，源项采用显式处理;在干湿边界处理上，为克服虚假速度问题，初始时刻在没有颗粒的地方也敷设一层很薄的颗粒，计算过程中限制颗粒流深度不小于此值。边界条件采用基于Riemann变量的特征边界条件。