适用于混合网格的改进雅可比迭代法及其应用

黄宇; 阎超; 袁武

doi:10.13700/j.bh.1001-5965.2015.0197

适用于混合网格的改进雅可比迭代法及其应用

doi: 10.13700/j.bh.1001-5965.2015.0197

黄宇¹,
阎超^1, ,,
袁武²

1.
北京航空航天大学航空科学与工程学院, 北京 100083
2. 中国科学院计算机网络信息中心超级计算中心, 北京 100190

详细信息

作者简介:
黄宇男,博士研究生。主要研究方向:计算流体力学、混合网格数值方法。Tel.:010-82338071 E-mail:huangyu@buaa.edu.cn;阎超男,博士,教授,博士生导师。主要研究方向:计算流体力学、空气动力学。Tel.:010-82317019 E-mail:yanchao@buaa.edu.cn

通讯作者:
阎超,Tel.:010-82317019 E-mail:yanchao@buaa.edu.cn

中图分类号: V221⁺.3;TB553
计量
- 文章访问数: 802
- HTML全文浏览量: 3
- PDF下载量: 698
- 被引次数: 0
出版历程
- 收稿日期: 2015-04-07
- 刊出日期: 2016-03-20

Improved Jacobi iterative method for hybrid grid and its application

HUANG Yu¹,
YAN Chao^{1
, ,},
YUAN Wu²

1.
School of Aeronautic Science and Engineering, Beijing University of Aeronautics and Astronautics, Beijing 100083, China
2. Super Computing Center, Computer Network Information Center, Chinese Academy of Sciences, Beijing 100190, China

摘要

摘要: LU-SGS因有较高的鲁棒性和小的内存需求而得到广泛应用,然而用于混合网格计算前需要进行网格排序和分组来实现算法并行;此外,LU-SGS格式收敛效率不高。针对这些缺点,本文提出了一种改进的适用于复杂混合网格的雅可比迭代方法,无需网格排序和分组就可实现算法的并行化,且有较快的收敛速度。该方法编程实现简单,易于采用OpenMP实现并行。算例研究表明,相比于LU-SGS格式,在各来流条件下,本文提出的方法收敛速度更快,鲁棒性好,并行和串行结果一致,且内存需求增加很少。
- 混合网格 /
- 隐式算法 /
- 并行计算 /
- OpenMP /
- 雅可比迭代 /
- 网格排序
Abstract: LU-SGS scheme is widely used today because of its robustness and cheap memory cost. However, the original LU-SGS shows less competitive convergence rate; in order to apply paralleled codes on hybrid unstructured grid, the grid reordering and regrouping procedure must be carried out beforehand. In this paper, an improved implicit method suitable for complex hybrid gird is developed to achieve fast convergence rate and to parallelize the algorithm without grid reordering and regrouping procedure. This method is simple for coding and easy to use OpenMP for code parallelization. The numerical results of Euler and viscous flows show that the method has a reliable performance, and it is able to achieve a significant efficiency improvement over implicit counterparts such as LU-SGS scheme with less requirement of extra memory, and parallel computation produce exactly the same result as serial case.
- hybrid grid /
- implicit method /
- parallel computation /
- OpenMP /
- Jacobi iteration /
- grid reorder

HTML全文

参考文献(22)

