We report three-dimensional simulation results of thermal mixing in square T-junction configuration at high Reynolds number. The square T-junction configuration with different channel sizes was studied experimentally at the Department of Mechanical Engineering of Mie University, Japan. The T-Junction was selected as a benchmark for thermal mixing in the ERCOFTAC Workshop held in EDF Chatou, France, 2011. RANS, URANS and SAS simulations were performed with CFD code using finite volume method The comparison shows excellent agreement between SAS and the data.