Frictional loss in turbulent regime is abnormally increased compared with in laminar regime. Thus the consideration of temperature rise across fluid film is significant in analysis and conventional isothermal theory loses its usefulness for performance prediction. This paper proposes to the three-dimensional thermohydrodynamic analysis of finite journal bearings operating under turbulent condition using two-equation model($kappa-varepsilon$ model) proposed by Hassid & Poreh. The equations are solved numerically by finite difference method. We make the analysis applicable even at large eccentricity when back flow of the lubricants occurs and axial flow is no longer ignored compared to circumferential flow.