(6.8.1)--5.8-网格生成与流场求解方法示例讲义.pdfVIP

《计算流体力学基础》讲义 第五章 网格与相应变换 5.8 网格生成与流场求解方法示例 本节主要介绍如何采用偏微分方程数值解法生成NACA0012 翼型O 型或C 型网格，并 采用速度势方程求解绕 NACA0012 翼型迎角为零时的低速不可压无黏流场。首先以 NACA0012 翼型为对象，讲解其网格生成方法。已知的条件为以下三方面。 一是对称翼型NACA0012 翼型的表面(包括上下表面点)坐标计算公式，可以表示为 4 3 2 yu ,l =0.6[-0.1015x + 0.2843x - 0.3576x - 0.1221x + 0.2969 x ] (A5.1) 二是网格生成的控制方程为Laplace 方程，其已经写成了计算平面中的变换形式，其中 系数α ，β ，γ 均由反变换系数组成  x 2 =+ y 2 ,  x x =+ y y ,  x 2 =+ y 2， J 2 (A5.2)         2 2 2  y  y  y  2 −2 + 2 0 (A5.3)    其中，  x 2 =+ y 2 ,  x x =+ y y ,  x 2 =+ y 2         三是流动控制方程速度势方程，可以表示为 2 2     + 0 (A5.4) x 2 y 2 即笛卡尔坐标系下的Laplace 方程。下面我们来简要讲解一下求解思路。 首先看看采用偏微分方程数值解法生成网格的大致步骤。第一步是设定边界条件。分为 两个部分。对于 O 型网格，远场一般设定为半径是一定倍数弦长的圆，例如可以设定以翼 型中心点为圆心，以15 倍弦长为半径的圆，作为远场边界，根据设置的离散点数，确定远 场网格点的(x ，y )坐标值。确定翼型表面点分布的原则是保证流动变化剧烈的区域分布密集， 变化平缓的区域分布稀疏。例如可以按照式(A5.5) 的关系式布置翼型表面点。 c x (1=+cos)， [0,2] (A5.5) 2 对于 C 型网格，远场边界可以设定成以翼型中心为圆心，以 15 倍弦长为半径的半圆， 加上15 倍弦长 30 倍弦长的矩形所组成的区域边界。确定翼型表面点分布的原则与上述O  型网格情况类似，保证流动变化剧烈的区域分布密集，变化平缓的区域分布稀疏。除了可以 使用式(A5.5) 的分布形式，也可以在x=0-0.5c 和x=0.5c-c 分段采用式(A5.6