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.
关键词
Keywords
references
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.