Stable and accurate numerical analytical method even at high Rayleigh numbers is desired, and in addition flexibility and economical efficiency are very important for numerical analysis. Unsteady flows of natural convention in a square cavity are investigated using the GSMAC (generalized and simplified marker and cell) finite-element method at high Rayleigh numbers from 10 {sup 6} to 10 {sup 8}. Validities of the multi-pass algorithm and BTD (balancing tensor diffusivity), which are highly accurate solutions for finite-element method, are investigated by introducing them to the GSMAC finite-element method. As the result, it was found that the multi-pass algorithm has little effect for the GSMAC finite-element method. The steady solutions at the Rayleigh numbers of 10 {sup 5} and 10 {sup 7} agreed well with the bench mark solutions and numerical solutions. Complex transient phenomena at Rayleigh numbers higher than 10 {sup 7} are successfully observed. 12 refs., 11 figs., 4 tabs.