The aim of this work is to study the distribution of the thermal field inhomogeneity in the melt along the crystal growth axis for different growth methods. Methodology. Within the framework of the SolidWorks Premium complex, in the Flow Simulation package for the constructed geometric models of the crystallization apparatus of the Bridgman, Czochralski and HDC (horizontal directional crystallization) methods, a curvilinear block-structured base mesh was formed. The conjugate problem was solved in the thermal unit and the water-cooled jacket to study the processes of heat-mass transfer and the distribution of the thermal field inhomogeneity in the melt along the crystal growth axis for different growth methods. Results. The results of computer simulation of the distribution of temperature profiles over the vertical and cross sections of the crystals grown by the Bridgman and Czochralski methods are presented. Pictures of the characteristic arrangement of facets are obtained. The obtained data fully agree with the known corresponding experimental results for yttrium-aluminum garnet crystals. Research implications. The analysis of the obtained results of computer simulation, carried out on the basis of consideration of the general theory of crystal growth mechanisms and morphological instability of the crystallization front for the three main methods of growth, confirmed the correctness of the obtained model pictures of the appearance of faceted growth forms.