A numerical study of heat island flows in an open domain: Stationary solutions