Natural convection between concentric spheres

A finite-difference solution for steady natural convective flow in a concentric spherical annulus with isothermal walls has been obtained. The stream function-vorticity formulation of the equations of motion for the unsteady axisymmetric flow is used; interest lying in the final steady solution. Forward differences are used for the time derivatives and second-order central differences for the space derivatives. The alternating direction implicit method is used for solution of the discretization equations. Local one-dimensional grid adaptation is used to resolve the steep gradients in some regions of the flow at large Rayleigh numbers. The break-up into multi-cellular flow is found at high Rayleigh numbers for air and water, and at significantly low Rayleigh numbers for liquid metals. Excellent agreement with previous experimental and numerical data is obtained.