基于无结构网格有限体积法的Roe隐格式探讨.pdfVIP

基于非结构网格有限体积法的Roe 隐格式初探 杨彬 河海大学，南京 (210098) E-mail ：tuohuangzhe119@163.com 摘 要：本文采用非结构网格有限体积法，尝试建立适合河流、湖泊和近海浅水水流流动的 Roe 隐式离散格式的数学模型。对于隐式离散方程组采用线性化的处理方法以后，用对称的 Gauss-Seidel(GS)法进行迭代求解，提高了求解效率。鉴于地形的处理对数学模型稳定性和 合理性的影响，将斜底模型的概念应用于模型的地形处理当中。运用模型对几种标准算例进 行验证模拟，并将模拟结果与相同条件下建立的 Roe 显式离散格式模型的计算结果进行分 析比较，通过计算表明本模型可以有效的缩短模拟计算时间，并且可以较好的处理间断解和 各种复杂流态的过渡。 关键词：二维浅水方程，无结构网格，Roe 隐格式，对称GS 法，地形处理 中图分类号：TV131.2 1. 引言 在实际计算不规则天然地形上的浅水流动时，为了贴切地模拟陆地边界，整个计算区域 的网格尺寸常相差较大，在利用显式格式时，时间步长往往受到所有网格稳定性限制中最小 者的控制，极大地降低了处理效率，特别是在计算历时较长的数值模拟中，计算机的实际耗 时往往让人难以接受，从而使数值模拟受到限制。鉴于隐式格式的高效率和较好的稳定性， 以及近年来非结构网格技术的发展，有限体积法的广泛应用，许多高性能的隐式新算法在计 算空气动力学领域中脱颖而出[1,2,3] 。然而在计算水力学中，目前利用无结构网格有限体积法， 基本上是用的显式格式。因此借鉴空气动力学中高性能的隐式格式，在浅水动力学中建立与 本领域相适应的高性能的隐式格式不但有利于恒定流或解的时间变率小、流动过程历时长等 一些实际问题的数值模拟，而且也能够有效的处理间断解。本文尝试在二维浅水方程的基础 上建立基于无结构网格有限体积法的 Roe 隐式离散格式数学模型，利用此模型对几种标准 算例进行验证模拟，并将计算结果与相同条件下建立的 Roe 显式离散格式模型的计算结果 进行分析比较。 2. 控制方程 [4] 基于圣维南方程的基本假设，本模型采用浅水方程组 : ∂U U E G ∂ ∂ ∂ +∇⋅F + + S (1) ∂t ∂t ∂x ∂y T 2 2 T 2 2 T [ ] [ ] 式中： U [h , hu , hv ] ; E hu , hu +gh / 2, huv ; G hu , huv , hv +gh / 2 ; S S +S T T u ,v分别为x ,y 方向的水深平均 [ ] [ ] 0 f 0, ghS 0 x , ghS 0 y + 0,−ghS fx ,−ghS fy ; 流速; t为时间;h为水深;g 为重力加速度; S 0 x , S 0 y 分别为x ,y 方向的底坡; S fx , S fy 分别为x ,y 方 向 的 摩 阻 坡 度 , 在 本 文 中 摩 阻 坡 度 可 用 曼 宁 公 式 计 算 , 即 : S fx n 2 u u 2 +v 2 / h 4 / 3 , S fy n 2 v u 2 +v 2 /