Approximate (variational) and exact (numerical) solutions of Bénard type problems with temperature dependent material properties