Se presentan resultados numéricos en dos dimensiones de flujos isotérmicos/térmicos de fluidos incompresibles y viscosos modelados matemáticamente por las ecuaciones de Navier-Stokes y por la aproximación de Boussinesq respectivamente en la formulación velocidad-vorticidad. Los flujos térmicos corresponden tanto a convección natural como a convección mixta. Los resultados se obtienen con un método simple el cual se basa en un proceso iterativo de punto fijo para resolver el sistema elíptico no lineal que resulta después de una discretización apropiada en el tiempo. El proceso iterativo lleva a la solución de problemas elípticos lineales, simétricos, bien condicionados y desacoplados para los cuales existen resolvedores eficientes sin importar la discretizacón espacial.