瞬间全溃溃坝波的传播反射和绕射的数值模拟.PDFVIP

　 辑第 15 卷第 1 期　　　　　 水 动 力 学 研 究 与 进 展 　　　 . , . 15, . 1 A Ser A V o l N o 　　2000 年 3 月　　　　　　　 　　　　　 . , 2000 JOU RNAL O F H YDROD YNAM ICS M ar 文章编号:(2000) 瞬间全溃溃坝波的传播、反射 和绕射的数值模拟 1 2 2 王嘉松 , 　　倪汉根 , 　　金　生 ( 1. 上海交通大学工程力学系, 上海 200030; 2. 大连理工大学土木系, 大连 116023) 　　摘　要: 　将二维浅水方程解耦并进行空间差分分裂, 应用通量限制的组合型 TVD 格式, 分 别求解分裂后的一维双曲型方程组, 针对瞬间全溃情况下溃坝波的运动问题, 对溃坝波的传播、反 射和绕射进行数值模拟, 揭示出溃坝波复杂的运动特性, 为溃坝灾害性分析和防灾减灾提供科学 依据。 关　键　词: 　浅水方程; 差分分裂; TVD 格式; 溃坝波; 反射; 绕射 中图分类号: 　TV 131. 4　　　文献标识码:A 1　前言 ( ) 溃坝波的运动属于带有自由表面的非恒定流, 由于存在间断波 涌波 而变得非常复杂, 常规的差分方法很难有效地处理这类问题。随着计算空气动力学的发展, 八十年代初中期提出 了各种高分辨率计算格式, 如 [ 1 ]、 [ 2 ] 格式, 对于求解欧拉方程, 解决分离涡、激波等复 TVD ENO 杂流场的计算问题效果十分显著, 具有很强的激波捕捉能力, 且在激波区不振荡。浅水方程在 性质上与欧拉方程类似, 这类格式为解决溃坝波的计算问题提供了思路。但上述高分辨率格式 是定义在一维双曲型守恒方程上, 对于二维双曲型方程, 只作形式上的推广, 至今没有严格的 [ 3 ] 理论基础 。因此本文将二维浅水方程进行空间的差分分裂, 分别应用 TVD 格式求解, 用以 模拟全溃溃坝波的传播、反射与绕射。 2　TVD 格式 针对下面的单变量双曲型守恒方程 u f (u ) + = 0 ( 1) t x 收稿日期: ( ) 基金项目: 　国家自然科学基金资助项目 ( ) 作者简介: 　王嘉松 1967～ , 男, 讲师, 博士。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 2 水　动　力　学　研　究　与　进　展　　　　　　　　　2000 年第 1 期 其差分格式为