We consider the problem of three-dimensional non-linear buoyant convection in ternary solidification. Under the limit of large far-field temperature, the convective flow is modeled to be in a rectangular cube composing of a horizontal liquid layer above a primary mushy layer, which itself is over a secondary mushy layer. We first apply linear stability analysis to calculate the conditions at the onset of motion. Next, we carry out weakly non-linear analyses to determine solutions in the form of hexagons and their possible stability and to obtain information about tendency for chimney formation. We find that if the flow is driven either from both mushy layers with equal critical conditions at the onset of motion or only by the primary mushy layer, then the flow can be in the form of a double-cell structure vertically with down-hexagons below or above up-hexagons. There is tendency for vertically oriented chimney formation at different horizontal locations in each mushy layer. For the cases where only the critical conditions at the onset of motion are equal in both mushy layers and depending on the values of the mush Rayleigh numbers, the flow can be subcritical (or supercritical) in both mushy layers or mixed subcritical in one layer and supercritical in another layer.