Numerical study of the Boussinesq model of natural convection in an annular space: having a horizontal axis bounded by circular and elliptical isothermal cylinders