2011
年第
35
卷
第
2
期
中国石油大学学报(自然科学版)
Joumal of China University of Petroleum
Vo
l. 35
No.2
Apr. 2011
文章编号:
1673-5005 (2011
)02
-0
131
-0
6
自动生成四边形网格的方法及其在数值
模拟中的应用
修荣荣
1
,
2
徐明海
2
黄善波
2
(1.中国石油大学期刊社,山东东营
257061;
2.
中国石油大学储运与建筑工程学院,山东青岛
266555)
摘要:采用间接方法生成四边形网格,首先利用改进的两点前沿推进法把计算区域剖分成三角形网格,然后采用插
点和细分的技术生成单元全部是四边形的网格,通过边互换、删点和局部插点技术进一步光滑平
j
帜,得到适用于数
值计算的网格。iftJ分结果表明,该方法能够在任意二维平面区域内自动生成全四边形网格,并能生成光滑过渡的局
部加密网格和贴体性较好的边界层网格。该方法具有算法简单,计算量少的特点。利用所生成的网格对计算传热
学中的典型算例-方腔自然对流进行求解,计算结果与基准解吻合,网格质量能够满足数值分析计算的要求。
关键词:二维区域;全四边形网格;三角形合并;变密度网格;数值模拟
中图分类号
:TK
124; TP 39
1.
41
文献标志码:
A doi: 10.
3969/j.
issn.
1673-5005.2011.02.023
Automatic generation method
of
quadrilateral meshes
and
its application to numerical simulation
XIU
Rong-rong
1
•
2
,
XU
Ming-hai
2
,
HUANG
Shan-ho
2
(1. Periodical Office
of
China University
of
Petroleum ,
Dongying
257061 ,
ChiTUL;
2. College
of
Storage & Transportation
and
Architectural Engineering
in
ChiTUL
University
of
Petroleum , Qingdao 266555 ,
ChiTUL)
Abstract:
Quadrilateral meshes on arbitrary 2D domain were generated
by
using indirect method. Firstly , triangle meshes
were generated by improved two-point advancing method. Then
by
means of inserting nodes and subdivison technique , the
mesh elements were composed with quadrilaterals completely. The edge interconversion
, node elimation and inserting node
locally were used to smooth meshes
, and the meshes suitable
to
numerical calculation were obtained. The proposed method is
characterized by its simplicity and ease of implementation. Mesh generation results demonstrate the robustness of the method.
The method can also generate smooth graded quadrilateral mesh and body-fitted boundary layer meshes. The generated me-
shes were used to solve natural convection in a
squ
缸
e
cavity. The solution results agree well with the reference results. The
quality of the generated meshes can be satisfied with the requirements of finite element numerical simulation.
Key
words:
two-dimensional domain; all-quadrilateral mesh; triangle merging; graded mesh; numerical simulation
在二维区域生成非结构化网格的方法中,三角
形网格的自动生成算法最为成熟,但有限元计算表
明,在相同网格步长下,三角形单元的计算精度不如
四边形单元
[l]
。
Aziz
等
[1)
研究表明,有限容积法数
值计算结果的精度与单元形状关系不大,与网格步
长关系则极为显著。在同样网格节点分布的情况
下,三角形单元的数目是四边形单元的两倍,因此,
收稿日期
:2010-04-20
最终形成的线性代数方程组的阶数也是两倍。这
样,无论是存储还是求解,都是不经济的。因此,为
了在一定计算工作量的条件下尽量提高精度,提出
了采用非结构化四边形网格计算思想,进而导致了
四边形网格生成的技术问题。生成四边形网格的方
法很多
[2
-8)按照单元生成方法大致分为两类:直接
法和间接法。直接法[矶
8)
包括几何分割法和前沿
基金项目:国家自然科学基金项目
(50904077)
;中央高校基本科研业务费专项资金资助项目(
IOCX050
Il
A)
作者简介:修荣荣(1
977-
),女(汉族)
,山东莱阳人,硕士,主要从事
CDF/NHT
技术研究。