Aiming at high accuracy and high resolution computation for large eddy simulation of compressible turbulence
a qusi-6th order central-WENO(weighted essentially non-oscillatory)hybrid scheme with adaptively controlled numerical dissipation is developed by introducing a local pressure gradient sensor and the upper and lower bounds of the weighting function. The mathematical characteristics of the developed hybrid scheme are theoretically analyzed with Fourier method. The 1D shock/density wave interaction problem and the 3D large eddy simulation of compressible isotropic turbulence are numerically analyzed to verify the validity and effectiveness of the scheme. The results show that the quasi-6th order central-WENO hybrid scheme has lower numerical dissipative error than the 5th order WENO scheme and can hold the high resolution for both the discontinuities and physical fluctuations in the flow fields. The decay curves of the turbulent statistics obtained by large eddy simulation with the central-WENO hybrid scheme coincide well with the available direct numerical simulation
and the -5/3 power law of the theoretical energy spectrum is successfully captured. The specific bounds of the weighting function applicable to the large eddy simulation of compressible flow are calibrated. This work lays a solid algorithm foundation for high accuracy large eddy simulation of compressible turbulence in fluid machinery.
关键词
Keywords
references
JIANG G S, SHU C W. Efficient implementation of weighted ENO Schemes [J]. Journal of Computational Physics, 1996, 126(1): 202-228.
MARTIN M P, TAYLOR E M, WU M, et al. A bandwidth-optimized WENO scheme for the effective direct numerical simulation of compressible turbulence [J]. Journal of Computational Physics, 2006, 220(1): 270-289.
KIM D, KWON J H. A high-order accurate hybrid scheme using a central flux scheme and a WENO scheme for compressible flowfield analysis [J]. Journal of Computational Physics, 2005, 210(2): 554-583.
BERMEJO-MORENO I, BODART J, LARSSON J, et al. Solving the compressible Navier-Stokes equations on up to 1.97 million cores and 4.1 trillion grid points [C]∥Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis. New York, USA: ACM, 2013: 1-10.
SAMTANEY R, PULLIN D I, KOSOVIC B. Direct numerical simulation of decaying compressible turbulence and shocklet statistics [J]. Physics of Fluids, 2001, 13(5): 1415-1430.
LI Xinliang, FU Dexun, MA Yanwen. Direct numerical simulation of compressible homogeneous isotropic turbulence [J]. Science in China: Series A, 2002, 32(8): 716-724.
LI Hu, ZHANG Shuhai. Direct numerical simulation of compressible isotropic decaying turbulence [J]. Chinese Journal of Theoretical and Applied Mechanics, 2012, 44(4): 673-686.
ERLEBACHER G, HUSSAINI M Y, SPEZIALE C G, et al. Toward the large-eddy simulation of compressible turbulent flows [J]. Journal of Fluid Mechanics, 1992, 238: 155-185.
KOTOV D V, YEE H C, WRAY A A, et al. Numerical dissipation control in high order shock-capturing schemes for LES of low speed flows [J]. Journal of Computational Physics, 2016, 307(1): 189-202.
NICOUD F, DUCROS F. Subgrid-scale stress modelling based on the square of the velocity gradient tensor [J]. Flow, Turbulence and Combustion, 1999, 62(3): 183-200.
LIOU M S. A sequel to AUSM: part II AUSM+-up for all speeds [J]. Journal of Computational Physics, 2006, 214(1): 137-170.
JAMESON A, SCHMIDT W, TURKEL E. Numerical solutions of the Euler equations by finite volume methods using Runge-Kutta time-stepping schemes: 81-1259 [R]. Reston, VA, USA: AIAA, 1981.
SHU C W, OSHER S. Efficient implementation of essentially non-oscillatory shock-capturing schemes [J]. Journal of Computational Physics, 1988, 77(2): 439-471.
Large Eddy Simulation of Flow Resistance and Heat Transfer Characteristics of Dual-Swirl Cooling on Turbine Blade Leading Edge
Research on the Thermal Mixing Mechanism of Fluids in a Circular T-Junction under Multi-Force Field Coupling Conditions
Large Eddy Simulation of Turbulent Flow and Heat Transfer in a Square Duct with Rayleigh-Benard Convection on the Cross Section
Study on the Special Vortical Structure Development During Flow Separation Control by VGJs in Boundary Layer of Flat Plate with Adverse Pressure Gradient
Large Eddy Simulation for Control of Boundary Layer Separation with Dimples
Related Author
LIU Xuebin
SONG Liming
YU Boyang
TAO Zhi
LI Jun
WANG Qiuwang
LI Anqi
SU Bo
Related Institution
MOE Key Laboratory of Thermal Fluid Science and Engineering, Xi'an Jiaotong University
Faculty of Infrastructure Engineering, Dalian University of Technology
Department of Fluid Machinery and Engineering, Xi'an Jiaotong University
College of Naval Architecture and Ocean Engineering, Naval University of Engineering
School of Energy and Power Engineering, Xi’an Jiaotong University