Laminar natural convection in an enclosure divided by a complete vertical conducting partition is studied numerically by using a stationary two dimensional model. Two opposite vertical sides of the enclosure were at constant different temperatures. The working fluid was glycerol. The heat conducting partition was positioned in the middle of the enclosure. The governing equations with variable density and thermophysical characteristics of the fluid were used. An appropriate finite volume numerical procedure based on SIMPLE, uniform collocated grid and deferred correction central difference scheme for convective terms has been developed. Numerical procedure was verified by comparing the results with numerical bench mark solutions for a cavity without partition filled with air and also by comparing them with the available experimental results for a partitioned cavity filled with glycerol. The developed procedure takes into account the nonlinear dependence of the thermophysical characteristics of the glycerol (especially viscosity) so that a better agreement with the experimental results was achieved.