Abstract
Recently it is widely recognized by LES and higher order upwind scheme without turbulence models that 3D calculation is needed for a flow around a square cylinder to consider a 3D structure of turbulence and its decay. However k-ε with a 2D calculation has possibility to represent turbulence energy balance because they are involved in it with a mathematical form. In this study modified production eddy viscosity model was applied to overcome difficulties of standard one for a flow around a stationary cylinder. This modified model was also applied to flows around vibrating cylinders in the along and across wind direction to clarify its effectiveness and limitation.