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 conceitoURANS (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.