中国空气动力研究与发展中心 计算空气动力研究所,绵阳,621000
张健(1989-), 男, 湖北钟祥人, 硕士, 研究方向:计算空气动力学.E-mail:j-zhang@foxmail.com
邱名(1984-), 男, 贵州织金人, 博士, 研究方向:叶轮机气体动力学.E-mail:qiu_ming_abc@163.com
纸质出版:2019
移动端阅览
张健, 唐静, 邱名, 邓有奇, 龚小权. 叶轮机械全环非定常大规模并行模拟程序设计[J]. 空气动力学学报, 2019, 37(4): 546-554.
ZHANG Jian, TANG Jing, QIU Ming, et al. Massive parallel program design for unsteady full annulus turbomachine simulation[J]. 2019, 37(4): 546-554.
张健, 唐静, 邱名, 邓有奇, 龚小权. 叶轮机械全环非定常大规模并行模拟程序设计[J]. 空气动力学学报, 2019, 37(4): 546-554. DOI: 10.7638/kqdlxxb-2017.0065.
ZHANG Jian, TANG Jing, QIU Ming, et al. Massive parallel program design for unsteady full annulus turbomachine simulation[J]. 2019, 37(4): 546-554. DOI: 10.7638/kqdlxxb-2017.0065.
为了更加精确地模拟叶轮机械内部复杂的三维非定常流动,开展了全环网格模拟方法探索,提出了针对滑移面的守恒变量插值和高效的大规模并行计算解决方案。基于千万亿次高性能计算平台,自主开发了非定常数值模拟程序,采用单级压气机、1.5级涡轮和3.5级高负荷压气机三个算例进行了程序的验证与确认,以及并行效率测试工作。模拟采用的最大网格数为2亿,计算核心数为4096个。结果表明:程序计算得到的整体气动性能与试验数据吻合良好;程序所采用的全环非定常计算模型能保证交界面处物理量的连续,削弱非物理熵增,得到更加可信的流场解;所采用的METIS的网格分区方法以及MPI并行策略使程序具有良好的负载均衡和并行效率;同时程序对超大规模的复杂问题具备良好的可扩展性和适应性,能够满足工程实际问题的流动分析和叶轮机械的精细流动研究要求。
In order to get an accurate simulation of the complex three-dimensional unsteady flows in turbomachine
an in-housed unsteady turbomachinery simulation program is developed based on a petaflops high performance computing platform. The full annulus grid simulation method has been developed. The conservative variable interpolation method for sliding interfaces and the high efficient massive parallel computation scheme are discussed in detail. A single stage compressor
a 1.5 stage turbine
and a 3.5 stage high-loaded compressor were simulated to verify and validate the program. The parallel efficiency of the program was also tested. The maximum grid number of the testcases is approximately 200 million
and the number of CPU cores is up to 4096. The results show that aerodynamic performance computed by this program agrees well with experimental data. The full unsteady computation method adopted in this paper guarantees variables' continuity at the sliding interface
weakens nonphysical entropy generation
and yields a reasonable flowfield. The grid partition methods by METIS and the parallel program scheme of MPI contribute to the program with a good load balance and a high parallel efficiency. This program has good expansibility and adaptability to the complex problem of large scale and satisfies the need of complex application flow analysis and exquisite turbomachiney flow filed study.
DENTON J D. The calculation of 3-D viscous flow through multistage turbomachines[J]. Journal of Turbomachinery, 1992, 114:18-26.
RAI M. Three-dimensional Navier-Stokes simulations of turbine rotor-stator interaction. I-Methodology[J]. AIAA Journal of Propulsion and Power, 1989, 5(3):305-311.
HE L, NING W. Efficient approach for analysis of unsteady viscous flows in turbomachines[J]. AIAA Journal, 1998, 36(11):2005-2012.
CHEN J P, BRILEY W R. A parallel flow solver for unsteady multiple blade row turbomachinery simulations: GT-2001-348[R]. New Orleans: ASME Turbo Expo, 2001.
CHEN J P, HATHAWAY M D, HERRICK G P. Pre-stall behavior of a transonic axial compressor stage via time-accurate numerical simulation: 2007-GT-29627[R]. Montreal: ASME Turbo Expo, 2007.
SHAO F, QI L, LUO J Q, et al. Numerical simulation of full scale unsteady flow in high pressure turbine[J]. Aeroengine, 2014, 40(2):31-37. (in Chinese)
FU H, LIAO J, YANG J, et al. The Sunway TaihuLight supercomputer:system and applications[J]. Science China Information Sciences, 2016, 59:1-16.
LIAO X K, XIAO L Q, YANG C Q, et al. MilkyWay-2 supercomputer:system and application[J]. Frontiers of Computer Science Selected Publications from Chinese Universities, 2014, 8(3):345-356.
SLOTNICK J, KHODADOUST A, ALONSO J, et al. CFD vision 2030 study: a path to revolutionary computational aerosciences[R]. NASA CR-2014-218178.
ZHANG Y B, DENG Y Q, WU X J, et al. Drag prediction of DLR-F6 using MFlow unstructured mesh solver[J]. Acta Aerodynamica Sinica, 2011, 29(2):163-169. (in Chinese)
YOON S, JAMESON A. Lower-upper symmetric gauss-seidel method for the Euler and Naiver-Stokes equations[J]. AIAA Journal, 1988, 26(9):1025-1026.
ROE P L. Approximate Riemann solvers, parameter vectors, and difference schemes[J]. J Computational Physics, 1981, 43:357-372.
VENKATAKRISHNAN V. On the accuracy oflimiters and convergence to steady state solutions. AIAA-93-0880[R]. Reston AIAA, 1993.
SPALART P R, ALLMARAS S R. A one equation turbulence model for aerodynamic flows. AIAA-92-0439[R]. Reston AIAA, 1992.
KARYPIS G, KUMAR V. A fast and high quality multilevel scheme for partitioning irregular graphs[J]. SIAM Journal on Scientific Computing, 1998, 20(1):359-392
BONET J. An alternating digital tree (ADT) algorithm for 3D geometric searching and insersection problems[J]. International Journal for Numerical Methods in Engineering, 1991, 31:1-17.
TANG Z G, LI B, ZHENG M, et al. Store seperation simulation using overset unstructured grid[J]. Acta Aerodynamica Sinica, 2009, 27(5):592-596 (in Chinese)
REID L, MOORE R D. Performance of a single-stage axial-flow transonic compressor with rotor and stator aspect ratios of 1.19 and 1.26, respectively, and with design pressure ratio of 1.82[R]. NASA TP-1338, 1978.
ZHANG J, TANG Z G, DENG Y Q, et al. Three-dimensional numerical simulations for transonic axial compressors using a mixing plane method[J]. Acta Aerodynamica Sinica, 2015, 33(5):631-635. (in Chinese)
HERRICK G P, HATHAWAY M D, CHEN J P. Unsteady full annulus simulations of a transonic axial compressor stage: AIAA-2009-1059[R]. Orlando: 47th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, 2009.
LIU J, BAI X Z, WANG W. High accurate method for transferring flow information among meshs and its validation[J]. Acta Aerodynamica Sinica, 2010, 28(1):1-6. (in Chinese)
HEMBERA M, LOOS A, KAU H P. Comparing a 40 Mio gridpoints full-annulus computation with a 7 Mio gridpoints nonlinear harmonic computation: AIAA-2009-930[R]. Orlando: 47th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, 2009.
GALLUS H E, WALRAEVENS R E. Stator-Rotor-Stator interaction in an axial flow turbine and its influence on loss mechanism. AGARD Report No. CP-571[R]. Derby: AGARD PEP 85th Symposium, 1995.
0
浏览量
0
下载量
0
CSCD
关联资源
相关文章
相关作者
相关机构
京公网安备11010802024621