基于三角形网格求解二维浅水方程的改进的HLL方法

该文采用HLL格式,在三角形网格下采用有限体积离散建立了求解二维浅水方程的数值模型。基于三角形网格,底坡源项采用简单的斜底模型离散,为保证计算格式的和谐性,对经典的HLL格式计算的数值通量中的静水压力项进行了修正,并证明了修正后的HLL格式具备和谐性。算例证明本文提出的方法是和谐性的并具有良好的间断捕捉能力与稳定性。     

基于修正Roe格式的有限体积法求解二维浅水方程

应用迎风有限体积法,对具有复杂地形的二维浅水方程进行数值模拟.基于非结构化网格,采用Roe格式的近似Riemann解计算非粘性界面通量.底坡源项采用简单的斜底模型离散,从而保证了地形的离散精度.为了保证该计算格式的压力项和底坡源项的和谐性,对经典的Roe格式计算数值通量中的静水压力项进行了修正,并证明了修正后的Roe格式具备和谐性.通过与超临界流倾斜水跃和有激波混合流算例理论解结果比较,验证了该模型具有良好的间断捕捉能力、计算稳定性和守恒性.     

求解非平底一维浅水方程的 KFVS 格式

以Boltzmann方程为基础，本文提出了求解非平底浅水方程的KFVS（Kinetic Flux Vector Splitting)格式。因通量计算中考虑了底坡项的作用，计算方法和谐，即在非平底静水计算中始终能保持流速为零，水位为常数。数值试验表明，无论是计算恒定流还是非恒定流，该方法具有健全性和高分辨率等特点。     

求解二维浅水流动方程的Godunov格式

采用网格变换和Strang算子分裂,以准确Riemann解为基础。建立了求解非平底浅水流动方程的Godunov格式,用“水位方程法（Water Level Formulation)”求Riemann解,结合中心差分和Riemann解离散底坡项,保证了计算格式的和谐性,经验证,方法健全。通用,且分辨率高。     

非结构化网格下求解二维浅水方程的和谐Roe格式

为保证底坡源项和重力梯度项的平衡离散,采用把底坡源项分解为两不同部分单独处理的方法,在非结构化网格上建立了求解带复杂地形的二维浅水方程数值模型.采用Roe格式的计算界面通量,隐式求解摩擦力源项以增加格式的稳定性,并给出计算格式在非结构化网格上满足和谐性条件的证明.通过实例验证了此格式是和谐的,并具有良好的间断捕捉能力和稳定性.     

基于和谐性离散格式求解带源项的二维浅水方程

本文采用基于HLL近似Riemann解的Godunov格式,在三角形网格下,建立了二维浅水方程的有限体积数值离散模型.底坡源项采用与静压梯度项统一的矩阵散度处理方式,并证明了格式具备和谐性.通过算例验证了格式的有效性.     

带源项浅水方程的高阶格式研究

该文以Roe格式的近似Riemann解为基础,构造了具有二阶时空精度的MUSCL型格式来模拟非平底、变河宽浅水流动,并保证了计算结果的和谐性.通过几个经典算例的计算,比较了不同重构变量对不同流态的影响,验证了本文构造的高阶方法具有和谐、通用性好、分辨率高等优点.     

二维浅水方程的高阶松弛格式求解

利用松弛方法,将二维浅水方程转化为松弛方程组,并用逐维五阶WENO重构和显隐式Runge-Kutta方法对松弛方程组的空间和时间方向进行离散,建立了求解二维浅水方程的五阶松弛格式。WENO重构方法的引入既提高了格式的精度,又可保证格式是无振荡的。应用该格式对圆柱溃坝等问题进行了数值模拟,计算结果与用其它方法所得结果吻合,表明了方法的有效性。     

计算二维浅水方程中静水压力项与底坡项的积分平衡法

考虑静水压力项随水深的非线性变化,推导出计算二维浅水方程中静水压力项和底坡项的积分平衡法,以消除静水条件下的虚假流动。采用三角形网格的有限体积法模拟浅水流动,在以斜底模型逼近实际地形和控制体内变量线性分布的条件下,准确积分得出沿控制体界面的静水压力项和控制体上底坡项的表达式,从理论上证明了在静水条件下两者能够达到平衡。过丘恒定流算例计算结果表明,该算法能够使水流保持静水状态,不产生虚假流动。     

具有复杂计算域和地形的二维浅水流动数值模拟

