非结构混合网格自适应并行技术

doi:10.7527/S1000-6893.2019.23202

流体力学与飞行力学

本期目录 | 过刊浏览 | 高级检索

前一篇 | 后一篇

非结构混合网格自适应并行技术

唐静^1,2, 张健¹, 李彬¹, 崔鹏程¹, 周乃春¹

1. 中国空气动力研究与发展中心计算空气动力研究所, 绵阳 621000;
2. 西北工业大学航空学院, 西安 710072

收稿日期:2019-06-10 修回日期:2019-06-24 出版日期:2020-01-15 发布日期:2019-07-02
通讯作者: 李彬 E-mail:leebin2008@hotmail.com

Parallel algorithms for unstructured hybrid mesh adaptation

TANG Jing^1,2, ZHANG Jian¹, LI Bin¹, CUI Pengcheng¹, ZHOU Naichun¹

1. Computational Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China;
2. School of Aeronautics, Northwestern Polytechnical University, Xi'an 710072, China

Received:2019-06-10 Revised:2019-06-24 Online:2020-01-15 Published:2019-07-02

摘要/Abstract

摘要： 计算流体力学（CFD）模拟实际工程问题所采用的网格规模可达千万量级，并行技术是减少计算时间的有效方法。耦合流场信息的网格自适应技术能有效动态优化计算网格，被NASA视为一项亟待发展的CFD关键技术。混合网格自适应系统包含网格分布优化、表面网格投影和空间网格匹配等关键技术。针对以上3项关键技术分别建立了高效的并行算法。首先，提出了"先唯一后同一"的两步法策略实现了网格单元分布优化过程的并行相容性；其次，基于局部曲面拟合思想，实现了曲面重构和新增物理网格点投影的完全并行；再次，提出了空间网格匹配技术的半并行算法，快速解决了网格单元交错问题。为了提高后续流场计算的并行效率，发展了基于并行重分区-网格数据迁移方法的动态负载平衡技术，并采用圆柱激波流场自适应模拟对动态负载平衡技术进行初步验证。最后，采用三角翼自适应加密测试了自适应系统的并行效率。结果表明，建立的混合网格自适应系统并行效率较高，且相比流场求解耗费总时间的比例低于1%。

关键词: 混合网格, 网格自适应, 并行技术, 并行相容性, 动态负载平衡, 并行效率

Abstract: The mesh elements can reach ten-millions for real engineering cases to simulate their flow field with CFD. The parallel technique is an effective way to reduce the computing time. Mesh adaptation coupled with flow information can dynamically optimize computing mesh and has been considered by NASA as a key technique to be urgently developed. Elements distribution optimization, surface points projection, and inner points matching are three essential techniques for hybrid mesh adaptation. In this study, parallel algorithms for these three techniques are established. First, the two-stage strategy, which is called uniqueness-then-identity, is proposed to obtain the parallel consistence for the process of elements distribution optimization. Next, based on the local surface fitting method, the entire parallel algorithm is developed for surface reconstruction and new physical points projection. Then, the half-parallel algorithm for inner points movement is devised to efficiently eliminate the overlapping of mesh elements. To increase the parallel efficiency of the sequent flow field calculation, the dynamic load balance method, based on parallel partitioning and mesh transferring methods, is developed, and the flow around a cylinder is simulated to verify the dynamic load balance method. Finally, mesh refinement with adaptation for a delta-wing was adopted to test the parallel efficiency of the adaptation system. Numerical results show that the parallel efficiency of the hybrid mesh adaptation system is high enough and the total time cost is less than 1% of the flow solver.

Key words: hybrid mesh, mesh adaptation, parallel algorithms, parallel consistence, dynamic load balance, parallel efficiency

中图分类号:

V211.3

唐静, 张健, 李彬, 崔鹏程, 周乃春. 非结构混合网格自适应并行技术[J]. 航空学报, 2020, 41(1): 123202-123202.

