NACA0012翼型O型网格生成及不可压无粘流场数值求解(0°攻角、2°攻角)

问题描述

        采用偏微分方程数值解法生成NACA0012翼型O型网格，并采用速度势方程求解绕NACA0012翼型迎角为零度和2°时的低速不可压无粘流场。

NACA0012翼型（修形）表面坐标为：

（42）

O型网格生成

用于网格生成的椭圆型偏微分方程        

        贴体网格的生成可以视为计算空间均匀网格向物理空间实际网格映射的过程。Laplace方程为椭圆型偏微分方程，计算域内点值由边界值唯一确定，因此可用于贴体网格的生成。应用椭圆型偏微分方程生成内点网格具有二阶可导性，且当边界有斜率不连续点时，内点网格分布仍然是光滑变化的。式（43-44）给出了经过雅可比变换后的Laplace方程，方程中0项已被略去。

		（43）
		（44）

        式中(x,y)对应物理空间网格点坐标，(ξ, η)对应计算空间网格点坐标。α β γ计算方法由式（45-47）给出。

		（45）
		（46）
		（47）

    

离散差分格式和数值求解方法

        采用中心差分格式对式（43-47）进行离散，离散格式较为简单，这里仅给出混合偏导项的离散格式，如式（48-49）所示，式中 ∆ξ=1，∆η=1。

		（48）
		（49）

       对离散后的差分方程整理可得简单迭代格式，如式（50-51）所示。式（52-54）给出了 α β γ  的差分格式。

		（50）
		（51）
		（52）
		（53）
		（54）

边界条件与初始条件

        网格生成中求解偏微分方程的初始条件与边界条件（壁面及物理空间计算域外边界网格点坐标）需根据翼型及流场特点给定，以在所需位置获得加密网格。对于二维翼型绕流，在无激波的条件下，其前缘及后缘流场参数变化相对剧烈。因此在给定壁面网格坐标时，其横坐标可由式（55）确定，以在前缘和尾缘附近获得更加密集的网格点分布。起始点位于翼型后缘点（c，0），θ 按顺时针方向式均匀递增。式中，c为翼型弦长，物理空间（0，0）坐标落在翼型前缘点，计算空间（0，0）坐标落在翼型后缘点。将式（55）代入翼型方程式（42）即可获得壁面网格点的纵坐标。

（55）

        物理空间计算域可设定为半径为50倍翼型弦长的圆，式（56）给出了其曲线方程，物理空间计算域外边界网格点自（0，c/2+50c）顺时针均匀分布在该曲线上。

（56）

       以上，确定了壁面网格点与外边界网格点的次序。按照一一对应关系，插值出每对壁面/外边界网格点之间的网格点坐标，作为物理空间内网格点坐标的初始值。采用周期性边界条件可得到割缝上 （ i=0）所有网格点的迭代计算公式。至此，边界条件与初始条件全部确定。

计算结果与分析

       因为网格求解中，迭代步数主要影响最终获得网格点的坐标，不会直接影响最终的流场求解，因此无需进行过多次数的迭代，对于本文算例，令收敛条件 eps<1e-4

        疏、密两网套格视图如下所示。