具有指定边界网格间距和正交性的BFC网格生成方法 Research on generating method of the BFC gird with the adjustable boundary mesh intervals and orthogonality.pdfVIP

第3期 水利 水运工程学报 No．3 ANDENGINEERING 2009年9月HYDRO·SCIENCE Sep．2009 具有指定边界网格间距和正交性的BFC网格生成方法 蒋光彪1’2，何永森1，舒适1 (1．湘潭大学数学与计算科学学院，湖南湘潭411105；2．中南林业科技大学流变力学与材料工程研究 所，湖南长沙410004) 摘要：边界拟合曲线坐标(BFC)网格生成的质量直接影响流场的数值模拟精度和收敛速度，而如何选取生成 BFC网格的Possion方程右端调节因子尸和Q，则关系到生成的BFC网格质量．本文提出了一种可以同时实现 控制网格正交性和控制边界网格间距的BFC网格生成的新方法．该方法可任意控制生成的网格边界间距大小， 同时，生成的BFC网格还具有边界及内部较好的正交性．应用实例的计算结果表明，该方法能够对复杂边界的 单连通域或多连通域生成较理想的BFC网格，网格分布能够适应流动参数的变化． 关 键 词：BFC网格；Poisson方程；调节因子；正交性；网格间距 中图分类号：0242 文献标识码：A 紊流数值计算中常采用有限差分方法，在流场的数值模拟中，求解域复杂的几何边界给数值模拟带来了 很大的困难．如何用差分法处理复杂几何区域的紊流数值模拟是许多学者研究的课题．边界拟合坐标(BFC) 网格生成质量对于流场的数值计算精度具有重要意义．目前生成BFC网格最常用的方法是求解偏微分方程 组的Thompson微分方程法…．微分方程网格生成方法是网格生成方法中的一类经典方法．这类方法利用微 分方程的解析性质，如调和函数的光滑性、变换中的正交不变性等，通过偏微分方程组将物理空间坐标和计 算空间坐标联系起来．根据用来生成贴体网格微分方程类型的不同，微分方程网格生成方法又分为椭圆型方 程方法、双曲型方程方法和抛物型方程方法．其中，常用的方法是椭圆型方程方法¨1．微分方程网格生成方 法所生成的网格较代数网格光滑、合理、通用性强，同时微分方程网格生成方法可以生成贴体坐标网格，且对 流场中物理量梯度大处任意加密．本文采用椭圆型微分方程数值网格生成方法． 1 BFC网格生成概述 BFC网格生成的主要思想是基于一种BFC坐标变换，将物理平面内复杂的几何区域变换到计算平面内 的规则几何区域，以便用有限差分方法求解控制方程．其相应的数学问题就是求解一个带第一类Dirichlet边 界条件的边值问题．设物理面的坐标变量和计算面的坐标变量的关系为：戈：戈(f，7／)，r=r(孝，叼)；f=亭(菇， r)，叼=田(x，r)，物理面和计算面用以下椭圆型方程相关联： 亭；，+孝。：P(f，叩) 刀，，+叩。=Q(亭，，7) (1) 当计算面内是正交直线网格(坐标为(孝，t7))时，为了求物理面上的对应点(坐标为(茹，，，))，要通过从计 算面向物理面的反变换，沿边界形状自动地生成适合于边界的曲线网格，进而可控制全部求解域内的曲线坐 标网格的配置．反变换偏微分方程如下： 收稿日期：2008一10—27 l10676031) 基金项目：国家自然科学基金资助项目(5937521 作者简介：蒋光彪(1975一)，男，湖南永州人，讲师，博士研究生，主要从事计算流体力学的研究 63．oorU E—mail：jiang～gb@1 万方数据 68 水 利 水 运 工 程 学 报 2009年9月 a算舒一2,8互翱+yx叶可=一．，2(Pxf+Q戈竹) ，^、 oty妊一2#y钿+yy可可=一J2(Pyf+Q)，可) Q=g(亭，田)可控制网格的配置，称为调节因子，作用是调节物理平面内曲线网格的形状及疏密程度． 求解反