A 2-D Lattice Boltzmann Method（LBM） coupled with a Sub-Grid Stress（SGS） model is proposed and validated by flows around a non-submerged spur dike in a channel.And then the LBM-SGS model is further applied to flows in a pumping-station forebay.Shallow water equations are numerically solved by the LBM and the turbulence can be taken into account and modeled efficiently by the Large Eddy Simulation（LES） model.The bounce-back scheme of the non-equilibrium part of the distribution function is used at the inlet boundary,the normal gradient of the distribution function is set as zero at the outlet boundary and the bounce-back scheme is applied to the solid wall to ensure non-slip boundary conditions.Firstly,the model successfully predicts the flow characteristics around a spur dike,such as circulating flow,velocity and water depth distributions.The results are verified by the experimental data and compared to the results obtained by conventional Smagoringsky Model（SM） of LES.Finally,the LBM-SGS model is used to further predict the flow characteristics in a forebay,such as secondary flow and water level.The comparisons show that the model scheme has the capacity to simulate complex flows in shallow water with reasonable accuracy and reliability.
A lattice Boltzmann model combined with curvilinear coordinate is proposed for lid-driven cavity three-dimensional (3D) flows. For particle velocity distribution, the particle collision process is performed in physical domain, and the particle streaming process is carried out in the corresponding computational domain, which is transferred from the physical domain using interpolation method. For the interpolation calculation, a second-order upwind interpolation method is adopted on internal lattice nodes in flow fields while a second-order central interpolation algorithm is employed at neighbor-boundary lattice nodes. Then the above-mentioned model and algorithms are used to numerically simulate the 3D flows in the lid-driven cavity at Reynolds numbers of 100, 400 and 1000 on non-uniform meshes. Various vortices on the x-y, y-z and x-z symmetrical planes are successfully predicted, and their changes in position with the Reynolds number increasing are obtained. The velocity profiles of u component along the vertical centerline and w component along the horizontal centerline are both in good agreement with the data in literature and the calculated results on uniform meshes. Besides, the velocity vector distributions on various cross sections in lid-driven cavity predicted on non-uniform meshes are compared with those simulated on uniform meshes and those in the literature. All the comparisons and validations show that the 3D lattice Boltzmann model and all the numerical algorithms on non-uniform meshes are accurate and reliable to predict effectively flow fields.
An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.