TANG Jing, ZHANG Jian, LI Bin, CUI Pengcheng, ZHOU Naichun. Parallel algorithms for unstructured hybrid mesh adaptation[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2020, 41(1): 123202-123202.

参考文献 32

[1]	BAKER T J. Mesh generation:Art or science[J]. Progress in Aerospace Sciences, 2005, 41(1):29-63.
[2]	SLOTNIK J, KHODADOUST A, ALONSO J, et al. CFD vision 2030 study:A path to revolutionary computational aeroscience:NASA/CR-2014-218178[R]. Washinton,D.C.:NASA, 2014.
[3]	孙旭, 张家忠, 黄科峰. 基于弹簧近似的非结构化网格自适应处理方法[J]. 西安交通大学学报, 2010, 44(9):104-108. SUN X, ZHANG J Z, HUANG K F. Spring analogy-based adaptive method for unstructured grids[J]. Journal of Xi'an Jiaotong University, 2010, 44(9):104-108(in Chinese).
[4]	GREENE P, SCHOFILD S, NOURGALIEV R. Dynamic mesh adaptation for front evolution using discontinuous Galerkin based weighted condition number relaxation[J]. Journal of Computational Physics, 2017, 335:664-687.
[5]	唐志共, 陈浩, 毕林, 等. 自适应笛卡尔网格超声速黏性流动数值模拟[J]. 航空学报, 2018, 39(5):121697. TANG Z G, CHEN H, BI L, et al. Numerical simulation of supersonic viscous flow based on adaptive Cartesian grid[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(5):121697(in Chinese).
[6]	MARCUM D, ALAUZET F. 3D metric-aligned and orthogonal solution adaptive mesh generation[J]. Procedia Engineering, 2017, 203(1):78-90.
[7]	SENGUTTUVAN V, CHALASANI S, LUKE E A, et al. Adaptive mesh refinement using general elements:AIAA-2005-0927[R]. Reston, VA:AIAA, 2005.
[8]	张扬, 张来平, 赫新, 等. 基于自适应混合网格的脱体涡模拟[J]. 航空学报, 2016, 37(12):3605-3614. ZHANG Y, ZHANG L P, HE X, et al. Detached eddy simulation based on adaptive hybrid grids[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(12):3605-3614(in Chinese).
[9]	唐静, 崔鹏程, 贾洪印, 等. 非结构混合网格鲁棒自适应技术[J]. 航空学报, 2019, 40(10):122894. TANG J, CUI P C, JIA H Y, et al. Robust adaptation techniques for unstructured hybrid mesh[J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(10):122894(in Chinese).
[10]	TANG J, LI B, CHEN J, et al. Large scale parallel computing for fluid dynamics on unstructured grid[C]//15th International Symposium on Parallel and Distributed Computting, 2016.
[11]	陈刚,王磊,陆忠华,等. 万核级并行飞机气动模拟软件CCFD研制[J]. 华中科技大学学报(自然科学版),2011,39(增刊1):99-101. CHEN G, WANG L, LU Z H, et al. Development of ten-thousand-core parallel software CCFD for aircraft aerodynamics simulation[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2011, 39(Suppl. 1):99-101(in Chinese).
[12]	YANG C, XUE W, FU H, et al. 10M-core scalable fully-implicit solver for nonhydrostatic atmospheric dynamics[C]//Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, 2016.
[13]	XU C, DENG X, ZHANG L, et al. Parallelizing a high-order CFD software for 3D, multi-block, structural grids on the TianHe-1A supercomputer[C]//International Supercomputing Conference, 2016:26-39.
[14]	KAVOUKLIS C, KALLINDERIS Y. Parallel adaptation of general three-dimensional hybrid meshes[J]. Journal of Computational Physics, 2010, 229:3454-3473.
[15]	VARADARAJAN R, HWANG I. An efficient dynamic load balancing algorithm for adaptive mesh refinement[C]//ACM Symposium on Applied Computing, 1970:464-472.
[16]	PARK M A, DARMOFAL D L. Parallel anisotropic tetrahedral adaptation:AIAA-2008-0917[R]. Reston, VA:AIAA, 2008.
[17]	BOILLAT J, BRUCE F, KROPF P. A dynamic load-balancing algorithm for molecular dynamics simulation on multi-processor systems[J]. Journal of Computational Physics, 1991, 96:1-14.
[18]	CYBENKO G. Dynamic load balancing for distributed memory multi-processors[J]. Journal of Parallel and Distributed Computing, 1989, 7:279-301.
[19]	LEPAGE C Y, STCYR A, HABASHI W G. Parallel unstructured mesh adaptation on distributed memory systems:AIAA-2004-2532[R]. Reston, VA:AIAA, 2004.
[20]	PARK Y M, KWON O J. Unsteady flow computations using a 3-D parallel unstructured dynamic mesh adaptation algorithm:AIAA-2001-0865[R]. Reston, VA:AIAA, 2001.
[21]	LIAN Y Y, HSU K H, SHAO Y L, et al. Parallel adaptive mesh-refining scheme on a three-dimensional unstructured tetrahedral mesh and its applications[J]. Computer Physics Communications, 2006, 175:721-737.
[22]	ANTEPARA O, LEHMKUHL O, CHIVA J, et al. Parallel adaptive mesh refinement simulation of the flow around a square cylinder at Re=22000[J]. Procedia Engineering, 2013, 61:246-250.
[23]	OLIKER L, BISWAS R, GABOW H N. Parallel tetrahedral mesh adaptation with dynamic load balancing[J]. Parallel Computing, 2000, 26(12):1583-1608.
[24]	CAVALLO P A. Further extension and validation of a parallel unstructured mesh adaptation package:AIAA-2005-0924[R]. Reston, VA:AIAA, 2005.
[25]	ALAUZET F, LI X, SEOL E S, et al. Parallel anisotropic 3D mesh adaptation by mesh modification[J]. Engineering with Computers, 2006, 21:247-258.
[26]	DIGONNET H, COUPEZ T, LAURE P, et al. Massively parallel anisotropic mesh adaptation[J]. International Journal of High Performance Computing Applications, 2017, 33(1):3-24.
[27]	SCHLOEGEL K, KARYPIS G, KUMAR V. Wave front diffusion and lmsr:Algorithms for dynamic repartitioning of adaptive meshes[J]. IEEE Transactions on Parallel and Distributed Systems, 2001, 12(5):451-466.
[28]	唐静, 邓有奇, 马明生, 等. 飞翼气动优化中参数化和网格变形技术研究[J]. 航空学报, 2015, 36(5):1480-1490. TANG J, DENG Y Q, MA M S, et al. Parametrization and grid deformation techniques for flying-wing aerodynamic optimization[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(5):1480-1490(in Chinese).
[29]	CHEN J T, ZHANG Y B, ZHOU N C, et al. Numerical investigations of the high-lift configuration with MFlow solver[J]. Journal of Aircraft, 2015, 52(4):1051-1062.
[30]	ROE P L. Approximate Riemann solvers, parameter vectors, and difference schemes[J]. Journal of Computational Physics, 1981, 43:357-372.
[31]	KIM J S, KWON O J. Improvement on block LU-SGS scheme for unstructured mesh Navier-Stokes computations:AIAA-2002-1061[R]. Reston, VA:AIAA, 2002.
[32]	GONG X Q, CHEN J T, ZHOU N C, et al. The effects of turbulence model corrections on drag prediction of NASA common research model:AIAA-2014-4371[R]. Reston, VA:AIAA, 2014.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

主管单位：中国科学技术协会主办单位：中国航空学会北京航空航天大学

非结构混合网格自适应并行技术

Parallel algorithms for unstructured hybrid mesh adaptation

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献 32

相关文章 15

编辑推荐

Metrics

本文评价

[1]	王年华, 常兴华, 赵钟, 张来平. 非结构CFD软件MPI+OpenMP混合并行及超大规模非定常并行计算的应用[J]. 航空学报, 2020, 41(10): 123859-123859.
[2]	唐静, 崔鹏程, 贾洪印, 李彬, 李欢. 非结构混合网格鲁棒自适应技术[J]. 航空学报, 2019, 40(10): 122894-122894.
[3]	赖剑奇, 李桦, 张冉, 常青. 多GPU并行可压缩流求解器及其性能分析[J]. 航空学报, 2018, 39(9): 121944-121953.
[4]	张培红, 张耀冰, 周桂宇, 陈江涛, 邓有奇. 面向混合网格高精度阻力预测的熵修正方法[J]. 航空学报, 2018, 39(9): 122019-122030.
[5]	张培红, 张耀冰, 周桂宇, 陈江涛, 邓有奇. 面向非结构混合网格高精度阻力预测的梯度求解方法[J]. 航空学报, 2018, 39(1): 121415-121415.
[6]	孙俊磊, 王和平, 周洲, 雷珊. 螺旋桨滑流对菱形翼布局无人机气动的影响[J]. 航空学报, 2018, 39(1): 121431-121431.
[7]	王科雷, 祝小平, 周洲, 王红波. 低雷诺数分布式螺旋桨滑流气动影响[J]. 航空学报, 2016, 37(9): 2669-2678.
[8]	张扬, 张来平, 赫新, 邓小刚. 基于自适应混合网格的脱体涡模拟[J]. 航空学报, 2016, 37(12): 3605-3614.
[9]	张健, 邓有奇, 李彬, 张耀冰. 一种适用于三维混合网格的GMRES加速收敛新方法[J]. 航空学报, 2016, 37(11): 3226-3235.
[10]	崔鹏程, 邓有奇, 唐静, 李彬. 基于伴随方程的网格自适应及误差修正[J]. 航空学报, 2016, 37(10): 2992-3002.
[11]	张扬, 张来平, 赫新, 邓小刚. 基于非结构/混合网格的脱体涡模拟算法[J]. 航空学报, 2015, 36(9): 2900-2910.
[12]	黄飞, 苗文博, 程晓丽, 沈清. 一种DSMC方法的并行策略[J]. 航空学报, 2014, 35(4): 968-974.
[13]	康忠良, 闫超. 适用于混合网格的约束最小二乘重构方法[J]. 航空学报, 2012, 33(9): 1598-1605.
[14]	陈科;曹义华;安克文;李栋. 应用混合网格分析复杂积冰翼型气动性能[J]. 航空学报, 2007, 28(增): 87-91.
[15]	陈德华;吴文华;王晋军. 战斗机弹射救生系统数值模拟研究[J]. 航空学报, 2006, 27(5): 778-783.