We present a numerical study of steady convection in a two-dimensional mushy layer during the directional solidification of a binary mixture. The calculations reveal the internal structure of strongly nonlinear states featuring upflow which has been focused into solid-free regions known as chimneys. The mushy layer is modelled as a porous medium whose permeability is a function of the local solid fraction. The mushy layer is coupled to a chimney that is modelled as a narrow vertical channel where lubrication scalings are used to simplify the Navier–Stokes equations. We use these methods to exhibit solutions which give the detailed structure of the temperature, solute, flow and solid-fraction fields within the mushy layer. A key finding of the numerics is that there are two distinct chimney solutions at low Rayleigh numbers, presumably corresponding to stable and unstable portions of a subcritical solution branch. We also explore the relationship between convective solutions with and without chimneys.