(6.1.1)--5.1-引言讲义传感器与检测技术.pdfVIP

《计算流体力学基础》讲义 第五章 网格与相应变换 5.1 引言 本节课我们学习第五章，网格与相应变换，本章的主要内容包括空间变换和网格生成两 大部分。在空间变换部分我们将介绍变换系数与雅可比行列式，以及强守恒型流动方程的推 导。在网格生成部分，将介绍用于有限差分法的网格生成和用于有限体积法的网格生成两种。 最后介绍用于复杂外形的网格生成技术，包括自适应网格、混合网格、分区网格等。 图A5.1.1 本章知识路线图 首先，我们先明确几个概念。第四章所讨论的有限差分法是要求在离散的节点上进行计 算的，那么流场计算域中分布的所有离散点以某种形式关联成一个整体，便构成了数值计算 所需要的网格。组成网格的这些离散点便称为网格节点或者网格点。确定网格节点位置和坐 标的过程，我们称为网格生成。 由第四章内容可以知道，标准的有限差分法要求在均匀网格上进行，如例u 对x 的二阶 x y 导数差分算子要求 必须相等，同理 也必须相等。那么对于非均匀网格怎么运用有限 差分法呢？此时必须将非均匀网格变换成为一个均匀的、正交的网格，伴随这个变换，流动 控制方程必须作相应变换，以便可以在这样的均匀正交网格上使用。因此，网格变换对有限 差分方法来说一般是必须的，因为实际工程问题中很少有均匀网格的情形。 2 u −2u +u  u i+1, j i , j i−1, j ( 2 )i, j 2 x (x ) 图A5.1.2 均匀网格上的有限差分法 - 1 - 《计算流体力学基础》讲义 例如，假设要计算绕图5.1 中翼型的流动，我们将翼型放在一个直角坐标均匀网格中。 我们可以清楚地看到这样的网格存在以下问题：一是有部分网格点落在翼型的内部，但它们 完全是在流场之外，显然这些网格点无法用于描述流动特性；二是很少或者说几乎没有网格 点恰好落在翼型的表面，而翼型表面是决定流场的重要边界条件，因此翼型表面信息必须明 确地由网格点表示出来，显然，这样的网格无法描述翼型边界。所以，可以得出结论：图中 的直角坐标均匀网格并不适合求解这样的流场。 图 5.1 在矩形网格内的翼型 这里图5.2(a)给出了一个合适的网格，它是一个非均匀的由曲线组成的网格，紧贴在翼  型的表面，环绕着翼型，我们称其所在的x-y 平面为物理平面。图中给出了新的坐标轴 和 ，在新坐标系下，翼型表面是一条坐标线， 等于常数，这种坐标系叫做贴体坐标系。   对于这样的曲线网格，网格点自然地落在翼型表面上，网格线并不是互相垂直的，也不是均 匀分布的。传统的差分方法很难在这样的网格上应用，必须将物理平面上的曲线网格变换到  ， 坐标系中的互相垂直的网格上，即如图 b 所示的正交网格，这种网格线互相垂直的   网格所在的 - 平面称为计算平面。所进行的变换必须保证：图b 计算平面和图a 物理平 面内的各点是一一对应的。例如，物理平面内的点a，b，c 对应于计算平面内的点a，b，c，   计算平面内的网格点是均匀分布的，间距分别为 ， 。偏微分控制方程在计算平面上 采用有限差分法来求解，然后通过网格点的一一对应直接转化到物理平面上。注意，在计算