Numerical calculation from the full differential equations for convection in an unbounded region is expensive, and hence approximate solutions are important. For small values of the Rayleigh number Ra, perturbation methods are appropriate. At large values of Ra thermal boundary layers are formed, and boundary layer theory is the obvious method of investigation. This approach forms the subject of much of this chapter. We follow to a large extent the discussion by Cheng (1985a).