A series of two-dimensional numerical models of magmatism and mantle convection in small planetary bodies are presented to discuss how the small size of the Moon exerts control over its mantle evolution. Mantle convection is modeled by a solid-state convection of internally heated materials with temperature-dependent viscosity. Magmatism is modeled as upward permeable flow of basaltic magma generated by decompression melting of upwelling mantle materials. Migration of the generated magma causes compaction/expansion of the coexisting matrix. The volume change of matrix and the buoyancy of magma induced by magmatism enhance the upwelling flow and hence the magmatism itself, when the Rayleigh number exceeds a threshold that is around the critical Rayleigh number for the onset of thermal convection. This positive feedback makes magmatism vigorous and episodic and causes efficient cooling and stirring of the mantle in large planets like Venus and the Earth. In small planetary bodies where the Rayleigh number is lower than the threshold, however, magmatism occurs more continuously with a characteristic time of several hundred million years, cooling and convective stirring of the mantle caused by magmatism are less important, and a compositionally stratified structure develops in the mantle. These features of magmatism and mantle structure fit in with the observations on the Moon. The positive feedback is a key to understanding why the evolution of the lunar mantle is so different from that of lager planets like the Earth and Venus.