[1]	阎超,于剑, 徐晶磊,等.CFD模拟方法的发展成就与展望[J].力学进展,2011,41(5):562-589. YAN C,YU J,XU J L,et al.On the achievements and prospects for the methods of computational fluid dynamics[J].Advances in Mechanics,2011,41(5):562-589(in Chinese).
[2]	VENKATAKRISHNAN V, MAVRIPLIS D J.Implicit solvers for unstructured meshes[J].Journal of Computational Physics,1993,105(1):83-91.
[3]	KNIGHT D D. A fully implicit Navier-Stokes algorithm using an unstructured grid and flux difference splitting:AIAA-1993-0875[R].Reston:AIAA,1993.
[4]	JAMESON A, YOON S.Lower upper implicit scheme with multiple grids for the Euler equations[J].AIAA Journal,1987,25(7):929-935.
[5]	MAVRIPLIS D J. On convergence acceleration techniques for unstructured meshes:AIAA-1998-2966[R].Reston:AIAA, 1998.
[6]	SHAROV D, LUO H,BAUM J D,et al.Implementation of unstructured grid GMRES+LU-SGS method on shared-memory,cache-based parallel computers:AIAA-2000-0927[R].Reston:AIAA,2000.
[7]	CHEN R F, WANG Z J.Fast,block lower-upper symmetric Gauss Seidel scheme for arbitrary grids[J].AIAA Journal,2000,38(12):2238-2245.
[8]	TOMARO R F, STRANG W Z,SANKAR L N.An implicit algorithm for solving time dependent flows on unstructured grids:AIAA-1997-0333[R].Reston:AIAA,1997.
[9]	SHAROV D, NAKAHASHI K.Reordering of 3-D hybrid unstructured grids for vectorized LU-SGS Navier-Stokes computations:AIAA-1997-2102[R].Reston:AIAA,1997.
[10]	LUO H, BAUM J D,LÖHNER R.A fast,matrix-free implicit method for compressible flows on unstructured grids[J].Journal of Computational Physics,1998,146(2):664-690.
[11]	SAAD Y.GMRES:A generalized minimum residual algorithm for solving nosymmetric linear systems[J].SIAM Journal on Scientific and Statistical Computing,1986,7(3):856-869.
[12]	ZHANG L P, CHANG X H,DUAN X P,et al.A block LU-SGS implicit unsteady incompressible flow solver on hybrid dynamic grids for 2D external bio-fluid simulations[J].Computers & Fluids,2009,38(2):290-308.
[13]	LUCAS P,BIJL H, VAN ZUIJLEN A H.Efficient unsteady high Reynolds number flow computations on unstructured grids[J].Computers & Fluids,2010,39(2):271-282.
[14]	WANG G, JIANG Y W,YE Z Y.An improved LU-SGS implicit scheme for high Reynolds number flow computations on hybrid unstructured mesh[J].Chinese Journal of Aeronautics,2012,25(1):33-41.
[15]	蒋跃文,叶正寅, 王刚.基于非结构网格的高效求解方法研究[J].计算力学学报,2012,29(2):217-223. JIANG Y W,YE Z Y,WANG G.Efficient solution of Euler/NS equations on unstructured grids[J].Chinese Journal of Computational Mechanics,2012,29(2):217-223(in Chinese).
[16]	BLAZEK J. Computational fluid dynamics:Principles and applications[M].Oxford:Elsevier,2001:16-17.
[17]	ROE P L. Approximate riemann solvers,parameter vectors and difference schemes[J].Journal of Computational Physics,1981,43(2):357-372.
[18]	BARTH T J. A 3-D upwind Euler solver for unstructured meshes:AIAA-1991-1548[R].Reston:AIAA,1991.
[19]	VENKATAKRISHNAN V. Convergence to steady state solutions of the Euler equations on unstructured grids with limiters[J].Journal of Computational Physics,1995,118(1):120-130.
[20]	SPALART P R, ALLMARAS S R.A one-equation turbulence model for aerodynamic flows:AIAA-1992-0439[R].Reston:AIAA,1992.
[21]	LEE-RAUSCH E M, FRINK N T,MAVRIPLIS D J,et al.Transonic drag prediction on a DLR-F6 transport configuration using unstru-ctured grid solvers:AIAA-2004-0554[R].Reston:AIAA,2004.
[22]	BENAY R,CHANETZ B, MANGIN B,et al.Shock wave transitional boundary-layer interactions in hypersonic flow[J].AIAA Journal,2006,44(6):1243-1254.