随着大容量、高速度的计算机发展,人们采用数值模型来研究河口水动力及泥沙输运问题越来越普遍,研究对象也从一维问题发展到二维问题以及三维问题。非结构网格有限体积法(FVM)物理意义清晰,能处理复杂的计算区域,逐渐被应用在水力学及河流动力学领域中。本文在广泛研究现有各类水沙数值模型基础之上,采用有限体积法(FVM)对方程进行离散,构建基于非结构网格在笛卡尔坐标系统下平面二维水沙数值模型。针对天然河道岸线的曲折和出水岛屿的特点,传统笛卡尔坐标系统下的矩形网格在边界上模拟常常达不到很高的精度,本文以阵面推进法(Advancing Front)为基础,采用节点间距函数来控制网格间距,形成了一套离散任意平面区域的三角形非结构网格自动生成方法,在边界上模拟能达到很高精度。本文中有限体积法,是基于Casulli[4](1992)隐式-显式分布模型。在对方程离散过程中,引入弥散修正项,成功地实现了重力项与运动项统一的格心格式。在对流项处理上,借鉴计算空气动力学的成果,采用TVD格式为基础的高阶对流格式,具有很高的计算精度(二阶精度),能控制间断解处非物理的震荡现象。数值模型计算程序采用Visual Fortran语言编制。最后利用一维驻波理论算例及前人做的床面演变实验资料对二维水沙数值模型进行了率定和验证,并且对长江口水动力进行了模拟,数值模拟结果与实测结果吻合的较好,模拟精度比较高,基本上能反映水沙运动的基本规律,模型是可靠的。