求解具有复杂地形二维浅水方程的修正HLL格式

采用非结构网格的有限体积方法,对具有复杂地形的二维浅水方程进行数值模拟.采用HLL近似Riemann解计算界面数值通量,基于三角形网格,底坡源项采用简单的斜底模型离散,保证了地形的离散精度,摩阻源项采用全隐方式求解以保证格式的稳定性.采用多维重构及多维限制器的方法获得空间二阶精度的格式,时间离散采用三阶Runge-Kutta法以获得高阶时间精度.为保证格式的和谐性,对经典的HLL格式计算的数值通量中的静水压力项进行了修正.数值计算的结果验证此格式具有良好的高精度捕捉间断的能力,可以应用到地形复杂的二维浅水问题计算中去.     

和谐WAF格式在带干河床浅水波方程中的应用

采用二阶精度和谐的加权平均通量（WAF）格式对浅水波方程中的间断解和干湿边界问题进行了研究．用基于HLL求解器的WAF格式近似数值通量,中心差分格式离散地形源项,在干湿边界处重新定义数值河床,保证格式是和谐和扩展和谐的．最后模拟了带激波和干湿边界的浅水问题,精确地捕捉到了激波和干湿边界面,验证了该格式具有高分辨率、无振荡以及正确模拟干湿边界的能力．     

澳门周围水域潮汐流动的数值模拟

采用有限体积法求解二维浅水流动方程,建立了一种在非规则结构化网格上数值求解平面二维浅水流动的数值方法,其中应用Roe格式计算数值对流通量,并通过MUSCL途径插值界面状态,以对水位插值代替水深插值来消除地形变化对格式的影响;对源项采用分步法及半隐式求解,防止了源项引起的计算不稳定;处理动边界采用人工小水深的方法.应用于澳门水域潮流的模拟,计算结果与实测吻合较好.     

澳门周围水域潮汐流动的数值模拟

采用有限体积法求解二维浅水流动方程,建立了一种在非规则结构化网格上数值求解平面二维浅水流动的数值方法,其中应用Roe格式计算数值对流通量,并通过MUSCL途径插值界面状态,以对水位插值代替水深插值来消除地形变化对格式的影响;对源项采用分步法及半隐式求解,防止了源项引起的计算不稳定;处理动边界采用人工小水深的方法．应用于澳门水域潮流的模拟,计算结果与实测吻合较好．     

高精度格式溃坝波模拟

采用高精度WENO和MWNEO格式对二维浅水方程变量进行数值重构,研究分析高精度格式捕捉溃坝波的能力和分辨率.对浅水方程,采用Roe-Riemann格式和波传算法进行数值离散,并与以上2种重构方法两两结合成4种计算方法.采用Fortran语言分别对其进行编程计算.计算结果显示,高精度格式在捕捉溃坝瞬时基波延续推进的过程较好.在4种计算方法计算结果中,MWENO型波传格式模拟结果最好,捕捉溃坝波能力强,分辨率很高.     

复杂地形上坝堤溃决洪水演进的非结构有限体积模型

针对复杂地形上坝堤溃决洪水演进的二维数值模拟,综合考虑底坡项、摩阻项和干湿界面对水流的影响,基于浅水二维水动力学理论建立了非结构有限体积模型.围绕传统底坡项处理技术应用于复杂地形时的不足,通过计算底坡项对物质通量和动量通量的影响,并进行相应的分解,提出了集成底坡项的近似Riemann求解器.分析并采用改进的隐式格式解决了摩阻项可能引起的刚性问题.通过求解干湿单元界面处"镜像湿单元"Riemann问题,结合变量重构技术,有效处理了动边界问题.算例研究结果表明,模型稳定健全,和谐性好,具有较好的推广应用价值.     

高精度WENO和MWENO格式溃坝波模拟

采用高精度WENO和MWENO格式对二维浅水方程变量进行数值重构,研究分析高精度格式捕捉溃坝波的能力和分辨率.对浅水方程,采用Roe-Riemann格式和波传算法进行数值离散,并与以上2种重构方法两两结合成4种计算方法.采用Fortran语言分别对其进行编程计算.计算结果显示,高精度格式在捕捉溃坝瞬时基波延续推进的过程较好.在4种计算方法计算结果中,MWENO型波传格式模拟结果最好,捕捉溃坝波能力强,分辨率很高.     

求解带源项双曲守恒律方程的IGVC格式

