Velocity–vorticity RANS turbulence modeling by boundary element method

Turbulent flow over various geometries is studied numerically. Incompressible set of Navier–Stokes equations is considered and solved by boundary domain integral method (BDIM). Governing equations are written in velocity–vorticity form. Turbulence models used are based on eddy-viscosity hypothesis. Integral form of equations, discretization and the solution algorithm are presented. The algorithm is tested with two separate test cases. The first is the turbulent channel flow for two different Reynolds numbers: Re τ = 180 and Re τ = 395 . Results show very good agreement with corresponding DNS data. The second test case is the flow over backward facing step for Reynolds number Re h = 5000 , which shows good agreement with literature data on mean reattachment length.