In turbulent simulation by the finite element method, when using the wall function and dealing with the tangential direction velocity of nodes on the body boundary as unknown variables, the tangential direction cannot be determined uniquely at the corner of nodes on the body. In this study, we considered that the tangential direction at the corner of the body will change according to the flow field. We evaluated the tangential direction of the velocity at the corner node by using of local velocities and inclinations of adjacent element edges along the body. By this way, we simulated a flow around the two-dimensional square cylinder, and computed aerodynamic forces acting to the square cylinder which are enough close to those of existing studies.