研究带源项双曲守恒律方程的IGVC(改进的群速度控制)格式。为了控制数值方法带来的非物理振荡,引入IGVC格式,考察其TV性质,对于源项我们采用分裂方法处理,该方法用于求解一维浅水方程,数值结果是令人满意的。     

用动力学的流矢量分裂法计算溃坝波问题

采用动力学的流矢量分裂法模拟浅水长波方程. 利用气体动力学方程与浅水长波方程 之间的比拟关系, 并结合TVD差分格式给出浅水长波方程的数值计算方法. 计算了溃坝波问 题. 数值结果与精确解的比较表明, 该方法能很好地模拟浅水长波方程中的间断解问题.     

非定常对流扩散方程的高阶差分格式

对流扩散方程主要包含对流和扩散两项。在数值计算中,方程中的扩散项一般采用具有优良物理特性和计算精度的中心差分离散格式,而对对流项的处理就稍显困难,处理不当便会产生数值震荡或数值弥散,给数值计算带来困难。针对对流扩散方程,通过引入指数变换将对流扩散方程变为扩散方程,避免对对流项的直接处理。利用四阶紧致差分格式,首先建立三类特殊方程的高精度差分格式,在此基础上建立一维非定常含源对流扩散方程的高阶格式,并进行稳定性分析,所得格式精度高且绝对稳定。数值算例表明了该格式的有效性。     

浅水方程Roe型格式的平衡性

采用类似对流项特征向量分解形式分裂底坡源项,证明了浅水方程中底坡源项和压力项的平衡性,有效避免了多变地形下伪流现象的产生.运用该模式对水下地形多变、岛屿星罗棋布的浙江省温岭龙门海域进行了潮流运动模拟,模拟计算结果与水文测验数据拟合程度较好.     

GAMMA格式及DATE方法在SIMPLE算法中的应用

为了提高N-S方程数值计算的精度并解决压力-速度的耦合问题,将GAMMA格式、DATE方法与SIMPLE算法相结合,得到了一种高精度非结构同位网格下的SIMPLE算法.该算法用GAMMA格式离散对流项,用DATE方法处理方程求解过程中压力-速度的耦合问题,可以弥补传统动量插值法不能处理强非线性压力场的不足,保证对流项离散的稳定性及高精度特性.     

一维浅水波方程有限体积流通量限制方法的数值研究

本文讨论有限体积方法求解一维浅水波方程的流通量限制方法.该方法将高阶数值流通量与低阶数值流通量作适应性组合,在解的光滑区域表现为高阶精度,在解的不光滑区域表现为低阶精度从而抑制数值振荡.数值实验表明该方法能更高分辨率地求解一维非线性浅水波方程组.     

一类求解浅水波方程的基本无振荡熵稳定格式

针对浅水波方程,提出了一类低耗散基本无振荡熵稳定格式.在Roe型熵稳定通量中添加熵守恒格式的熵数值黏性绝对值的量来抵消解在跨过激波时所产生的熵增,从而抑制伪振荡;并且,利用通量限制器函数构造出相应的高分辨率熵稳定格式.利用新格式模拟一维和二维经典问题,数值结果表明,该格式具有低耗散、高分辨率、基本无振荡性等特点,是求解浅水波方程较为理想的方法.     

激波捕捉格式求解复杂床面浅水方程的源项处理

为了解决激波捕捉格式应用于浅水模拟所带来的源项,尤其是底坡项的处理这一难点问题,建立了一种边界条件给定简单的方法。采取构造单元内部流动变量分布的手段,避开了对底坡项的直接处理,综合了常用的Godunov型通量差分裂格式(flux difference splitting,FDS)和LeVeque波传播方法的思想。通过与"人工渗透法"相结合,该方法被应用于存在移动的干湿边界的情况。对两个有解析解和一个有实验观测数据的算例进行了模拟,模拟结果与解析解或观测数据吻合良好。该方法可有效处理复杂地形变化的二维浅水流动。     

激波捕捉格式求解复杂床面浅水方程的源项处理

为了解决激波捕捉格式应用于浅水模拟所带来的源项,尤其是底坡项的处理这一难点问题,建立了一种边界条件给定简单的方法。采取构造单元内部流动变量分布的手段,避开了对底坡项的直接处理,综合了常用的Godunov型通量差分裂格式(flux difference splitting,FDS)和LeVeque波传播方法的思想。通过与“人工渗透法”相结合,该方法被应用于存在移动的干湿边界的情况。对2个有解析解和1个有实验观测数据的算例进行了模拟,模拟结果与解析解或观测数据吻合良好。该方法可有效处理复杂地形变化的二维浅水流动。