结构网格下求解非规则区域流场的一种差分法.pdfVIP

第 42卷 第 4期 中山大学学报 (自然科学版) Vo1．42 No．4 2003年 7月 ACTA SCIENTIARUM NATURAUUM UNIVERS1TATIS SUNYATSENI Ju1． 20o3 结构网格下求解非规则区域流场的一种差分法 吴春亮，詹杰民 (中山大学应用力学与工程系，广东广州 510275) 摘 要：对于粘性不可压缩流场的计算，提出一种标识边界的差分方法，对不规则区域的流场边界做标记，应 用SIMPLE算法在直角坐标系的结构网格下离散求解N—S方程；对圆柱绕流，不规则腔体流等问题进行数值模 拟，取得了较好的效果。 关键词 ：结构网格；非规则区域；差分法 中图分类号：0357 文献标识码：A 文章编号：0529-6579(2003)04-0009．04 由于N—s方程中速度和压力通过连续性方程 可以方便给出，这一类型的网格广泛应用于有限差 的耦合，直接求解是 困难 的，针对这一 困难， 分法中。显然 ，经过边界拟合坐标变换的网格也是 Spading和Patankar提出了SIMPLE算法，该方法在 结构网格。一种网格称为非结构网格 ，其单元编 不可压缩流场的计算中得到了广泛的应用，特别是 号 、结点编号都是随机的，如二维空间中的三角形 应用 结构 的规 则 网格求解 规则 区域 的流场； 网格。这一类型的网格是非有限差分型的，多用于 SIMPEL 算法中的速度压力校正思想后来应用到边 无粘流场的计算，如有限体积法。 界拟合坐标网格中，可以求解复杂区域的流场，然 结构网格有很多优点，形成的算法简单 ；当计 而由于求解区域经过边界拟合坐标变换后，相应地 算的流场区域为不规则区域时，产生适当的结构网 作变换后的流体力学方程也变得较为复杂，应用 格是困难的，即便是采用边界拟合坐标变换 ，也有 SIMPEL 算法并不方便。 出现网格奇异性的时候。Zhan针对这一情况 曾提 对于不规则区域的流场计算，考虑到 SMIPEL 出一种混网格系统 ，即采用两套网格系统，成功计 算法的特点，即在已知速度边界或是压力边界的情 算了弯曲渠道的流动[3】。即采用非结构网格，可以 况下，压力校正不用提边界条件…，这也意味着在 消除这些困难。然而，采用非结构网格，程序的编 这种情况下，可以采用交错网格，先构造出压力校 制较复杂，系统的存储量也较大，特别是计算界面 正的计算网格点，然后根据压力计算网格构造出与 通量的时候要应用邻近网格的信息，对于单元的邻 之交错的速度计算差分基点；对于构成不规则边界 近网格基点的搜索算法也是复杂的。虽然，利用本 的网格点，进行标记，在计算过程中自动识别。在 文的方法，在边界处理方面有一定的复杂性，但考 文中可以看到，这样处理仍然保留差分格式的精 虑到在计算区域的大部分地方，求解过程计算量 度；另外 ，有限差分和有限体积法中常用的分区计 小，精度高，数值稳定性好，值得发展一套高效 算以及局部加密的思想在这里仍然适用，这一方面 的，适合于求解非规则区域流场的计算方法． 的工作有待以后的研究探讨。从算例的结果来看， 在直角坐标系下，流体力学方程具有原始变量 该方法是有效的，可以发展为一种处理复杂边界形 的形式，差分离散和应用 SIMPEL 算法求解是十分 状的通用方法。 方便的；为了处理不规则 的边界，在直角坐标系 下，采用结构网格 ，通常要进行网格的调整，即插 1 计算