A comprehensive mathematical model of the traveling heater method is applied to understand interactions of flow, segregation, and stability for the crystal growth by the traveling heater method. We demonstrate that the formation of a secondary, counter-rotating vortex near the growth interface limits the transport of segregated solvent, leading to large concentration gradients and supercooling of the liquid near the interface. The secondary vortex arises from the same mechanism responsible for the formation of lee waves in atmospheric flows, and its spatial form scales with the Brunt-Vaisala frequency. Significantly, this supercooled layer of liquid arises from the lateral transport of solute, which is predominantly affected by the secondary flow vortex, rather than the axial diffusion of solute, as assumed in the classical derivation of constitutional supercooling by Mullins and Sekerka. Thus, supercooling in a model cadmium zinc telluride system occurs at growth rates of order 1 mm/day, nearly an order of magnitude smaller than expected from classical arguments. Paradoxically, the traditional strategy to alleviate constitutional supercooling during crystal growth, namely increasing the temperature gradient at the interface, is expected to strengthen the lee-wave vortex and accentuate the onset of instability. Other approaches will be needed to overcome growth limits in the traveling heater method. (C) 2015 Elsevier B.V. All rights reserved.