(6.7.1)--5.7-基于椭圆型方程的贴体网格生成方法讲义-zuixin.pdfVIP

《计算流体力学基础》讲义 第五章 网格与相应变换 5.7 基于椭圆型方程的贴体网格生成方法 本节介绍基于椭圆型方程的贴体网格生成方法。首先我们看一个简单问题的贴体网格坐 标系统。考虑图5.6 所示的扩张管道的网格系统。该扩张管道的壁面是曲线de，简单的直角 网格不再适用这样的外形。因此，需要曲线网格。将图5.6a 中曲线网格的上下边界de 和fg 作为坐标线，也就是使坐标线与物理边界完全重合。然后，我们需要将此曲线网格变换成计 算平面中的直角网格。为此，我们令ys f (x) 为图(5.6a)中上边界的坐标，通过下面的转换 （，） 即可得到一个 平面内的直角网格，即 ξ= x (5.65) (a)物理平面 (b)计算平面 图 5.6 扩张管道的贴体网格系统 y  = (5.66) y s 通过这个变换关系式，可以清楚地看出边界上的 c，d，e 点在计算平面中的坐标均为  1，也即在计算平面中，点c，d，e 有相同的η 坐标。因此，物理平面上边界上所有的 点都在计算平面内 1的水平坐标线上。这就在计算平面上形成了一个均匀的直角网格， 完成了该扩张管道的曲线网格向计算平面的变换。 - 1 - 《计算流体力学基础》讲义 刚才的图5.6 只是给出了一个简单的贴体网格坐标系统的例子。这里，图5.7 给出了一 个稍微复杂的翼型周围贴体曲线坐标系统。坐标线 1 常数在翼型表面上，是网格的内 边界Γ 。坐标线  常数给出的是网格外边界Γ2 。内边界Γ1 的形状和位置由其所在 1 2 的翼型表面所决定。外边界Γ2 的形状和位置原则上是任意的，一般根据实际情况确定一个 较大的区域作为计算域外边界。由内边界Γ1 出发到外边界Γ2 为止的网格线是ξ 为常数的 线。可以发现，η 为常数的线都是围绕翼型的封闭曲线，像是被拉伸的圆，这样的网格叫做 翼型的O 型网格。如果η 为常数的网格线一直向右侧下游延展，除了内边界Γ1 外，并不是 完全封闭的，呈C 字型，这样的网格叫做C 型网格。我们接着分析图5.7 中的O 型网格变 换。对于物理平面中的O 型翼型网格，如何将其变换成如图5.7b 所示的计算平面上的均匀 网格呢？ (a)物理平面 (b)计算平面 图 5.7 用椭圆型方程生成贴体网格实例 考察图 5.7a，对于一个外边界Γ2 的形状已经确定的给定翼型来说，内边界Γ1 上各点 物理坐标(x，y) 即为翼型边界点坐标，是已知的，外边界Γ2 上各点的物理坐标(x，y)也是已 知的，那么该问题本质上就是一个已知所有边界条件的边值问题。从数学物理方程的知识可 以知道，椭圆型偏微分方程所描述的问题即为边值问题。因此，我们可以考虑采用椭圆型偏 微分方程作为变换关系式。 一个最简单的椭圆型方程便是Laplace 方程： - 2 - 《计算流体力学基础》讲义