基于复杂网格处理高精度数值积分技术.docVIP

基于复杂网格处理高精度数值积分技术 　　摘要：目前，对飞机表面进行网格划分并进行流体计算以优化整体结构的手段被广泛运用，然而由于外形复杂，采用笛卡尔网格划分时会出现带曲边或曲面的复杂网格，若直接进行数值积分则精度有限。该文提出了对复杂网格进行映射变换处理后再进行数值积分，并且数值实验表明采用该处理方法大大提高了精度。 　　关键词：网格划分；数值积分；复杂网格处理 　　中图分类号：TP393 文献标识码：A 文章编号：1009-3044（2015）29-0179-02 　　High-accuracy Numerical Integration Technology Based on Processing Complex Grids 　　LIN Yi-miao 　　（School of Mathematics and System Sciences， Beihang University， Beijing 100191， China） 　　Abstract： Currently， the method of aircraft surface meshing and computational fluid dynamics to optimize the overall structure is widely used， however， due to complex shape， there exist complex grids with curved edge or surface when using Descartes grid. The accuracy is limited if doing the numerical integration directly. This paper proposes that processing complex grids with mapping transformation before numerical integration， and the numerical experiment shows that using the method of processing complex grids can greatly improve the accuracy. 　　Key words： meshing； numerical integration； processing complex grids 　　随着数字化风洞的不断研究，使用计算机模拟飞机在风洞中的表现，大大降低实验成本来优化飞机结构，将飞机表面进行网格划分并且进行流体力学计算是在所难免的手段。 　　然而，由于飞机外形结构复杂，采用笛卡尔网格对飞机表面进行网格划分时，所出现的数值积分区域就不是标准的积分区域了。例如二维翼型网格划分后会出现带曲边边界的三角形和四边形网格，三维翼型网格划分后会出现带曲面的四面体网格。 　　目前，一般的数值积分技术中，对于带曲边边界的网格是直接线段化处理，然而这样做精度有限，而如果将网格加密处理，精度提高有限却会增加计算负担和降低计算效率，所以在本文将提出一种保留曲边边界进行处理并且提高了数值积分精度的方法。 　　1 复杂网格介绍 　　网格划分属于流体力学计算中的预处理部分，是将计算区域划分为较小的、不重叠的子域或单元（网格），从而得到有限单元及节点。而单元的数目决定了CFD的计算精度，网格的细密程度决定了数值解的精度和计算时间。 　　在采用笛卡尔网格进行划分时，由于飞机外形复杂，例如对飞机机翼表面进行网格划分时，会出现一边为曲边的四边形网格（或三角形网格）和两边为曲边的四边形网格（或三角形网格），而如今随着不断研究，飞机外形也是日趋复杂化，相应地对飞机表面进行网格划分所得到的网格也是日渐复杂。 　　当然，若是网格划分足够密，每个网格足够小，那么曲边在一定程度上也就相当于直边了，所以在目前的数值积分技术中，都是采用直接将曲边看成直线进行处理，但如此做法必然会影响到最后的数值积分精度，所以本文提出一种基于复杂网格处理的数值积分技术，来提高数值积分的精度。 　　2 复杂网格处理 　　现假设在二维平面[?2]上有一个带曲边边界的四边形网格，要求解函数[fx，y]在其积分区域[Ω]上的数值积分，一般情况下，直接将曲边看成直线边进行处理，利用上面介绍的数值积分方法（Simpson或Gauss方法等）求数值积分得到所要的结果，下面来介绍一下我们的处理方法。 　　首先建立一般化的直角坐标，坐标平面为[xOy]平面，四边形网格的四个顶点坐标所对应的向量依次为[γ1，γ2，γ3，γ4]，其中向量[γ2，γ4]所对应的顶点连接的正是网格的曲边。现有另一坐标系[O-ξη]