An explanation of the return flow onset behavior in a long horizontal CVD reactor is proposed on the basis of combined; analytical and numerical approach. This approach allows one to obtain the condition of vortex onset in the form, that includes not only Grashof and Reynolds numbers, but also the temperature difference at the reactor top and bottom walls. Also some three-dimensional effects are investigated, which accompany the return flow formation in a horizontal CVD reactor.