Calculation of lid-driven cavity flow using GAO-YONG turbulence equations

A two-dimensional stationary incompressible lid-driven cavity flow was numerically simulated employing staggered grid arrangement, incompressible SIMPLER method and QUICK scheme to solve GAO-YONG turbulence equations. The vorticity contour in the cavity and the velocity changing curve of cavity's central horizontal section line were acquired and compared with directly numerical simulation(DNS) result qualitatively and the consistency was very good. It shows that the mechanical energy equation corresponding to the scales of nonlinearity in the turbulence equations is very important for the turbulence equations to reasonably simulate the complex distribution of eddy viscosity coefficients and the inversion phenomenon of turbulence energy in the lid-driven cavity flow. Farther the GAO-YONG turbulence equations being able to gain real viscosity field in complex flow was proved.