A scalable parallel solver is developed to simulate the Earth’s core convection. With the help from the “multiphysics” data structure and the restricted additive Schwarz preconditioning in PETSc the iterative solution of the linear solver converges rapidly at every time-step. The solver gains nearly 20 times speedup compared to a previous solver using least-squares polynomial preconditioning in Aztec. We show the efficiency and effectiveness of our new solver by giving numerical results obtained on a BlueGene/L supercomputer with thousands of processor cores.