Three-dimensional simulations of gas flow and chemical reactions are made. The epitaxy model is newly developed. It is found that the gas flow rate Q is a key parameter which controls the spatial distribution of the growth rate above the susceptor. Here, Q is proportional to pressure and flow velocity. At low Q, the peak position of the growth rate shifts in the upstream direction. At an optimum Q, the growth rate distribution is spatially uniform. At high Q, the distribution is a split one with respect to the reactor symmetry plane. This is due to the peaking of gas flow velocity around the symmetry plane, which is driven by a pair of large return-flows around tapered gas entrance region.Introduction Gallium nitride (GaN) has been confirmed as a key material for achieving violet laser diodes (LDs) through recent splendid development [1,2]. Metalorganic vapor phase epitaxy (MOVPE) is a leading technology for this purpose. This technology is, however, not a satisfactory level for high power commercial LDs from the standpoint of crystal quality itself and the overall epitaxy process control. Numerical calculation is one of the powerful tools for analyzing the physical and chemical mechanism and finding the pathway to an optimal process condition. In this paper, three-dimensional simulations of gas flow and chemical reactions are made in a horizontal reactor for GaN MOVPE. Attention is paid to the effect of the gas flow velocity and gas pressure on the spatial distribution of both streamlines and the surface reaction rate.