本文采用非结构化网格的有限体积方法,对具有复杂计算域和地形的二维河道进行数值模拟。采用Roe格式的近似Riemann解计算界面通量,对地形变化源项采用特征分解平衡界面通量以保证格式的和谐性,对摩擦力源项采用隐式或半隐式求解以增加格式的稳定性。通过与混合流算例理论解比较,验证了格式是和谐性的并具有良好的间断捕捉能力与稳定性。应用此方法对地形复杂,存在多处险滩和深潭的某河段进行了实际模拟,数值计算结果和模型试验结果吻合较好。     

非结构网格浅水方程隐式解法

基于非结构网格模拟河道复杂边界的研究域,应用算子分裂法剖分浅水方程,采用有限体积法离散,建立隐式格式的通量计算,构造一种模拟二维水流运动的隐式算法。该方法具有计算稳定的特点,并应用于长江南京段分汊河道的流场模拟,计算结果与实测结果吻合良好,验证了算法的可行性和有效性。     

A TVD-WAF-based hybrid finite volume and finite difference scheme for nonlinearly dispersive wave equations

A total variation diminishing-weighted average flux(TVD-WAF)-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing equations were rewritten in the conservative form and then discretized on a uniform grid. The finite volume method was used to discretize the flux term while the remaining terms were approximated with the finite difference method. The second-order TVD-WAF method was employed in conjunction with the Harten-Lax-van Leer(HLL) Riemann solver to calculate the numerical flux, and the variables at the cell interface for the local Riemann problem were reconstructed via the fourthorder monotone upstream-centered scheme for conservation laws(MUSCL). The time marching scheme based on the third-order TVD RungeKutta method was used to obtain numerical solutions. The model was validated through a series of numerical tests, in which wave breaking and a moving shoreline were treated. The good agreement between the computed results, documented analytical solutions, and experimental data demonstrates the correct discretization of the governing equations and high accuracy of the proposed scheme, and also conforms the advantages of the proposed shock-capturing scheme for the enhanced version of the Boussinesq model, including the convenience in the treatment of wave breaking and moving shorelines and without the need for a numerical filter.     

Two-dimensional shallow water equations with porosity and their numerical scheme on unstructured grids

In this study, porosity was introduced into two-dimensional shallow water equations to reflect the effects of obstructions, leading to the modification of the expressions for the flux and source terms. An extra porosity source term appears in the momentum equation. The numerical model of the shallow water equations with porosity is presented with the finite volume method on unstructured grids and the modified Roe-type approximate Riemann solver. The source terms of the bed slope and porosity are both decomposed in the characteristic direction so that the numerical scheme can exactly satisfy the conservative property. The present model was tested with a dam break with discontinuous porosity and a flash flood in the Toce River Valley. The results show that the model can simulate the influence of obstructions, and the numerical scheme can maintain the flux balance at the interface with high efficiency and resolution.     

移动边界的有限体积KFVS方法在一维溃坝波计算中的应用

本文从Boltzmaun方程出发,采用有限体积KFVS（Kinetic flux-vector splitting）方法,对微观粒子的速度分布进行组合,给出了模拟河道洪水演进时干湿交界处移动边界的一种新方法。并成功地模拟了大坝下游为干河道的一维溃坝模型,通过和精确解的比较得出,此方法在处理河道移动边界问题上是成功和有效的。     

基于FD/MLSFD的混合算法计算 二维浅水方程

将基于最小二乘法的无网格有限差分方法(MLSFD方法)和传统有限差分方法(FDM)相结 合,生成了一种混合算法。研究并解决了差分网格/无网格混合算法的藕合问题,将混合算法运用于 求解二维浅水方程,通过计算高雷诺数(Re =105)条件下的圆柱绕流模型揭示了该混合算法在进行二 维浅水模型的数值模拟方面具有较高的精确度。     

2-D FLOW AND SEDIMENT SIMULATION FOR THE FLOOD REGULATION OF A RESERVOIR WITH WATER INTAKE

In this paper, the flood regulation by operating the downstream sluice gates for a reservoir with a water intake is studied. The two-dimensional depth-averaged flow equations are solved by the boundary fitted finite volume method (FVM) based on MacCormsck prediction-correction scheme. The bed deformation caused by both the bed load and incoming suspended sediment is determined in a coupled way. The model is used to simulate the practical flood regulation operation of a reservoir. The results have been compared with the physical experiment.