Poucos estudos investigam, em conjunto, os efeitos t?rmicos com os fen?menos associados aos processos iniciais da turbul?ncia. Em geral, estas quest?es s?o estudadas separadamente. Assim, o prop?sito deste trabalho ? usar o M?todo de Fronteira Imersa (MFI) para estudar escoamentos Newtonianos bidimensionais incompress?veis em torno de corpos aquecidos com temperatura constante ? superf?cie, considerando a transfer?ncia de calor por convec??o mista em conjunto com os processos iniciais da turbul?ncia (n?meros de Reynolds moderados) no intervalo [1; 5000]. Para as simula??es num?ricas, s?o utilizados os modelos de turbul?ncia de Smagorinsky e Spalart-Allmaras. O primeiro modelo utiliza a metodologia Large Eddy Simulation (LES) e baseia-se na hip?tese de equil?brio local para as pequenas escalas associadas ? hip?tese de Boussinesq, de tal forma que a energia injetada no espectro da turbul?ncia equilibra a energia dissipada pelos efeitos viscosos. O segundo modelo, utiliza o conceito URANS (Unsteady Reynolds Averaged Navier-Stokes Equations), utilizando apenas uma equa??o de transporte para a viscosidade turbulenta, sendo calibrada nas camadas de gradiente de press?o.