This contribution describes preliminary results dealing with computation of the turbulent flow around a square cylinder by 2D implicit large eddy simulation. The solid wall boundary conditions were described by the immersed boundary method.