This paper is the first part of a two-part study, and it presented numerical research of double-diffusive natural convection within an annulus area, situated in two horizontal confocal elliptic solids charged by a Newtonian fluid. The elliptical coordinates were used to transform the physical domain into a rectangular one. To resolve the governing equations and the boundary layer conditions, a calculator code based on the finite volume approach was developed. The details of the influences of thermal Rayleigh number on heat and mass transfer were investigated. Results obtained were compared with those existing in other reasearch works. K E Y W O R D S double-diffusive, elliptical cylinders, natural convection, Rayleigh number, stream function formulation