A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusio...A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.展开更多
Numerical instability may occur when simulating high Reynolds number flows by the lattice Boltzmann method(LBM).The multiple-relaxation-time(MRT)model of the LBM can improve the accuracy and stability,but is still sub...Numerical instability may occur when simulating high Reynolds number flows by the lattice Boltzmann method(LBM).The multiple-relaxation-time(MRT)model of the LBM can improve the accuracy and stability,but is still subject to numerical instability when simulating flows with large single-grid Reynolds number(Reynolds number/grid number).The viscosity counteracting approach proposed recently is a method of enhancing the stability of the LBM.However,its effectiveness was only verified in the single-relaxation-time model of the LBM(SRT-LBM).This paper aims to propose the viscosity counteracting approach for the multiple-relaxationtime model(MRT-LBM)and analyze its numerical characteristics.The verification is conducted by simulating some benchmark cases:the two-dimensional(2D)lid-driven cavity flow,Poiseuille flow,Taylor-Green vortex flow and Couette flow,and threedimensional(3D)rectangular jet.Qualitative and Quantitative comparisons show that the viscosity counteracting approach for the MRT-LBMhas better accuracy and stability than that for the SRT-LBM.展开更多
文摘A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.
基金supported by the National Natural Science Foundation of China(NSFC,Grant Numbers 10572106,10872153 and 11172219)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20130141110013).
文摘Numerical instability may occur when simulating high Reynolds number flows by the lattice Boltzmann method(LBM).The multiple-relaxation-time(MRT)model of the LBM can improve the accuracy and stability,but is still subject to numerical instability when simulating flows with large single-grid Reynolds number(Reynolds number/grid number).The viscosity counteracting approach proposed recently is a method of enhancing the stability of the LBM.However,its effectiveness was only verified in the single-relaxation-time model of the LBM(SRT-LBM).This paper aims to propose the viscosity counteracting approach for the multiple-relaxationtime model(MRT-LBM)and analyze its numerical characteristics.The verification is conducted by simulating some benchmark cases:the two-dimensional(2D)lid-driven cavity flow,Poiseuille flow,Taylor-Green vortex flow and Couette flow,and threedimensional(3D)rectangular jet.Qualitative and Quantitative comparisons show that the viscosity counteracting approach for the MRT-LBMhas better accuracy and stability than that for the SRT-LBM.