一种处理不规则边界的二维切削网格方法研究.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：093491 一种处理不规则边界的二维切削网格方法 1 2 3 1 1 菅振 ，孙东亮 ，吴志根 ，李增耀 ，陶文铨 （1.西安交通大学 动力工程及多相流国家重点实验室 能源与动力工程学院，西安 710049； 2.华北电力大学 可再生能源学院，北京 102206； 3. School of Engineering，University of Warwick，United Kingdom） (Tel:029 Email: wqtao@mail.xjtu.edu.cn) 摘要：切削网格方法(Cut-cell)是一种处理复杂边界的网格生成方法，本文从切削网格分类、网格合并 技术和切削界面相关物理量求解等角度对方法进行了介绍。最后，将该网格生成技术用于计算小雷诺 数下二维圆柱绕流问题，并与实验结果做了比较，验证方法和所开发程序的可靠性。 关键词：SIMPLE 算法 正交“背景网格” 切削网格 速度插值。 0 前言 对于复杂区域生成高质量的网格，一直是数值模拟中十分重要的研究方向。为了解 [1] [2] 决这一问题，多种网格生成方法相继产生。其中非结构化网格 和适体网格 是两类十 分重要的网格生成方法。但这两类方法的实施相对复杂。 切削网格方法是另一类网格生成方法，其本质上是在正交“背景网格”的基础上， 对固体边界处网格进行修正。本文从控制方程、网格划分、网格合并、控制方程离散、 切削界面插值和切削网格固体边界梯度求解等方面介绍了一种二维切削网格的生成方 法，并对圆柱绕流问题进行了实验与模拟的结果对比分析。 1 控制方程 计算针对二维定常流动、不可压缩、常物性流体展开，采用守恒形式的控制方程可 表示为： ∂(ρu) ∂(ρv) 连续性方程 + 0 （1） ∂x ∂y ∂ ∂ ∂p ∂ ∂u ∂ ∂u uu + vu =− + + +S u −动量方程 (ρ ) (ρ ) (ν ) (ν ) u （2 ） ∂x ∂y ∂x ∂x ∂x ∂y ∂y ∂ ∂ ∂p ∂ ∂v ∂ ∂v v − uv + vv =− + + +S 动量方程 (ρ ) (ρ ) (ν ) (ν ) v （3 ） ∂x ∂y ∂x ∂x ∂x ∂y ∂y 其中，S , S 分别为u 、v 动量方程的源项。 u v 致谢：本文得到国家自然科学基金的资助（No.，No.） 2 切削网格方法 切削网格方法的主要思想是，在一套简单的正交“背景网格”的基础上，对不规则 固体边界区域的网格根据固体边界的形状进行修正。从而在整体上网格划分相对简单的 前提下，又使固体边界处的网格可以得到更好的逼近，以提高计算精度。 2.1 网格种类