Recently there was a flurry of activities in studying the effects of phase transitions in the Earth's mantle. From petrological and geophysical considerations, phase-transitions would also play an important role in venusian dynamics. The basic differences between the two planets are the surface boundary conditions, both thermally and mechanically. In this vein we have studied time-dependent mantle convection with multiple phase transitions and depth-dependent thermal expansivity (alpha is approximately rho(exp -6)), based on high-pressure and temperature measurements. Both the olivine-spinel and spinel-perovskite transitions were simulated by introducing an effective thermal expansivity, as described. Used together with the extended Boussinesq Approximation this method serves as a powerful tool to examine the effects of phase transitions on convection at relatively low computational costs.