非结构网格中WENO型有限体积法溃坝模拟.pdfVIP

非结构网格中WENO 型有限体积法溃坝模拟 王伟 河海大学环境科学与工程学院，南京 210098 Email ：wangwe2005@126.com 摘 要：本文采用非结构网格中WENO 格式构造方法，在Roe 的通量差分分裂的基础上， 对二维浅水方程中的物理量h,u,v 进行数值重构，并与Upwind 格式和ENO 格式进行对比， 结合二维溃坝模型对三者分别进行编程模拟，模拟结果表明，WNEO 格式不但精度高，分 辨率高，而且具有较强的捕捉激波的能力，明显优越于Upwind 格式和ENO 格式。 关键词： 非结构网格；ENO 格式；WENO 格式；溃坝模拟 中图分类号：0175.2 1 引言 有限体积法 finite volume method，简记为FVM 是八十年代以来发展起来的一种新型的 微分方程的离散方法，具有独特的优点，目前已成为偏微分方程问题和计算流体力学问题的 数值计算、数值模拟中一个重要的方法。 然而在非结构网格有限体积法 FVM 中,数值重构是个重要环节，直接决定数值方法的 [1] [2] 空间精度，分辨率。目前较为常见的重构方法有MUSCL 格式 、迎风 Upwind 格式 、Quick 格式等。本文采用WENO 格式对二维浅水方程中物理量进行数值重构，结合溃坝模型进行 数值模拟，并与Upwind 格式和ENO 格式进行对比分析。 2 二维浅水方程 Euler 型 及离散 2.1 二维浅水方程 Euler 型 二维浅水方程守恒形式为[3] U + f U + g U 0 1 t x y 其中，U 为速度矢量，f ，g 分别为x 和y 方向上的通量。 h ? hu ? ? hv ? ? ? ? ? ? 2 2 ? ? ? U hu ，f hu +gh / 2 ，g ? huv ? ? ? ? ? ? ? ? huv ? ? 2 2 ? hv hv +gh / 2 ? ? ? ? ? ? 2.2 空间离散 本文采用有限体积法 FVM 离散方程，在计算过程中，我们用的是三角形网格，而且构 造格式是格子中心式。 记控制单元为V ，并在其上对方程 1 积分得[4] ? ? U dxdy + f U + g U dxdy 0 ∫V t ∫V ? x y ? 由格林公式可得 AU + F ?ndl 0 2 t ∫?V 其中A 控制元V 的有向面积，F f , g ，n n n 为单元边的外法向量。 x , y 将式 2 线性积分化得其半离散格式 - 1 -