Periodicity Analysis and Benchmark Solution of Lid-driven Cavity at Re=15 000

Abstract

Abstract: We study periodic cavity flow at Re=15 000, with a high quality non-uniform mesh obtained by boundary-fitted coordinated method and extra quadtree encryption method based on structural grid. Fourier power spectra analysis and phase-space trajectory were applied to understand periodic evolution and kinetic energy distribution of the flow field. Fundamental frequency and corresponding period were determined. It was found that during periodic evolution process of the flow field, the primary vortex maintained a low frequency of kinetic energy, while high-order vortices near boundaries and corners showed high frequencies of kinetic energy. Benchmark solutions of two characteristic states in period evolution were obtained.

Key words: lid-driven cavity, periodic evolution, boundary-fitted coordinated method, energy distribution, benchmark solution

CLC Number:

O357.1

ZHAN Jiemin, BIAN Shiqi, LUO Yingying, HU Wenqing, GONG Yejun. Periodicity Analysis and Benchmark Solution of Lid-driven Cavity at Re=15 000[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2020, 37(4): 413-421.

References

[1] ERTURK E. Discussions on driven cavity flow[J]. International Journal for Numerical Methods in Fluids, 2009, 60(3):275-294.
[2] GHIA U, GHIA K N, SHIN C T. High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method[J]. Journal of Computational Physics, 1982, 48(3):387-411.
[3] ERTURK E, CORKE T C, GOKCOL C. Numerical solutions of 2-D steady incompressible driven cavity flow at high Reynolds numbers[J]. International Journal for Numerical Methods in Fluids, 2005, 48(7):747-774.
[4] ERTURK E, GOKCOL C, DURSUN B, et al. Fourth order compact formulation of Navier-Stokes equations and driven cavity flow at high Reynolds numbers[J]. International Journal for Numerical Methods in Fluids, 2006, 50(4):421-436.
[5] ABDELMIGID T A, SAQR K M, KOTB M A, et al. Revisiting the lid-driven cavity flow problem:Review and new steady state benchmarking results using GPU accelerated code[J]. Alexandria Engineering Journal, 2016, 56(1).
[6] BRUNEAU C H, SAAD M. The 2D lid-driven cavity problem revisited[J]. Computers & Fluids, 2006, 35(3):326-348.
[7] HOU S, ZOU Q, CHEN S, et al. Simulation of cavity flow by the lattice Boltzmann method[J]. Journal of Computational Physics, 1995, 118(2):329-347.
[8] GUO Z, SHI B, WANG N. Lattice BGK model for incompressible Navier-Stokes equation[M]. Academic Press Professional, Inc, 2000.
[9] FRAGA Filho C A D, CHACALTANA J T A, PINTO W J N. Meshless Lagrangian SPH method applied to isothermal lid-driven cavity flow at low-Re numbers[J]. Computational Particle Mechanics, 2018(11):1-9.
[10] YU P X, TIAN Z F. An upwind compact difference scheme for solving the stream function-velocity formulation of the unsteady incompressible Navier-Stokes equation[J]. Computers & Mathematics with Applications, 2018, 75(9):3224-3243.
[11] YANG X C, SHANG Y Q. Two-level subgrid stabilized methods for Navier-Stokes equations at high Reynolds numbers[J]. Chinese Journal of Computational Physics, 2017, 34(6):657-665.
[12] LI L X. Parallel finite element computation of impressible viscous flows based on block preconditioning strategy[J]. Chinese Journal of Computational Physics, 2018, 35(2):151-160.
[13] BRUNEAU C H, JOURON C. An efficient scheme for solving steady incompressible Navier-Stokes equations[J]. Journal of Computational Physics, 1990, 89(2):389-413.
[14] PENG Y F, SHIAU Y H, HWANG R R. Transition in a 2-D lid-driven cavity flow[J]. Computers & Fluids, 2003, 32(3):337-352.
[15] HUSER A, BIRINGEN S. Calculation of two-dimensional shear-driven cavity flows at high Reynolds numbers[J]. International Journal for Numerical Methods in Fluids, 1992, 14(9):1087-1109.
[16] BAO Y, YE M X, LUO J H. An efficient parallel direct method for turbulent thermal convection[J]. Chinese Journal of Computational Physics, 2017, 34(6):651-656.
[17] SUN J C, DU P, LI P S, et al. Natural convection heat transfer in a porous medium with partially thermally active walls:Lattice Boltzmann method[J]. Chinese Journal of Computational Physics, 2017, 34(5):583-592.
[18] ZHAN J M, LI Y S. Numerical simulation of unsteady flow in natural and thermosolutal convection using boundary-fitted coordinate system[J]. International Journal of Numerical Methods for Heat & Fluid Flow, 2003, 13(8):1031-1056.
[19] ZHAN J M, LUO Y Y, LI Y S. A high accuracy hybrid method for two-dimensional Navier-Stokes equations[J]. Applied Mathematical Modelling, 2008, 32(5):873-888.
[20] LAL R, LI Z. Sensitivity analysis of a mesh refinement method using the numerical solutions of 2D lid-driven cavity flow[J]. Journal of Mathematical Chemistry, 2015, 53(3):844-867.
[21] XIA Z H, SHI Y P, CAI Q D, et al. Multiple states in turbulent plane Couette flow with spanwise rotation[J]. Journal of Fluid Mechanics, 2018, 837:477-490.
[22] HUANG Y H, XIA Z H, WAN M P, et al. Numerical investigation of plane Couette flow with weak spanwise rotation[J]. Science China:Physics, Mechanics & Astronomy, 2019, 62:044711.
[23] XIE Y C, XIA K Q. Turbulent thermal convection over rough plates with varying roughness geometries[J]. Journal of Fluid Mechanics, 2017, 825:573-599.