In this work a technique of implementing boundary conditions for Lattice Boltzmann method is given. In particular it focuses on the description of the wall-turbulence interaction. To this end a robust multidimensional WENO (Weighted Essential Non Oscillatory) interpolation technique is adopted, using the large scale strain tensor components. The subgrid-scale model is a conventional Smagorinsky model. This technique is tested on several classic problems, like flows past on flate plate, channel flows and flows past transverse thin plate. Numerical experiments show good efficiency and numerical stability.