To effectively design a large furnace for producing large-size AlN crystals, a fully coupled compressible flow solver was developed to study the sublimation and mass transport processes in AlN crystal growth. Compressible effect, buoyancy effects, flow coupling between aluminum gas and nitrogen gas, and Stefan effect are included. Two sets of experimental data were used to validate the present solver. Simulation results showed that the distributions of Al and N-2 partial pressures are opposite along the axial direction due to constant total pressure and Stefan effect, with the Al and nitrogen partial pressures being highest at the source and seed crystals positions, respectively. The distributions of species inside the growth chamber are obviously two-dimensional, which can curve a flat crystal surface. Simulation results also showed that AlN crystal growth rate can be increased by reducing total pressure or by increasing seed temperature or by increasing source-seed temperature difference. High nitrogen pressure causes decrease in growth rate, but it is beneficial for obtaining uniform growth rate in the radial direction. Results of simulation also showed that there is an optimized temperature difference (40 degrees C) in the present furnace for obtaining good homogeneity of growth rate. (C) 2011 Elsevier B.V. All rights reserved.